diff ( y2 , x , 5 ) = y1 ;
diff ( y1 , x , 1 ) = m1 * y2  ;
!
Digits := 32;
max_terms := 30;
!
x_start := 0.0;
x_end := 5.0;
array_y1_init[1] := exact_soln_y1(x_start);
array_y2_init[1] := exact_soln_y2(x_start);
array_y2_init[2] := exact_soln_y2p(x_start);
array_y2_init[3] := exact_soln_y2pp(x_start);
array_y2_init[4] := exact_soln_y2ppp(x_start);
array_y2_init[5] := exact_soln_y2pppp(x_start);
glob_h := 0.00001;
glob_look_poles := true;
glob_max_iter := 20;
!
exact_soln_y1 := proc(x)
	cos(x);
end;
exact_soln_y2 := proc(x)
	sin(x);
end;
exact_soln_y2p := proc(x)
	cos(x);
end;
exact_soln_y2pp := proc(x)
 -sin(x);
end;
exact_soln_y2ppp := proc(x)
	-cos(x);
end;
exact_soln_y2pppp := proc(x)
   sin(x);
end;