A new finite element (FE) scheme, based on the decomposition of a second order differential equation into a set of first order symplectic (Hamiltonian) equations, is presented and tested on a one-dimensional, driven Sturm-Liouville problem. Error analysis ...
1996