We report alternating group explicit (AGE) iteration method to solve the tri-diagonal linear system derived from a new finite difference discretization of sixth order accuracy of the two point singular boundary value problem
,
,
and 2 subject to boundary conditions
,
, where A and B are finite constants. We also discuss Newton-AGE iteration method for the sixth order numerical solution of two point non-linear boundary value problem. The proof for the convergence of the AGE iteration method when the coefficient matrix is real and unsymmetric is discussed. Numerical results are presented to illustrate the proposed iterative methods.