function [coef] = fund(data)
[NN, NT]=size(data);
if NT == 1;
    data=data';
    [NN, NT]=size(data);
end

for nn=1:NN;
H=fft(data(nn,:));
coef(nn,1)=abs(H(2)/(NT/2));
end
