remove the hardcoded steps, add display of a0 values

search_O1.m

@@ -34,7 +34,8 @@ fmean=f(1)+0.5; % Take the mean value of each 1-Hz sub-band
 g=[-0.1:0.01:0.1]; % asini grid (+/-10% uncertainty, may search larger uncertainty)
 a0=(1+g)*a; % asini values being searched
 
-for m=1:21
+for m=1:length(g)
+    disp(sprintf('a0=%g',a0(m)))
     tic
     % Generate Bessel-weighted matched filter
     [b,h] = binary_fstat(f,fmean,P,a0(m));