Fortran - De Matlab a Fortran

clear;

clc;

display('Proyecto Primer Bimestre Luis Jimenez, Diego Bustamante');

t = input('Ingrese la malla temporal: ');

x = input('Ingrese la malla espacial: ');

tm = input('Ingrese el tiempo máximo: ');

L = 0.1;

dx = L / x;

dt = tm / t;

k = 50;

Ti = 300;

Tinf = 400;

h = 200;

alfa = 15e-6;

Fo =  (alfa * dt) / (dx ^ 2);

 if Fo > 0.5

    display('El Fo calculado tiene un valor superior a 0.5, se procederá a calcular el tiempo máximo.');

    dt = (0.5 * (dx ^ 2)) / alfa;

    display('El tiempo máximo es: ');

    display(dt * t);

    Fo =  (alfa * dt) / (dx ^ 2);

 end;

Bi = h * dx / k;

%Inicializamos los valores de la matriz.

Ti1 = zeros(x + 1, t);

Ti1(:, 1) = Ti * ones(x + 1, 1);

for j=1:t

    for i=1:x+1

        if i==1

            %Ecuación 2

            Ti1(1,j+1)=2*Fo*(Ti1(2,j)+(Bi*Tinf))+((1-(2*Fo)-(2*Bi*Fo))*Ti1(1,j));

        end

        if i==x+1

            %Ecuación 3

            Ti1(i,j+1)=(2*Fo*Ti1(i-1,j))+((1-(2*Fo))*Ti1(i,j));

        end

        if i > 1 && i <= x

            %Ecuación 1

            Ti1(i,j+1)=Fo*(Ti1(i-1,j) +Ti1(i+1,j))+(1-2*Fo)*Ti1(i,j);

        end

    end

end

eje_x = (0:dx:L);

t_1 = (0:dt:tm)

figure(1)

plot(eje_x, Ti1(:, 2), eje_x, Ti1(:, round(t / 10)), eje_x, Ti1(:, round(t/5)), eje_x, Ti1(:, round(t/2)), eje_x, Ti1(:, t))

xlabel('Malla Espacial')

ylabel('Temperatura')

title('temperatura vs tiempo')

figure(2)

plot(t_1, Ti1(1, :), t_1, Ti1(x / 2, :),t_1, Ti1(x, :))

xlabel('Malla Espacial')

ylabel('Temperatura')

title('temperatura vs tiempo')

figure(3)

contourf(Ti1(1:2, :));

ylabel('Malla Espacial')

xlabel('Temperatura')

figure(4)

contourf(Ti1(x/2:x/2+1,:))

ylabel('Malla Espacial')

xlabel('Temperatura')

figure(5)

contourf(Ti1(x:x+1,:))

ylabel('Malla Espacial')

xlabel('Temperatura')