Метод сеток для решения ДУ в частных производных

нестационарное уравнение теплопроводности:
Явная и неявная схемы
Задача Дирихле для двумерного уравнения

распределение не проводя эксперимент?

gi(i)*u(i))+tau*f(i);
end
u1(1)=be0;
u1(N1)=be1;
Plot(x,u1);
u=u1;
t=t+tau;
end;
…

i=2:N
a(i)=
b(i)=
c(i)=
d(i)= …
end
ks(1)=-c(1)/b(1); et(1)=d(1)/b(1);
for i=2:N1
z=b(i)+a(i)*ks(i-1);
ks(i)=-c(i)/z; et(i)=(d(i)-a(i)*et(i-1))/z;
end;
u1(N1)=be1;
For i=N:-1:1
u1(i)=ks(i)*u1(i+1)+et(i);
End;
Plot(x,u1);
u=u1;
t=t+tau;
end;

b=… c=… d=… e=…
up= a*u(i-1,j)+b*u(i+1,j)+c*u(i,j-1)+d*u(i,j+1)+e;
u1(i,j)=wr*up+(1-wr)*u(i.j);
end;end;

d=… e=…
up= a*u(i-1,j)+b*u(i+1,j)+c*u(i,j-1)+d*u(i,j+1)+e;
If abs(u(i,j)-up)>d then d= abs(u(i,j)-up);
u(i,j)=wr*up+(1-wr)*u(i.j);

решаем систему

Для решаем систему

end; end; (ij)
ks(1)=…; et(1)=…;
for i=2:N1 for j=2:M1
a=…b=…c=…d=…;
z=b+a*ks(j-1);
ks(j)=-c/z; et(j)=(d-a*et(j-1))/z;
end;
u(I,M1)=et(M1);
For j=M:-1:1
u(i,j)=ks(j)*u1(i,j+1)+et(j+1);
End;