\(\begin{align*} \vec{a}\cdot (x,y,z) = 0 \implies a_1x+a_2y+a_3z =0 &&(1)\end{align*}\)
\(\begin{align*} \vec{b}\cdot (x,y,z) = 0 \implies b_1x+b_2y+b_3z =0 && (2)\end{align*}\)
\(\begin{align*}\dfrac{x}{(a_2b_3 - a_3b_2)} &= \dfrac{y}{(a_3b_1 - a_1b_3)}\end{align*}\)
If we repeat the same procedure, instead multiplying equation \((1)\) by \(b_1\) and equation \((2)\) by \(a_1\), we may eliminate \(x\) to get
\[\dfrac{y}{(a_3b_1 - a_1b_3)} = \dfrac{z}{(a_1b_2 - a_2b_1)} \]
Therefore, we have
\[\frac{x}{(a_2b_3 - a_3b_2)}=\frac{y}{(a_3b_1 - a_1b_3)} = \frac{z}{(a_1b_2 - a_2b_1)}=c\]
for some \(c\in\mathbb{R}\).
Finally, we see that
\[\begin{align*}(x, y, z) &= (c(a_2b_3 - a_3b_2), c(a_3b_1 - a_1b_3), c(a_1b_2 - a_2b_1)) \\ &= c(a_2b_3 - a_3b_2, a_3b_1 - a_1b_3, a_1b_2 - a_2b_1)\end{align*}\]
Since any scalar multiple of \((x,y,z)\) will be orthogonal to both \(\vec{a}\) and \(\vec{b}\) (you can verify this by taking the dot product), we may simplify the vector above by letting \(c=1\).
In fact, by definition, the cross product of \(\vec{a}=(a_1,a_2,a_3)\) and \(\vec{b}=(b_1,b_2,b_3)\) is exactly the vector given above when \(c=1\).