Conservative numerical methods for \(\ddot x=f(x)\). (English) Zbl 0561.65056

Two numerical methods are developed for solving the initial value problem \(\ddot x=f(x)\), \(x(0)=\alpha\), \(\dot x(0)=\beta\). Both methods conserve the same total energy (kinetic plus potential) as does this differential equation, i.e. \({1/2}\dot x(t)^ 2 + \phi(x(t))\). Here \(\phi(x)\) is any function such that \(-d\phi /dx=f(x)\). Two examples describe the situation.
Reviewer: M.Bartušek


65L05 Numerical methods for initial value problems involving ordinary differential equations
34A34 Nonlinear ordinary differential equations and systems
Full Text: DOI


[1] Feynman, R.P.; Leighton, R.B.; Sands, M., The Feynman lectures on physics, (1963), Addison-Wesley Reading, Massachusetts · Zbl 0131.38703
[2] Greenspan, D., Discrete numerical methods in physics and engineering, (1974), Academic Press New York · Zbl 0288.65001
[3] {\scD. Greenspan}, Technical Report No. 205, Dept. of Math., Univ. of Texas at Arlington, Arlington, Texas.
[4] Labudde, R.A.; Greenspan, D., Numer. math., 25, 323-346, (1976)
[5] {\scR. A. Labudde and D. Greenspan}, Int. J. Math. and Math. Sci., in press.
[6] Langdon, A.B., J. comput. phys., 12, 247-268, (1973)
[7] Popov, Y.P.; Samarskii, A.A., USSR comput. math. math. phys. (engl. transl.), 10, 265-273, (1970)
[8] Potter, D., Computational physics, (1973), Wiley New York · Zbl 0314.70008
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.