Computation of Optical Fiber Mode Field Distribution and Propagation Constants

This program computes the mode field distribution and propagation constants of optical fibers. Users can input specific mode orders to generate corresponding mode patterns for different fiber configurations. The core implementation involves solving the Helmholtz wave equation using a finite-difference method or transfer matrix approach, with key functions handling refractive index profile generation and boundary value problem resolution. Simply execute main.py while adjusting parameters like core radius, refractive indices, and wavelength to analyze structures with up to 5 layered dielectric interfaces. Note that accurate calculation requires careful consideration of material dispersion models, perfectly matched layer (PML) boundary conditions, and numerical convergence criteria - particularly important for applications in telecommunications (e.g., dispersion compensation) and medical technology (e.g., fiber-optic sensors). The code structure separates preprocessing (parameter validation), solver routines (matrix eigenvalue computation), and postprocessing (field visualization via matplotlib).