The Numerical Integration method for solving two point Boundary value problem can be used as an optional method to solve Linear two point boundary value problem. The method is applied for the Dirichlet boundary conditions by adjusting a small positive value and the step size. Like the finite difference method it uses the tridiagonal matrix and the solution of the matrix can be obtained by Thomas Algorithm.