ATOMFT ver. 3.11, copyright(C) 1979-93, Y. F. Chang. ----------------------------------------------------- COPTION DOUBLE,DUMP=1 DIFF ( Y2 , X , 5 ) = Y1 DIFF ( Y1 , X , 1 ) = - DIFF (Y2 , X,4) $ C The ODE, system of ODE, entered is linear. C MSTIFF has been set to 10. $ MPRINT = 2 NSTEPS = 1000 START = 0.1d0 END = 0.7d0 DLTXPT = 0.1d0 Y1(1) = ( COS(START) ) Y2(1) = ( SIN(START) + 10.0d0 * START + 10.0d0) Y2(2) = ( COS(START) + 10.0d0) Y2(3) = ( -SIN(START)) Y2(4) = ( -COS(START)) Y2(5) = ( SIN(START)) WRITE(LIST,120) START,END,DLTXPT,Y1(1),Y2(1),Y2(2),Y2(3),Y2(4) C ,Y2(5) 120 FORMAT(8F16.10) $ $ ATOMFT completed