Number of solutions of a linear equation $AX=B$ I have a question and a proposed solution - Please tell me if I'm correct.
Problem: Prove that if $A$ and $B$ are real matrices and the system of equations $AX=B$ has more than one solution, then it has infinitely many.
Solution: Assume that the system of equations $AX=B$ has more than one solution. This means that the reduced row echelon matrix form of the solution to the equation $AX=B$ has at least one free variable, because if all of the variables were pivot variables, then we would be left with a set of unique values for the variables and hence one solution. Therefore, with at least one free variable, we have that the variable(s) can vary over the real numbers. Hence, there are infinitely many solutions.
Thank you!
 A: Another, more abstract way, to see it is this:  suppose the system has at least two distinct solutions.  Then there exist $X_1$ and $X_2$ such that  $AX_1 = AX_2 = B$ with $X_1 \ne X_2$.   Thus $A(X_1 - X_2) = 0$, while $X_1 - X_2 \ne 0$.  Let $C$ be the set of non-vanishing columns of $X_1 - X_2$.  $C \ne \phi$ since $X_1 - X_2 \ne 0$.  Let $D$ be any matrix formed by selecting its columns from $span(C)$.  Then $AD = 0$, so $A(X_1 + D) = AX_1 + AD = AX_1 = B$; but clearly the number of such matrices $D$ is infinite, whence the number of matrices $X_1 + D$ is infinite as well.  QED and Cheers.
A: The "infinitely many" comes from the fact that there are infinitely many scalars (presumably, as you are probably working over the real or complex numbers; if you consider the possibility of working over a finite field, then the result is obviously not true). You can see this without playing with columns and echelon forms, as follows.
Let $x_0,x_1$ be two different solutions: $Ax_0=B=Ax_1$. Then putting $x_\lambda=\lambda x_1+(1-\lambda)x_0$ one has by linearity $Ax_\lambda=\lambda Ax_1+(1-\lambda)Ax_0=\lambda B+(1-\lambda)B=B$, so $x_\lambda$ is a solution as well, and there are infinitely many choices for $\lambda$. And these are all different, since $x_\lambda=x_\mu$ implies $0=x_\lambda-x_\mu=(\lambda-\mu)(x_1-x_0)$, which given that $x_0\neq x_1$ can only hold if $\lambda-\mu$.
