In this paper, we present algorithms to solve a complex system of partial integro-differential equations (PIDE's) of parabolic type. The system is motivated by applications in finance where the solution of the system gives the price of European options in a regime-switching jump diffusion model....