%off nat;
on rounded;


XX:=MAT((X(1)),(X(2)));


*** X declared operator 

      [X(1)]
XX := [    ]
      [X(2)]



S:=mat((a,b),(c,d));


     [A  B]
S := [    ]
     [C  D]


S*XX;


[X(2)*B + X(1)*A]
[               ]
[X(2)*D + X(1)*C]


S*XX - XX;


[X(2)*B + X(1)*A - X(1)]
[                      ]
[X(2)*D - X(2) + X(1)*C]



bye;