On Fri, Oct 04, 2013 at 10:23:11AM -0700, Aaron Meurer wrote:
>    In sympy/core/add.py:
>  > @@ -705,13 +705,13 @@ def _eval_as_leading_term(self, x):
>  >          else:
>  >              plain = self.func(*[s for s, _ in self.extract_leading_order(x)])
>  >              rv = factor_terms(plain, fraction=False)
>  > -            rv_fraction = factor_terms(rv, fraction=True)
>  > +            rv_simplify = rv.simplify()
> 
>    Is there a more specific simplification algorithm that could be applied
>    here?

Yes, I think - something trigonometry-related.

Yet, it's clear that you can adopt the failing
example to use bessel functions (for example).