program numerical1

      real x,t,dt,f
      open(unit=13,file="data.d")
      dt  = 0.25

      x = 1

      do 10 t=0,5,dt
         write(13,*) t,x,exp(-t)
         f = -x
         x = x + f*dt
 10   continue

      end