diff ( y , x , 1 ) = arcsin(sqrt(0.1 * x + 0.2));
!
Digits := 32;
max_terms := 30;
!
x_start := 0.0;
x_end := 0.5 ;
diff(y,0,exact_soln_y(x_start));
glob_look_poles := true;
glob_max_iter := 1000000;
glob_display_interval := 0.1;
glob_max_minutes := 10;
glob_max_h := 0.0001;
#
# Not Given = 0
# Real = 1
# Complex = 2 
# No Pole = 3
# Impossible Eq = 4
#
glob_type_given_pole := 3;
#  Real Part 
array_given_rad_poles[1,1] := 0.0;  
# Imag Part 
array_given_rad_poles[1,2] := 1.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(10.0 * (0.1 * x + 0.2) * arcsin(sqrt ( 0.1 * x + 0.2)) + 5.0 * sqrt( 0.1 * x +
0.2) * sqrt( 0.8 - 0.1 * x) - 5.0 * arcsin(sqrt( 0.1 * x + 0.2)));
end;