Numerical procedure to integrate non-linear partial differential equations where the dependent variables are expanded in Fourier series. The spectral simulations are particularly suited if extremely high accuracy is required as in laminar to turbulent transition simulations.