diff ( y , x , 1 ) = m1 * 2.0 / x / x / x ;
!
Digits := 32;
max_terms := 20;
!
x_start := -1.0;
x_end := -0.7;
diff(y,0,exact_soln_y(x_start));
glob_look_poles := true;
glob_max_iter := 100;
!
exact_soln_y := proc(x)
return(1.0/x/x);
end;