Mathematica 9 is now available

Wolfram Library Archive

Courseware Demos MathSource Technical Notes
All Collections Articles Books Conference Proceedings

Solving the Equations - uxx - euyy = f(x,y,u) by an O(h^4) Finite Difference Method

E. Lungu
T. Motsumi
T. Stys
Journal / Anthology

Numerical Methods for Partial Differential Equations
Year: 2000
Volume: 16
Issue: 4
Page range: 395-407

The semi-linear equation - uxx - euyy = f(x,y,u) with Dirichlet boundary conditions is solved by an O(h^4) finite different method, which has local truncation error O(h^2) at the mesh points neighboring the boundary and O(h^4) at most interior mesh points. It is proved that the finite difference method is O(h^4) 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 ht method. To illustrate the method, the equation of twisting a springy rod is solved.

*Applied Mathematics > Numerical Methods
*Mathematics > Calculus and Analysis > Differential Equations