【1階微分　オイラー法にて】


****** main program ******
      integer i
      real*8 t,y,tini,tfin,dt
      real*8 f
      
      f(y,t)=t
      write(*,*)'dt='
      read(*,*)dt
      tini=0d0
      write(*,*)'t='
      read(*,*)tfin
c     yini=1d0
      write(*,*)'  t               ','          y  '

******* subroutine *******
c start euler method      
         t=0.0
         y=1.0
      do 20 i=1,10000
         t=tini+dt*dble(i)
         y=y+f(y,t)*dt
         write(*,*)t,y
      if(t.ge.tfin) stop
   20 continue
      end