用差分法求解偏微分方程

function u = peEllip5（nxminxmaxxnyminymaxy）
format long;
hx = （maxx-minx）/（nx-1）;
hy = （maxy-miny）/（ny-1）;
u0 = zeros（nxny）;
for j=1:ny
    u0（j1） = EllIni2Uxl（minxminy+（j-1）*hy）;
    u0（jnx） = EllIni2Uxr（maxxminy+（j-1）*hy）;
end
for j=1:nx
    u0（1j） = EllIni2Uyl（minx+（j-1）*hxminy）;
    u0（nyj） = EllIni2Uyr（minx+（j-1）*hxmaxy）;
end 

A = -4*eye（（nx-2）*（ny-2）（nx-2）*（ny-2））;
b = zeros（（nx-2）*（ny-2）1）;
for i=1:（nx-2）*（ny-2）
    if mod（inx-2） == 1
        if i==1
            A（12） = 1;
            A（1nx-1） = 1;
            b（1） = - u0（12） - u0（21）;
        else
            if i == （ny-3）*（nx-2）+1
                A（ii+1） = 1;
                A（ii-nx+2） = 1;
                b（i） = - u0（ny-11） - u0（ny2）;
            else
                A（ii+1） = 1;
                A（ii-nx+2） = 1;
                A（ii+nx-2） = 1;
                b（i） = - u0（floor（i/（nx-2））+21）;
            end
        end
    else
        if mod（inx-2） == 0
            if i == nx-2
                A（ii-1） = 1;
                A（ii+nx-2） = 1;
                b（i） = - u0（1nx-1） - u0（2nx）;
            else