updated Saha equation code, fixed numerical errors