원통형 코일의 인덕턴스(R=반지름, L=길이, N=권수)

μ0 = 4*pi*10^(-7);
D=2*R;
x=D/L;
k=1/sqrt(1+(L/(2*R))^2);
k2 = k²;
Kn = (4 / (3 * pi * sqrt(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*pi*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 /(pi*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);