try to base the small perturbation to L in steady on tol instead of float.eps, since the latter constantly gives warnings about poorly conditioned problem.
use warnings.warn instead of print calls