Commutation matrix (memoised)
commutation_matrix0(...)
...
'r' integer; row dimension.
'c' integer; column dimension.
Use empty argument to clear the cache.