Discrete and Continuous Models and Applied Computational ScienceDiscrete and Continuous Models and Applied Computational Science2658-46702658-7149Peoples' Friendship University of Russia named after Patrice Lumumba (RUDN University)3032710.22363/2658-4670-2022-30-1-62-78Research ArticleFinite-difference methods for solving 1D Poisson problemNdayisengaSerge<p>Student of Department of Applied Probability and Informatics</p>1032195775@rudn.ruhttps://orcid.org/0000-0002-9297-9839SevastianovLeonid A.<p>Doctor of Physical and Mathematical Sciences, Professor of Department of Applied Probability and Informatics of Peoples’ Friendship University of Russia (RUDN University), Leading Researcher of Bogoliubov Laboratory of Theoretical Physics, JINR</p>sevastianov-la@rudn.ruhttps://orcid.org/0000-0002-1856-4643LovetskiyKonstantin P.<p>Candidate of Physical and Mathematical Sciences, Associate Professor of Department of Applied Probability and Informatics</p>lovetskiy-kp@rudn.ruhttps://orcid.org/0000-0002-3645-1060Peoples’ Friendship University of Russia (RUDN University)Bogoliubov Laboratory of Theoretical Physics Joint Institute for Nuclear Research01042022301627825022022Copyright © 2022, Ndayisenga S., Sevastianov L.A., Lovetskiy K.P.2022<p style="text-align: justify;">The paper discusses the formulation and analysis of methods for solving the one-dimensional Poisson equation based on finite-difference approximations - an important and very useful tool for the numerical study of differential equations. In fact, this is a classical approximation method based on the expansion of the solution in a Taylor series, based on which the recent progress of theoretical and practical studies allowed increasing the accuracy, stability, and convergence of methods for solving differential equations. Some of the features of this analysis include interesting extensions to classical numerical analysis of initial and boundary value problems. In the first part, a numerical method for solving the one-dimensional Poisson equation is presented, which reduces to solving a system of linear algebraic equations (SLAE) with a banded symmetric positive definite matrix. The well-known tridiagonal matrix algorithm, also known as the Thomas algorithm, is used to solve the SLAEs. The second part presents a solution method based on an analytical representation of the exact inverse matrix of a discretized version of the Poisson equation. Expressions for inverse matrices essentially depend on the types of boundary conditions in the original setting. Variants of inverse matrices for the Poisson equation with different boundary conditions at the ends of the interval under study are presented - the Dirichlet conditions at both ends of the interval, the Dirichlet conditions at one of the ends and Neumann conditions at the other. In all three cases, the coefficients of the inverse matrices are easily found and the algorithm for solving the problem is practically reduced to multiplying the matrix by the vector of the right-hand side.</p>1D Poisson equationfinite difference methodtridiagonal matrix inversionThomas algorithmGaussian elimination1D уравнение Пуассонаметод конечных разностейобращение трехдиагональной матрицыалгоритм Томасаисключение Гаусса[A. N. Tikhonov and A. A. Samarskii, Equations of Mathematical Physics [Uravneniya matematicheskoy fiziki], 7th ed. Moscow: Moscow State University, Nauka, 2004, in Russian.][D. A. Yakovlev, V. G. Chigrinov, and H. S. Kwok, Modeling and optimization of LCD optical performance. New York: Wiley, 2015. DOI: 10.1002/9781118706749.][L. N. Trefethen, Approximation theory and approximation practice. Philadelphia: SIAM - Society for Industrial and Applied Mathematics, 2019.][M. Planitz et al., Numerical Recipes: The Art of Scientific Computing, 3rd ed. New York: Cambridge University Press, 2007.][N. N. Kalitkin and P. V. Koryakin, “Numerical methods [Chislennyye metody],” in Methods of Mathematical Physics, 1st ed. Moscow: Academia, 2013, vol. 2, in Russian.][A. A. Abramov and V. B. Andreyev, “On the application of the method of successive substitution to the determination of periodic solutions of differential and difference equations,” USSR Computational Mathematics and Mathematical Physics, vol. 3, no. 2, pp. 498-504, 1963. DOI: 10.1016/0041-5553(63)90034-x.][A. A. Samarskiy and A. V. Gulin, Numerical Methods of Mathematical Physics [Chislennyye metody matematicheskoy fiziki]. Moscow: Scientific world, 2003, in Russian.][L. H. Thomas, Elliptic problems in linear difference equations over a network. New York: Waston Sci. Comput. Lab. Rept., Columbia University, 1949.][S. B. Gueye, K. Talla, and C. Mbow, “Generalization of the exact solution of 1D Poisson equation with robin boundary conditions, using the finite difference method,” Journal of Electromagnetic Analysis and Applications, vol. 6, no. 12, pp. 372-381, 2014. DOI: 10.4236/jemaa.2014.612038.][S. B. Gueye, K. Talla, and C. Mbow, “Solution of 1D Poisson equation with Neumann-Dirichlet and Dirichlet-Neumann boundary conditions, using the finite difference method,” Journal of Electromagnetic Analysis and Applications, vol. 6, no. 10, pp. 309-318, 2014. DOI: 10.4236/jemaa.2014.610031.][S. B. Gueye, “The exact formulation of the inverse of the tridiagonal matrix for solving the 1D Poisson equation with the finite difference method,” Journal of Electromagnetic Analysis and Applications, vol. 6, no. 10, pp. 303-308, 2014. DOI: 10.4236/jemaa.2014.610031.][N. N. Kalitkin and E. A. Alshina, “Numerical Methods [Chislennyye metody],” in Numerical analysis. Moscow: Academia, 2013, vol. 1, in Russian.][A. Amosov, Y. Dubinsky, and N. Kopchenova, Computational Methods [Vychislitel’nyye metody], 4th ed. St. Petersburg: Lan’, 2021, in Russian.]