Three variables (infinite number of solutions)

> e:=array(1.. 3,1.. 4,[[4,-3,10,27],[-2,1,-4,-15],[-1,1,-3,-6]]);

> glinsys(e);

Multiply row 1 by 1/4

> e1:=linalg[mulrow](e,1,1/4);

Add to row 2 2 times row 1

> e2:=linalg[addrow](e1,1,2,2);

> glinsys(e2);

Add to row 3 1 times row 1

> e3:=linalg[addrow](e2,1,3,1);

> glinsys(e3);

Multiply row 2 by -2

> e4:=linalg[mulrow](e3,2,-2);

Add to row 1 3/4 times row 2

> e5:=linalg[addrow](e4,2,1,3/4);

> glinsys(e5);

Add to row 3 -1/4 times row 2

> e6:=linalg[addrow](e5,2,3,-1/4);

> glinsys(e6);

> geSlideShow([e,e1,e2,e3,e4,e5,e6]);

>