Wavelength < 0.1λ: Wire Antenna Gain and Matching Coil at 828 kHz

c = 2.99792458*10^8;
λ = c/frequency;
Rr = (40*pi^2)*(height/λ)^2;
diameter = diametermm/1000;
Za = 60 * (ln(4 * height / diameter) - 1);
θ = 2 * pi * height / λ;
Xc = Za / tan(θ);
Xa = Za * tan(θ);
XL = Xc - Xa;
ω = 2 * pi * frequency;
L = XL/ω;
Ca = 1/(ω*Xc);
RL = XL/Q;
μ0 = 4*pi*10^(-7);
δ = sqrt(2*ρ/(ω*μ0));
radius = diameter /2;
Rdc = ρ * height / (π * radius²);
Rac = ρ * height / (π * (radius² - (radius - δ)²));
δmicron = δ * 10⁶;
η = Rr/(Rr+Rac+RL)*100;
heff = height/2;
Vo = Q*heff/2*Eradio;
Rp = Q^2*(Rr+Rac+RL);
print(λ);
print(Rr);
print(Za);
print(θ);
print(Xc);
print(Xa);
print(XL);
print(L*10^6);
print(RL);
print(η);
print(heff);
print(δmicron);
print(Rdc);
print(Rac);
print(ρ);
print(Ca*10^12);
print(Vo);
print(Rp);