Du kannst hier ja die ersten (n-1) Variablen beliebig wählen a1, a2 etc und dann für
xn = -a1 -a2 -.... -an-1 einsetzen
Der Lösungsvektor x sieht dann so aus
x = ( a1 , a2 , ...., an-1 , -a1 -a2 -.... -an-1 )
oder eben so
x = a1*(1,0,0,...,0,-1) + a2*(0,1,0,0,...,0,-1) + ....+ an-1*(0,0,0,...,1,-1).
Und die darin vorkommenden Vektoren erzeugen den ganzen Lösungsraum
und sind lin. unabh., also bilden sie eine Basis.