Recently I came across with respect to this post of mine hyperbolic solution to the cubic equation for one real root given by
t=−2√p3sinh(13sinh−1(3q2p√3p))
Intuitively I sought to find the related definitely integral,
I=∫1−113sinh−1(3√32(1−t2))dt
Unfortunately, there was no closed form solution. However, the Integral is amazingly near √2.
I=0.8285267994716327,I2+1=1.4142633998
To investigated more, I tried a heuristic expansion of the integral into Egyptian fractions. Although it gets problematic after the 4th term,
The first four terms are,
I2+1=1+12−112−1416
Here the denominators can be given by,
an=n∑k=0nCk(2n−2kq)n−kqk,q=√2
(Likewise, the denominators in the expansion for √2 are related to Pell numbers, which makes me believe that my integral too is somewhat related to the numbers an.) Therefore, I am finding either a closed form or possibly a fast converging infinite series solution to the integral, just any of these. Thanks for any help.
For t=sinz and applying integration by parts, I get another, somewhat simpler, indefinite integral,
sinz3sinh−1(3√32cos2z)+2√3∫sin2zcoszdz√27cos4z+4
Then again I am stuck. Moreover, this expression ensures that my definite integral is an improper one.
A closed solution in terms of incomplete elliptic integrals with complex arguments is, as given by a user in the comments section,
49(9+2√3i)[F(sin−1√331(9+2√3i);131(23−12√3i))−E(sin−1√331(9+2√3i);131(23−12√3i))]
However, I am still wondering how to transform this into a real number, especially the an connection of the integral is fascinating my mind.
Answer
We have from symmetry that I=23∫10sinh−1[3√32(1−x2)]dx
So we define
f(a)=∫10sinh−1[a(1−x2)]dx
Then we recall that
sinh−1(x)=x2F1(12,12;32;−x2)=∑n≥0(−1)n(1/2)2n(3/2)nx2n+1n!
so
sinh−1[a(1−x2)]=a(1−x2)2F1(12,12;32;−a2(1−x2)2)=∑n≥0(−1)na2n+1n!(1/2)2n(3/2)n(1−x2)2n+1
so
f(a)=∑n≥0(−1)na2n+1n!(1/2)2n(3/2)n∫10(1−x2)2n+1dx
For this integral, we use x=sin(t):
jn=∫10(1−x2)2n+1dx=∫π/20cos(t)4n+3dt
I leave it as a challenge to you to show that ∫π/20sin(t)acos(t)bdt=Γ(a+12)Γ(b+12)2Γ(a+b2+1)
choosing b=4n+3, a=0 we have
jn=Γ(1/2)Γ(2n+2)2Γ(2n+5/2)
Then defining
tn=(1/2)2n(3/2)njn
we have tn+1tn=(n+12)2(n+1)(n+74)(n+54)
Which gives f(a)=a3F2(12,12,1;74,54;−a2)
And since I=23f(3√3/2) we have (assuming I've made no mistakes),
I=√33F2(12,12,1;74,54;−274)
No comments:
Post a Comment