Fix root of negative values in steadyState
Set up the initial conditions correctly in parallel