This paper illustrates the use of certain orthogonal polynomials (Chebyshev polynomials) as trial functions in a variational iterative approach for solving 6th order boundary value problems (BVPs). The proposed method is validated through some numerical examples without any form of discretization or perturbation. The resulting numerical evidences show that the method is effective and accurate. All calculations are analyzed using maple18 software
