Comparison between the performance of the second and first degree stationary iterative techniques is performed on two representations to two point boundary value problems of the second and fourth orders. The numerical treatment of the Fredholm integral equation representation for the second and fourth order boundary value problems (BVP) has illustrated the effective use of their integral representation. The finite difference method with the same accuracy is employed to construct linear systems from the differential or the equivalent Fredholm form. Second degree linear stationary iterative methods are extensively used. The second degree Gauss-siedle (S.G.S) method is presented, by measuring the asymptotic rates of conver-gence of the sequence; we are able to determine when the Gauss-Seidle second-degree iterative method is superior to its cor-responding first-degree one. Two numerical examples are considered, one of them is of the second degree BVP and the other is fourth order BVP. All calculations are done with the help of computer algebra system (MATHEMATICA 10.2).