townini 12.3

clear all;
close all;

k=0.0001; % krok czasowy
h=1; % krok odległościowy
D=1000; % współczynnik dyfuzji dla aluminium
lambda=D*k/(h^2);
N=40;
M=10000;
T=zeros(N,M);
T(1,:)=200;

for i=2:N-1
for j=1:M
T(i,j+1)=T(i,j)+lambda*(T(i+1,j)-2*T(i,j)+T(i-1,j));
T(N,j+1)=T(N,j)+lambda*(2*T(N-1,j)-2*T(N,j));
end
end
surf(T);