For Lebesgue integrals, I know how to approximate $\int_{A} f(x)dx$ for some measurable set $A\subset \mathbb{R}^{n}$, I just decompose $A$ into measurable subsets $A_{i}$ and then make some finite "Riemann sum" $\sum f(x_{i})\mu (A_{i})$ where $x_{i}$ is a well-chosen point in $A_i$.
For integration over forms, I guess that I cannot anymore do this, since now I need to deal with orientation and stuff. Basically, I cannot just treat $dx\wedge dy \wedge dz$ as $dxdydz$. How do I proceed then to obtain a numerical approximation of $\int_{A} f(x,y,z) dx\wedge dy \wedge dz$ for some $A\subset \mathbb{R}^{3}$ ?