Solving the equation - uxx - ∈Uyy = f(x, y, u) by an O(h4) finite difference method

E. Lungu, T. Motsumi, Tadeusz Styś

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

The semi-linear equation - uxx - ∈uyy = f(x,y,u) with Dirichlet boundary conditions is solved by an O(h4) finite difference method, which has local truncation error O(h2) at the mesh points neighboring the boundary and O(h4) at most interior mesh points. It is proved that the finite difference method is O(h4) uniformly convergent as h → 0. The method is considered in the form of a system of algebraic equations with a nine diagonal sparse matrix. The system of algebraic equations is solved by an implicit iterative method combined with Gauss elimination. A Mathematica module is designed for the purpose of testing and using the method. To illustrate the method, the equation of twisting a springy rod is solved.

Original languageEnglish
Pages (from-to)395-407
Number of pages13
JournalNumerical Methods for Partial Differential Equations
Volume16
Issue number4
Publication statusPublished - Jul 2000

Fingerprint

Finite difference method
Difference Method
Finite Difference
Algebraic Equation
Iterative methods
Linear equations
Mesh
Semilinear Equations
Truncation Error
Mathematica
Implicit Method
Diagonal matrix
Boundary conditions
Sparse matrix
Dirichlet Boundary Conditions
Gauss
Elimination
Testing
Interior
Iteration

All Science Journal Classification (ASJC) codes

  • Analysis
  • Applied Mathematics
  • Computational Mathematics

Cite this

@article{742daedf885048db988efacaef9e0594,
title = "Solving the equation - uxx - ∈Uyy = f(x, y, u) by an O(h4) finite difference method",
abstract = "The semi-linear equation - uxx - ∈uyy = f(x,y,u) with Dirichlet boundary conditions is solved by an O(h4) finite difference method, which has local truncation error O(h2) at the mesh points neighboring the boundary and O(h4) at most interior mesh points. It is proved that the finite difference method is O(h4) uniformly convergent as h → 0. The method is considered in the form of a system of algebraic equations with a nine diagonal sparse matrix. The system of algebraic equations is solved by an implicit iterative method combined with Gauss elimination. A Mathematica module is designed for the purpose of testing and using the method. To illustrate the method, the equation of twisting a springy rod is solved.",
author = "E. Lungu and T. Motsumi and Tadeusz Styś",
year = "2000",
month = "7",
language = "English",
volume = "16",
pages = "395--407",
journal = "Numerical Methods for Partial Differential Equations",
issn = "0749-159X",
publisher = "John Wiley and Sons Inc.",
number = "4",

}

Solving the equation - uxx - ∈Uyy = f(x, y, u) by an O(h4) finite difference method. / Lungu, E.; Motsumi, T.; Styś, Tadeusz.

In: Numerical Methods for Partial Differential Equations, Vol. 16, No. 4, 07.2000, p. 395-407.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Solving the equation - uxx - ∈Uyy = f(x, y, u) by an O(h4) finite difference method

AU - Lungu, E.

AU - Motsumi, T.

AU - Styś, Tadeusz

PY - 2000/7

Y1 - 2000/7

N2 - The semi-linear equation - uxx - ∈uyy = f(x,y,u) with Dirichlet boundary conditions is solved by an O(h4) finite difference method, which has local truncation error O(h2) at the mesh points neighboring the boundary and O(h4) at most interior mesh points. It is proved that the finite difference method is O(h4) uniformly convergent as h → 0. The method is considered in the form of a system of algebraic equations with a nine diagonal sparse matrix. The system of algebraic equations is solved by an implicit iterative method combined with Gauss elimination. A Mathematica module is designed for the purpose of testing and using the method. To illustrate the method, the equation of twisting a springy rod is solved.

AB - The semi-linear equation - uxx - ∈uyy = f(x,y,u) with Dirichlet boundary conditions is solved by an O(h4) finite difference method, which has local truncation error O(h2) at the mesh points neighboring the boundary and O(h4) at most interior mesh points. It is proved that the finite difference method is O(h4) uniformly convergent as h → 0. The method is considered in the form of a system of algebraic equations with a nine diagonal sparse matrix. The system of algebraic equations is solved by an implicit iterative method combined with Gauss elimination. A Mathematica module is designed for the purpose of testing and using the method. To illustrate the method, the equation of twisting a springy rod is solved.

UR - http://www.scopus.com/inward/record.url?scp=0041688077&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0041688077&partnerID=8YFLogxK

M3 - Article

AN - SCOPUS:0041688077

VL - 16

SP - 395

EP - 407

JO - Numerical Methods for Partial Differential Equations

JF - Numerical Methods for Partial Differential Equations

SN - 0749-159X

IS - 4

ER -