Galerkin Method with Modified Shifted Lucas Polynomials for Solving the 2D Poisson Equation

This study looks at how to solve the two-dimensional Poisson equation, a math problem common in physics and engineering. We focus on spectral methods, which are good at solving problems with smooth solutions. We introduce a spectral Galerkin method that uses tensor products of modified shifted Lucas polynomials. These polynomials haven’t been used this way before. By adding a factor of x(1-x) to the Lucas polynomials, our method automatically meets certain boundary conditions, which makes it easier to use while keeping its accuracy. Our goal is to create and test this method for solving the Poisson equation on a square. We create fast algorithms for putting together matrices and study how well the method converges using math and computer experiments. The tests show that our method has similar convergence rates to other methods like Chebyshev and Legendre. The errors go down exponentially for smooth source terms. The method is efficient and has good conditioning, which suggests that Lucas polynomials could be a good alternative to regular polynomials in spectral methods. This research could lead to using Lucas polynomial-based spectral methods for more general problems.