Corrected bruttone indexing

dt_psd=1/(2*Nyq);

dt = dt_psd;

N = fix(Tseg/dt);

f_supamin = 49

f_supamax = 1550

f_supamin = 49 % Hz

f_supamax = 1550 % Hz

fr_int_width = 100 % Hz (-1)

n_fr_intervals = (f_supamax-f_supamin-1)/fr_int_width;

Porb = 68023.92;

porb_unc = 0.02; % s Killestein

porb_min = Porb - 3*porb_unc; % s

% a_step = a_max;

tol = 0.9

for fifo = 1:15

    f_min = f_supamin + (100*(fifo-1))

    f_max = f_supamax - 1400 + (100*(fifo-1))

for fifo = 1:n_fr_intervals

    f_min = f_supamin + (fr_int_width*(fifo-1))

    f_max = f_supamax - fr_int_width*(n_fr_intervals - fifo)

    MJDREF = zeros(numfiles,1);

    f_gr = (f_min:f_step:f_max).';

    lefgr = length(f_gr);

    %fprintf("The workspace will need %.5f KB or %.5f MB or %.5f GB of memory", inKB, inMB, inGB)

    timerapp = 24.0*summina/518260001575;

    bruttone(1+(fifo-1)*(lefgr):(fifo)*(lefgr)) = bru;    

    bruttone((1+(fifo-1)*(fr_int_width/f_step)):((fifo*fr_int_width+1)/f_step)+1) = bru;    

    % disp(bestpar);

    % save('C:\Users\Filippo\Desktop\XMM_Jxxx\risultelli.mat');

    disp(['I am about to save around ',num2str(round(inGB)),' GB of data (+ headers), it should take roughly ',num2str(ceil(timerapp)),' minutes.']);