Add new command line options: newton solver tol and "damping override"
Don't try to do Mallinson solution for zero damping!

we get a division by zero
Add various bits of code for alternative damping constant recalculation methods
When doing parameter sweeps write command to a script file
Small improvements to python scripts