COPTION DOUBLE,DUMP=1 DIFF ( Y , X , 3 ) = -1.0 * DIFF ( Y , X , 1 ) $ $ MPRINT = 2 NSTEPS = 1000 START = -1.0d0 END = 1.0d0 DLTXPT = 0.1d0 Y(1) = ( - DCOS(START)) Y(2) = ( DSIN(START)) Y(3) = ( DCOS(START)) WRITE(LIST,120) START,END,DLTXPT,X(1),Y(1),Y(2),Y(3) 120 FORMAT(8F16.10) $ $ C DIGITS := 32; C MAX_TERMS := 40; C ! C X_START := -1.0; C X_END := 1.0 ; C DIFF(Y,0,EXACT_SOLN_Y(X_START)); C DIFF(Y,1,EXACT_SOLN_YP(X_START)); C DIFF(Y,2,EXACT_SOLN_YPP(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 EXACT_SOLN_Y := PROC(X) C RETURN( - COS(X)); C END; C EXACT_SOLN_YP := PROC(X) C RETURN( SIN(X)); C END; C EXACT_SOLN_YPP := PROC(X) C RETURN( COS(X)); C END;