# linear matrix equation

I want to solve the following equation for a matrix solution X:

'A*X*B == C'
# in general, they are rectangular with dimensions:
# A m-by-n, B p-by-q, and hence, X n-by-p, C m-by-q
# for simplicity, assume: they are all n-by-n square matrices


Are there any tools in SAGE (or Numpy/Scipy, Fortran, ...) for solving equations of this type?