diff ( y , x , 1 ) =  arctan ( x ) ;
!
Digits := 32;
max_terms := 30;
!
x_start := -1.0;
x_end := 5.00 ;
diff(y,0,exact_soln_y(x_start));
glob_h := 0.00001 ;
glob_look_poles := true;
glob_max_iter := 100;
!
exact_soln_y := proc(x)
2.0 +	x * arctan(x) - log(x * x + 1.0)/2.0
end;