In this study, we construct a bivariate market model combining the copula function with the affine GARCH(p,q) process used to describe the marginal dynamics of the log price. We then provide a numerical procedure for pricing European spread option contracts. To assess the accuracy of our approach we present a comparison with the Monte Carlo simulation method.