Not to mention complex numerical optimization code that mixes closed-form approximations and something like Newton.

Without guides as to why a particular hairy expression is a good idea as a first estimate, the code is pretty much unreadable. (E.g. is it setting derivatives to zero, using a polynomial approximation, or something else?)