This paper develops a fast method for the computation of option prices for models whose characteristic function is time-consuming to compute due to the need to solve ordinary differential equations or difference equations numerically, which is the case for a wide class of models of stocks,...