¡Descarga metodos computacionales y más Apuntes en PDF de Matemáticas solo en Docsity!
PROGRAM EULER_sistemaEDO IMPLICIT NONE REAL, DIMENSION (1000)::T,S,Y REAL:: H,F1,F2,A,B INTEGER::N,I OPEN(20, FILE='enfermedad.TXT',STATUS='UNKNOWN',ACTION='WRITE') WRITE(,10,ADVANCE='NO')'INGRESE LIMITE INFERIOR a:' READ,A WRITE(,10,ADVANCE='NO')'INGRESE LIMITE SUPERIOR b:' READ,B WRITE(,10,ADVANCE='NO')'NUMERO DE INTERACIONES N:' READ,N 10 FORMAT(1X,A) !PASO 2 T(0)=A !t= S(0)=700 !CONDICION INICIAL S(0)= Y(0)=1 !CONDICION INICAL Y(0)= H=(B-A)/FLOAT(N) !PASO 3 !WRITE(,200) T(0),S(0),Y(0) ! PASO 4 DO I=0,N- T(I+1)=T(I)+H S(I+1)=S(I)+HF1(T(I),S(I),Y(I)) Y(I+1)=Y(I)+HF2(T(I),S(I),Y(I)) WRITE(20,200) T(I),S(I),Y(I) !PASO 6 END DO 200 FORMAT(3X,3(2X,F20.8)) END PROGRAM EULER_sistemaEDO !************************************ FUNCTION F1(T,S,Y) REAL::T,S,Y,F A=0. F1=-ASY+0T! LA DERIVADA RETURN END FUNCTION F2(T,S,Y) REAL::T,S,Y,F A=0. B=0. F2=ASY-BY+0T! LA DERIVADA RETURN END
clc clear all A=load('enfermedad.TXT'); t=A(:,1); S=A(:,2); Y=A(:,3); subplot(2,2,1); plot(t,S,'r.-') xlabel('tiempo (dias)') ylabel('Suceptible a enfermedad') subplot(2,2,2); plot(t,Y,'k.-') xlabel('tiempo (dias)') ylabel('con enfermedad') subplot(2,2,3); plot(S,Y,'m.-') xlabel('Suceptible a enfermedad') ylabel('con enfermedad') subplot(2,2,4); plot(t,S,t,Y) xlabel('tiempo (dias)') ylabel('S Y Y')
WRITE(20,200) T(I),Y(I)
!PASO 6
END DO
200 FORMAT(3X,3(2X,F20.8))
END PROGRAM RK
FUNCTION F(T,Y)
REAL::T,Y,F
A=0.
F=10EXP(-(T-2)2/(2(0.075*2)))-0.6Y! LA DERIVADA
RETURN
END
PROGRAM RK
IMPLICIT NONE
REAL, DIMENSION (1000)::T,Y
REAL:: H,F,A,B,K1,K2,K3,K
INTEGER::N,I
OPEN(20, FILE='enfermedad1.TXT',STATUS='UNKNOWN',ACTION='WRITE') WRITE(,10,ADVANCE='NO')'INGRESE LIMITE INFERIOR a:' READ,A WRITE(,10,ADVANCE='NO')'INGRESE LIMITE SUPERIOR b:' READ,B WRITE(,10,ADVANCE='NO')'NUMERO DE INTERACIONES N:' READ,N 10 FORMAT(1X,A) !PASO 2 T(0)=A !t= Y(0)=0.5 !CONDICION INICAL Y(0)=
H=(B-A)/FLOAT(N)
!PASO 3
!WRITE(*,200) T(0),S(0),Y(0)
! PASO 4
DO I=0,N-
T(I+1)=T(I)+H
K1=H*F(T(I),Y(I))
K2=HF(T(I)+0.5H,Y(I)+0.5*K1)
K3=HF(T(I)+0.5H,Y(I)+0.5*K2)
K4=H*F(T(I)+H,Y(I)+K3)
Y(I+1)=Y(I)+((K1+2*(K2+K3)+K4)/6.0)
WRITE(20,200) T(I),Y(I)
!PASO 6
END DO
200 FORMAT(3X,3(2X,F20.8))
END PROGRAM RK
FUNCTION F(T,Y)
REAL::T,Y,F
A=0.
F=10EXP(-(T-2)2/(2(0.075*2)))-0.6Y! LA DERIVADA
RETURN
END
