No need to calculate the reduced row echolon form manually, you may use rref.
But as rank(M)=rank(ME)=2<3 we already know, that we have an infinite number of solutions. You may google "Rouché-Capelli theorem".
Strange enough symbolic evaluation will yield wrong result!!?? Probably the result of numerical inaccuracy due to the awful automatic switch to float mode.
