Numerical solution of multi-point boundary value problems in Chebyshev series theory of the method