diff ( y , x , 1 ) = expt( 2.0 * x + 1.0 , 3.0 * x + 2.0 ) * ( 3.0 * ln( 2.0 * x + 1.0 )+ ( 2.0 * ( 3.0 * x + 2.0 ) ) / ( 2.0 * x + 1.0) ) ; ! # Comment 1 Digits := 32; max_terms := 30; # Comment 2 ! # Comment 3 x_start := 0.1; x_end := 1.0 ; diff(y,0,exact_soln_y(x_start)); glob_look_poles := true; glob_max_iter := 1000000; # Comment 4 ! # Comment 5 exact_soln_y := proc(x) # Comment 6 return(expt(2.0*x+1.0,3.0*x+2.0)); # Comment 7 end;