Sistemato pezzo che ci avrebbe messo ~88 ore, ora ci mette 0.05 secondi.

x = zeros(M,N);

toc

tic

aux1 = M*N;

aux2 = C(1:aux1);

x = reshape(aux2,[M,N]);

toc

tic

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% Lasciarlo così per circa 13E+6 bin temporali richiede ~88 ore (su questo

%% computer): trovare modo più sensato

for m = 1:M

    for j = 1:N

        pippo = t>=tm(m,j)-dt/2 & t<tm(m,j)+dt/2;

        Cmj=C(pippo); clear pippo

        x(m,j)=sum(Cmj);

    end

end

%for m = 1:M

%    for j = 1:N

%        pippo = t>=tm(m,j)-dt/2 & t<tm(m,j)+dt/2;

%        Cmj=C(pippo); clear pippo

%        x(m,j)=sum(Cmj);

%    end

%end

toc

% Try s* and check \nu_s range

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