A finite element method involving collocation method with Cubic B-splines as basis functions have been developed to solve sixth order boundary value problems. The third order, fourth order, fifth order and sixth order derivatives for the dependent variable are approximated by the central differences of second order derivatives. The basis functions are redefined into a new set of basis functions which in number match with the number of collocated points selected in the space variable domain. The proposed method is tested on six linear and five non-linear boundary value problems. The solution of a non-linear boundary value problem has been obtained as the limit of a sequence of solutions of linear boundary value problems generated by quasi-linearization technique. Numerical results obtained by the present method are in good agreement with the exact solutions available in the literature.