aplicaciÓn de las ecuaciones diferenciales ordinarias.docx

7/26/2019 APLICACIN DE LAS ECUACIONES DIFERENCIALES ORDINARIAS.docx

1/5

APLICACIN DE LAS ECUACIONES DIFERENCIALES ORDINARIAS (EDO)

PARA EL EJEMPLO 1:

1. MOSTRAR GRFICAMENTECDIGOPROGRAM EULER !SOLUCION DE EDO PARA UN SISTEMA DE DOS ECUACIONESCONSIDERANDO !Y(I+1)=Y(I)+H*G(Y) Y (I+1)=(I)+H*G(Y) IMPLICIT NONE REAL DIMENSION("###)::Y$ REAL::FA%CDNH INTEGER::I OPEN(UNIT=FILE='EULER.TT'STATUS='UNNON') PRINT*'METODO DE EULER EN LA SOLUCION DE SISTEMAS DE ED'

PRINT*'INGRESE LOS ALORES ETREMOS A Y % DEL TIEMPO' READ*A% PRINT*'INGRESE EL NUMERO DE INTERALOS N' READ*N PRINT*'INGRESE LAS CONDICIONES INICIALES C ' READ*CD (1)=A !PRIMER ALOR DE T Y(1)=C !1 CONDICION INICIAL $(1)=C !& CONDICION INICIALH=#.1

PRINT*'ALORES DE LA FUNCION ' RITE(#)(1)Y(1)$(1) 1# FORMAT(1F1#."F1#."F1#.) DO I=1N (I+1)=(#)+I*H Y(I+1)= Y(I)+H*F((I)Y(I)) $(I+1)=,-*EP(,(I)),&*(I)+&

RITE(*)(I+1)Y(I+1)$(I+1) END DO

STOP

END

FUNCTION F(Y) IMPLICIT NONE REAL::FY F=,&*,Y RETURN END


2/5

GRFICO

# #.1 #.& #.- #. #." #.

,1.&

,1

,#./

,#.

,#.

,#.&

#

C023 T4356

&. MODIFICAR EL ALGORITMO ANTERIOR PARA DIFERENTES 7H8

PARA EL EJEMPLO &:

1. MOSTRAR GRFICAMENTECDIGO

PROGRAM EULER !SOLUCION DE EDO PARA UN SISTEMA DE DOS ECUACIONESCONSIDERANDO !Y(I+1)=Y(I)+H*G(Y) Y (I+1)=(I)+H*G(Y) IMPLICIT NONE REAL DIMENSION("###)::Y$ REAL::FA%CDNH INTEGER::I OPEN(UNIT=FILE='EULER.TT'STATUS='UNNON') PRINT*'METODO DE EULER EN LA SOLUCION DE SISTEMAS DE ED' PRINT*'INGRESE LOS ALORES ETREMOS A Y % DEL TIEMPO' READ*A%

PRINT*'INGRESE EL NUMERO DE INTERALOS N' READ*N PRINT*'INGRESE LAS CONDICIONES INICIALES C ' READ*CD (1)=A !PRIMER ALOR DE T Y(1)=C !1 CONDICION INICIAL $(1)=C !& CONDICION INICIALH=#.1

PRINT*'ALORES DE LA FUNCION ' RITE(#)(1)Y(1)$(1)

1# FORMAT(1F1#."F1#."F1#.) DO I=1N (I+1)=(#)+I*H


3/5

Y(I+1)= Y(I)+H*F((I)Y(I)) $(I+1)=EP(,(I))+#."*(SIN((I)),COS((I)))

RITE(*)(I+1)Y(I+1)$(I+1) END DO

STOP END

FUNCTION F(Y) IMPLICIT NONE REAL::FY F=,&*,Y RETURN ENDGRFICO

# #.1 #.& #.- #. #." #. #.9 #./ #.

,#.-

,#.&

,#.1

#

#.1

#.&

#.-

#.

#."

#.

C023 T4356

&. MODIFICAR EL ALGORITMO ANTERIOR PARA DIFERENTES 7H8CDIGOPROGRAM PRIMORDIALIMPLICIT NONEREAL::A%CD

INTEGER::NCHARACTER(1)::OP

1 PRINT*'CALCULO DE LAS RAICES DE LA FUNCION'!APLICACION DE LOS METODOSPRINT*',,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,'PRINT*'(1) APLICAR METODO DE EULER'

PRINT*'(#) PARA SALIR DEL PROGRAMA'PRINT*',,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,'& PRINT*'SELECCIONE UNA OPCION ;#,1


4/5

CASE('#')CASE('1') 55 EULER(>?@)

PRINT*'EL PROGRAMA HA FINALI$ADO'CASE DEFAULTGOTO &END SELECTGOTO 1END PROGRAM PRIMORDIAL

SU%ROUTINE EULER(A%CDN) !SOLUCION DE EDO PARA UN SISTEMA DE DOS ECUACIONESCONSIDERANDO !Y(I+1)=Y(I)+H*G(Y) Y (I+1)=(I)+H*G(Y) REAL DIMENSION("###)::Y$ REAL::FA%CDH

INTEGER::IN OPEN(UNIT=FILE='EULER [email protected]'STATUS='UNNON') PRINT*'METODO DE EULER EN LA SOLUCION DE SISTEMAS DE ED' PRINT*'INGRESE LOS ALORES ETREMOS A Y % DEL TIEMPO' READ*A% PRINT*'INGRESE EL NUMERO DE INTERALOS N' READ*N PRINT*'INGRESE LAS CONDICIONES INICIALES C ' READ*CD (1)=A !PRIMER ALOR DE T Y(1)=C !1 CONDICION INICIAL $(1)=C !& CONDICION INICIALH=(%,A)BN

PRINT*'ALORES DE LA FUNCION ' RITE(#)(1)Y(1)$(1) 1# FORMAT(1F1#."F1#."F1#.) DO I=1N (I+1)=(#)+I*H Y(I+1)= Y(I)+H*F((I)Y(I)) $(I+1)=EP(,(I))+#."*(SIN((I)),COS((I))) PRINT*(I+1)Y(I+1)$(I+1) RITE(#)(I+1)Y(I+1)$(I+1)

END DO

RETURN END

FUNCTION F(Y) REAL::FY F=,&*,Y RETURN END

GRFICO


5/5

# #." 1 1." & &." - -." ."

,9

,

,"

,

,-

,&

,1#

1

&

ECUACIN DE EULER

ANALITICO

UNODOS

TRES

CUATRO

CINCO

SEIS

ITERACIONES

F (ITERACIONES)

PARA EL EJEMPLO -:

1. MOSTRAR GRFICAMENTE&. MODIFICAR EL ALGORITMO ANTERIOR PARA DIFERENTES 7H8

PARA EL EJEMPLO :


PARA EL EJEMPLO ":


aplicaciÓn de las ecuaciones diferenciales ordinarias.docx

Documents