diff( y , x , 1 ) = cos ( x ) / ( 2.0 * x + 1.0 ) - 2.0 * sin ( x ) / ( 2.0 * x + 1.0 ) / ( 2.0 * x + 1.0 ) ;
!
Digits := 32;
max_terms := 30;
!
x_start := 0.1;
x_end := 1.0 ;
diff(y,0,exact_soln_y(x_start));
glob_look_poles := true;
#
# Not Given = 0
# Real = 1
# Complex = 2 
# No Pole = 3
# Impossible Eq = 4
# Some Pole = 5
#
glob_type_given_pole := 1;
#  Real Part 
array_given_rad_poles[1,1] := -0.5;  
# Imag Part 
array_given_rad_poles[1,2] := 0.0;
# Order 
array_given_ord_poles[1,1] := 1.0; 
# Not Used
array_given_ord_poles[1,2] := 0.0;  

!
exact_soln_y := proc(x)
return(sin(x)/(2*x+1));
end;