function psu=pspec(u)
n=length(u);
ffu=fft(u-mean(u));
psu=ffu.*conj(ffu)/n/n;
%