In [1]:
from time import clock
from sympy import *
from __future__ import division
x = Symbol('x')
. . .
In [2]:
def asec(x):
return acos(1/x)
. . .
In [3]:
Charlwood = dict()
Charlwood[ 1] = (asin(x)*log(x),-2*sqrt(1-x**2)+atanh(sqrt(1-x**2))-x*asin(x)*(1-log(x))+sqrt(1-x**2)*log(x))
Charlwood[ 2] = (x*asin(x)/sqrt(1-x**2),x-sqrt(1-x**2)*asin(x))
Charlwood[ 3] = (asin(sqrt(x+1)-sqrt(x)),((sqrt(x)+3*sqrt(1+x))*sqrt(-x+sqrt(x)*sqrt(1+x)))/(4*sqrt(2))-(Rational(3,8)+x)*asin(sqrt(x)-sqrt(1+x)))
Charlwood[ 4] = (log(1+x*sqrt(1+x**2)),-2*x+sqrt(2*(1+sqrt(5)))*atan(sqrt(-2+sqrt(5))*(x+sqrt(1+x**2)))-sqrt(2*(-1+sqrt(5)))*atanh(sqrt(2+sqrt(5))*(x+sqrt(1+x**2)))+x*log(1+x*sqrt(1+x**2)))
Charlwood[ 5] = (cos(x)**2/sqrt(cos(x)**4+cos(x)**2+1),x/3+Rational(1,3)*atan((cos(x)*(1+cos(x)**2)*sin(x))/(1+cos(x)**2*sqrt(1+cos(x)**2+cos(x)**4))))
Charlwood[ 6] = (tan(x)*sqrt(1+tan(x)**4),(-Rational(1,2))*asinh(tan(x)**2)-atanh((1-tan(x)**2)/(sqrt(2)*sqrt(1+tan(x)**4)))/sqrt(2)+Rational(1,2)*sqrt(1+tan(x)**4))
Charlwood[ 7] = (tan(x)/sqrt(1+sec(x)**3),(-Rational(2,3))*atanh(sqrt(1+sec(x)**3)))
Charlwood[ 8] = (sqrt(tan(x)**2+2*tan(x)+2),asinh(1+tan(x))+sqrt(Rational(1,2)*(1+sqrt(5)))*atan((-sqrt(-1+sqrt(5))+sqrt(1+sqrt(5))*tan(x))/(sqrt(2)*sqrt(2+tan(x)*(2+tan(x)))))-sqrt((1/2)*(-1+sqrt(5)))*atanh((sqrt(1+sqrt(5))+sqrt(-1+sqrt(5))*tan(x))/(sqrt(2)*sqrt(2+tan(x)*(2+tan(x))))))
Charlwood[ 9] = (sin(x)*atan(sqrt(sec(x)-1)),Rational(1,2)*atan(sqrt(-1+sec(x)))-atan(sqrt(-1+sec(x)))*cos(x)+Rational(1,2)*cos(x)*sqrt(-1+sec(x)))
Charlwood[10] = (x**3*exp(asin(x))/sqrt(1-x**2),Rational(1,10)*exp(asin(x))*(3*x+x**3-3*sqrt(1-x**2)-3*x**2*sqrt(1-x**2)))
Charlwood[11] = ((x*log(1+x**2)*log(x+sqrt(1+x**2)))/sqrt(1+x**2),4*x-2*atan(x)-2*sqrt(1+x**2)*log(x+sqrt(1+x**2))+log(1+x**2)*(-x+sqrt(1+x**2)*log(x+sqrt(1+x**2))))
Charlwood[12] = (atan(x+sqrt(1-x**2)),-(asin(x)/2)+Rational(1,4)*sqrt(3)*atan((-1+sqrt(3)*x)/sqrt(1-x**2))+Rational(1,4)*sqrt(3)*atan((1+sqrt(3)*x)/sqrt(1-x**2))-Rational(1,4)*sqrt(3)*atan((-1+2*x**2)/sqrt(3))+x*atan(x+sqrt(1-x**2))-Rational(1,4)*atanh(x*sqrt(1-x**2))-Rational(1,8)*log(1-x**2+x**4))
Charlwood[13] = (x*atan(x+sqrt(1-x**2))/sqrt(1-x**2),-(asin(x)/2)+Rational(1,4)*sqrt(3)*atan((-1+sqrt(3)*x)/sqrt(1-x**2))+Rational(1,4)*sqrt(3)*atan((1+sqrt(3)*x)/sqrt(1-x**2))-Rational(1,4)*sqrt(3)*atan((-1+2*x**2)/sqrt(3))-sqrt(1-x**2)*atan(x+sqrt(1-x**2))+Rational(1,4)*atanh(x*sqrt(1-x**2))+Rational(1,8)*log(1-x**2+x**4))
Charlwood[14] = (asin(x)/(1+sqrt(1-x**2)),((-1+sqrt(1-x**2))*asin(x))/x+asin(x)**2/2-log(1+sqrt(1-x**2)))
Charlwood[15] = (log(x+sqrt(1+x**2))/(1-x**2)**Rational(3,2),(-Rational(1,2))*asin(x**2)+(x*log(x+sqrt(1+x**2)))/sqrt(1-x**2))
Charlwood[16] = (asin(x)/(1+x**2)**Rational(3,2),(x*asin(x))/sqrt(1+x**2)-asin(x**2)/2)
Charlwood[17] = (log(x+sqrt(x**2-1))/(1+x**2)**Rational(3,2),(-Rational(1,2))*acosh(x**2)+(x*log(x+sqrt(-1+x**2)))/sqrt(1+x**2))
Charlwood[18] = (log(x)/(x**2*sqrt(x**2-1)),sqrt(-1+x**2)/x-atanh(x/sqrt(-1+x**2))+(sqrt(-1+x**2)*log(x))/x)
Charlwood[19] = (sqrt(1+x**3)/x,(2*sqrt(1+x**3))/3-Rational(2,3)*atanh(sqrt(1+x**3)))
Charlwood[20] = (x*log(x+sqrt(x**2-1))/sqrt(x**2-1),-x+sqrt(-1+x**2)*log(x+sqrt(-1+x**2)))
Charlwood[21] = (x**3*(asin(x)/sqrt(1-x**4)),Rational(1,4)*x*sqrt(1+x**2)-Rational(1,2)*sqrt(1-x**4)*asin(x)+asinh(x)/4)
Charlwood[22] = (x**3*(asec(x)/sqrt(x**4-1)),-(sqrt(-1+x**4)/(2*sqrt(1-1/x**2)*x))+Rational(1,2)*sqrt(-1+x**4)*asec(x)+Rational(1,2)*atanh((sqrt(1-1/x**2)*x)/sqrt(-1+x**4)))
Charlwood[23] = (x*atan(x)*log(x+sqrt(1+x**2))/sqrt(1+x**2),(-x)*atan(x)+Rational(1,2)*log(1+x**2)+sqrt(1+x**2)*atan(x)*log(x+sqrt(1+x**2))-Rational(1,2)*log(x+sqrt(1+x**2))**2)
Charlwood[24] = (x*log(1+sqrt(1-x**2))/sqrt(1-x**2),sqrt(1-x**2)-(1+sqrt(1-x**2))*log(1+sqrt(1-x**2)))
Charlwood[25] = (x*log(x+sqrt(1+x**2))/sqrt(1+x**2),-x+sqrt(1+x**2)*log(x+sqrt(1+x**2)))
Charlwood[26] = (x*log(x+sqrt(1-x**2))/sqrt(1-x**2),sqrt(1-x**2)+atanh(sqrt(2)*x)/sqrt(2)-atanh(sqrt(2)*sqrt(1-x**2))/sqrt(2)-sqrt(1-x**2)*log(x+sqrt(1-x**2)))
Charlwood[27] = (log(x)/(x**2*sqrt(1-x**2)),-(sqrt(1-x**2)/x)-asin(x)-(sqrt(1-x**2)*log(x))/x)
Charlwood[28] = (x*atan(x)/sqrt(1+x**2),-asinh(x)+sqrt(1+x**2)*atan(x))
Charlwood[29] = (atan(x)/(x**2*sqrt(1-x**2)),-((sqrt(1-x**2)*atan(x))/x)-atanh(sqrt(1-x**2))+sqrt(2)*atanh(sqrt(1-x**2)/sqrt(2)))
Charlwood[30] = (x*atan(x)/sqrt(1-x**2),-asin(x)-sqrt(1-x**2)*atan(x)+sqrt(2)*atan((sqrt(2)*x)/sqrt(1-x**2)))
Charlwood[31] = (atan(x)/(x**2*sqrt(1+x**2)),-((sqrt(1+x**2)*atan(x))/x)-atanh(sqrt(1+x**2)))
Charlwood[32] = (asin(x)/(x**2*sqrt(1-x**2)),-((sqrt(1-x**2)*asin(x))/x)+log(x))
Charlwood[33] = (x*log(x)/sqrt(x**2-1),-sqrt(-1+x**2)+atan(sqrt(-1+x**2))+sqrt(-1+x**2)*log(x))
Charlwood[34] = (log(x)/(x**2*sqrt(1+x**2)),-(sqrt(1+x**2)/x)+asinh(x)-(sqrt(1+x**2)*log(x))/x)
Charlwood[35] = (x*asec(x)/sqrt(x**2-1),sqrt(x**2-1)*asec(x)-(sqrt(1-1/x**2)*x*log(x))/sqrt(x**2-1))
Charlwood[36] = (x*log(x)/sqrt(1+x**2),-sqrt(1+x**2)+atanh(sqrt(1+x**2))+sqrt(1+x**2)*log(x))
Charlwood[37] = (sin(x)/(1+sin(x)**2),-(atanh(cos(x)/sqrt(2))/sqrt(2)))
Charlwood[38] = ((1+x**2)/((1-x**2)*sqrt(1+x**4)),(1/sqrt(2))*atanh(sqrt(2)*(x/sqrt(1+x**4))))
Charlwood[39] = ((1-x**2)/((1+x**2)*sqrt(1+x**4)),atan((sqrt(2)*x)/sqrt(1+x**4))/sqrt(2))
Charlwood[40] = (log(sin(x))/(1+sin(x)),-x-atanh(cos(x))-(cos(x)*log(sin(x)))/(1+sin(x)))
Charlwood[41] = (log(sin(x))*sqrt(1+sin(x)),(4*cos(x))/sqrt(1+sin(x))-(2*cos(x)*log(sin(x)))/sqrt(1+sin(x))-4*atanh(cos(x)/sqrt(1+sin(x))))
Charlwood[42] = (sec(x)/sqrt(sec(x)**4-1),-(atanh((cos(x)*cot(x)*sqrt(sec(x)**4-1))/sqrt(2))/sqrt(2)))
Charlwood[43] = (tan(x)/sqrt(1+tan(x)**4),-(atanh((1-tan(x)**2)/(sqrt(2)*sqrt(1+tan(x)**4)))/(2*sqrt(2))))
Charlwood[44] = (sin(x)/sqrt(1-sin(x)**6),atanh((sqrt(3)*cos(x)*(1+sin(x)**2))/(2*sqrt(1-sin(x)**6)))/(2*sqrt(3)))
Charlwood[45] = (sqrt(sqrt(sec(x)+1)-sqrt(sec(x)-1)),sqrt(2)*(sqrt(-1+sqrt(2))*atan((sqrt(-2+2*sqrt(2))*(-sqrt(2)-sqrt(-1+sec(x))+sqrt(1+sec(x))))/(2*sqrt(-sqrt(-1+sec(x))+sqrt(1+sec(x)))))-sqrt(1+sqrt(2))*atan((sqrt(2+2*sqrt(2))*(-sqrt(2)-sqrt(-1+sec(x))+sqrt(1+sec(x))))/(2*sqrt(-sqrt(-1+sec(x))+sqrt(1+sec(x)))))-sqrt(1+sqrt(2))*atanh((sqrt(-2+2*sqrt(2))*sqrt(-sqrt(-1+sec(x))+sqrt(1+sec(x))))/(sqrt(2)-sqrt(-1+sec(x))+sqrt(1+sec(x))))+sqrt(-1+sqrt(2))*atanh((sqrt(2+2*sqrt(2))*sqrt(-sqrt(-1+sec(x))+sqrt(1+sec(x))))/(sqrt(2)-sqrt(-1+sec(x))+sqrt(1+sec(x)))))*cot(x)*sqrt(-1+sec(x))*sqrt(1+sec(x)))
Charlwood[46] = (x*log(x**2+1)*atan(x)**2,(-x)*atan(x)*(-3+log(1+x**2))+Rational(1,4)*(-6+log(1+x**2))*log(1+x**2)+Rational(1,2)*atan(x)**2*(-3-x**2+(1+x**2)*log(1+x**2)))
Charlwood[47] = (atan(x*sqrt(1+x**2)),Rational(1,2)*atan(sqrt(1+x**2)/x**2)+x*atan(x*sqrt(1+x**2))+Rational(1,2)*sqrt(3)*atanh((sqrt(3)*sqrt(1+x**2))/(2+x**2)))
Charlwood[48] = (atan(sqrt(x+1)-sqrt(x)),sqrt(x)/2+(1+x)*atan(sqrt(1+x)-sqrt(x)))
Charlwood[49] = (asin(x/sqrt(1-x**2)),x*asin(x/sqrt(1-x**2))+atan(sqrt(1-2*x**2)))
Charlwood[50] = (atan(x*sqrt(1-x**2)),x*atan(x*sqrt(1-x**2))-sqrt(Rational(1,2)*(1+sqrt(5)))*atan(sqrt(Rational(1,2)*(1+sqrt(5)))*sqrt(1-x**2))+sqrt(Rational(1,2)*(-1+sqrt(5)))*atanh(sqrt(Rational(1,2)*(-1+sqrt(5)))*sqrt(1-x**2)))
. . .
In [ ]:
def test(suite, problem):
integrand, antideriv = suite[problem]
print [problem], "Charlwood"
print "integrand :", integrand
print "antideriv :", antideriv
try:
start = clock()
solution = integrate(integrand, x)
duration = clock() - start
except RuntimeError, e:
print "integrator: RuntimeError:", e
print
return
except KeyboardInterrupt:
print "time :", round(clock() - start)
print "integrator: Timed out! \n"
return
except AttributeError, e:
print "integrator: AttributeError:", e
print
return
except ValueError, e:
print "integrator: ValueError:", e
print "comment : Additional constraints? Assume? \n"
return
print "sympy :", solution
str_solution = str(solution)
diffsolution = solution.diff(x)
matheq = 0 == simplify((integrand - diffsolution).expand())
if len(str_solution) > 2*len(str(antideriv)):
print "comment : Solution is messy."
if 'Integral' in str_solution:
print "comment : Unevaluated or not an elementary solution!"
else:
if 'I' in str_solution:
print "comment : Not a real solution!"
if 'Piecewise' in str_solution:
print "comment : Not a continous solution!"
if not matheq:
print "comment : Diffs not back to zero. (Check use of abs.)"
print "time :", round(duration)
if duration > 60 :
print "You are wasting my time!"
print
def test_all(suite):
start = clock()
for problem in suite.iteritems():
test(suite, problem[0])
duration = clock() - start
print "total time:", round(duration)
. . .
In [5]:
test(Charlwood, 1)
[1] Charlwood integrand : log(x)*asin(x) antideriv : -x*(-log(x) + 1)*asin(x) + sqrt(-x**2 + 1)*log(x) - 2*sqrt(-x**2 + 1) + atanh(sqrt(-x**2 + 1)) sympy : Integral(log(x)*asin(x), x) comment : Unevaluated or not an elementary solution! time : 9.0
. . .
In [6]:
test(Charlwood, 2)
[2] Charlwood integrand : x*asin(x)/sqrt(-x**2 + 1) antideriv : x - sqrt(-x**2 + 1)*asin(x) sympy : x - sqrt(-x**2 + 1)*asin(x) time : 0.0
. . .
In [7]:
test(Charlwood, 3)
[3] Charlwood integrand : -asin(sqrt(x) - sqrt(x + 1)) antideriv : sqrt(2)*(sqrt(x) + 3*sqrt(x + 1))*sqrt(sqrt(x)*sqrt(x + 1) - x)/8 - (x + 3/8)*asin(sqrt(x) - sqrt(x + 1)) sympy : Integral(-asin(sqrt(x) - sqrt(x + 1)), x) comment : Unevaluated or not an elementary solution! time : 8.0
. . .
In [8]:
test(Charlwood, 4)
[4] Charlwood integrand : log(x*sqrt(x**2 + 1) + 1) antideriv : x*log(x*sqrt(x**2 + 1) + 1) - 2*x + sqrt(2 + 2*sqrt(5))*atan(sqrt(-2 + sqrt(5))*(x + sqrt(x**2 + 1))) - sqrt(-2 + 2*sqrt(5))*atanh(sqrt(2 + sqrt(5))*(x + sqrt(x**2 + 1))) time : 75.0 integrator: Timed out!
. . .
In [9]:
test(Charlwood, 5)
[5] Charlwood integrand : cos(x)**2/sqrt(cos(x)**4 + cos(x)**2 + 1) antideriv : x/3 + atan((cos(x)**2 + 1)*sin(x)*cos(x)/(sqrt(cos(x)**4 + cos(x)**2 + 1)*cos(x)**2 + 1))/3 time : 60.0 integrator: Timed out!
. . .
In [10]:
test(Charlwood, 6)
[6] Charlwood integrand : sqrt(tan(x)**4 + 1)*tan(x) antideriv : sqrt(tan(x)**4 + 1)/2 - asinh(tan(x)**2)/2 - sqrt(2)*atanh(sqrt(2)*(-tan(x)**2 + 1)/(2*sqrt(tan(x)**4 + 1)))/2 sympy : Integral(sqrt(tan(x)**4 + 1)*tan(x), x) comment : Unevaluated or not an elementary solution! time : 1.0
. . .
In [11]:
test(Charlwood, 7)
[7] Charlwood integrand : tan(x)/sqrt(sec(x)**3 + 1) antideriv : -2*atanh(sqrt(sec(x)**3 + 1))/3 sympy : Integral(tan(x)/sqrt(sec(x)**3 + 1), x) comment : Unevaluated or not an elementary solution! time : 1.0
. . .
In [12]:
test(Charlwood, 8)
[8] Charlwood integrand : sqrt(tan(x)**2 + 2*tan(x) + 2) antideriv : asinh(tan(x) + 1) + sqrt(1/2 + sqrt(5)/2)*atan(sqrt(2)*(sqrt(1 + sqrt(5))*tan(x) - sqrt(-1 + sqrt(5)))/(2*sqrt((tan(x) + 2)*tan(x) + 2))) - sqrt(-0.5 + 0.5*sqrt(5))*atanh(sqrt(2)*(sqrt(-1 + sqrt(5))*tan(x) + sqrt(1 + sqrt(5)))/(2*sqrt((tan(x) + 2)*tan(x) + 2))) sympy : Integral(sqrt(tan(x)**2 + 2*tan(x) + 2), x) comment : Unevaluated or not an elementary solution! time : 1.0
. . .
In [13]:
test(Charlwood, 9)
[9] Charlwood integrand : sin(x)*atan(sqrt(sec(x) - 1)) antideriv : sqrt(sec(x) - 1)*cos(x)/2 - cos(x)*atan(sqrt(sec(x) - 1)) + atan(sqrt(sec(x) - 1))/2 time : 62.0 integrator: Timed out!
. . .
In [14]:
test(Charlwood, 10)
[10] Charlwood integrand : x**3*exp(asin(x))/sqrt(-x**2 + 1) antideriv : (x**3 - 3*x**2*sqrt(-x**2 + 1) + 3*x - 3*sqrt(-x**2 + 1))*exp(asin(x))/10 sympy : x**3*exp(asin(x))/10 - 3*x**2*sqrt(-x**2 + 1)*exp(asin(x))/10 + 3*x*exp(asin(x))/10 - 3*sqrt(-x**2 + 1)*exp(asin(x))/10 time : 5.0
. . .
In [15]:
test(Charlwood, 11)
[11] Charlwood integrand : x*log(x + sqrt(x**2 + 1))*log(x**2 + 1)/sqrt(x**2 + 1) antideriv : 4*x + (-x + sqrt(x**2 + 1)*log(x + sqrt(x**2 + 1)))*log(x**2 + 1) - 2*sqrt(x**2 + 1)*log(x + sqrt(x**2 + 1)) - 2*atan(x) time : 61.0 integrator: Timed out!
. . .
In [16]:
test(Charlwood, 12)
[12] Charlwood integrand : atan(x + sqrt(-x**2 + 1)) antideriv : x*atan(x + sqrt(-x**2 + 1)) - log(x**4 - x**2 + 1)/8 - asin(x)/2 - sqrt(3)*atan(sqrt(3)*(2*x**2 - 1)/3)/4 + sqrt(3)*atan((sqrt(3)*x - 1)/sqrt(-x**2 + 1))/4 + sqrt(3)*atan((sqrt(3)*x + 1)/sqrt(-x**2 + 1))/4 - atanh(x*sqrt(-x**2 + 1))/4 time : 70.0 integrator: Timed out!
. . .
In [17]:
test(Charlwood, 13)
[13] Charlwood integrand : x*atan(x + sqrt(-x**2 + 1))/sqrt(-x**2 + 1) antideriv : -sqrt(-x**2 + 1)*atan(x + sqrt(-x**2 + 1)) + log(x**4 - x**2 + 1)/8 - asin(x)/2 - sqrt(3)*atan(sqrt(3)*(2*x**2 - 1)/3)/4 + sqrt(3)*atan((sqrt(3)*x - 1)/sqrt(-x**2 + 1))/4 + sqrt(3)*atan((sqrt(3)*x + 1)/sqrt(-x**2 + 1))/4 + atanh(x*sqrt(-x**2 + 1))/4 time : 63.0 integrator: Timed out!
. . .
In [20]:
test(Charlwood, 14)
[14] Charlwood integrand : asin(x)/(sqrt(-x**2 + 1) + 1) antideriv : -log(sqrt(-x**2 + 1) + 1) + asin(x)**2/2 + (sqrt(-x**2 + 1) - 1)*asin(x)/x sympy : Integral(asin(x)/(sqrt(-x**2 + 1) + 1), x) comment : Unevaluated or not an elementary solution! time : 0.0
. . .
In [21]:
test(Charlwood, 15)
[15] Charlwood integrand : log(x + sqrt(x**2 + 1))/(-x**2 + 1)**(3/2) antideriv : x*log(x + sqrt(x**2 + 1))/sqrt(-x**2 + 1) - asin(x**2)/2 time : 63.0 integrator: Timed out!
. . .
In [22]:
test(Charlwood, 16)
[16] Charlwood integrand : asin(x)/(x**2 + 1)**(3/2) antideriv : x*asin(x)/sqrt(x**2 + 1) - asin(x**2)/2 sympy : Integral(asin(x)/(x**2 + 1)**(3/2), x) comment : Unevaluated or not an elementary solution! time : 48.0
. . .
In [23]:
test(Charlwood, 17)
[17] Charlwood integrand : log(x + sqrt(x**2 - 1))/(x**2 + 1)**(3/2) antideriv : x*log(x + sqrt(x**2 - 1))/sqrt(x**2 + 1) - acosh(x**2)/2 time : 61.0 integrator: Timed out!
. . .
In [24]:
test(Charlwood, 18)
[18] Charlwood integrand : log(x)/(x**2*sqrt(x**2 - 1)) antideriv : -atanh(x/sqrt(x**2 - 1)) + sqrt(x**2 - 1)*log(x)/x + sqrt(x**2 - 1)/x sympy : Integral(log(x)/(x**2*sqrt(x**2 - 1)), x) comment : Unevaluated or not an elementary solution! time : 41.0
. . .
In [25]:
test(Charlwood, 19)
[19] Charlwood integrand : sqrt(x**3 + 1)/x antideriv : 2*sqrt(x**3 + 1)/3 - 2*atanh(sqrt(x**3 + 1))/3 sympy : 2*x**(3/2)/(3*sqrt(1 + x**(-3))) - 2*asinh(x**(-3/2))/3 + 2/(3*x**(3/2)*sqrt(1 + x**(-3))) comment : Diffs not back to zero. (Check use of abs.) time : 1.0
. . .
In [26]:
test(Charlwood, 20)
[20] Charlwood integrand : x*log(x + sqrt(x**2 - 1))/sqrt(x**2 - 1) antideriv : -x + sqrt(x**2 - 1)*log(x + sqrt(x**2 - 1)) time : 59.0 integrator: Timed out!
. . .
In [27]:
test(Charlwood, 21)
[21] Charlwood integrand : x**3*asin(x)/sqrt(-x**4 + 1) antideriv : x*sqrt(x**2 + 1)/4 - sqrt(-x**4 + 1)*asin(x)/2 + asinh(x)/4 sympy : Integral(x**3*asin(x)/sqrt(-x**4 + 1), x) comment : Unevaluated or not an elementary solution! time : 37.0
. . .
In [28]:
test(Charlwood, 22)
[22] Charlwood integrand : x**3*acos(1/x)/sqrt(x**4 - 1) antideriv : sqrt(x**4 - 1)*acos(1/x)/2 + atanh(x*sqrt(1 - 1/x**2)/sqrt(x**4 - 1))/2 - sqrt(x**4 - 1)/(2*x*sqrt(1 - 1/x**2)) time : 62.0 integrator: Timed out!
. . .
In [29]:
test(Charlwood, 23)
[23] Charlwood integrand : x*log(x + sqrt(x**2 + 1))*atan(x)/sqrt(x**2 + 1) antideriv : -x*atan(x) + sqrt(x**2 + 1)*log(x + sqrt(x**2 + 1))*atan(x) - log(x + sqrt(x**2 + 1))**2/2 + log(x**2 + 1)/2 time : 64.0 integrator: Timed out!
. . .
In [30]:
test(Charlwood, 24)
[24] Charlwood integrand : x*log(sqrt(-x**2 + 1) + 1)/sqrt(-x**2 + 1) antideriv : sqrt(-x**2 + 1) - (sqrt(-x**2 + 1) + 1)*log(sqrt(-x**2 + 1) + 1) time : 61.0 integrator: Timed out!
. . .
In [31]:
test(Charlwood, 25)
[25] Charlwood integrand : x*log(x + sqrt(x**2 + 1))/sqrt(x**2 + 1) antideriv : -x + sqrt(x**2 + 1)*log(x + sqrt(x**2 + 1)) time : 64.0 integrator: Timed out!
. . .
In [32]:
test(Charlwood, 26)
[26] Charlwood integrand : x*log(x + sqrt(-x**2 + 1))/sqrt(-x**2 + 1) antideriv : -sqrt(-x**2 + 1)*log(x + sqrt(-x**2 + 1)) + sqrt(-x**2 + 1) + sqrt(2)*atanh(sqrt(2)*x)/2 - sqrt(2)*atanh(sqrt(2)*sqrt(-x**2 + 1))/2 time : 66.0 integrator: Timed out!
. . .
In [33]:
test(Charlwood, 27)
[27] Charlwood integrand : log(x)/(x**2*sqrt(-x**2 + 1)) antideriv : -asin(x) - sqrt(-x**2 + 1)*log(x)/x - sqrt(-x**2 + 1)/x sympy : Integral(log(x)/(x**2*sqrt(-x**2 + 1)), x) comment : Unevaluated or not an elementary solution! time : 41.0
. . .
In [34]:
test(Charlwood, 28)
[28] Charlwood integrand : x*atan(x)/sqrt(x**2 + 1) antideriv : sqrt(x**2 + 1)*atan(x) - asinh(x) sympy : x**2*atan(x)/sqrt(x**2 + 1) - asinh(x) + atan(x)/sqrt(x**2 + 1) time : 6.0
. . .
In [35]:
test(Charlwood, 29)
[29] Charlwood integrand : atan(x)/(x**2*sqrt(-x**2 + 1)) antideriv : sqrt(2)*atanh(sqrt(2)*sqrt(-x**2 + 1)/2) - atanh(sqrt(-x**2 + 1)) - sqrt(-x**2 + 1)*atan(x)/x sympy : Integral(atan(x)/(x**2*sqrt(-x**2 + 1)), x) comment : Unevaluated or not an elementary solution! time : 25.0
. . .
In [36]:
test(Charlwood, 30)
[30] Charlwood integrand : x*atan(x)/sqrt(-x**2 + 1) antideriv : -sqrt(-x**2 + 1)*atan(x) - asin(x) + sqrt(2)*atan(sqrt(2)*x/sqrt(-x**2 + 1)) time : 63.0 integrator: Timed out!
. . .
In [37]:
test(Charlwood, 31)
[31] Charlwood integrand : atan(x)/(x**2*sqrt(x**2 + 1)) antideriv : -atanh(sqrt(x**2 + 1)) - sqrt(x**2 + 1)*atan(x)/x time : 67.0 integrator: Timed out!
. . .
In [38]:
test(Charlwood, 32)
[32] Charlwood integrand : asin(x)/(x**2*sqrt(-x**2 + 1)) antideriv : log(x) - sqrt(-x**2 + 1)*asin(x)/x sympy : Integral(asin(x)/(x**2*sqrt(-x**2 + 1)), x) comment : Unevaluated or not an elementary solution! time : 0.0
. . .
In [39]:
test(Charlwood, 33)
[33] Charlwood integrand : x*log(x)/sqrt(x**2 - 1) antideriv : sqrt(x**2 - 1)*log(x) - sqrt(x**2 - 1) + atan(sqrt(x**2 - 1)) sympy : Integral(x*log(x)/sqrt(x**2 - 1), x) comment : Unevaluated or not an elementary solution! time : 10.0
. . .
In [40]:
test(Charlwood, 34)
[34] Charlwood integrand : log(x)/(x**2*sqrt(x**2 + 1)) antideriv : asinh(x) - sqrt(x**2 + 1)*log(x)/x - sqrt(x**2 + 1)/x sympy : Integral(log(x)/(x**2*sqrt(x**2 + 1)), x) comment : Unevaluated or not an elementary solution! time : 44.0
. . .
In [41]:
test(Charlwood, 35)
[35] Charlwood integrand : x*acos(1/x)/sqrt(x**2 - 1) antideriv : -x*sqrt(1 - 1/x**2)*log(x)/sqrt(x**2 - 1) + sqrt(x**2 - 1)*acos(1/x) time : 68.0 integrator: Timed out!
. . .
In [42]:
test(Charlwood, 36)
[36] Charlwood integrand : x*log(x)/sqrt(x**2 + 1) antideriv : sqrt(x**2 + 1)*log(x) - sqrt(x**2 + 1) + atanh(sqrt(x**2 + 1)) sympy : Integral(x*log(x)/sqrt(x**2 + 1), x) comment : Unevaluated or not an elementary solution! time : 10.0
. . .
In [43]:
test(Charlwood, 37)
[37] Charlwood integrand : sin(x)/(sin(x)**2 + 1) antideriv : -sqrt(2)*atanh(sqrt(2)*cos(x)/2)/2 time : 59.0 integrator: Timed out!
. . .
In [44]:
test(Charlwood, 38)
[38] Charlwood integrand : (x**2 + 1)/((-x**2 + 1)*sqrt(x**4 + 1)) antideriv : sqrt(2)*atanh(sqrt(2)*x/sqrt(x**4 + 1))/2 sympy : Integral((x**2 + 1)/((-x**2 + 1)*sqrt(x**4 + 1)), x) comment : Unevaluated or not an elementary solution! time : 1.0
. . .
In [45]:
test(Charlwood, 39)
[39] Charlwood integrand : (-x**2 + 1)/((x**2 + 1)*sqrt(x**4 + 1)) antideriv : sqrt(2)*atan(sqrt(2)*x/sqrt(x**4 + 1))/2 sympy : Integral((-x**2 + 1)/((x**2 + 1)*sqrt(x**4 + 1)), x) comment : Unevaluated or not an elementary solution! time : 1.0
. . .
In [46]:
test(Charlwood, 40)
[40] Charlwood integrand : log(sin(x))/(sin(x) + 1) antideriv : -x - atanh(cos(x)) - log(sin(x))*cos(x)/(sin(x) + 1) time : 73.0 integrator: Timed out!
. . .
In [47]:
test(Charlwood, 41)
[41] Charlwood integrand : sqrt(sin(x) + 1)*log(sin(x)) antideriv : -4*atanh(cos(x)/sqrt(sin(x) + 1)) - 2*log(sin(x))*cos(x)/sqrt(sin(x) + 1) + 4*cos(x)/sqrt(sin(x) + 1) sympy : Integral(sqrt(sin(x) + 1)*log(sin(x)), x) comment : Unevaluated or not an elementary solution! time : 35.0
. . .
In [48]:
test(Charlwood, 42)
[42] Charlwood integrand : sec(x)/sqrt(sec(x)**4 - 1) antideriv : -sqrt(2)*atanh(sqrt(2)*sqrt(sec(x)**4 - 1)*cos(x)*cot(x)/2)/2 sympy : Integral(sec(x)/sqrt(sec(x)**4 - 1), x) comment : Unevaluated or not an elementary solution! time : 0.0
. . .
In [49]:
test(Charlwood, 43)
[43] Charlwood integrand : tan(x)/sqrt(tan(x)**4 + 1) antideriv : -sqrt(2)*atanh(sqrt(2)*(-tan(x)**2 + 1)/(2*sqrt(tan(x)**4 + 1)))/4 sympy : Integral(tan(x)/sqrt(tan(x)**4 + 1), x) comment : Unevaluated or not an elementary solution! time : 0.0
. . .
In [50]:
test(Charlwood, 44)
[44] Charlwood integrand : sin(x)/sqrt(-sin(x)**6 + 1) antideriv : sqrt(3)*atanh(sqrt(3)*(sin(x)**2 + 1)*cos(x)/(2*sqrt(-sin(x)**6 + 1)))/6 time : 60.0 integrator: Timed out!
. . .
In [51]:
test(Charlwood, 45)
[45] Charlwood integrand : sqrt(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1)) antideriv : sqrt(2)*sqrt(sec(x) - 1)*sqrt(sec(x) + 1)*(sqrt(-1 + sqrt(2))*atan(sqrt(-2 + 2*sqrt(2))*(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1) - sqrt(2))/(2*sqrt(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1)))) - sqrt(1 + sqrt(2))*atan(sqrt(2 + 2*sqrt(2))*(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1) - sqrt(2))/(2*sqrt(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1)))) - sqrt(1 + sqrt(2))*atanh(sqrt(-2 + 2*sqrt(2))*sqrt(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1))/(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1) + sqrt(2))) + sqrt(-1 + sqrt(2))*atanh(sqrt(2 + 2*sqrt(2))*sqrt(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1))/(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1) + sqrt(2))))*cot(x) sympy : Integral(sqrt(-sqrt(sec(x) - 1) + sqrt(sec(x) + 1)), x) comment : Unevaluated or not an elementary solution! time : 17.0
. . .
In [52]:
test(Charlwood, 46)
[46] Charlwood integrand : x*log(x**2 + 1)*atan(x)**2 antideriv : -x*(log(x**2 + 1) - 3)*atan(x) + (log(x**2 + 1)/4 - 3/2)*log(x**2 + 1) + (-x**2 + (x**2 + 1)*log(x**2 + 1) - 3)*atan(x)**2/2 sympy : x**2*log(x**2 + 1)*atan(x)**2/2 - x**2*atan(x)**2/2 - x*log(x**2 + 1)*atan(x) + 3*x*atan(x) + log(x**2 + 1)**2/4 + log(x**2 + 1)*atan(x)**2/2 - 3*log(x**2 + 1)/2 - 3*atan(x)**2/2 time : 19.0
. . .
In [53]:
test(Charlwood, 47)
[47] Charlwood integrand : atan(x*sqrt(x**2 + 1)) antideriv : x*atan(x*sqrt(x**2 + 1)) + atan(sqrt(x**2 + 1)/x**2)/2 + sqrt(3)*atanh(sqrt(3)*sqrt(x**2 + 1)/(x**2 + 2))/2 time : 64.0 integrator: Timed out!
. . .
In [54]:
test(Charlwood, 48)
[48] Charlwood integrand : -atan(sqrt(x) - sqrt(x + 1)) antideriv : sqrt(x)/2 - (x + 1)*atan(sqrt(x) - sqrt(x + 1)) time : 64.0 integrator: Timed out!
. . .
In [55]:
test(Charlwood, 49)
[49] Charlwood integrand : asin(x/sqrt(-x**2 + 1)) antideriv : x*asin(x/sqrt(-x**2 + 1)) + atan(sqrt(-2*x**2 + 1)) time : 67.0 integrator: Timed out!
. . .
In [56]:
test(Charlwood, 50)
[50] Charlwood integrand : atan(x*sqrt(-x**2 + 1)) antideriv : x*atan(x*sqrt(-x**2 + 1)) - sqrt(1/2 + sqrt(5)/2)*atan(sqrt(1/2 + sqrt(5)/2)*sqrt(-x**2 + 1)) + sqrt(-1/2 + sqrt(5)/2)*atanh(sqrt(-1/2 + sqrt(5)/2)*sqrt(-x**2 + 1)) time : 67.0 integrator: Timed out!
. . .