Matrix eqn

Added solve_discrete_lyapunov and solve_discrete_riccati.
Another possible operator precedence thing. Should we have `G0 = B * (R_hat \ B')`?