Let $\Omega\subset \mathbb R^2$ closed and bounded. If we make the substitution $x=g(u,v)$ and $y=h(u,v)$, Why $$\iint_\Omega f(x,y)dxdy=\iint_{\Sigma}f(g(u,v),h(u,v))|J|dudv$$ where $|J|$ is the jacobian of $(u,v)\longmapsto (g(u,v),h(u,v))$ and $\Sigma=\{(u,v)\mid (x(u,v),y(u,v))\in \Omega \}$. It's in fact the Jacobian that I don't understand. Indeed,
$$dx=\frac{\partial g}{\partial u}du+\frac{\partial g}{\partial v}dv$$ $$dy=\frac{\partial h}{\partial u}du+\frac{\partial h}{\partial v}dv$$ and thus, using the convention that $dudu=dvdv=0$ (because it's very small), I get $$dxdy= \frac{\partial g}{\partial u}\frac{\partial h}{\partial v}dudv+\frac{\partial g}{\partial v}\frac{\partial h}{\partial u}dvdu=\det\begin{pmatrix}\frac{\partial g}{\partial u}&-\frac{\partial g}{\partial v}\\\frac{\partial h}{\partial u}&\frac{\partial h}{\partial v}\end{pmatrix}dudv,$$ which is not exactly the Jacobian (but very near to be). Is there any explanation ?