Einheitsmatrix:
#(* Beispiel *) with(linalg): m1=array(identity,1..4,1..4);
Diagonalmatrix:
#(* Beispiel *) with(linalg): m2=diag(a,b,c,d);
Variable Matrix:
#(* Beispiel *) m4=matrix([[a1,b1,c1],[a2,b2,c2],[a3,b3,c3]]);