clear;
load whistle_data.mat

for n = 1:5000
  A = trapz(t,y.*sin(2*pi*t*n));
  B = trapz(t,y.*cos(2*pi*t*n));
  C(n) = sqrt(A^2 +B^2)*2;
 
end
loglog([1:5000],C)