r sqrt(x^2+y^2+z^2+eps) "Radial coordinate"
rh sqrt(x^2+y^2+eps) "Radial coordinate in xy-plane"
dr 1e-2*b[1/m] "Refractive index smoothing function"
cos_theta z/r ""
sin_theta rh/r ""
pr px*sin_theta*cos_phi+py*sin_theta*sin_phi+pz*cos_theta "Wave vector, r-component"
p_theta px*cos_theta*cos_phi+py*cos_theta*sin_phi-pz*sin_theta "Wave vector, theta-component"
p_phi -px*sin_phi+py*cos_phi "Wave vector, phi-component"
n_r_cloak b/(b-a) "Refractive index in cloak, r-component"
n_theta_cloak b/(b-a)*(r-a)/r "Refractive index in cloak, theta-component"
n_phi_cloak n_theta_cloak "Refractive index in cloak, phi-component"
cos_phi x/rh ""
sin_phi y/rh ""
n_r n_air+(n_r_cloak-n_air)*flc2hs(b[1/m]-r[1/m],dr) "Refractive index, r-component"
n_theta n_air+(n_theta_cloak-n_air)*flc2hs(b[1/m]-r[1/m],dr) "Refractive index, theta-component"
n_phi n_air+(n_phi_cloak-n_air)*flc2hs(b[1/m]-r[1/m],dr) "Refractive index, phi-component"
H_photon c_const*sqrt(pr^2/n_r^2+p_theta^2/n_theta^2+p_phi^2/n_phi^2+eps) "Dispersion relation"