Jump diffusion models are being used more and more often in financial applications. Consisting of a Brownian motion (with drift) and a jump component, such models have a number of parameters that have to be set at some level. Maximum Likelihood Estimation (MLE) turns out to be suitable for this task, however it is computationally demanding. For a complicated likelihood function it is seldom possible to find derivatives. The global maximum of a likelihood function defined for a jump diffusion model can however, be obtained by numerical methods. I chose to use the Bound Optimization BY Quadratic Approximation (BOBYQA) method which happened to be effective in this case. However, results of Maximum Likelihood Estimation (MLE) proved to be hard to interpret.