Nagaoka Coefficient Ks and Inductance L of a Square Loop Antenna

μ0 = 4*pi*10^(-7);
k0 = p/w;
k1 = w/p;
k2 = k1^2;
k3 = k1^3;
ksr1 = sqrt(k2+1);
ksr2 = sqrt(k2+2);
Ks = (2/pi)*((k1-1/k1)*ln(ksr1/(ksr2+1)) + k1*ln((ksr1+1)/k1) + 2*atan(k1/ksr1) +
(1/k1)*(ksr1 - ksr2 - ln(sqrt(2) + 1) + (1/3)*(k3 + (k2 + 2)^(3/2) - 2*(k2 + 1)^(3/2) + sqrt(2) - 1)));
L = Ks * μ0 * p^2 * N^2 / w;
LmicroH = L * 10^6;
S = p^2;
print(k0);
print(Ks);
print(LmicroH);
print(S);
ω = 2*pi*frequency;
XL = ω*L;
C = 1/(ω*XL);
CpF = C *10^12;
print(XL);
print(CpF);
LengthWire = 4*p*N;
print(LengthWire);
λ = (3*10^8)/frequency;
Rr = (8/3)*pi^3*120*pi*(N*S/λ^2)^2;
print(Rr);
ρ = 1.68*10^(-8);
δ = 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*p^2/λ;
print(heff);
Q = XL/Rac;
print(Q);
print(ρ);
Vo = Q*heff*Eradio;
print(Vo);
Rp = (Rr+Rac)*Q^2;
print(Rp);