The integration by substitution technique is dervied from the following statement: $$\int _{a}^{b}f(\varphi (x))\varphi '(x)\,dx=\int _{\varphi (a)}^{\varphi (b)}f(u)\,du$$ Now almost all the explanations I saw about the actual technique start playing around with the Leibniz notation and $dx$, $du$ to get to the new form, e.g stating $du=\varphi '(x)\,dx$ from ${\frac {du}{dx}}=\varphi '(x)$.
Now the problem is that $dx$ and $du$ are said to have no actual meaning, at least in elementary calculus, and yet we're playing around with them as if they were actual variables forming actual quotients. As far as I am aware this wouldn't be a problem in nonstandard analysis, but how should one make sense of it in the elementary settings? How can we apply this technique without "heuristic" methods (which actually seems to be used almost all the time for different techniques)?
Note: My problem is actually applying the method, not the proof, as there are many proofs available not using the Leibniz's notation.