五点差分格式

function L=Laplace_N(x0,y0,Lx,Ly,h) I=Lx/h;J=Ly/h; for(i=1:I+1) u(1,i)=feval('Bound_Con',x0+(i-1)*h,y0); u(J+1,i)=feval('Bound_Con',x0+(i-1)*h,y0+Ly);