c Runge Kutta for first order differential equations c PROGRAM Euler IMPLICIT none c c declarations c nsteps:number of steps, tstep:length of steps, y: initial position c REAL*8 t, y, tstep INTEGER j, nsteps nsteps=10 tstep=0.5 y=1.0 c c open file OPEN(6, FILE='eulerf.dat') WRITE (6,*) 0, y c c do loop nsteps of Euler algorithm DO 60 j = 1, nsteps t=j*tstep y=y-tstep*y WRITE (6,*) t, y 60 CONTINUE c CLOSE(6) STOP END