bug in dealing with <fname>_steadystate.m functions that modify parameters
The changes to M_.params introduced in steadystate.m functions are not taken into accounts by calling functions that obtain M as argument and not as global variable.
Need to add M_.params = evalin('base','M_.params'); after each call to _steadystate.m function