obsv_mat — observability matrix
[O]=obsv_mat(A,C) [O]=obsv_mat(sl)
obsv_mat returns the observability matrix:
O=[C; CA; CA^2;...; CA^(n-1) ]
