Approximating a double integral with prisms