% wcisk
% N i m
Ms=input('Podaj ms: ')
k=2;
mi=0.1;
dw=0.035;
dz=0.053;
DZ=0.064;
DW=0.053;
l=0.036;

p=(2*k*Ms)/(mi*pi*dz^2*l)% 5.6 MPa

%Lame
Ew=100000000000;
Ez=210000000000;
ni_w=0.35;
ni_z=0.3;
wl=dz*p*((1/Ew)*((dz^2+dw^2)/(dz^2-dw^2)+ni_w)+(1/Ez)*((DZ^2+DW^2)/(DZ^2-DW^2)+ni_z)) %% 17 mikrom
a=0.6;
Rzw=3.2*10^-6;
Rzz=6.3*10^-6;

%calkowity
disp('calkowity w mikrometrach')
w_prim=(wl+2*a*(Rzw+Rzz))*1000000