Die Bedingung für U ist ja:
x1+x2-x3+x4=0
Also kannst du 3 Variablen frei wählen, etwa
x1=r und x2=s und x3=t und musst dann nur für
x4= r+s-t nehmen.
Dann sind alle Lösungen von der Art
(r;s;t;r+s-t) =
r*(1;0;0;1) + s*(0;1;0;1)+t*(0;0;1;-1)
also wäre eine Basis für U
(1;0;0;1) , (0;1;0;1) , (0;0;1;-1)
Die kannst du ja nach Gram-Schmidt zu einer
Orthonormalbasis umbauen.