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