I have a sort-of sane behaving saha equation function.
Revised my math to take the electron number density as input.