圆柱线圈的电感（R=半径，L=长度，N=匝数）

μ0 = 4*π*10⁻⁷;
D=2*R;
x=D/L;
k=1/√(1+(L/(2*R))²);
k2=k^2;
Kn=(4/(3*π*√(1-k2))*(((1-k2)/k2)*ellipticK(k)-((1-2*k2)/k2)*ellipticE(k)-k));
S = π*R²;
Lcoil=Kn*ur*μ0*S*N²/L*10⁶;
print(D);
print(x);
print(Kn);
print(Lcoil);
print(S);
ω = 2*π*frequency;
LcoilH = Lcoil/10^6;
XL = ω*LcoilH ;
C = 1/(ω*XL);
CpF = C *10^12;
print(XL);
print(CpF);
LengthWire = pi*D*N;
print(LengthWire);
λ=(299792458)/frequency;
Rr = (8/3)*pi^3*120*pi*(N*S/λ^2)^2;
print(Rr);
δ = sqrt(2*ρ/(ω*μ0));
radius = radius_mm/1000;
Rdc = ρ*LengthWire /(π*radius^2);
Rac = ρ*LengthWire /(π*(radius^2-(radius-δ)^2));
Diameter_mm = 2*radius_mm;
δmicron = δ*10^6;
print(Diameter_mm );
print(Rdc);
print(δmicron);
print(Rac);
η = Rr/(Rr+Rac)*100;
print(η);
heff = 2*pi*N*S/λ;
print(heff);
Q = XL/Rac;
print(Q);
print(ρ);
Vo = Q*heff*Eradio;
Rp = Q^2*(Rr+Rac);
print(Vo);
print(Rp);