$$

\left\{

\begin{array}{l}

x + y + az + at = a \\

2x + y – a^2z + 2at = 1 \\

4x + 3y + a^2 z + 4a^2 t = 1

\end{array}

\right.

$$

Is there a way to predict the number of solutions $(x, y, z, t)$, depending on the value of the parameter $a$, without solving the system explicitely?

## Best Answer

Yes, by reducing the augmented matrix of the system to row echelon form. For your system, the augmented matrix is $$ \left[ \begin{array}{cccc|c} 1 & 1 & a & a & a \\ 2 & 1 & -a^2 & 2a & 1 \\ 4 & 3 & a^2 & 4a^2 & 1 \end{array}\right] $$ and the row echelon form $$ \left[ \begin{array}{cccc|c} 1 & 1 & a & a & a \\ 0 & 1 & a^2+2a & 0 & 2a-1 \\ 0 & 0 & 2(a^2-a) & 4(a^2-a) & -2a \end{array}\right]. $$ First thing to look for is consistency, so that the last row does not correspond to $0=$non-zero. This guarantees existence of at least one solution. After that, the dimension of the solution space = the number of free variables = $n - \rm{rank}(A)$, where $n$ is the number of variables of the system, $A$ is the augmented matrix, whose rank is given by the number of non-zero rows in the row echelon form. (Of course if you've calculated the row echelon form, you're only 30 seconds from the explicit solution. If you need to know the dimension of the solution space, I don't think there's any alternative: you need the rank of $A$.)