A Higher Order Finite Difference Method in Reservoir Stimulation