diff ( y , x , 3 ) = m1 * diff ( y , x , 1 ) ;
!
Digits := 32;
max_terms := 30;
!
x_start := 0.1;
x_end := 5.0 ;
array_y_init[1] := exact_soln_y(x_start);
array_y_init[2] := exact_soln_yp(x_start);
array_y_init[3] := exact_soln_ypp(x_start);
glob_h := 0.001 ;
glob_look_poles := true;
glob_max_iter := 10000;
!
exact_soln_y := proc(x)
		-cos(x);
end;
exact_soln_yp := proc(x)
		sin(x)
end;
exact_soln_ypp := proc(x)
		cos(x)
end;
exact_soln_yppp := proc(x)
		-sin(x)
end;