Lecture 14 - Implicit methods for ODEs (Hermite integration) in Mathematica

Implicit methods (Trapezium rule, Hermite rule, Lotkin rule) based on a 2-point Taylor expansion are time-symmetric and conserve energy and symplectic structure when used to integrate Ordinary Differential Equations.