Several bug fixes

1. In negbin_ll() and gen_yule_ll(), implement built-in log functions
for each term to prevent -inf.
2. Rewrite gen_yule_ll() to be consistent with notation in
gen_yule_solver().
3. In gen_yule_solver(), implement maximum number of iterations to
prevent the optimization process from running out of hand.
Fix bug in iteration