home bbs files messages ]

Forums before death by AOL, social media and spammers... "We can't have nice things"

   sci.math.symbolic      Symbolic algebra discussion      10,432 messages   

[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]

   Message 9,816 of 10,432   
   clicliclic@freenet.de to Albert Rich   
   Re: Need a Rubi rule for (x*exp(x))/(1-e   
   21 Feb 18 14:52:02   
   
   Albert Rich schrieb:   
   >   
   > On Monday, February 19, 2018 at 7:26:06 AM UTC-10, Nasser M. Abbasi wrote:   
   > > Rubi 4.14.6  can integrate   
   > >   
   > >         (x Exp[x])/(1 - Exp[x])     (1)   
   > >   
   > > but not   
   > >   
   > >         (x Exp[x])/(1 - Exp[-x])    (2)   
   > >   
   > > I was just wondering why that is. Could not figure what changes   
   > > in rules needed. Since the result are very similar (in terms   
   > > of PolyLog, which Rubi supprts). For example, for (1) it   
   > > gives   
   > >   
   > > In[59]:= rubi = Int[(x Exp[x])/(1 - Exp[x]), x]   
   > > Out[59]= -x Log[1 - E^x] - PolyLog[2, E^x]   
   > >   
   > > But for (2) it does not evaluate it. Mathematica gives for (2)   
   > >   
   > > In[60]:= mma = Integrate[(x Exp[x])/(1 - Exp[-x]), x]   
   > > Out[60]= E^x (-1 + x) + x Log[1 - E^x] + PolyLog[2, E^x]   
   > >   
   > > thanks   
   > > --Nasser   
   >   
   > Thank you for reporting this embarrassing deficiency in Rubi. The just   
   > released version 4.14.7 includes a new rule which algebraically   
   > expands your example integrand into x*E^x+x/(1-E^(-x)). This sum is   
   > then integrated in 7 steps using existing rules.   
   >   
   > Among other enhancements, Rubi 4.14.7 can find elementary   
   > antiderivatives for integrands of the form x/((c+d*x^3)*sqrt(a+b*x^3))   
   > when b*c-4*a*d=0, or b*c+8*a*d=0, or b*c-2*(5+3*sqrt(3))*a*d=0.  These   
   > pseudo-elliptic integrals are discussed in the sci.math.symbolic   
   > thread "More teething help".   
   >   
   > I am still trying to find optimal antiderivatives for when   
   > b*c-2*(5-3*sqrt(3))*a*d=0. The best so far has 3 arctan(h) terms, but   
   > involves the imaginary unit. Maybe some system or person can do   
   > better...   
   >   
      
   Oops, are you perhaps still looking for some of these:   
      
   INT(x/((x^3 + 6*SQRT(3) + 10)*SQRT(x^3 + 1)), x)   
    = SQRT(2)*3^(3/4)*(2 - SQRT(3))/36*   
   (2*ATANH(12^(1/4)*(2*x - SQRT(3) - 1)/(2*SQRT(x^3 + 1)))   
    + ATANH(12^(1/4)*(SQRT(3) - 1)*(x + 1)/(2*SQRT(x^3 + 1))))   
    + 12^(1/4)*(2 - SQRT(3))/36*   
   (2*ATAN(12^(1/4)*(SQRT(3) - 1)*SQRT(x^3 + 1)/6)   
    - 3*ATAN(12^(1/4)*(SQRT(3) + 1)*(x + 1)/(2*SQRT(x^3 + 1))))   
      
   INT(x/((x^3 - 6*SQRT(3) + 10)*SQRT(x^3 + 1)), x)   
    = 12^(1/4)*(2 + SQRT(3))/36*   
   (2*ATANH(12^(1/4)*(SQRT(3) + 1)*SQRT(x^3 + 1)/6)   
    - 3*ATANH(12^(1/4)*(SQRT(3) - 1)*(x + 1)/(2*SQRT(x^3 + 1))))   
    + SQRT(2)*3^(3/4)*(2 + SQRT(3))/36*   
   (2*ATAN(12^(1/4)*(2*x + SQRT(3) - 1)/(2*SQRT(x^3 + 1)))   
    - ATAN(12^(1/4)*(SQRT(3) + 1)*(x + 1)/(2*SQRT(x^3 + 1))))   
      
   INT(x/((x^3 - 6*SQRT(3) - 10)*SQRT(x^3 - 1)), x)   
    = 12^(1/4)*(2 - SQRT(3))/36*   
   (2*ATANH(12^(1/4)*(SQRT(3) - 1)*SQRT(x^3 - 1)/6)   
    - 3*ATANH(12^(1/4)*(SQRT(3) + 1)*(x - 1)/(2*SQRT(x^3 - 1))))   
    + SQRT(2)*3^(3/4)*(2 - SQRT(3))/36*   
   (2*ATAN(12^(1/4)*(2*x + SQRT(3) + 1)/(2*SQRT(x^3 - 1)))   
    + ATAN(12^(1/4)*(SQRT(3) - 1)*(x - 1)/(2*SQRT(x^3 - 1))))   
      
   INT(x/((x^3 + 6*SQRT(3) - 10)*SQRT(x^3 - 1)), x)   
    = SQRT(2)*3^(3/4)*(2 + SQRT(3))/36*   
   (2*ATANH(12^(1/4)*(2*x - SQRT(3) + 1)/(2*SQRT(x^3 - 1)))   
    - ATANH(12^(1/4)*(SQRT(3) + 1)*(x - 1)/(2*SQRT(x^3 - 1))))   
    + 12^(1/4)*(2 + SQRT(3))/36*   
   (2*ATAN(12^(1/4)*(SQRT(3) + 1)*SQRT(x^3 - 1)/6)   
    - 3*ATAN(12^(1/4)*(SQRT(3) - 1)*(x - 1)/(2*SQRT(x^3 - 1))))   
      
   ? To flip the sign of SQRT(3), one must systematically replace 3^(1/4)   
   by #i*3^(1/4), SQRT(3) by -SQRT(3), and 3^(3/4) by -#i*3^(3/4) in both   
   the integrand and the evaluation. Risch's disciples will know why.   
      
   Or are you trying to exorcise jumps with imaginary units?   
      
   Martin.   
      
   --- SoupGate-Win32 v1.05   
    * Origin: you cannot sedate... all the things you hate (1:229/2)   

[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]


(c) 1994,  bbs@darkrealms.ca