
Approximate the volume under a 3D surface z = f(x, y) using a grid of vertical rectangular Riemann prisms. Drag the sample probe or domain boundary handle, switch between Midpoint, Lower Corner, and Upper Corner rules, and increase subdivisions N from 2 to 20 to watch the approximation error shrink toward zero.