COPTION DOUBLE,DUMP=1
      DIFF ( Y , X , 1 ) = - 2.5 / (( X - 6.0 ) ** 3.5) 
   $
   $
      MPRINT = 2
      NSTEPS = 1000
      START = 3.0d0
      END = 2.0d0
      DLTXPT = 0.1d0
      Y(1) = ((START - 6.0) ** 2.5 )
      WRITE(LIST,120) START,END,DLTXPT,Y(1)
  120 FORMAT(8F16.10)
   $
   $
C DIGITS := 32;
C MAX_TERMS := 20;
C !
C X_START := 3.0;
C X_END := 2.0;
C DIFF(Y,0,EXACT_SOLN_Y(X_START));
C GLOB_LOOK_POLES := TRUE;
C # NOT GIVEN = 0
C # NO POLE = 3
C # POLE = 4
C GLOB_TYPE_GIVEN_POLE := 4;
C #  REAL PART 
C ARRAY_GIVEN_RAD_POLES[1,1] := 6.0;  
C # IMAG PART 
C ARRAY_GIVEN_RAD_POLES[1,2] := 0.0;
C # ORDER 
C ARRAY_GIVEN_ORD_POLES[1,1] := 3.5; 
C # NOT USED
C ARRAY_GIVEN_ORD_POLES[1,2] := 0.0;  
C !
C EXACT_SOLN_Y := PROC(X)
C RETURN( EXPT(X - 6.0,2.5 ));
C END;
C 
C