COPTION DOUBLE,DUMP=1 DIFF ( Y , X , 1 ) = (SIN(0.2 * X + 0.3) ** 2.0) $ $ MPRINT = 2 NSTEPS = 1000 START = 0.1d0 END = 0.2d0 DLTXPT = 0.01d0 Y(1) = (-2.5d0 C *SIN(0.2d0 C *START+0.3d0) C *COS(0.2d0 C *START+0.3d0) C +0.5d0 C *START+0.75d0) 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 := 0.2 ; C DIFF(Y,0,EXACT_SOLN_Y(X_START)); C GLOB_LOOK_POLES := TRUE; C ! C EXACT_SOLN_Y := PROC(X) C RETURN(-2.5000000000000000000000000000000*SIN(0.20000000000000000000000000000000*X+0.30000000000000000000000000000000)*COS(0.20000000000000000000000000000000*X+0.30000000000000000000000000000000)+0.50000000000000000000000000000000*X+0.75000000000000000000000000000000); C END;