Recently, the valuation of variable annuity products has become a hot topic in actuarial science. In this paper, we use the Fourier cosine series expansion (COS) method to value the guaranteed minimum death benefit (GMDB) products. We first express the value of GMDB by the discounted density function approach, then we use the COS method to approximate the valuation Equations. When the distribution of the time-until-death random variable is approximated by a combination of exponential distributions and the price of the fund is modeled by an exponential Lévy process, explicit equations for the cosine coefficients are given. Some numerical experiments are also made to illustrate the efficiency of our method.