In this paper, we propose a partial differential equation formulation for the value of an option when the underlying asset price is described by a discrete-time GARCH process. Our numerical approach involves a spectral Fourier-Chebyshev interpolation. Numerical illustrations are provided, and the results are compared with other available valuation methods. Our numerical procedure converges exponentially fast and allows for the efficient computation of option prices, achieving a high level of precision in a few seconds of computing time.