In this paper, we present a computationally efficient technique based on the Method of Lines for the approximation of the Bermudan option values via the associated partial differential equations. The method of lines converts the Black Scholes partial differential equation to a system of ordinary differential equations. The solution of the system of ordinary differential equations so obtained only requires spatial discretization and avoids discretization in time. Additionally, the exact solution of the ordinary differential equations can be obtained efficiently using the exponential matrix operation, making the method computationally attractive and straightforward to implement. An essential advantage of the proposed approach is that the associated Greeks can be computed with minimal additional computations. We illustrate, through numerical experiments, the efficacy of the proposed method in pricing and computation of the sensitivities for a European call, cash-or-nothing, powered option, and Bermudan put option.