Recent development in numerical simulations of molecular beam epitaxial equation