COPTION DOUBLE,DUMP=1 DIFF ( Y , X , 1 ) = - 5.0 / ( X - 6.0 ) C / ( X - 6.0 ) / ( X - 6.0) / ( X- 6.0) C / ( X - 6.0) / ( X - 6.0) $ $ MPRINT = 2 NSTEPS = 1000 START = 3.0d0 END = 2.0d0 DLTXPT = 0.1d0 Y(1) = (1.0d0/ ( START - 6.0d0 ) / ( START - 6.0d0 ) C / ( START - 6.0d0 ) C / ( START - 6.0d0 )/ ( START - 6.0d0 )) WRITE(LIST,120) START,END,DLTXPT,Y(1) 120 FORMAT(8F16.10) $ $ C DIGITS := 64; C MAX_TERMS := 30; C ! C X_START := 3.0; C X_END := 2.0; C DIFF(Y,0,EXACT_SOLN_Y(X_START)); C GLOB_LOOK_POLES := TRUE; 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 GLOB_TYPE_GIVEN_POLE := 1; C # REAL PART C ARRAY_GIVEN_RAD_POLES[1,1] := 6.0; C # IMAG PART C #ARRAY_GIVEN_RAD_POLES[1,2] := 0.0; C # ORDER C #ARRAY_GIVEN_ORD_POLES[1,1] := 6.0; C # NOT USED C #ARRAY_GIVEN_ORD_POLES[1,2] := 0.0; C ! C EXACT_SOLN_Y := PROC(X) C RETURN(1.0/ ( X - 6.0 ) / ( X - 6.0 )/ ( X - 6.0 )/ ( X - 6.0 )/ ( X - 6.0 )) ; C END; C C