In this dissertation, we consider the problem of simulation of stochastic differential equations driven by pure jump Levy processes with infinite jump activity. Examples include, the class of stochastic differential equations driven by stable and tempered stable Levy processes, which are suited for modeling of a wide range of heavy tail phenomena. We replace the small jump part of the driving Levy process by a suitable Brownian motion, as proposed by Asmussen and Rosinski, which results in a jump-diffusion equation. We obtain Lp [the space of measurable functions with a finite p-norm], for p greater than or equal to 2, and weak error estimates for the error resulting from this step. Combining this with numerical schemes for jump diffusion e...