In this talk, I will present our recent results on calculations of the neutrinoless double beta decay matrix elements for two isotopes 76Ge and 150Nd with different approaches. 76Ge is one of the isotopes widely used in various double beta decay experiments such as GERDA, and I will present our recent results for the comparison among different methods: QRPA, shell model, etc. 150Nd is another important candidate with the largest phase space factor, I will give our recent results with deformed QRPA calculations accounting isospin symmetry restoration. These results show that 150Nd is possibly the one with the shortest half-lives for neutrinoless double beta decay among all candidates.