def _eval_as_leading_term(self, x): from diofant import factor_terms expr = self.func(*[t.as_leading_term(x) for t in self.args]).removeO() if not expr: # simple leading term analysis gave us 0 but we have to send # back a term, so compute the leading term (via series) return self.compute_leading_term(x) elif not expr.is_Add: return expr else: plain = expr.func(*[s for s, _ in expr.extract_leading_order(x)]) rv = factor_terms(plain, fraction=False) rv_simplify = rv.simplify() # if it simplifies to an x-free expression, return that; # tests don't fail if we don't but it seems nicer to do this if x not in rv_simplify.free_symbols: if rv_simplify.is_zero and plain.is_zero is not True: return (expr - plain)._eval_as_leading_term(x) return rv_simplify return rv
def simp_hyp(expr): return factor_terms(expand_mul(expr)).rewrite(sin)
def simp_hyp(expr): return factor_terms(expand_mul(expr)).rewrite(sin)