PyPI does not support this highlighting
Solvers: rename the gamma step size to step
Avoid complex casting to real when yn starts real and becomes complex