where \(N^{\prime }_{n,k}= N_{n,k}/2\pi a^2\)  is dimensionless. Propagation constants required in the calculations are determined using the well-known characteristic equation for the scalar modes of a step-index fiber[1]}, [2]}. All numerical calculations are done in MatLab.
