Changed grid

% end

% Tseg=512; %segment length in seconds

f_gr = (500:1/Tseg:600).';

porb_gr=(16265:(51.5151515151515):17965).';

a_gr = (0.25:3e-3:0.4).';

MJDREF = 58107.1236342593;

tasc_gr = (((58106.9994:3.532451105175235e-4:58107.1095) - MJDREF).*86400).';

f_min = min(f_gr)

f_max = max(f_gr)

porb_min = min(porb_gr)

porb_max = max(porb_gr)

a_min = min(a_gr)

a_max = max(a_gr)

tasc_min = min(tasc_gr)

tasc_max = max(tasc_gr)

% f_gr = (500:1/Tseg:600).';

% porb_gr=(16265:(51.5151515151515):17965).';

% a_gr = (0.25:3e-3:0.4).';

% MJDREF = 58107.1236342593;

% tasc_gr = (((58106.9994:3.532451105175235e-4:58107.1095) - MJDREF).*86400).';

f_min = 550

f_max = 600

porb_min = 17115.3216 % s

porb_max = 17115.3216; porb_gr = (porb_max).'; % s

a_fact = porb_max/(2.0*pi*299792458);

a_min = a_fact*256000*0.04

a_max = a_fact*280000*0.72

tasc_min = (5.8107009477e+04 - MJDREF).*86400 - porb_max/2

tasc_max = (5.8107009477e+04 - MJDREF).*86400 + porb_max/2

% % Parameters for Sco X-1

%  % s 

% % porb_step = (porb_max-porb_min)/ceil((porb_max-porb_min)/stepparino)

% porb_gr = (porb_min:porb_step:porb_max).';

%%

% % stepparino = delasenc(porb_min,f_max,toff,Tseg,M,tol)

delasenc(porb_min,f_max,toff,Tseg,M,tol)

stepparino = delasenc(porb_min,f_max,toff,Tseg,M,tol)

% delasenc(porb_min,f_max,toff,Tseg,M,tol)

1.0/(2.0*f_max*tol)

%%

% a_step = (a_max-a_min)/ceil((a_max-a_min)/stepparino)

a_step = (a_max-a_min)/ceil((a_max-a_min)/stepparino)

% % a_step = 1.0/(2.0*f_max*0.9)

%%

deltasc(a_max,porb_min,f_max,toff,Tseg,M,tol)

% deltasc(a_max,porb_min,f_max,toff,Tseg,M,tol)

0.1025*porb_min/(pi*f_max*a_max*tol)

% % stepparino = max(deltasc(a_max,porb_min,f_max,toff,Tseg,M,tol),0.1025*porb_min/(pi*f_max*a_max*tol))

stepparino = deltasc(a_max,porb_min,f_max,toff,Tseg,M,tol)

%%

% % tasc_step = (tasc_max-tasc_min)/ceil((tasc_max-tasc_min)/stepparino)

% % a_gr = (a_min:a_step:a_max).';

% % tasc_gr = (tasc_min:tasc_step:tasc_max).';

tasc_step = (tasc_max-tasc_min)/ceil((tasc_max-tasc_min)/stepparino)

a_gr = (a_min:a_step:a_max).';

tasc_gr = (tasc_min:tasc_step:tasc_max).';

trupar = [5.924214682724856e+02,1.71155216592e+04,0.343356,5.8107009477e+04];

% f_min = min(f_gr)