diff ( y , x , 1 ) = expt(sin(0.2 * x + 0.3) , 2.0); ! Digits := 32; max_terms := 40; ! x_start := c(0.1); x_end := c(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); return(c(-2.5000000000000000000000000000000)*sin(c(0.20000000000000000000000000000000)*c(x)+c(0.30000000000000000000000000000000))*cos(c(0.20000000000000000000000000000000)*c(x)+c(0.30000000000000000000000000000000))+c(0.50000000000000000000000000000000)*c(x)+c(0.75000000000000000000000000000000)); end;