Numerical approximation of nonlinear BVPs by means of BVMs