Most of file has no units for given parameters. For example see steady_state_PFR.m line 34. It is also recommended to give references for problem solved here. Thank you for this great library