diff ( y , x , 1 ) = m1 * 5.0 / ( x - 6.0 ) / ( x - 6.0 ) / ( x - 6.0) / ( x- 6.0) / ( x - 6.0) / ( x - 6.0) ;
!
Digits := 64;
max_terms := 30;
!
x_start := 2.0;
x_end := 3.0;
diff(y,0,exact_soln_y(x_start));
glob_look_poles := true;
glob_max_iter := 100000;
glob_max_h := 0.0001;
# Not Given = 0
# No Pole = 3
# Pole = 4
glob_type_given_pole := 4;
array_given_rad_poles[1][1] = 6.0;
array_given_rad_poles[1][2] = 0.0;
array_given_ord_poles[1][1] = 6.0;
array_given_ord_poles[1][2] = 0.0;
#  Real Part 
array_given_rad_poles[1,1] := 6.0;  
# Imag Part 
array_given_rad_poles[1,2] := 0.0;
# Order 
array_given_ord_poles[1,1] := 6.0; 
# Not Used
array_given_ord_poles[1,2] := 0.0;  
!
exact_soln_y := proc(x)
return(1.0/ ( x - 6.0 ) / ( x - 6.0 )/ ( x - 6.0 )/ ( x - 6.0 )/ ( x - 6.0 )) ;
end;