Model and simulation:procedure OP;(r,t)/ot=-Lp:8F/8P (r,1), Ou;(r,t)/Ot=-L8F/8u;(r,1), (i=x,y以 Initial State Unpoled State Poled State T>>T;E=0 T=0;E=0.0 pole T=0;E=0.0 T=T+△T T=T+△T E=E+AE E=E+AE Model and simulation: procedure ∂Pi (r,t)/∂t= − L P·δF/ δ Pi (r, t), ∂ ui (r,t)/∂t= − L u·δF/ δ ui (r, t), (i=x,y)  Initial State T>>T c; E=0 Unpoled State T=0; E=0.0 T=T+ ∆ T ZFC pole E = E + ∆ E T=T+ ∆ T E = E + ∆ E Poled State T=0; E=0.0