Many of the different numerical techniques in the partial differential equations framework for solving option pricing problems have employed only standard second-order discretization schemes. A higher-order discretization has the advantage of producing low size matrix systems for computing...