We propose a new method to solve the boundary value problem for a class of second order linear ordinary differential equations, which has a non-negative solution. The method applies the maximum entropy principle to approximating the solution numerically. The theoretical analysis and numerical examples show that our method is convergent
