COPTION DOUBLE,DUMP=1
C File ex2_2.ato
C Example 2-2.  Illustrates simple ODEINP file
C Y.F. Chang, ATOMFT User Manual, Version 3.11
C Partial results:
C        Step number     0 at     T =  0.00000000D+00
C         Y  1.00000000000000D+00  0.00000000000000D+00
C   
C        Step number     1 at     T =  3.90600000D-01
C         Y  1.55090677298022D+00  3.34940599980569D+00
C   
C        Step number     2 at     T =  9.88800000D-01
C         Y  2.10504054608518D+01  1.93239147009624D+02
C   
C        Step number     3 at     T =  1.10000000D+00
C         Y  8.77740601591450D+01  1.64472417712645D+03
C
C First Painleve transcendent
      DIFF(Y,T,2) = 6*Y*Y + T
   $
   $
      MSTIFF = 1
      OPEN(1,FILE='DATA')
C Read integration interval and initial conditions.
      READ(1,*) 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 /)
   $
   $