>> b=[1 3 5 -4 0;1 3 2 -2 1;1 -2 1 -1 -1;1 -4 1 1 -1];
>> rref(b)
ans =
1.0000 0 0 0 0.5000
0 1.0000 0 0 0.5000
0 0 1.0000 0 0
0 0 0 1.0000 0.5000