diff ( y , x , 1 ) =  sinh (2.0 * x + 3.0) ;
!
Digits := 32;
max_terms := 30;
!
x_start := 0.0;
x_end := 10.0 ;
diff(y,0,exact_soln_y(x_start));
glob_look_poles := true;
glob_max_iter := 10;
!
exact_soln_y := proc(x)
return(cosh(2.0*x+3.0)/2.0);
end;