A novel approach to the numerical solution of boundary value problems on infinite intervals