LAUSR

In this paper, we present a method based on the combination of the finite difference and the meshless techniques for solving 2D generalized telegraph equations in single and multi-connected domains. The three-layer Crank–Nicolson time-stepping scheme transforms the equation into a sequence of second-order elliptic partial differential equations with mixed derivatives and variable coefficients. The approximate solution on each time layer is sought as series over radial basis functions based approximations. These are constructed in such a way that they satisfy the homogeneous boundary conditions. Thus, any their linear combination satisfies the homogeneous boundary condition too. The coefficients of the series are enforced to satisfy the governing equation in the solution domain. Several numerical examples are presented to demonstrate stability and accuracy.