Skip to content

Instantly share code, notes, and snippets.

@asmeurer
Created June 5, 2011 05:31
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save asmeurer/1008685 to your computer and use it in GitHub Desktop.
Save asmeurer/1008685 to your computer and use it in GitHub Desktop.
32-bit Python risch_integrate(logarithms)
In [1]: for i in range(100):
...: print risch_integrate(Add(*(log(x**j) for j in range(i))), x)
...:
0
0
-x + x*log(x)
-3*x + 3*x*log(x**2)/2
-6*x + 3*x*log(x**2)
-10*x + 5*x*log(x**4)/2
-15*x + 15*x*log(x**4)/4
-21*x + 21*x*log(x**4)/4
-28*x + 7*x*log(x**4)
-36*x + 9*x*log(x**4)
-45*x + 45*x*log(x**4)/4
-55*x + 55*x*log(x**4)/4
-66*x + 33*x*log(x**4)/2
-78*x + 39*x*log(x**4)/2
-91*x + 91*x*log(x**4)/4
-105*x + 105*x*log(x**4)/4
-120*x + 30*x*log(x**4)
-136*x + 34*x*log(x**4)
-153*x + 153*x*log(x**4)/4
-171*x + 171*x*log(x**4)/4
-190*x + 95*x*log(x**4)/2
-210*x + 105*x*log(x**4)/2
-231*x + 231*x*log(x**4)/4
-253*x + 253*x*log(x**13)/13
-276*x + 276*x*log(x**13)/13
-300*x + 300*x*log(x**13)/13
-325*x + 25*x*log(x**13)
-351*x + 27*x*log(x**13)
-378*x + 378*x*log(x**13)/13
-406*x + 406*x*log(x**13)/13
-435*x + 435*x*log(x**13)/13
-465*x + 465*x*log(x**13)/13
-496*x + 496*x*log(x**13)/13
-528*x + 528*x*log(x**13)/13
-561*x + 561*x*log(x**13)/13
-595*x + 595*x*log(x**13)/13
-630*x + 630*x*log(x**13)/13
-666*x + 37*x*log(x**36)/2
-703*x + 703*x*log(x**36)/36
-741*x + 247*x*log(x**36)/12
-780*x + 65*x*log(x**36)/3
-820*x + 205*x*log(x**36)/9
-861*x + 287*x*log(x**36)/12
-903*x + 301*x*log(x**36)/12
-946*x + 86*x*log(x**11)
-990*x + 90*x*log(x**11)
-1035*x + 1035*x*log(x**11)/11
-1081*x + 1081*x*log(x**11)/11
-1128*x + 1128*x*log(x**11)/11
-1176*x + 1176*x*log(x**11)/11
-1225*x + 1225*x*log(x**11)/11
-1275*x + 1275*x*log(x**11)/11
-1326*x + 1326*x*log(x**11)/11
-1378*x + 1378*x*log(x**11)/11
-1431*x + 1431*x*log(x**11)/11
-1485*x + 135*x*log(x**11)
-1540*x + 140*x*log(x**11)
-1596*x + 1596*x*log(x**11)/11
-1653*x + 1653*x*log(x**11)/11
-1711*x + 1711*x*log(x**11)/11
-1770*x + 1770*x*log(x**11)/11
-1830*x + 1830*x*log(x**11)/11
-1891*x + 1891*x*log(x**11)/11
-1953*x + 1953*x*log(x**11)/11
-2016*x + 2016*x*log(x**11)/11
-2080*x + 2080*x*log(x**11)/11
-2145*x + 195*x*log(x**11)
-2211*x + 201*x*log(x**11)
-2278*x + 2278*x*log(x**11)/11
-2346*x + 2346*x*log(x**11)/11
-2415*x + 2415*x*log(x**11)/11
-2485*x + 2485*x*log(x**11)/11
-2556*x + 2556*x*log(x**11)/11
-2628*x + 2628*x*log(x**11)/11
-2701*x + 2701*x*log(x**11)/11
-2775*x + 2775*x*log(x**11)/11
-2850*x + 2850*x*log(x**11)/11
-2926*x + 266*x*log(x**11)
-3003*x + 39*x*log(x**77)
-3081*x + 3081*x*log(x**77)/77
-3160*x + 3160*x*log(x**77)/77
-3240*x + 3240*x*log(x**77)/77
-3321*x + 3321*x*log(x**77)/77
-3403*x + 3403*x*log(x**77)/77
-3486*x + 498*x*log(x**77)/11
-3570*x + 510*x*log(x**77)/11
-3655*x + 3655*x*log(x**77)/77
-3741*x + 3741*x*log(x**11)/11
-3828*x + 348*x*log(x**11)
-3916*x + 356*x*log(x**11)
-4005*x + 4005*x*log(x**11)/11
-4095*x + 4095*x*log(x**11)/11
-4186*x + 4186*x*log(x**11)/11
-4278*x + 4278*x*log(x**11)/11
-4371*x + 4371*x*log(x**11)/11
-4465*x + 4465*x*log(x**11)/11
-4560*x + 4560*x*log(x**11)/11
-4656*x + 4656*x*log(x**11)/11
-4753*x + 4753*x*log(x**11)/11
-4851*x + 441*x*log(x**11)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment