We develop a highly accurate numerical method for pricing discrete double barrier options under the Black–Scholes (BS) model. To this aim, the BS partial differential equation is discretized in space by the parabolic finite element method, which is based on a variational formulation and thus...