円筒コイルのインダクタンス(R=半径, L=長さ,N＝巻数)

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