Then a simple analytical formula is obtained, which expresses the constitutive relation between polarization P and loading electric field E of ferroelectric materials.