Abstract
The direct solution of two inhomogeneous ordinary differential equations provides a simple method of calculating accurate non-adiabatic matrix elements when the basic adiabatic wavefunctions are obtained on a grid by finite difference techniques.