By using four particle cluster model Hamiltonian along with third and fourth order phonon anharmonic interaction terms for KDP-type crystals, expressions for renormalized soft mode frequency and inverse dielectric susceptibility have been obtained. The method of double time temperature dependent Green’s function has been used for calculation. Fitting the values of model parameters, the temperature dependence of soft mode frequency and inverse dielectric susceptibility has been evaluated. The inverse dielectric susceptibility was observed increases linearly with temperature, for KDP crystal, in its paraelectric phase. Theoretical results are compared with experimental results of Kim et al [Physics Review B73, 134114 (2006) and Current Applied Physics 9(2009) 1307]. A good agreement has been found between the reported and present results.