Added quasi-logarithmic grid for visual inspection

ted = (0:dt:Tseg);

tm = (0.5*dt:dt:(Tseg-0.5*dt));

F1=gpuArray((0:N-1)./(N*dt)).';

F0 = F1(2:N/2+1);

F02 = F0(end/2:end);

F0 = F1(1:N/2+1);

% F02 = F0(end/2:end);

nblog = 200;

Flog = logspace(log10(F1(2)/2),log10(F0(end)+0.5/Tseg),nblog);

for i = 1:length(Flog)-1

[~,minind] = min(abs(F0(:)-0.5/512-(Flog(i))));

Fqlog(i) = F0(minind)-0.5/Tseg;

end

Fqlog(nblog) = F0(end)+0.5/Tseg;

Fqlog = unique(Fqlog);

Fqlmid(1:length(Fqlog)-1) = 0.5*(Fqlog(2:end)+Fqlog(1:end-1));

Fqlbin(1:length(Fqlog)-1) = Fqlog(2:end)-Fqlog(1:end-1);

for i = 1:length(mario)

    filoc = [pwd,'/',char(mario(i))]

    % file = fits.openFile(filoc)

    xtemp = gpuArray(xtemp);

    Y0 = abs(fft(xtemp));

    Y0 = 2.0.*(Y0(1:N/2+1).^2)./Y0(1);

    % Y02 = Y0(end/2:end);

    for nino = 1:length(Fqlbin)

        Yqlbin = sum(Y0((Fqlog(nino)<=F0)&&(F0<Fqlog(nino+1))));

    end

    Yqlbin(:) = Yqlbin(:)./Fqlbin(:);

    semilogx(Flmid,Yqlbin,''+'')

    pause

end

% length(xtemp(xtemp == 1))/N

% mean(xtemp)

% round(length(xtemp(xtemp == 0)))/N

% ptot0 = length(xtemp(xtemp == 0))/N;

% mi = -log(ptot0);

ptot0 = length(xtemp(xtemp == 0))/N;

mi = -log(ptot0);

% creps = 1 - (length(xtemp(xtemp == 1))/N)/(mi*ptot0);

% ncross = 8;

% crn = 8;

% ncr = 8;

% epscr = 1 - (length(xtemp(xtemp == 1))/N)/(mi*ptot0);

ncr = 8;

epscr = 1 - (length(xtemp(xtemp == 1))/N)/(mi*ptot0);

% clear creps ncross crn

pcr = 1- (1-epscr)^(1/ncr);

expec1 = 1+ epscr + (3*(ncr-1)/(2*ncr))*epscr^2;