C File ex2_1.ato
C Example 2-1.  Illustrates simple ODEINP file
C Y.F. Chang, ATOMFT User Manual, Version 3.11
C Partial results:
C        Step number     0 at     T =  0.00000E+00
C         Y  1.00000E+00  0.00000E+00
C    This solution will be limited by the machine roundoff.
C
C        Step number     1 at     T =  6.30000E-01
C         Y  3.00802E+00  1.03304E+01
C
C        Step number     2 at     T =  1.10000E+00
C         Y  8.77740E+01  1.64472E+03
C
C First Painleve transcendent
      DIFF(Y,T,2) = 6*Y*Y + T
   $
   $
      MSTIFF = 1
      OPEN(1,FILE='DATA',STATUS='OLD')
C Read integration interval and initial conditions.
      READ(1,110) START,END,Y(1),Y(2)
  110 FORMAT(4F10.3)
      WRITE(LIST,120) START,END,Y(1),Y(2)
  120 FORMAT(' Solve first Painleve transcendent' /
     A ' interval:          ',2F10.3 /
     A ' initial conditions:',2F10.3 /)
   $
   $