Vol 30, No 2 (2022)

We consider Cauchy problem for ordinary differential equation with solution possessing a sequence of multiple poles. We propose the generalized reciprocal function method. It reduces calculation of a multiple pole to retrieval of a simple zero of accordingly chosen function. Advantages of this approach are illustrated by numerical examples. We propose two representative test problems which constitute interest for verification of other numerical methods for problems with poles.

The classical pseudospectral collocation method based on the expansion of the solution in a basis of Chebyshev polynomials is considered. A new approach to constructing systems of linear algebraic equations for solving ordinary differential equations with variable coefficients and with initial (and/or boundary) conditions makes possible a significant simplification of the structure of matrices, reducing it to a diagonal form. The solution of the system is reduced to multiplying the matrix of values of the Chebyshev polynomials on the selected collocation grid by the vector of values of the function describing the given derivative at the collocation points. The subsequent multiplication of the obtained vector by the two-diagonal spectral matrix, ‘inverse’ with respect to the Chebyshev differentiation matrix, yields all the expansion coefficients of the sought solution except for the first one. This first coefficient is determined at the second stage based on a given initial (and/or boundary) condition. The novelty of the approach is to first select a class (set) of functions that satisfy the differential equation, using a stable and computationally simple method of interpolation (collocation) of the derivative of the future solution. Then the coefficients (except for the first one) of the expansion of the future solution are determined in terms of the calculated expansion coefficients of the derivative using the integration matrix. Finally, from this set of solutions only those that correspond to the given initial conditions are selected.

The model of adiabatic waveguide modes (AWMs) in a smoothly irregular integrated optical waveguide is studied. The model explicitly takes into account the dependence on the rapidly varying transverse coordinate and on the slowly varying horizontal coordinates. Equations are formulated for the strengths of the AWM fields in the approximations of zero and first order of smallness. The contributions of the first order of smallness introduce depolarization and complex values characteristic of leaky modes into the expressions of the AWM electromagnetic fields. A stable method is proposed for calculating the vertical distribution of the electromagnetic field of guided modes in regular multilayer waveguides, including those with a variable number of layers. A stable method for solving a nonlinear equation in partial derivatives of the first order (dispersion equation) for the thickness profile of a smoothly irregular integrated optical waveguide in models of adiabatic waveguide modes of zero and first orders of smallness is described. Stable regularized methods for calculating the AWM field strengths depending on vertical and horizontal coordinates are described. Within the framework of the listed matrix models, the same methods and algorithms for the approximate solution of problems arising in these models are used. Verification of approximate solutions of models of adiabatic waveguide modes of the first and zero orders is proposed; we compare them with the results obtained by other authors in the study of more crude models.