Two point Gauss quadrature formula has been used to calculate the integral values. Owing to the nonlinearity of the systems
of equations an iterative scheme is required to solve the nonlinear algebraic matrix system. The system is linearized by
incorporating known function f , h, h, which are assumed to be known values of the functions f , h, h. An initial guess of zero istaken at each node point. After applying the given boundary conditions, remaining system of equations has been solved
using Gauss-elimination method. This gives us new values of unknowns. This process continues till the absolute differences
of two successive iterate value of unknowns is less than the accuracy of 104.