I think I've an error in the function "fiber".
The 2 first lines
nx = round((sum(r)+side)/dx);
ny = round((sum(r)+side)/dy);
must be replaced by
nx = round((max(r)+side)/dx);
ny = round((max(r)+side)/dy);
as r is the radius of the layer and not its thickness.
Nice job, and looks to be very useful.
I'am interested by scalar solutions, but I see only semi-vectorial or full-vectorial solutions. Where is the scalar method hidden
Does a tutorial exist which may describe the exact role and definition of each variable describing the guide? My application is for scalar modes of optical fiber.
Thank you for your help.
Hi, thanks for posting this package online. I wonder that if you can tell me what references used for the formulation in 'svmode'. I am using the book 'Introduction to Optical Waveguide Analysis:Solving Maxwell's Equations and the Schrodinger Equation' by Kenji Kawano. It seems the equations used here for aw, ae, an and as are partially different from that in the book.