Local solution method for numerical solving of the wave propagation problem Original Research Article

A new method for numerical solving of boundary problem for ordinary differential equations with slowly varying coefficients which is aimed at better representation of solutions in the regions of their rapid oscillations or exponential increase (decrease) is proposed. It is based on the approximation of the solution sought for in the form of a superposition of certain polynomial-exponential basic functions. The method is studied for the Helmholtz equation in comparison with the standard finite difference and finite element methods. The numerical tests have shown the convergence of the method proposed. In comparison with the standard methods the same accuracy is obtained on substantially coarser mesh. This advantage becomes more pronounced, if the solution varies very rapidly.