sylv — Sylvester equation.
sylv(A,B,C,flag)
X = sylv(A,B,C,'c') computes X, solution
of the "continuous time" Sylvester equation
A*X+X*B=C
X=sylv(A,B,C,'d') computes X, solution
of the "discrete time" Sylvester equation
A*X*B-X=C