The general matrix for a boost in an arbitrary direction, is given by:
$$\Lambda (\vec{v})= \begin{pmatrix} \gamma & -\gamma\beta_x & -\gamma\beta_y & -\gamma\beta_z\\ -\gamma\beta_x & 1+(1-\gamma)\frac{\beta_x^2}{\beta^2} & (\gamma-1)\frac{\beta_x\beta_y}{\beta^2} & (\gamma-1)\frac{\beta_x\beta_z}{\beta^2}\\ -\gamma\beta_y & (\gamma-1)\frac{\beta_x\beta_y}{\beta^2} & 1+(\gamma-1)\frac{\beta_y^2}{\beta^2} & (\gamma-1)\frac{\beta_y\beta_z}{\beta^2}\\ -\gamma\beta_z & (\gamma-1)\frac{\beta_x\beta_z}{\beta^2} & (\gamma-1)\frac{\beta_y\beta_z}{\beta^2} & 1+ (\gamma-1)\frac{\beta_z^2}{\beta^2} \end{pmatrix} $$
And the boosts in the cartesian directions, are given by
$$ \Lambda_x = \begin{pmatrix} \gamma_x & -\gamma_x\beta_x & 0 & 0\\ -\gamma_x & \gamma_x & 0 & 0 \\ 0 & 0 &1&0\\ 0 & 0& 0 & 1\\ \end{pmatrix}, $$ $$ \Lambda_y = \begin{pmatrix} \gamma_y & 0 & -\gamma_y\beta_y & 0 \\ 0 & 1 & 0 & 0\\ -\gamma_y\beta_y & 0 & \gamma_y & 0 \\ 0 & 0 & 0 &1 \end{pmatrix},$$ $$ \Lambda_z = \begin{pmatrix} \gamma_z & 0 & 0 & -\gamma_z\beta_z \\ 0 & 1 & 0 & 0\\ 0 & 0 & 1 & 0\\ -\gamma_z\beta_z & 0 & 0 & \gamma_z \end{pmatrix}, $$
Where $\gamma_i = \sqrt{1-\beta_i^2}$, and $\beta_i = \frac{v_i}{c} $ for $i=x,y,z$.
From what I have searched, in order to have a continuous transformation, we can obtain the general matrix by applying successive boosts and rotations, but for me is not clear exactly how to combine this to get the general boost matrix. What I have tried is: $\Lambda_xR_x\Lambda_yR_z\Lambda_x$, this is; the first boost in the $x$ direction, then rotate around the $z$-axes counterclockwise by 90 degrees to position the $x$-axis onto the $y$-axis, boost in the $y$ direction, rotate about the $x$-axis by 90 degrees counterclockwise to position the actual $y$-axis on the $z$-axis and finally boost in the $z$-direction. However, this doesn't give the general boost. Since rotations do not commute, I think there must be a concrete order.
On the other hand; I have looked at the way the boost in cartesian $x$-direction is obtained through the Lie Theory of infinitesimal rotations, but I would like to obtain the general matrix with this procedure. If possible, I would appreciate any other method to obtain the general boost in a formal way.