P. 733, line 19

From Bill Goodwine's Wiki
Jump to: navigation, search

The program should be:

      program systemeuler

c     This is a sample FORTRAN program that solves the differential
c     equation
c
c     x' = y, y' = (1 - x^2)y - x
c     x(0) = 0.02, y(0) = 0.0
c     
c     using Euler's method.

      double precision x(2),t,dt
      double precision copy(2)

      open(unit=13,file="systemfortran.d")

      dt = 0.001
      x(1) = 0.02
      x(2) = 0.0

      do 10 t = 0, 20, dt
         write(13,*) t,x(1),x(2)
         copy(1) = x(1)
         copy(2) = x(2)
         x(1) = x(1) + (copy(2))*dt
         x(2) = x(2) + ((1.0 - copy(1)**2)*copy(2) - copy(1))*dt
 10   continue
      stop
      end


Return to errata.