function u = peEllip5(nx,minx,maxx,ny,miny,maxy)
format long;
hx = (maxx-minx)/(nx-1);
hy = (maxy-miny)/(ny-1);
u0 = zeros(nx,ny);
for j=1:ny
u0(j,1) = EllIni2Uxl(minx,miny+(j-1)*hy);
u0(j,nx) = EllIni2Uxr(maxx,miny+(j-1)*hy);
end
for j=1:nx
u0(1,j) = EllIni2Uyl(minx+(j-1)*hx,miny);
u0(ny,j) = EllIni2Uyr(minx+(j-1)*hx,maxy);
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(i,nx-2) == 1
if i==1
A(1,2) = 1;
A(1,nx-1) = 1;
b(1) = – u0(1,2) – u0(2,1);
else
if i == (ny-3)*(nx-2)+1
A(i,i+1) = 1;
A(i,i-nx+2) = 1;
b(i) = – u0(ny-1,1) – u0(ny,2);
else</