diff ( y2 , x , 1 ) = diff ( y1, x , 1) ;
diff ( y1 , x , 1 ) = sin ( x ) ;
!
Digits := 32;
max_terms := 30;
!
x_start := 0.1;
x_end := 5.0;
array_y1_init[1] := exact_soln_y1(x_start);
array_y2_init[1] := exact_soln_y2(x_start);
glob_h := 0.00001 ;
glob_max_iter := 20;
!
exact_soln_y1 := proc(x) 
  -cos(x);
end;
exact_soln_y2 := proc(x) 
 -cos(x);
end;