COPTION DOUBLE,DUMP=1
      DIFF ( Y2 , X , 1 ) = - Y1 + 1.0 
      DIFF ( Y1 , X , 1 ) = Y2 - 1.0
   $
   $
      MPRINT = 2
      NSTEPS = 1000
      START = 0.1d0 
      END = 0.5d0
      DLTXPT = 0.01d0
      Y1(1) = (1.0d0 + DCOS(START))
      Y2(1) = (1.0d0 - DSIN(START))
      WRITE(LIST,120) START,END,DLTXPT,Y1(1),Y2(1)
  120 FORMAT(8F16.10)
   $
   $
C DIGITS := 16;
C MAX_TERMS := 40;
C !
C X_START := 0.1;
C X_END := 0.5;
C DIFF(Y1,0,EXACT_SOLN_Y1(X_START));
C DIFF(Y2,0,EXACT_SOLN_Y2(X_START));
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(1.0 +  COS(X));
C END;
C EXACT_SOLN_Y2 := PROC(X)
C RETURN(1.0 - SIN(X));
C END;
C