Numerical integration schemes of time-domain Maxwell equations in optical dispersive and nonlinear media are considered. In this context, we study quantitatively the impact of numerical dispersion on a typical nonlinear parametric conversion process, comparing the widely used finite difference (FDTD) approach and pseudo-spectral (PSTD) methods. Our results show that, unless using very dense grid, only fourth-order PSTD gives results in good quantitative agreement with standard coupled-mode theory