Only check proper ni

            % Y1=Y1(1:round(length(Y1)/2));

            % toc

            % tic

            if n ~= 1 && n ~= length(f_gr)

                cond = F1>=0.5*(nizero+f_gr(n-1)) & F1<0.5*(nizero+f_gr(n+1));

            elseif n == 1

                cond = F1>=(nizero) & F1<0.5*(nizero+f_gr(n+1));

            else

                cond = F1>=0.5*(nizero+f_gr(n-1)) & F1<=nizero;

            end

            Y1=Y1(cond); clear F1

            % if n ~= 1 && n ~= length(f_gr)

            %     cond = F1>=0.5*(nizero+f_gr(n-1)) & F1<0.5*(nizero+f_gr(n+1));

            % elseif n == 1

            %     cond = F1>=(nizero) & F1<0.5*(nizero+f_gr(n+1));

            % else

            %     cond = F1>=0.5*(nizero+f_gr(n-1)) & F1<=nizero;

            % end

            [nnfreq,nnind] = min(abs(F1-nizero));

            % cond = F1>=F1(nnind-1) & F1<=F1(nnind+1);

            Y1=Y1(nnind); clear F1

            %Calcolo della detection statistic

            % Lambda(i,n,m)=sum(abs(Y1).^2)/sum(xtemp(:)); %CREDO (oppure prendono la potenza massima?)