diff ( y , x , 1 ) = sinh(sqrt(0.1 * x + 0.2));
!
Digits := 32;
max_terms := 40;
!
x_start := 2.0;
x_end := 3.0 ;
diff(y,0,exact_soln_y(x_start));
glob_look_poles := true;
glob_max_h := 0.0000001;
#
# Not Given = 0
# Real = 1
# Complex = 2 
# No Pole = 3
# Impossible Eq = 4
#
glob_type_given_pole := 1;
#  Real Part 
array_given_rad_poles[1,1] := -2.0;  
# Imag Part 
array_given_rad_poles[1,2] := 0.0;
# Order 
array_given_ord_poles[1,1] := 0.5; 
# Not Used
array_given_ord_poles[1,2] := 0.0;  

!
exact_soln_y := proc(x)
return(20.0 * sqrt(0.1 * x + 0.2) * cosh( sqrt(0.1 * x + 0.2)) - 20.0 * sinh( sqrt(0.1 * x + 0.2))); 
end;