ATOMFT ver. 3.11, copyright(C) 1979-93, Y. F. Chang. ----------------------------------------------------- COPTION DOUBLE,DUMP=1 DIFF ( Y , X , 3 ) = SIN(X) $ $ MPRINT = 2 NSTEPS = 1000 START = 0.1d0 END = 1.5d0 DLTXPT = 0.1d0 Y(1) = (1.0d0 + DCOS(START)) Y(2) = (-DSIN(START)) Y(3) = (-DCOS(START)) WRITE(LIST,120) START,END,DLTXPT,Y(1),Y(2),Y(3) 120 FORMAT(8F16.10) $ $ ATOMFT completed