Got it...
I was making it harder than it was...
I solve for one of the variables, say a...
a=b+2c-d
Then I plug this value into L.
[b+2c-d, b,c,d]
I could set up in aug matrix, but I can pull the constants b, c, d right out to get the kernel basis.
Then my vectors are...