The vibrational multi-reference coupled cluster method is developed to calculate the vibrational excitation energies of polyatomic molecules. The method is implemented on ozone and formaldehyde molecules and the results are compared with full vibrational configuration interaction (FVCI) method. A good agreement is found between the vibrational multi-reference coupled cluster method and converged FVCI method for lower lying vibrational states.