The authors analyze and implement a numerical scheme for the efficient numerical solution of Fokker Planck equations for Markov process with jumps. The discontinuous Galerkin discretization in time and a wavelet discretization in space are applied. An error analysis on the approximation is presented. The considered results by numerical experiments are illustrated. Moreover, applications to purely discontinuous Lévy processes arising in finance are given.