diff ( y , x , 1 ) = x * x ;
!
Digits := 32;
max_terms := 30;
!
x_start := 0.1;
x_end := 10.0 ;
array_y_init[1] := exact_soln_y(x_start);
glob_h := 0.00001 ;
glob_look_poles := true;
glob_max_iter := 20;
!
exact_soln_y := proc(x)
   x * x * x / 3.0;
end;