diff(y2,x,1) = y1 - 2.0; diff(y1,x,1) = diff(y2,x,5); ! Digits := 32; max_terms := 30; ! x_start := 0.5; x_end := 10.0; diff(y1,0,exact_soln_y1(x_start)); diff(y2,0,exact_soln_y2(x_start)); diff(y2,1,exact_soln_y2p(x_start)); diff(y2,2,exact_soln_y2pp(x_start)); diff(y2,3,exact_soln_y2ppp(x_start)); diff(y2,4,exact_soln_y2pppp(x_start)); glob_h := 0.00001 ; glob_look_poles := true; glob_max_iter := 10; glob_subiter_method := 3; ! exact_soln_y1 := proc(x) 2.0 + sin(x); end; exact_soln_y2 := proc(x) 2.0 - cos(x); end; exact_soln_y2p := proc(x) sin(x); end; exact_soln_y2pp := proc(x) cos(x); end; exact_soln_y2ppp := proc(x) -sin(x); end; exact_soln_y2pppp := proc(x) -cos(x); end;