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

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

mu   = 
EI   = 

omega = zeros(12,1);

for n=1:12
   kn = 
   omega(n) = 
   fprintf('omega(%2d) = %8.1f\n',n,omega(n,1));
end