simplify(expression,{relations}) works for non-function variables so I suppose I could write out each function instance explicitly.

uh, you didn't really answer the question. using linalg, define two matrices using "matrix". How do you take the tensor product of these in maple 9? I've looked online and through help but am not finding it as quick as I should be as usual. Maybe Maple 10 is better.