%
%  one_third_octave_frequencies.m  ver 1.0  October 19, 2012
%
function[fl,fc,fu,imax]=one_third_octave_frequencies()
%
fc(1)=0.02;   % 0.02 Hz needed for ISO 2631 calculation
fc(2)=0.025;  
fc(3)=0.0315;  
fc(4)=0.04;
fc(5)=0.05;
fc(6)=0.063;
fc(7)=0.08;
fc(8)=0.1;
fc(9)=0.125;
fc(10)=0.16;
%
i=1;
while(1)
    fc(i+10)=10*fc(i);
    if(fc(i+10)>=20000)
        break;
    end
    i=i+1;
end
%
imax=length(fc);
%
fl=fc/(2^(1/6));
%
for i=1:imax-1
    fu(i)=fl(i+1);
end
%
fu(imax)=fc(imax)*(2^(1/6));