riemann_prisms3d function

Riemann rectangular prisms over a planar region