That's a bug. integrate() shouldn't raise an exception. Aaron Meurer
On Wed, Feb 24, 2016 at 7:00 AM, Abhishek Verma <[email protected]> wrote: > Hi everyone, > I am Abhishek Verma, and I will be applying for GSoC this year. > I have trouble while calculating the Indefinite Integration. > > I have a Expr=X^(log(x^log(x))) ,while Solving this by integrate() function > i get the result - >>>> integrate(x**log(x**log(x)),x) > ⌠ > ⎮ 3 > ⎮ log (x) > ⎮ ℯ dx > ⌡ > > > But for Expr=X^(log(x^log(x^log(x)))) > > > >>>> integrate(x**log(x**log(x**log(x))),x) > Traceback (most recent call last): > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 365, in from_expr > poly = self._rebuild_expr(expr, mapping) > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 359, in _rebuild_expr > return _rebuild(sympify(expr)) > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 351, in _rebuild > return reduce(add, list(map(_rebuild, expr.args))) > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 353, in _rebuild > return reduce(mul, list(map(_rebuild, expr.args))) > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 351, in _rebuild > return reduce(add, list(map(_rebuild, expr.args))) > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 353, in _rebuild > return reduce(mul, list(map(_rebuild, expr.args))) > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 351, in _rebuild > return reduce(add, list(map(_rebuild, expr.args))) > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 353, in _rebuild > return reduce(mul, list(map(_rebuild, expr.args))) > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 357, in _rebuild > return domain.convert(expr) > File > "/usr/local/lib/python3.4/dist-packages/sympy/polys/domains/domain.py", line > 146, in convert > raise CoercionFailed("can't convert %s of type %s to %s" % (element, > type(element), self)) > sympy.polys.polyerrors.CoercionFailed: can't convert _x0**_x1 of type <class > 'sympy.core.power.Pow'> to > QQ[_A0,_A1,_A2,_A3,_A4,_A5,_A6,_A7,_A8,_A9,_A10,_A11,_A12,_A13,_A14,_A15,_A16,_A17,_A18,_A19,_A20,_A21,_A22,_A23,_A24,_A25,_A26,_A27,_A28,_A29,_A30,_A31,_A32,_A33,_A34,_A35,_A36,_A37,_A38,_A39,_A40,_A41,_A42,_A43,_A44,_A45,_A46,_A47,_A48,_A49,_A50,_A51,_A52,_A53,_A54,_A55,_A56,_A57,_A58,_A59,_A60,_A61,_A62,_A63,_A64,_A65,_A66,_A67,_A68,_A69,_A70,_A71,_A72,_A73,_A74,_A75,_A76,_A77,_A78,_A79,_A80,_A81,_A82,_A83,_A84,_A85,_A86,_A87,_A88,_A89,_A90,_A91,_A92,_A93,_A94,_A95,_A96,_A97,_A98,_A99,_A100,_A101,_A102,_A103,_A104,_A105,_A106,_A107,_A108,_A109,_A110,_A111,_A112,_A113,_A114,_A115,_A116,_A117,_A118,_A119,_B0] > > During handling of the above exception, another exception occurred: > > Traceback (most recent call last): > File "<stdin>", line 1, in <module> > File > "/usr/local/lib/python3.4/dist-packages/sympy/utilities/decorator.py", line > 35, in threaded_func > return func(expr, *args, **kwargs) > File > "/usr/local/lib/python3.4/dist-packages/sympy/integrals/integrals.py", line > 1232, in integrate > risch=risch, manual=manual) > File > "/usr/local/lib/python3.4/dist-packages/sympy/integrals/integrals.py", line > 487, in doit > conds=conds) > File > "/usr/local/lib/python3.4/dist-packages/sympy/integrals/integrals.py", line > 862, in _eval_integral > h = heurisch_wrapper(g, x, hints=[]) > File "/usr/local/lib/python3.4/dist-packages/sympy/integrals/heurisch.py", > line 128, in heurisch_wrapper > unnecessary_permutations) > File "/usr/local/lib/python3.4/dist-packages/sympy/integrals/heurisch.py", > line 566, in heurisch > solution = _integrate('Q') > File "/usr/local/lib/python3.4/dist-packages/sympy/integrals/heurisch.py", > line 555, in _integrate > numer = ring.from_expr(raw_numer) > File "/usr/local/lib/python3.4/dist-packages/sympy/polys/rings.py", line > 367, in from_expr > raise ValueError("expected an expression convertible to a polynomial in > %s, got %s" % (self, expr)) > ValueError: expected an expression convertible to a polynomial in Polynomial > ring in _x0, _x1, _x2, _x3, _x4, _x5, _x6 over > QQ[_A0,_A1,_A2,_A3,_A4,_A5,_A6,_A7,_A8,_A9,_A10,_A11,_A12,_A13,_A14,_A15,_A16,_A17,_A18,_A19,_A20,_A21,_A22,_A23,_A24,_A25,_A26,_A27,_A28,_A29,_A30,_A31,_A32,_A33,_A34,_A35,_A36,_A37,_A38,_A39,_A40,_A41,_A42,_A43,_A44,_A45,_A46,_A47,_A48,_A49,_A50,_A51,_A52,_A53,_A54,_A55,_A56,_A57,_A58,_A59,_A60,_A61,_A62,_A63,_A64,_A65,_A66,_A67,_A68,_A69,_A70,_A71,_A72,_A73,_A74,_A75,_A76,_A77,_A78,_A79,_A80,_A81,_A82,_A83,_A84,_A85,_A86,_A87,_A88,_A89,_A90,_A91,_A92,_A93,_A94,_A95,_A96,_A97,_A98,_A99,_A100,_A101,_A102,_A103,_A104,_A105,_A106,_A107,_A108,_A109,_A110,_A111,_A112,_A113,_A114,_A115,_A116,_A117,_A118,_A119,_B0] > with lex order, got _x0**5*_x5 - _x0**3*(2*_x0**_x2*_x2*(2*_A103*_x2*_x6 + > _A106*_x4*_x5 + _x4**2*_A107 + _A108*_x1*_x2 + _A110*_x4 + _A116*_x0*_x1 + > 2*_A12*_x1*_x6 + 2*_A19*_x4*_x6 + _A2 + _A20*_x2*_x3 + _A21*_x1*_x3 + > _A23*_x3*_x4 + _A26*_x0*_x3 + _A3*_x2 + _x5**2*_A32 + _A35*_x0*_x2 + > _A39*_x1 + _x2**2*_A41 + _A50*_x3 + 3*_x6**2*_A52 + 2*_A58*_x6 + _A6*_x0*_x4 > + _A62*_x0*_x5 + _x1**2*_A69 + _A72*_x2*_x5 + _A74*_x2*_x4 + 2*_A77*_x0*_x6 > + _x0**2*_A78 + _A8*_x5 + _x3**2*_A81 + 2*_A83*_x3*_x6 + _A89*_x0 + > 2*_A90*_x5*_x6 + _A95*_x1*_x4 + _A96*_x1*_x5 + _A99*_x3*_x5) + _A101*_x0*_x1 > + _A102*_x0*_x4 + _x6**2*_A103 + _A105*_x0*_x3 + _A108*_x1*_x6 + > _A109*_x0*_x5 + _A11*_x3*_x5 + _A14*_x3 + 2*_A15*_x2*_x3 + _A18*_x4 + > _A20*_x3*_x6 + 2*_A24*_x1*_x2 + 2*_A27*_x2*_x5 + _A3*_x6 + _A31*_x1*_x4 + > _x5**2*_A33 + _A35*_x0*_x6 + 2*_A36*_x2 + 3*_x2**2*_A38 + _A40*_x3*_x4 + > 2*_A41*_x2*_x6 + _A42*_x1*_x5 + _A61*_x1 + _A70*_x0 + _x0**2*_A71 + > _A72*_x5*_x6 + _A74*_x4*_x6 + _A75*_x5 + 2*_A76*_x0*_x2 + 2*_A80*_x2*_x4 + > _x1**2*_A82 + _A86*_x4*_x5 + _A87*_x1*_x3 + _x4**2*_A9 + _x3**2*_A93 + _A94 > + 2*_x2*(2*_A1*_x1*_x5 + _A10 + _A101*_x0*_x2 + _A108*_x2*_x6 + _x5**2*_A112 > + _A113*_x3*_x4 + _A116*_x0*_x6 + 2*_A117*_x0*_x1 + _x6**2*_A12 + 2*_A13*_x1 > + _A21*_x3*_x6 + _x4**2*_A22 + _x2**2*_A24 + _A31*_x2*_x4 + _A34*_x0*_x5 + > _A37*_x3 + _A39*_x6 + _A4*_x3*_x5 + _A42*_x2*_x5 + _A45*_x5 + _A47*_x0*_x4 + > 3*_x1**2*_A5 + _A51*_x0*_x3 + _x0**2*_A53 + _A56*_x0 + 2*_A59*_x1*_x4 + > _A60*_x4*_x5 + _A61*_x2 + _x3**2*_A64 + 2*_A66*_x1*_x3 + 2*_A69*_x1*_x6 + > 2*_A82*_x1*_x2 + _A87*_x2*_x3 + _A95*_x4*_x6 + _A96*_x5*_x6 + _A98*_x4) + > (_x1 + 2*_x2**2)*(_A100*_x0*_x3 + _A102*_x0*_x2 + _A106*_x5*_x6 + > 2*_A107*_x4*_x6 + _A110*_x6 + _A113*_x1*_x3 + _x5**2*_A115 + _A16*_x0*_x5 + > 2*_A17*_x3*_x4 + _A18*_x2 + _x6**2*_A19 + 2*_A22*_x1*_x4 + _A23*_x3*_x6 + > _A31*_x1*_x2 + _A40*_x2*_x3 + 2*_A46*_x0*_x4 + _A47*_x0*_x1 + _x0**2*_A48 + > _x3**2*_A55 + _x1**2*_A59 + _A6*_x0*_x6 + _A60*_x1*_x5 + _A63 + _A65*_x0 + > _A7*_x3*_x5 + _A74*_x2*_x6 + 2*_A79*_x4*_x5 + _x2**2*_A80 + _A84*_x3 + > _A85*_x5 + _A86*_x2*_x5 + 2*_A88*_x4 + 2*_A9*_x2*_x4 + 3*_x4**2*_A91 + > _A95*_x1*_x6 + _A98*_x1) + (_x0**_x1*_x1 + 2*_x0**_x1*_x2**2)*(_A100*_x0*_x4 > + 2*_A104*_x3*_x5 + _A105*_x0*_x2 + _A11*_x2*_x5 + _A113*_x1*_x4 + > 2*_A118*_x0*_x3 + _A119*_x0*_x5 + _A14*_x2 + _x2**2*_A15 + _x4**2*_A17 + > _A20*_x2*_x6 + _A21*_x1*_x6 + _A23*_x4*_x6 + _A26*_x0*_x6 + _A29 + _A30*_x5 > + _A37*_x1 + _A4*_x1*_x5 + _A40*_x2*_x4 + _x0**2*_A43 + _A50*_x6 + > _A51*_x0*_x1 + _x5**2*_A54 + 2*_A55*_x3*_x4 + _A57*_x0 + 2*_A64*_x1*_x3 + > _x1**2*_A66 + _A7*_x4*_x5 + 2*_A81*_x3*_x6 + _x6**2*_A83 + _A84*_x4 + > _A87*_x1*_x2 + 2*_A92*_x3 + 2*_A93*_x2*_x3 + 3*_x3**2*_A97 + _A99*_x5*_x6) + > (_x1*_x2*_x5 + 2*_x2**3*_x5 + _x4*_x5)*(_x1**2*_A1 + _x3**2*_A104 + > _A106*_x4*_x6 + _A109*_x0*_x2 + _A11*_x2*_x3 + _A111*_x0 + 2*_A112*_x1*_x5 + > _x0**2*_A114 + 2*_A115*_x4*_x5 + _A119*_x0*_x3 + _A16*_x0*_x4 + _x2**2*_A27 > + _A30*_x3 + 2*_A32*_x5*_x6 + 2*_A33*_x2*_x5 + _A34*_x0*_x1 + _A4*_x1*_x3 + > _A42*_x1*_x2 + 3*_x5**2*_A44 + _A45*_x1 + 2*_A49*_x5 + 2*_A54*_x3*_x5 + > _A60*_x1*_x4 + _A62*_x0*_x6 + 2*_A68*_x0*_x5 + _A7*_x3*_x4 + _A72*_x2*_x6 + > _A73 + _A75*_x2 + _x4**2*_A79 + _A8*_x6 + _A85*_x4 + _A86*_x2*_x4 + > _x6**2*_A90 + _A96*_x1*_x6 + _A99*_x3*_x6)) - _x0**2*(_x0**2*(_A100*_x3*_x4 > + _A101*_x1*_x2 + _A102*_x2*_x4 + _A105*_x2*_x3 + _A109*_x2*_x5 + _A111*_x5 > + 2*_A114*_x0*_x5 + _A116*_x1*_x6 + _x1**2*_A117 + _x3**2*_A118 + > _A119*_x3*_x5 + _A16*_x4*_x5 + 3*_x0**2*_A25 + _A26*_x3*_x6 + 2*_A28*_x0 + > _A34*_x1*_x5 + _A35*_x2*_x6 + 2*_A43*_x0*_x3 + _x4**2*_A46 + _A47*_x1*_x4 + > 2*_A48*_x0*_x4 + _A51*_x1*_x3 + 2*_A53*_x0*_x1 + _A56*_x1 + _A57*_x3 + > _A6*_x4*_x6 + _A62*_x5*_x6 + _A65*_x4 + _A67 + _x5**2*_A68 + _A70*_x2 + > 2*_A71*_x0*_x2 + _x2**2*_A76 + _x6**2*_A77 + 2*_A78*_x0*_x6 + _A89*_x6 + > _B0) + _x0*(-_A0 - _x1**2*_A1*_x5 - _A10*_x1 - _A100*_x0*_x3*_x4 - > _A101*_x0*_x1*_x2 - _A102*_x0*_x2*_x4 - _x6**2*_A103*_x2 - _x3**2*_A104*_x5 > - _A105*_x0*_x2*_x3 - _A106*_x4*_x5*_x6 - _x4**2*_A107*_x6 - > _A108*_x1*_x2*_x6 - _A109*_x0*_x2*_x5 - _A11*_x2*_x3*_x5 - _A110*_x4*_x6 - > _A111*_x0*_x5 - _x5**2*_A112*_x1 - _A113*_x1*_x3*_x4 - _x0**2*_A114*_x5 - > _x5**2*_A115*_x4 - _A116*_x0*_x1*_x6 - _x1**2*_A117*_x0 - _x3**2*_A118*_x0 - > _A119*_x0*_x3*_x5 - _x6**2*_A12*_x1 - _x1**2*_A13 - _A14*_x2*_x3 - > _x2**2*_A15*_x3 - _A16*_x0*_x4*_x5 - _x4**2*_A17*_x3 - _A18*_x2*_x4 - > _x6**2*_A19*_x4 - _A2*_x6 - _A20*_x2*_x3*_x6 - _A21*_x1*_x3*_x6 - > _x4**2*_A22*_x1 - _A23*_x3*_x4*_x6 - _x2**2*_A24*_x1 - _x0**3*_A25 - > _A26*_x0*_x3*_x6 - _x2**2*_A27*_x5 - _x0**2*_A28 - _A29*_x3 - _A3*_x2*_x6 - > _A30*_x3*_x5 - _A31*_x1*_x2*_x4 - _x5**2*_A32*_x6 - _x5**2*_A33*_x2 - > _A34*_x0*_x1*_x5 - _A35*_x0*_x2*_x6 - _x2**2*_A36 - _A37*_x1*_x3 - > _x2**3*_A38 - _A39*_x1*_x6 - _A4*_x1*_x3*_x5 - _A40*_x2*_x3*_x4 - > _x2**2*_A41*_x6 - _A42*_x1*_x2*_x5 - _x0**2*_A43*_x3 - _x5**3*_A44 - > _A45*_x1*_x5 - _x4**2*_A46*_x0 - _A47*_x0*_x1*_x4 - _x0**2*_A48*_x4 - > _x5**2*_A49 - _x1**3*_A5 - _A50*_x3*_x6 - _A51*_x0*_x1*_x3 - _x6**3*_A52 - > _x0**2*_A53*_x1 - _x5**2*_A54*_x3 - _x3**2*_A55*_x4 - _A56*_x0*_x1 - > _A57*_x0*_x3 - _x6**2*_A58 - _x1**2*_A59*_x4 - _A6*_x0*_x4*_x6 - > _A60*_x1*_x4*_x5 - _A61*_x1*_x2 - _A62*_x0*_x5*_x6 - _A63*_x4 - > _x3**2*_A64*_x1 - _A65*_x0*_x4 - _x1**2*_A66*_x3 - _A67*_x0 - > _x5**2*_A68*_x0 - _x1**2*_A69*_x6 - _A7*_x3*_x4*_x5 - _A70*_x0*_x2 - > _x0**2*_A71*_x2 - _A72*_x2*_x5*_x6 - _A73*_x5 - _A74*_x2*_x4*_x6 - > _A75*_x2*_x5 - _x2**2*_A76*_x0 - _x6**2*_A77*_x0 - _x0**2*_A78*_x6 - > _x4**2*_A79*_x5 - _A8*_x5*_x6 - _x2**2*_A80*_x4 - _x3**2*_A81*_x6 - > _x1**2*_A82*_x2 - _x6**2*_A83*_x3 - _A84*_x3*_x4 - _A85*_x4*_x5 - > _A86*_x2*_x4*_x5 - _A87*_x1*_x2*_x3 - _x4**2*_A88 - _A89*_x0*_x6 - > _x4**2*_A9*_x2 - _x6**2*_A90*_x5 - _x4**3*_A91 - _x3**2*_A92 - > _x3**2*_A93*_x2 - _A94*_x2 - _A95*_x1*_x4*_x6 - _A96*_x1*_x5*_x6 - > _x3**3*_A97 - _A98*_x1*_x4 - _A99*_x3*_x5*_x6)) > > > > I get this error message but as we know from Documentation that If integrate > is unable to compute an integral, it returns an unevaluated Integral object. > > What's happening ,unable to figure out .Can anyone tell me ???? > > -- > You received this message because you are subscribed to the Google Groups > "sympy" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to [email protected]. > To post to this group, send email to [email protected]. > Visit this group at https://groups.google.com/group/sympy. > To view this discussion on the web visit > https://groups.google.com/d/msgid/sympy/21f0b5ae-75ee-4c1a-b2b6-e62f6a196818%40googlegroups.com. > > For more options, visit https://groups.google.com/d/optout. -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To post to this group, send email to [email protected]. Visit this group at https://groups.google.com/group/sympy. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/CAKgW%3D6L1v-eykX3Wn1Mex%2BcH0e0kkfCWgaQSSWv_h-zuTMTSHg%40mail.gmail.com. For more options, visit https://groups.google.com/d/optout.
