Higher Order Spectral method