Morphotropic tetragonal and rhombohederal phase boundary(MPB) in PZT system was calculated and the effect of mechanical stress and electric field on the MPB was examined using phenomenological thermodynamic theory proposed by Devonshire. Mechanical stress and electric field was calculated to change free energies of phases and shift MPB. The stable composition range of the phase, of which spontaneous strain and spontaneous polarization aligned to the direction of tensile stress or electirc field, was expanded.