No, it will not work if you have a steady state file. Instead, I would change the tolerance parameter of the fsolvecalled on line 224 of prometheo_steadystate.m (see the Mathworks documentation for fsolve) and maybe the maximum number of iterations. Also you do not test the exit value of fsolve. So we do not know if the the solver has converged to a solution. It is good practice to test for the exitflag of fsolve and return from the steady state file with check equal to 1 (if exitflag is not positive).