Barel any change

        nismax(m,s) = f_max*a_max*(Omega_max^s)*singah(s);

    end

end

%SE T0 È IL MIDPOINT DI TUTTA L'OSSERVAZIONE, ALLORA NIMIN E NIMAX NON

%DIPENDONO DA M E SONO UGUALI PER TUTTI I SEGMENTI. CONTROLLARE!!

% Try s* and check \nu_s range

g_jj=((pi*T)^2)/3.*[1; (T^2)/60; (T^3)/1344; (T^4)/172800]; %eq. 22 M2015 + calcoli da eq. 21 M2015

mu_s=0.999; %massimo mismatch da scegliere

mu_s=0.001; %massimo mismatch sulla griglia coerente da scegliere

s_s=int16(4);

while(1)

    for m=1:M

    end

end

% Resampling the times over nism{m,s} templates

% ---- Andrà fatto per ogni possibile combinazione di nism a un dato m mi

% sa

%Fourier transform on original time-series --------------------------------

%per mantenere l'informazione di fase, non faccio il valore assoluto al quadrato della fft

    Y=Y(cond); 

    X=ifft(Y); %inverse-fourier transf.

    %SCRIVERE CODICE PER RESAMPLED TIME SERIES!!!!

    %X1 è la timeseries resempled (controllare che sia un vettore colonna)

    %X1 è la timeseries resempled (controllare che sia un vettore colonna)

    %zero-padding (metto gli zeri alla fine) -------------------------------

    X1 = [X1;zeros(L-length(F),1)]; %concateno alla timeseries downsapled and resampled gli zeri

    X1 = [X1;zeros(L-length(F),1)]; %concateno alla timeseries downsampled and resampled gli zeri

    %il secondo è un vettore colonna di zeri con lunghezza pari al numero di elementi che mi serve per tornare alla risoluzione originale

    %Fourier transform -----------------------------------------------------