Matrix eqn

Added solve_discrete_lyapunov and solve_discrete_riccati.
use `error = Base.maxabs(H1 - H0)`

it is more efficient