p 1 "1 for TE, 1/m^2 for TM"
l 1 "Radial mode number"
nu n_m+1/2 ""
xn "nu/m - Ai(l)/m*(nu/2)^(1/3) + dk(0)/(nu^(0/3) * (m^2 - 1)^((0+1)/2)) + dk(1)/(nu^(1/3) * (m^2 - 1)^((1+1)/2)) + dk(2)/(nu^(2/3) * (m^2 - 1)^((2+1)/2)) + dk(3)/(nu^(3/3) * (m^2 - 1)^((3+1)/2))" "Resonance size parameter"
omega_res "xn *c_const/(n_a*r_res)" "Angular frequency"
f_res omega_res/(2*pi) "Resonance frequency"
lam_res c_const/f_res "Resonance wavelength"