Fortran - llevar de matlab a fortran

clear all

close all

clc

dt=0.1

x(1)=100

y(1)=40

t(1)=1

alf=0.25

bet=-0.01

gam=-1.0

del=0.01

for n=1:1000 ;

x(n+1)=x(n)+dt*(alf*x(n)+bet*x(n)*y(n)) ;

y(n+1)=y(n)+dt*(gam*y(n)+del*x(n)*y(n)) ;

t(n+1)=t(n)+dt ;

end

figure()

hold on

plot(t,x)

plot(t,y,'r')

grid on

figure(2)

plot(x,y,'g')

%%%%%GRACIAS%%%%%

%% Si lo tienen para gfortran mejor%%%%

PROGRAM MAIN

   IMPLICIT NONE

   REAL, PARAMETER  :: dt=0.1, alf=0.25, bet=-0.01, gam=-1.0, del=0.01

   REAL             :: x, y, t

   INTEGER          :: i

     x = 100

     y = 40

     t = 1

     OPEN(UNIT=11,FILE="x_t.dat") ! CREA UN ARCHIVO DE TEXTO LLAMADO x_t.dat PARA GUARDAR LOS DATOS DEL PRIMER GRAFICO

     OPEN(UNIT=12,FILE="y_t.dat")

     OPEN(UNIT=13,FILE="x_y.dat")

     WRITE(11,*) x, t

     WRITE(12,*) y, t

     WRITE(13,*) x, y

     DO i=1,1000

        x = x + dt*(alf*x + bet*x*y)

        y = y + dt*(gam*y + del*x*y)

        t = t + dt

        WRITE(11,*) x, t

        WRITE(12,*) y, t

        WRITE(13,*) x, y

     END DO

     CALL SYSTEM ('xmgrace x_t.dat y_t.dat &') ! EN LA TERMINAL LE PIDE A GRACE QUE GRAFIQUE EL ARCHIVO x_t.dat Y y_t.dat EN LA MISMA FIGURA

     CALL SYSTEM ('xmgrace x_y.dat &')

END PROGRAM