a: matrix([1, 2], [3, 4]);
b: matrix([2, 4], [3, 1]);

a * b;
a / b;
a + b;
a - b;
a^3;
a^b;  /* won't work */
fullmapl("^", a, b);
sin(a);