COPTION DOUBLE,DUMP=1 DIFF ( Y , X , 1 ) = COSH(SQRT(0.1 * X + 0.2)) $ $ MPRINT = 2 NSTEPS = 1000 START = -1.9D0 END = -1.8D0 DLTXPT = 0.01D0 Y(1) = (20.0D0 * DSQRT(0.1D0 * START + 0.2D0) C * DSINH( DSQRT(0.1D0 * START + 0.2D0)) - 20.0D0 C * DCOSH( DSQRT(0.1D0 * START + 0.2D0))); WRITE(LIST,120) START,END,DLTXPT,Y(1) 120 FORMAT(8F16.10) $ $ C DIGITS := 32; C MAX_TERMS := 30; C ! C X_START := -1.9; C X_END := -1.8 ; C DIFF(Y,0,EXACT_SOLN_Y(X_START)); C GLOB_LOOK_POLES := TRUE; C GLOB_MAX_H := 0.01; C GLOB_UPPER_RATIO_LIMIT := 1.001; C GLOB_LOWER_RATIO_LIMIT := 0.999; C # GLOB_MIN_H := 0.000001; 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 := 1; C # REAL PART C ARRAY_GIVEN_RAD_POLES[1,1] := -2.0; C # IMAG PART C ARRAY_GIVEN_RAD_POLES[1,2] := 0.0; C # ORDER C ARRAY_GIVEN_ORD_POLES[1,1] := 0.5; C # NOT USED C ARRAY_GIVEN_ORD_POLES[1,2] := 0.0; C ! C EXACT_SOLN_Y := PROC(X) C RETURN(20.0 * SQRT(0.1 * X + 0.2) * SINH( SQRT(0.1 * X + 0.2)) - 20.0 * COSH( SQRT(0.1 * X + 0.2))); C END;