In this work we investigate the pricing of swing options in a model where the underlying asset follows a jump diffusion process. We focus on the derivation of the partial integro-differential equation (PIDE) which will be applied to swing contracts and construct a novel pay-off function from a tree-based pay-off matrix that can be used as initial condition in the PIDE formulation. For valuing swing type derivatives we develop a theta implicit-explicit finite difference scheme to discretize the PIDE using a Gaussian quadrature method for the integral part. Based on known results for the classical theta-method the existence and uniqueness of solution to the new implicit-explicit finite difference method is proven. Various numerical examples illustrate the usability of the proposed method and allow us to analyse the sensitivity of swing options with respect to model parameters. In particular the effects of number of exercise rights, jump intensities and dividend yields will be investigated in depth.
