Your $(2, 5, 2, 4)$ is not a solution: it doesn't satisfy the second equation.
Substituting $w$ and $x$ from the first two equations into the last two and simplifying gives you
$$ \eqalign{-{\frac {y \left( 120\,{y}^{2}-34\,yz-17\,{z}^{2} \right) }{ \left( 2\,z+12\,y \right) \left( 3\,z+10\,y \right) }}&=0\cr
{\frac {z \left( 110\,{y}^{2}-15\,yz-6\,{z}^{2} \right) }{ \left( 2\,z
+12\,y \right) \left( 3\,z+10\,y \right) }}&=0\cr}$$
Of course $y$ and $z$ can't be $0$, so that leaves you with
$$ \eqalign{ 120\,{y}^{2}-34\,yz-17\,{z}^{2} &= 0\cr
110\,{y}^{2}-15\,yz-6\,{z}^{2} &= 0\cr}$$
Each of these is the equation of two lines through the origin in the $yz$ plane,
and the slopes are all different, so the intersection is only $(0,0)$, which is not a solution of the original system. Thus your system has no solutions.
There are some problems here. You have computed, correctly,
$$\begin{vmatrix}
\lambda & 1 & 1\\
1& \lambda & 1\\
1 & 1 & \lambda
\end{vmatrix}=\lambda^3-3\lambda + 2.$$
However, this is not $0$ if and only if $\lambda = 0$. Indeed, substituting $\lambda = 0$ into this produces a determinant of $2$.
Further, if the determinant is non-zero, as you say, you'll get a unique solution. This means you can't use $z$ (or any of the variables) as a free parameter. The solution is unique, hence there must be one and only one value of $z$ that produces a solution.
Let's look at the polynomial $\lambda^3-3\lambda + 2$. We can solve this via the usual methods: guess a root, use the factor theorem, then solve the remaining quadratic. If we look at the matrix itself, it jumps out at me that $\lambda = 1$ must be a root, as the corresponding matrix would have all the same rows/columns, making the rows/columns linearly dependent. So, I expect $\lambda = 1$ to be a root (it is, check it!), and $\lambda - 1$ to be a factor. Using polynomial division,
$$\lambda^3 - 3\lambda + 2 = (\lambda - 1)(\lambda^2 + \lambda - 2).$$
Solving/factorising the quadratic by whatever means you prefer, you should then get
$$\lambda^3 - 3\lambda + 2 = (\lambda - 1)^2(\lambda - 2).$$
While this information will be important, it will not, by itself, help us solve the system in place. Cramer's Rule could be used when $\lambda \neq -2, 1$, but just because the system doesn't have a unique solution doesn't mean it can't be solved.
Instead, just set up an augmented matrix, as usual:
$$\left[\begin{array}{ccc|c}
\lambda & 1 & 1 & 1 \\
1 & \lambda & 1 & \mu \\
1 & 1 & \lambda & \mu^2
\end{array}\right].$$
We can now perform all the usual elementary row operations. Remember, $\lambda$ is just a scalar. Though, we must be careful when dividing a row by a scalar expression in terms of $\lambda$, in case that expression happens to be $0$. First, let's swap row $1$ and $3$, to get $1$ in the top left:
$$\left[\begin{array}{ccc|c}
1 & 1 & \lambda & \mu^2 \\
1 & \lambda & 1 & \mu \\
\lambda & 1 & 1 & 1
\end{array}\right].$$
Subtract row $1$ from $2$, and $\lambda$ lots of row $1$ from row $3$:
$$\left[\begin{array}{ccc|c}
1 & 1 & \lambda & \mu^2 \\
0 & \lambda - 1 & 1 - \lambda & \mu - \mu^2 \\
0 & 1 - \lambda & 1 - \lambda^2 & 1 - \lambda\mu^2
\end{array}\right]. \tag{1}$$
The next step I'd want to do is divide a row by $\lambda - 1$ or $1 - \lambda$. But, what if this scalar is $0$, i.e. if $\lambda = 1$? Then the matrix turns into:
$$\left[\begin{array}{ccc|c}
1 & 1 & 1 & \mu^2 \\
0 & 0 & 0 & \mu - \mu^2 \\
0 & 0 & 0 & 1 - \mu^2
\end{array}\right],$$
which can only have a solution if $1 - \mu^2 = 0$ and $\mu - \mu^2 = 0$. The first equation is satisfied if and only if $\mu = \pm 1$, while the second is satisfied if and only if $\mu = 0$ or $1$. Thus, $\mu = 1$ is the only possible value that produces solution. In this case, we get one non-zero row:
$$x + y + z = 1,$$
which we solve by making $y$ and $z$ both free parameters.
Moving back to $(1)$, we can now divide rows $2$ and $3$ by $1 - \lambda$, giving us,
$$\left[\begin{array}{ccc|c}
1 & 1 & \lambda & \mu^2 \\
0 & -1 & 1 & \frac{\mu - \mu^2}{1 - \lambda} \\
0 & 1 & \lambda + 1 & \frac{1 - \lambda\mu^2}{1 - \lambda}
\end{array}\right].$$
Add row $2$ to row $3$, then multiply row $2$ by $-1$ to get
$$\left[\begin{array}{ccc|c}
1 & 1 & \lambda & \mu^2 \\
0 & 1 & -1 & \frac{\mu^2 - \mu}{1 - \lambda} \\
0 & 0 & \lambda + 2 & \frac{1 - \lambda\mu^2 + \mu - \mu^2}{1 - \lambda}
\end{array}\right].$$
Again, we find ourselves wanting to divide by the scalar $\lambda + 2$, but we don't know if it's $0$. Try substituting $\lambda = -2$ into the matrix, and find the two values of $\mu$ that produce solutions. This time, when $\mu$ is one of those two values, there will only be a one-parameter solution.
Otherwise, $\lambda + 2 \neq 0$, and we get
$$\left[\begin{array}{ccc|c}
1 & 1 & \lambda & \mu^2 \\
0 & 1 & -1 & \frac{\mu^2 - \mu}{1 - \lambda} \\
0 & 0 & 1 & \frac{1 - \lambda\mu^2 + \mu - \mu^2}{(1 - \lambda)(\lambda + 2)}
\end{array}\right].$$
This gives us immediately our unique value of $z$ (assuming $\lambda \neq 1, -2$), and the other values can be found through back-substitution.
It's a long question for sure, but it can be done with all the elementary techniques.
Best Answer
Just for simplicity's sake, make a change of variables $$(a,b,c) = (2x,2y,2z)$$ We have the system of equations \begin{align} 16b + 4a^3 &= 3a+17 \\ 8c + 4b^3 &= 3b+9 \\ 4a + 4c^3 &= 3c+5 \\ \end{align} or \begin{align} 16(b-1) &= - 4a^3+3a+1 \\ 8(c-1) &= - 4b^3+3b+1 \\ 4(a-1) &=-4c^3+ 3c+1 \\ \end{align} or \begin{align} 16(b-1) &= - (a-1)(4a^2+4a+1) \tag{1}\\ 8(c-1) &= - (b-1)(4b^2+4b+1) \tag{2}\\ 4(a-1) &=-(c-1)(4c^2+4c+1) \tag{3}\\ \end{align}
Suppose $b >1$ then from $(1)$: $a<1$ $\implies$ from $(3)$: $c>1$ $\implies$ from $(2)$: $b>1$ $\implies$ contracdiction!
Same for $b <1$, we have also a contradiction.
For $b=1$, we have $c =1$ and $a =1$.
Conclusion: the system of equations has a unique solution $(a,b,c) = (1,1,1)$ or $$(x,y,z) = \left(\frac{1}{2},\frac{1}{2},\frac{1}{2} \right)$$