program Fg_01_04;

uses
   Utils;

var
   Lambda,Omega,B: Double;

begin
   Compile_Full_Dataset(325,1917,True,'North');

   Compute_Method(0,5,1,12.0,16.0,0.001,0.0,0.0,0.0,'North',mtHits,True,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG01');

   Compute_Method(0,5,1,13.5,15.0,0.001,0.0,0.0,0.0,'North',mtHits,True,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG03a');
   
   Compute_Method(0,5,1,13.5,15.0,0.001,0.0,6.0,0.1,'North',mtHits,True,True,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG03b');
   Compute_Method(0,5,1,Omega-0.005,Omega+0.005,0.0001,B-0.2,B+0.2,0.002,'North',mtHits,True,False,Lambda,Omega,B);
   Write_Results(Lambda,Omega,B,'North','');

   Write_Shifts(Lambda,Omega,B,'North','FG04',True);

   Compile_Strip(325,1917,True,0.0,10.0,'North_0_10');
   Compile_Strip(325,1917,True,15.0,25.0,'North_15_25');
   Compile_Strip(325,1917,True,30.0,90.0,'North_30_90');

   Compute_Method(0,5,1,12.0,16.0,0.001,0.0,0.0,0.0,'North_0_10',mtHits,True,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG02a');

   Compute_Method(0,5,1,12.0,16.0,0.001,0.0,0.0,0.0,'North_15_25',mtHits,True,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG02b');

   Compute_Method(0,5,1,12.0,16.0,0.001,0.0,0.0,0.0,'North_30_90',mtHits,True,False,Lambda,Omega,B);
   Compute_Omega_Vector;
   Write_Vector('FG02c');


end.