diff ( y1 , x , 1 ) = m1 * diff(y2,x,0)  ;
diff ( y2 , x , 1 ) = diff(y1,x,0) ;
!
Digits := 32;
max_terms := 30;
!
x_start := 0.1;
x_end := 10.0;
array_y1_init[1] := exact_soln_y1(x_start);
array_y2_init[1] := exact_soln_y2(x_start);
glob_h := 0.00001 ;
glob_look_poles := true;
glob_max_iter :=   10;
!
exact_soln_y1 := proc(x)
cos(x);
end;
exact_soln_y2 := proc(x) 
sin(x);
end;