MAINT: enable np.matrix inputs to solve_discrete_lyapunov

closes https://github.com/scipy/scipy/issues/4176