function?[thitan0mm]?=?regulation（xymthitan0）
%UNtitleD3?此處顯示有關此函數的摘要
%???此處顯示詳細說明
%%?定義初值
Nx=7;??%%領域的大小
n=（Nx-1）/2;
lamda=2;?%%?控制調制項的參數?1或者2
tao=10;%%τ的范圍是在（0.110）之間

%%?設置Ur的初值，即在（xy）處的梯度

k=1;%%k表示循環次數，k=1是初值
errf=1;
while?errf>1e-8&&k<5
for?kesi=（x-n+1）:（x+n-1）?%%x表示行，y表示列
????for?yita?=（y-n+1）:（y+n-1）
????????
????????p（kesiyita）=thitan0（1）+thitan0（2）*（x-kesi）+thitan0（3）*（y-yita）;
????????U_phi0（kesiyita）=2*（In（kesiyita）-cos（p（kesiyita）））*sin（p（kesiyita））-2*lamda*（thitan0（1）-p（kesiyita）