How does one prove the midpoint rule for triple integral?
That is,
$$\sum_i \sum_j \sum_k f(\hat{x_i}, \hat{y_j}, \hat{z_k}) \Delta x_i \Delta y_j \Delta z_k$$
where $\displaystyle (\hat{x_i}, \hat{y_i}, \hat{z_i})=\bigg( \frac{x_{i-1}+x_i}{2}, \frac{y_{j-1}+y_j}{2}, \frac{z_{k-1}+z_k}{2}\bigg)$
It seems that this formula relies on that "every midpoint of every subinterval is visited".
However, I wonder if this is some kind of permutation algebra proof or whether there's some other reasoning to it?
Or does one somehow know that to get the midpoints of all of some $\mathbb{R}^n$, then midpoints of each individual $\mathbb{R}$ are composed through a product rule?