Implementation of numerical Fourier method for second order Taylor schemes