Matrix eqn

Added solve_discrete_lyapunov and solve_discrete_riccati.
Cool, thanks. This conversation will disappear when you have pushed the updated file.