COPTION DOUBLE,DUMP=1 DIFF ( Y , X , 1 ) = SIN(0.3 * X + 0.1) - (0.1 * X + 0.2) $ $ MPRINT = 2 NSTEPS = 1000 START = 0.1d0 END = 5.0d0 DLTXPT = 0.1d0 Y(1) = (-0.05d0 * START * START - 0.2d0 * START C - DCOS(0.3d0 * START + 0.1d0) / 0.3d0) WRITE(LIST,120) START,END,DLTXPT,Y(1) 120 FORMAT(8F16.10) $ $ C DIGITS := 32; C MAX_TERMS := 40; C ! C X_START := 0.1; C X_END := 5.0 ; C DIFF(Y,0,EXACT_SOLN_Y(X_START)); C GLOB_LOOK_POLES := TRUE; C # C # NOT GIVEN = 0 C # REAL = 1 C # COMPLEX = 2 C # NO POLE = 3 C # IMPOSSIBLE EQ = 4 C # C GLOB_TYPE_GIVEN_POLE := 3; C C ! C EXACT_SOLN_Y := PROC(X) C RETURN(-0.05 * X * X - 0.2 * X - COS(0.3 * X + 0.1) / 0.3); C END;