clc;
length        = 82.00e-3;
height        =  1.90e-3;
width         = 20.00e-3;

e_modul       =  2.10e+11;
poissonratio  =       0.3;
density       =      7850;

mass = density*length*height*width;
mu   = mass / length;
EI   = e_modul*(width)*(height^3)/12;

omega = zeros(12,1);

for n=1:12
   kn = n*pi/length;
   omega(n) = sqrt(EI/mu)*kn*kn/(2*pi);
   fprintf('omega(%2d) = %8.1f\n',n,omega(n,1));
end