Three variables (no solution)
> d:=array(1.. 3,1.. 4,[[1,1,1,2],[1,2,2,4],[2,3,3,9]]);
> glinsys(d);
Add to row 2 -1 times row 1
> d1:=linalg[addrow](d,1,2,-1);
> glinsys(d1);
Add to row 3 -2 times row 1
> d2:=linalg[addrow](d1,1,3,-2);
> glinsys(d2);
Add to row 1 -1 times row 2
> d3:=linalg[addrow](d2,2,1,-1);
> glinsys(d3);
Add to row 3 -1 times row 2
> d4:=linalg[addrow](d3,2,3,-1);
> glinsys(d4);
> geSlideShow([d,d1,d2,d3]);
>