COPTION DOUBLE,DUMP=1 DIFF ( Y1 , X , 1 ) = - Y2 DIFF ( Y2 , X , 1 ) = Y1 $ $ MPRINT = 2 NSTEPS = 1000 START = 0.2D0 END = 0.8D0 DLTXPT = 0.01D0 Y1(1) = DCOS(START) Y2(1) = DSIN(START) MSTIFF = 1 ERRLIM = 1.0D-04 WRITE(LIST,120) START,END,DLTXPT,Y1(1),Y2(1) 120 FORMAT(8F16.10) $ $ C DIGITS := 64; C MAX_TERMS:=40; C ! C X_START := 0.2; C X_END := 0.8; C DIFF(Y1,0,EXACT_SOLN_Y1(X_START)); C DIFF(Y2,0,EXACT_SOLN_Y2(X_START)); C GLOB_MAX_H := 0.00001; C GLOB_LOOK_POLES := TRUE; C # GLOB_MAX_HOURS := 15; 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_Y1 := PROC(X) C RETURN(COS(X)); C END; C EXACT_SOLN_Y2 := PROC(X) C RETURN(SIN(X)); C END; C C C C