diff ( y , x , 1 ) = expt(sin(0.2 * x + 0.3) , 2.0); ! Digits := 32; max_terms := 40; ! x_start := 0.1; x_end := 0.2 ; diff(y,0,exact_soln_y(x_start)); glob_look_poles := true; ! exact_soln_y := proc(x) return(-2.5000000000000000000000000000000*sin(0.20000000000000000000000000000000*x+0.30000000000000000000000000000000)*cos(0.20000000000000000000000000000000*x+0.30000000000000000000000000000000)+0.50000000000000000000000000000000*x+0.75000000000000000000000000000000); end;