We apply interior penalty discontinuous Galerkin finite element method (dGFEM) for pricing of European and American options for Heston PDE. The advantages of dGFEM space discretization with Rannacher smoothing as time integrator with non-smooth initial and boundary conditions are illustrated in several numerical examples for European call, as well as butterfly spread, digital call and American put options. The convection dominated Heston PDE for vanishing volatility is efficiently solved utilizing the adaptive dGFEM algorithm. The linear complementary problem for the American option is solved using the norm preconditioned projected successive over relaxation (PSOR) method. Numerical experiments illustrate that dGFEM is an accurate and efficient method for pricing options.