Inductance of a Cylindrical Coil (R = radius, L = length, N = number of turns)

μ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);