clc;
clear;
//coeficientes
TD = 283.15;
TE = 313.15;
TINI = 298.15;
L=1;
alpha = 40*10^-6;
celulas = 200;
dt = 1;
dx = L/celulas;
S = (alpha*dt)/(dx*dx);
B=0.1;
Tnova(1) = TE;
Tnova(celulas) = TD;
Tvelha(1)= TE;
Tvelha(celulas)= TD;
t=0;
tmax=300;
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);
            //Taux=(2*S*B*(Tnova(i-1)+Tnova(i+1))+2*S*(1-B)*(Tvelha(i-1)+Tvelha(i+1))+(1-2*S*(1-B))*Tvelha(i))/(1+2*S*B);
            //disp(Taux)
            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
//disp(barra)
//disp(Tnova);
//a=get("current_axes")
//a.x_location="origin";
//a.y_location="origin";
scf(0)
plot (barra,Tnova)
xlabel("Celula", "fontsize", 3,"color", "blue");
ylabel("Temperatura (Kelvin)", "fontsize", 3, "color", "blue");
xs2png(0,'Grafico01 dt1 200cel B0.1 300segs ')//saida do grafico em png
write('TNOVA01 dt1 200cel B0.1 300segs .txt',Tnova)//saida do vetor TNOVA