kas = find(abs(fft(sin(0.5346*n),2048))>400);
ka = kas(1);
kbs = find(abs(fft(sin(0.5906*n),2048))>400);
kb = kbs(1);
kcs = find(abs(fft(sin(0.6535*n),2048))>400);
kc = kcs(1);
kds = find(abs(fft(sin(0.7217*n),2048))>400);
kd = kds(1);
kes = find(abs(fft(sin(0.9273*n),2048))>400);
ke = kes(1);
kfs = find(abs(fft(sin(1.0247*n),2048))>400);
kf = kfs(1);
kgs = find(abs(fft(sin(1.1328*n),2048))>400);
kg = kgs(1);
k = [ka kb kc kd ke kf kg];
omegak = 2*pi*k/2048;

% Results are:
% k =
%    175   193   214   236   303   335   370
% omegak =
%    0.5369    0.5921    0.6565    0.7240    0.9296    1.0278    1.1351
