Let's look at this in two dimensions first. After this you should know how to do it in three!
Suppose that you are working in the plane and have two linearly independent vectors $v$ and $w$. You want to make $v$ and $w$ orthogonal to each other in terms of the standard euclidean inner product. How can you do it?
Well notice that we can subtract a certain multiple of $w$ from $v$, let's call it $cw$ where $c$ is some constant so that $v - cw$ will be orthogonal to $w$. In other words, after some algebraic manipulation we find that $c$ must be equal to
$$\frac{\langle v,w \rangle }{\langle w,w \rangle}$$
simply by solving the equation $\langle v - cw, w \rangle = 0$ for $c$. Then now you will have an orthogonal basis for $\mathbb{R}^2$, namely the vectors
$$w \quad \text{and} \quad v - \frac{\langle v,w \rangle }{\langle w,w \rangle} w.$$
To find an orthonormal basis, you just need to divide through by the length of each of the vectors.
In $\mathbb{R}^3$ you just need to apply this process recursively as shown in the wikipedia link in the comments above. However you first need to check that your vectors are linearly independent! You can check this by calculating the determinant of the matrix whose columns are the vectors that you have stated in your question.
\begin{pmatrix} ... \end{pmatrix}
orbmatrix
instead ofmatrix
to get the parentheses or brackets. $\endgroup$