Ticket #8459: trac-8459.patch

sage/calculus/calculus.py

@@ -1425 +1425 @@
 
 #######################################################
 
-symtable = {'%pi':'pi', '%e': 'e', '%i':'I', '%gamma':'euler_gamma',
-            'li[2]':'polylog2', 'li[3]':'polylog3'}
+symtable = {'%pi':'pi', '%e': 'e', '%i':'I', '%gamma':'euler_gamma'}
 
 from sage.symbolic.pynac import register_symbol
 from sage.rings.infinity import infinity, minus_infinity
@@ -1475 +1474 @@
         x#0
         sage: a.sage()
         x != 0
+
+    TESTS:
+    
+    Trac #8459 fixed::
+
+        sage: maxima('3*li[2](u)+8*li[33](exp(u))').sage()
+        3*polylog(2, u) + 8*polylog(33, e^u)
+            
     """
     syms = dict(_syms)
 
@@ -1505 +1512 @@
                 syms[X[2:]] = function_factory(X[2:])
         s = s.replace("?%","")
 
+    regex=re.compile('li\[([0-9]+?)\]\(')
+    s = regex.sub('polylog(\\1,',s) 
     s = multiple_replace(symtable, s)
     s = s.replace("%","")
 

sage/calculus/tests.py

@@ -105 +105 @@
     sage: expand(integrate(log(1-x^2), x))
     x*log(-x^2 + 1) - 2*x - log(x - 1) + log(x + 1)
     sage: integrate(log(1-x^2)/x, x)
-    log(-x^2 + 1)*log(x) + 1/2*polylog2(-x^2 + 1)
+    log(-x^2 + 1)*log(x) + 1/2*polylog(2, -x^2 + 1)
     sage: integrate(exp(1-x^2),x)
     1/2*sqrt(pi)*e*erf(x)
     sage: integrate(sin(x^2),x)