clc;
clear;
//coeficientes
TD = 283.15;
TE = 313.15;
TINI = 298.15;
L=1;
alpha = 80*10^-6;
celulas = 200;
dt = 1;
dx = L/celulas;
S = (alpha*dt)/(dx*dx);
B=0.7;
Tnova(1) = TE;
Tnova(celulas) = TD;
Tvelha(1)= TE;
Tvelha(celulas)= TD;
t=0;
tmax=3600;
totalit=0;
itmax=10000
tolerancia=0.00001;
for i=2:(celulas-1)
    Tnova(i)=TINI;
    Tvelha(i)=TINI;
end;
while(t<tmax)
    while (totalit<itmax)
        erro=-1;
        for i=2:(celulas-1)
            Taux=(Tvelha(i)+S*(1-B)*(Tvelha(i-1)-2*Tvelha(i)+Tvelha(i+1))+S*B*Tnova(i-1)+S*B*Tnova(i+1))/(1+2*S*B);
            erro=max(erro,abs(Taux-Tnova(i)));
            Tnova(i) = Taux;
        end;
        if (erro<tolerancia) then
            break;
        end;
    totalit=totalit+1;
end;
   Tvelha= Tnova;
   t=t+dt;
end;
for i=1:celulas
    barra(i)=i;
end
scf(0)
plot (barra,Tnova)
xlabel("Celula", "fontsize", 3,"color", "blue");
ylabel("Temperatura (Kelvin)", "fontsize", 3, "color", "blue");
xs2png(0,'Grafico01 dt1 200cel B0.7 3600segs alpha80')//saida do grafico em png
write('TNOVA01 dt1 200cel B0.7 3600segs alpha80.txt',Tnova)//saida do vetor TNOVA