Use master function lyapunov_solver.m to call individual solvers
Closes #821 (closed) by switching to Schur-based method by default when iteration fails.
Closes #1077 (closed) by moving call to Matlab's dlyapchol
and to disclyap_fast
to master function lyapunov_solver