Numerical solution of Boundary value problems