COPTION DOUBLE,DUMP=1
      DIFF ( Y2 , X , 5 ) = Y1 
      DIFF ( Y1 , X , 1 ) = - DIFF (Y2 , X,4) 
   $
   $
      MPRINT = 2
      NSTEPS = 1000
      START = 0.1d0
      END = 0.7d0
      DLTXPT = 0.1d0
      Y1(1) = (	 COS(START) )
      Y2(1) = (	SIN(START) + 10.0d0 * START + 10.0d0)
      Y2(2) = (	COS(START) + 10.0d0)
      Y2(3) = ( -SIN(START))
      Y2(4) = (	-COS(START))
      Y2(5) = (   SIN(START))
      WRITE(LIST,120) START,END,DLTXPT,Y1(1),Y2(1),Y2(2),Y2(3),Y2(4)
     C ,Y2(5)
  120 FORMAT(8F16.10)
   $
   $
C DIGITS := 64;
C MAX_TERMS := 40;
C !
C X_START := 0.1;
C X_END := 0.7;
C #
C # TROUBLE ABOUT PI/4???
C #
C DIFF(Y1,0,EXACT_SOLN_Y1(X_START));
C DIFF(Y2,0,EXACT_SOLN_Y2(X_START));
C DIFF(Y2,1,EXACT_SOLN_Y2P(X_START));
C DIFF(Y2,2,EXACT_SOLN_Y2PP(X_START));
C DIFF(Y2,3,EXACT_SOLN_Y2PPP(X_START));
C DIFF(Y2,4,EXACT_SOLN_Y2PPPP(X_START));
C GLOB_LOOK_POLES := TRUE;
C GLOB_MAX_H := 0.1;
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) + 10.0 * X + 10.0);
C END;
C EXACT_SOLN_Y2P := PROC(X)
C RETURN(	COS(X) + 10.0);
C END;
C EXACT_SOLN_Y2PP := PROC(X)
C RETURN( -SIN(X));
C END;
C EXACT_SOLN_Y2PPP := PROC(X)
C RETURN(	-COS(X));
C END;
C EXACT_SOLN_Y2PPPP := PROC(X)
C RETURN(   SIN(X));
C END;
C