118
Views
7
CrossRef citations to date
0
Altmetric
Original Articles

A new L-stable Simpson-type rule for the diffusion equation

&
Pages 601-607 | Received 23 Aug 2004, Published online: 20 Aug 2006
 

Abstract

The well-known Simpson rule is an optimal two-step fourth-order method which is unconditionally unstable. In the present paper we describe a new L-stable version of the method. A suitable combination of the arithmetic average approximation with the explicit backward Euler formula provides a third-order approximation at the midpoint which, when plugged into the Simpson rule, gives a third-order L-stable scheme. The L-stable Simpson-type rule (LSIMP3) obtained is then employed to derive a third-order time integration scheme for the diffusion equation. Numerical illustrations are provided to compare the performance of the new LSIMP3 scheme with the Crank–Nicolson scheme. While the Crank–Nicolson scheme can produce unacceptable oscillations in the computed solution, the present LSIMP3 scheme can provide both stable and accurate approximations.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.