function?[p?e?v?x?y?k]=possion（hMxMykmaxep）??%?g-s迭代法解五點(diǎn)差分法問(wèn)題?
%kmax為最大迭代次數(shù)??
%mn為xy方向的網(wǎng)格數(shù)，例如（2-0）/0.01=200;?
%e為誤差，p為精確解?
syms?temp;??
v=zeros（Mx+1My+1）;?
u1=zeros（Mx+1My+1）;
u2=zeros（Mx+1My+1）;
u3=zeros（Mx+1My+1）;
f=zeros（Mx+1My+1）;
x=0+（0:Mx）*h;?
y=0+（0:My）*h;?
w=4/pi*atan（0.2144*Mx+0.4822）;
w=1;
for?i=Mx/4:Mx*3/4??????
v（i1）=1500;??????
end
t=zeros（My-1Mx-1）;
for?k=1:kmax
????for?i=2:Mx
????????for?j=2:My
????????????%temp=（1-w）*u（ij）+w*（h*h*f（ij）/4+（u（ij+1）+u（ij-1）+u（i+1j）+u（i-1j）））/4;
????????????temp=-h*h*f（ij）/4+（u（ij+1）+u（i