Placa plana en estado estable

Placa plana en estado estable
Resuelva el problema de la placa plana sin generación de calor con condiciones que se indican en la figura utilizando 10, 20, y 50 nodos distribuidos uniformemente enlas direcciones coordenadas. Estado estable.
Utilice un ε=0.1X10-5.
Haga sus cálculos para un factor de relajación λ=1 y compare con los de λ=1.2 y λ=1.5 graficando el error máximo en cadaiteración contra el numero de iteraciones.
Lx=Ly=10 cm
Código de MATLAB
clear;
clc;
N=input('\nIndique el número de nodos para el análisis N: ');
lamda=input('\nIndique el factor de relajacion: ');Tlxcero=75;
Tlxfinal=50;
Tlyfinal=100;
Lx=10;
n=N-1;
deltax=Lx/(n);
e=.00001;
error=1;
for i=1:N;
for j=1:N;
if i==1
T(j,i)=75;
elseif i==N
T(j,i)=50;
elseif j==N
T(j,i)=100;
else
T(j,i)=0;
end
end
endTold=T;
Tuno=T;
parae=0;
while error>e
parae=parae+1;
%barrido en X
for X=2:1:n;
for j=1:1:n;
for i=1:1:n;
if j==i
ax(j,i)=4;
elseif i==j+1
ax(j,i)=-1;
elseif i==j-1
ax(j,i)=-1;else
ax(j,i)=0;
end
end
end
ax(1,2)=-2;
for j=1:1:n;
if j==1
b(j,1)=T(j,X-1)+T(j,X+1);
elseif j==n
b(j,1)=T(j,X-1)+T(j,X+1)+T(j+1,X);
else
b(j,1)=T(j,X-1)+T(j,X+1);
end
endfor r=2:1:n;
ax(r,r-1)= ax(r,r-1)/ax(r-1,r-1);
ax(r,r)= ax(r,r)-(ax(r,r-1)*ax(r-1,r));
b(r,1)= b(r,1)-(ax(r,r-1)*b(r-1,1));
end
if b(n,1)==0;
C(n,1)=0;
else
C(n,1)=b(n,1)/ax(n,n);
endT(n,X)=C(n,1);
for j=1:1:n-1;
nume=b(n-j,1)-(C(n-j+1,1)*ax(n-j,n-j+1));
deno=ax(n-j,n-j);
C(n-j,1)=nume/deno;
T(n-j,X)=C(n-j,1);
end
end
%barrido en Y
for Y=1:1:n;
for j=1:1:n-1;
for i=1:n-1;if j==i
a(j,i)=4;
elseif i==j+1
a(j,i)=-1;
elseif i==j-1
a(j,i)=-1;
else
a(j,i)=0;
end
end
end
if Y==1
for j=1:1:n-1;
if j==1
b(j,1)=T(Y,j)+(2*T(Y+1,j+1));
elseif j==n-1b(j,1)=T(Y,j+2)+(2*T(Y+1,j+1));
else
b(j,1)=2*T(Y+1,j+1);
end
end
else
for j=1:1:n-1;
if j==1
b(j,1)=T(Y,j)+T(Y+1,j+1)+T(Y-1,j+1);
elseif j==n-1
b(j,1)=T(Y,j+2)+T(Y+1,j+1)+T(Y-1,j+1);
else...