M={{1,3,1,1,0,0},{2,1,3,0,1,0},{3,4,5,0,0,1}};
Print["М=",MatrixForm[M]];
M1={{0,0,0,0,0,0},{0,0,0,0,0,0},{0,0,0,0,0,0}};
M2={{0,0,0,0,0,0},{0,0,0,0,0,0},{0,0,0,0,0,0}};
M3={{0,0,0,0,0,0},{0,0,0,0,0,0},{0,0,0,0,0,0}};
M4={{0,0,0,0,0,0},{0,0,0,0,0,0},{0,0,0,0,0,0}};
M5={{0,0,0,0,0,0},{0,0,0,0,0,0},{0,0,0,0,0,0}};
M6={{0,0,0,0,0,0},{0,0,0,0,0,0},{0,0,0,0,0,0}};
(Subscript[M1, [[1]]])=(Subscript[M, [[1]]]);
(Subscript[M1, [[2]]])=-((Subscript[M, [[2,1]]])/(Subscript[M, [[1,1]]]))*(Subscript[M, [[1]]])+(Subscript[M, [[2]]]);
(Subscript[M1, [[3]]])=-((Subscript[M, [[3,1]]])/(Subscript[M, [[1,1]]]))*(Subscript[M, [[1]]])+(Subscript[M, [[3]]]);
(Subscript[M2, [[1]]])=(Subscript[M1, [[1]]]);
(Subscript[M2, [[2]]])=(Subscript[M1, [[2]]])/(Subscript[M, [[2,2]]]);
(Subscript[M2, [[3]]])=(Subscript[M1, [[3]]]);
(Subscript[M3, [[1]]])=(Subscript[M2, [[1]]]);
(Subscript[M3, [[2]]])=(Subscript[M2, [[2]]])/(Subscript[M2, [[2,2]]]);
(Subscript[M3, [[3]]])=(Subscript[M2, [[3]]]);
(Subscript[M4, [[1]]])=(Subscript[M3, [[1]]]);
(Subscript[M4, [[2]]])=(Subscript[M3, [[2]]]);
(Subscript[M4, [[3]]])=-(Subscript[M3, [[3,2]]])*(Subscript[M3, [[2]]])+(Subscript[M3, [[3]]]);
(Subscript[M5, [[1]]])=(Subscript[M4, [[1]]])-(Subscript[M4, [[3]]]);
(Subscript[M5, [[2]]])=-(Subscript[M4, [[2,3]]])*(Subscript[M4, [[3]]])+(Subscript[M4, [[2]]]);
(Subscript[M5, [[3]]])=(Subscript[M4, [[3]]]);
(Subscript[M6, [[1]]])=-(Subscript[M5, [[1,2]]])*(Subscript[M5, [[2]]])+(Subscript[M5, [[1]]]);
(Subscript[M6, [[2]]])=(Subscript[M5, [[2]]]);
(Subscript[M6, [[3]]])=(Subscript[M5, [[3]]]);
Print["М1=",MatrixForm[M1]];
Print["М2=",MatrixForm[M2]];
Print["М3=",MatrixForm[M3]];
Print["М4=",MatrixForm[M4]];
Print["М5=",MatrixForm[M5]];
Print["М6=",MatrixForm[M6]];