possion方程

function [p e v x y k]=possion（hMxMykmaxep）  % g-s迭代法解五点差分法问题 
%kmax为最大迭代次数  
%mn为xy方向的网格数，例如（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