lda0 1.55[um] "Operating wavelength" f0 c_const/lda0 "Operating frequency" height 4000[nm] "Height of waveguide" n_core_o sqrt(epsr_o) "Core refractive index, ordinary" n_core_eo sqrt(epsr_eo) "Core refractive index, extraordinary" eps_clad 2.05 "Cladding relative permittivity" k0_1 2*pi/lda0 "Wavenumber in vacuum" theta 0[deg] "x-y plane rotation" phi 45[deg] "y-z plane rotation" zeta 90[deg] "x-z plane rotation" epsr_eo 2.19 "Core relative permittivity, extraordinary" epsr_o 2.31 "Core relative permittivity, ordinary"