Status: Accepted
Owner: ----
Labels: Type-Defect Priority-Medium
New issue 1323 by Vinzent.Steinberg: integrate(1/sqrt(16 + 4*x**2), x) fails
http://code.google.com/p/sympy/issues/detail?id=1323
>>> integrate(1/sqrt(16 + 4*x**2), x)
Traceback (most recent call last):
File "<stdin>", line 2, in <module>
File "sympy/utilities/decorator.py", line 55, in threaded_decorator
return func(expr, *args, **kwargs)
File "sympy/integrals/integrals.py", line 381, in integrate
return integral.doit()
File "sympy/integrals/integrals.py", line 133, in doit
antideriv = self._eval_integral(function, x)
File "sympy/integrals/integrals.py", line 306, in _eval_integral
h = heurisch(g, x, hints=[])
File "sympy/integrals/risch.py", line 216, in heurisch
denom = reduce(lambda p, q: lcm(p, q, V), denoms)
File "sympy/integrals/risch.py", line 216, in <lambda>
denom = reduce(lambda p, q: lcm(p, q, V), denoms)
File "<string>", line 3, in lcm
File "sympy/polys/wrappers.py", line 23, in _map_basic
result = f(*_conv_args(n, args), **kwargs)
File "sympy/polys/algorithms.py", line 370, in poly_lcm
f, g = f.unify_with(g)
File "sympy/polys/polynomial.py", line 1525, in unify_with
other = cls(other, *symbols, **flags)
File "sympy/polys/polynomial.py", line 402, in __new__
terms = Poly._decompose(poly, *symbols)
File "sympy/polys/polynomial.py", line 545, in _decompose
raise PolynomialError("Can't decompose %s" % factor)
sympy.polys.polynomial.PolynomialError: Can't decompose (4 + _x3**2)**(1/2)
Result should be (according to integrals.wolfram.com):
(1/2) * asinh(x/2)
--
You received this message because you are listed in the owner
or CC fields of this issue, or because you starred this issue.
You may adjust your issue notification preferences at:
http://code.google.com/hosting/settings
--~--~---------~--~----~------------~-------~--~----~
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en
-~----------~----~----~----~------~----~------~--~---
