COPTION DOUBLE,DUMP=1 DIFF ( Y , X , 1 ) = SIN ( X ) $ $ MPRINT = 2 NSTEPS = 1000 START = -5.0d0 END = 5.0d0 DLTXPT = 0.1d0 Y(1) = (1.0d0 - DCOS(START)) WRITE(LIST,120) START,END,DLTXPT,Y(1) 120 FORMAT(8F16.10) $ $ C DIGITS := 32; C MAX_TERMS := 30; C ! C X_START := -5.0; C X_END := 5.0 ; C DIFF(Y,0,EXACT_SOLN_Y(X_START)); C GLOB_LOOK_POLES := TRUE; C GLOB_MAX_ITER := 10000; C GLOB_DISPLAY_INTERVAL := 0.1; C GLOB_MAX_MINUTES := 10; C ! C EXACT_SOLN_Y := PROC(X) C RETURN(1.0 - COS(X)); C END;