A finite element method involving collocation method with quintic B-splines as basis functions have been developed to solve sixth order boundary value problems. The sixth order and fifth order derivatives for the dependent variable are approximated by the central differences of fourth 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 several linear and 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 quasilinearization technique. Numerical results obtained by the present method are in good agreement with the exact solutions or numerical solutions available in the literature.