Numerical methods for large scale matrix equations with applications in LTI system model reduction