Abstract
The smoothing problem for two-point boundary value models is investigated. It is shown that by triangularizing the relevant hamiltonian system, a new smoothing algorithm can be derived which is twice as fast as that obtained by diagonalizing the hamiltonian.