ω = 3140;
gama1 = (1/0.018);
gama2 = 10;
mir1 = 5;
mir2 = 300;
mi0 = 4  π*(10^(-6));
μ1 = mir1*mi0;
μ2 = mir2*mi0;
a1 = (Sqrt[(2/(ω*gama1*μ1))]);
a2 = (Sqrt[(2/(ω*gama2*μ2))]);

D1 = 0.5;
D2 = 0.4;
k1 = 1.10;
k2 = 1.04;
l1 = 1.3;
l2 = 1;
alfa1 = 0.75;
alfa2 = 0.75;
F = 0.8;
d1 = D1 + a1;
d2 = D2 - a2;
R1 = k1*((π*d1)/(gama1*a1*l1))
R2 = k2*((π*d2)/(gama2*a2*l1))
L1 = μ1*((π*(d1^2)*alfa1)/(4*l1))
L2 = mi0*((π*(d2^2)*alfa2)/(4*l2))
M = mi0*((π*((d2^2)/4))/l2)*F
p = (Sqrt[(((ω^2)*(M^2))/((R2^2)*(ω^2) (L2^2)))])