The Collocation Method in the Numerical Solution of Boundary Value Problems for Neutral Functional Differential Equations. Part I: Convergence Results