Use RatioInf as initial value during optimization
Add simulateOptimizeWithConstraints and improve zetaOnS to wider range of arguments
Improve mOnS using RatioInf