High Order Difference Methods for Time Dependent PDE