COPTION DOUBLE,DUMP=1
      DIFF ( Y , X , 5 ) = - DIFF ( Y , X , 3) 
   $
   $
      MPRINT = 2
      NSTEPS = 1000
      START = 0.1d0
      END = 1.5d0
      DLTXPT = 0.1d0
      Y(1) = (DCOS(START))
      Y(2) = (-DSIN(START))
      Y(3) = (-DCOS(START))
      Y(4) = (DSIN(START))
      Y(5) = (DCOS(START))
      WRITE(LIST,120) START,END,DLTXPT,Y(1),Y(2),Y(3),Y(4),Y(5)
  120 FORMAT(8F16.10)
   $
   $
C DIGITS := 64;
C MAX_TERMS := 40;
C !
C X_START := 0.1;
C X_END := 1.5 ;
C #
C # TROUBLE ABOUT PI/2
C #
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 DIFF(Y,3,EXACT_SOLN_YPPP(X_START));
C DIFF(Y,4,EXACT_SOLN_YPPPP(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;
C EXACT_SOLN_YPPP := PROC(X)
C RETURN(SIN(X));
C END;
C EXACT_SOLN_YPPPP := PROC(X)
C RETURN(COS(X));
C END;