//SINGULAR Example 2.1.6



ring A = 0,(x,y,z),dp;
matrix M[2][3] = 1, x+y, z2,        //2x3 matrix
                 x,   0, xyz;
matrix N[3][3] = 1,2,3,4,5,6,7,8,9; //3x3 matrix

M;                             //lists all entries of M
print(N);                      //displays N as usual
print(M+M);                    //addition of matrices
print(x*N);
print(M*N);                    //multiplication of matrices
M[2,3];                        //access to single entry
M[2,3]=37;                     //change single entry
print(M);


LIB "matrix.lib"; LIB "inout.lib";

print(power(N,3));       //exponentiation of matrices
pmat(power((x+y+z)*N,3),15);//show first 15 terms of entries
matrix K = concat(M,N);  //concatenation
print(K);
ideal(M);                //converts matrix to ideal
print(unitmat(5));       //5x5 unit matrix
intmat I[2][3]=1,2,3,4,5,6;
I;