numerical and analytical jacobian agree
analytical and numerical jacobian nearly agree

with mu for all nuclei variable, including me_n in m_star, still have
incorrect Y_e on convergence