diff (x1,t,1) = 4.0 * x2 - 2.0 * diff (x2,t ,1) - 2.0 * x1; diff (x2,t,2) = 3.0 * diff(x2,t,1) - 2.0 * x2 - diff(x1,t,2) - diff (x1,t,1) + x1; ! Digits := 32; max_terms:=30; ! t_start := 1.5; t_end := 1.6; diff(x1,0,exact_soln_x1(t_start)); diff(x1,1,exact_soln_x1p(t_start)); diff(x2,0,exact_soln_x2(t_start)); diff(x2,1,exact_soln_x2p(t_start)); glob_max_iter := 1000000000; glob_display_interval := 0.1; glob_max_minutes := 2; glob_desired_digits_correct := 10; glob_look_poles := true; ! exact_soln_x1 := proc(t) local c1,c2,c3; c1 := 1.0; c2 := 0.0002; c3 := 0.0003; return(2.0 * c1 + 6.0 * c3 * exp(-t)); end; exact_soln_x1p := proc(t) local c1,c2,c3; c1 := 1.0; c2 := 0.0002; c3 := 0.0003; return( - 6.0 * c3 * exp(-t)); end; exact_soln_x2 := proc(t) local c1,c2,c3; c1 := 1.0; c2 := 0.0002; c3 := 0.0003; return(c1 + c2 * exp(2.0 * t) + c3 * exp(-t)); end; exact_soln_x2p := proc(t) local c1,c2,c3; c1 := 1.0; c2 := 0.0002; c3 := 0.0003; return( 2.0 * c2 * exp(2.0 * t) - c3 * exp(-t)); end;