n=100;
f=600+200*rand(n,1);
r=0.02+0.004*rand(n,1);
emod=(185+15*rand(n,1))*1e9;
l=1+0.05*rand(n,1);
inert=pi*r.^4/4;