This is an attempt to get the system of equations for the calculation of the coefficients aᵢ of the Taylor polynomial developed around 0.
We want the coefficients aᵢ of this Taylor expansion:
\( {^xa}=\sum_{n=0}^{\infty}{a_n .x^n} \)
They should match this equation:
\( {\color{Blue} {^{x+1}a} \)\( \,=\, \)\( {\color{Red} a^{^xa}} \)
Let's start with the left side. We compose the Taylor expansion (on the first equation) with (x+1):
\( {\color{Blue} {{^{(x+1)}a}=\sum_{n=0}^{\infty}{a_n.(x+1)^n}} \)
We use the binomial theorem
\(
\left.\begin{matrix}
(1+x)^n = \sum_{k=0}^n {n \choose k}x^k
\\
{n \choose k} = \frac{n!}{k!\,(n-k)!}
\end{matrix}\right\}
{^{(x+1)}a}=\sum_{n=0}^{\infty}{a_n.\sum_{k=0}^n {n \choose k}x^k}
\)
Each term on the summation is
\(
a_0.(1+x)^0 = a_0. \left({0 \choose 0}x^0 \right)\\
a_1.(1+x)^1 = a_1. \left({1 \choose 0}x^0 + {1 \choose 1}x^1\right)\\
a_2.(1+x)^2 = a_2. \left({2 \choose 0}x^0 + {2 \choose 1}x^1 + {2 \choose 2}x^2\right)\\
\vdots \\
a_n.(1+x)^n = a_n. \left({n \choose 0}x^0 + {n \choose 1}x^1 + {n \choose 2}x^2 + \cdots + {n \choose {n-1}}x^{n-1} + {n \choose n}x^n\right) \)
In resume:
\( {\color{Blue} {{^{(x+1)}a}=\sum_{n=0}^{\infty}\left({\sum_{m=n}^{\infty}{{m \choose n}.a_m}\right).x^n \)
\( ^{x+1}a=( a_{0} + a_{1} + a_{2} + a_{3} + a_{4} + a_{5} + a_{6} + a_{7} + a_{8} +a_{9} + a_{10} + ...)+\\
(a_{1} + 2a_{2} + 3a_{3} + 4a_{4} + 5a_{5} + 6a_{6} + 7a_{7} + 8a_{8} + 9a_{9} + 10a_{10} + ...)x +\\
(a_{2} + 3a_{3} + 6a_{4} + 10a_{5} + 15a_{6} + 21a_{7} + 28a_{8} + 36a_{9} + 45a_{10} + ...)x^2 +\\
(a_{3} + 4a_{4} + 10a_{5} + 20a_{6} + 35a_{7} + 56a_{8} + 84a_{9} + 120a_{10} + ...)x^3... \)
It makes a super nice linear system, with determinant=1.
I bet that it is well studied and known.
His inverse matrix is the same matrix, but each element has alternate sign.
Here is the inverse of the coefficients matrix, generated with PariGP
The inverse can be directly generated with m=matrix(Size, Size,n,i,(-1)^(n+i)*binomial(i-1,n-1)); if you want to test it with Sheldonison kneser.gp, use Size=60, and get the coefficients with ai=matrix(60,1,f,c,polcoeff(xsexp,f-1)).
Now, let's deal with the right side of the first equation:
This is the Taylor expansion for a^x:
\( a^x\,=\, \sum_{n=0}^{\infty}\frac{ln(a)^n}{n!}.x^n \)
Let's compose with the Taylor polynomial at the start of the post:
\( {\color{Red} {a^{^xa}\,=\, \sum_{n=0}^{\infty}\frac{ln(a)^n}{n!}.{\left(\sum_{m=0}^{\infty}{a_m .x^m}\right)}^n \)
(We should apply the multinomial theorem)
I made an expansion of a^x with 12 terms, replaced each x with a polynomial with 12 terms:
...and threw it to PariGp to expand. I got this 24Mb text file
http://s000.tinyupload.com/index.php?fil...4013545746
Since the series are truncated, I made some guess about what each term means:
\( a^{^xa}=
\left(\sum_{n=0}^{...}{\frac{lna^n}{n!}} \right) \,+\,
\left(a_1.\sum_{n=0}^{...}{\frac{lna^{n+1}}{n!}} \right).x \,+\,
\left(a_2 .\sum_{n=0}^{...}{\frac{lna^{n+1}}{n!}} \,+\,
a_1^2.\sum_{n=0}^{...}{\frac{lna^{2+n}}{2.n!}} \right).x^2\,+\, \\
\left(a_3 .\sum_{n=0}^{...}{\frac{lna^{n+1}}{n!}} \,+\,
a_1.a_2.\sum_{n=0}^{...}{\frac{lna^{2+n}}{n!}} \,+\,
a_1^3 .\sum_{n=0}^{...}{\frac{lna^{3+n}}{6.n!}} \right).x^3
\,+\,
\left(a_4 .\sum_{n=0}^{...}{\frac{lna^{n+1}}{n!}} \,+\,
a_2^2 .\sum_{n=0}^{...}{\frac{lna^{2+n}}{2.n!}} \,+\,
a_3.a_1 .\sum_{n=0}^{...}{\frac{lna^{2+n}}{n!}} \,+\,
a_2.a_1^2.\sum_{n=0}^{...}{\frac{lna^{3+n}}{2.n!}} \,+\,
a_1^4 .\sum_{n=0}^{...}{\frac{lna^{4+n}}{24.n!}}
\right).x^4 \,+\, ... \)
\( {\color{Red} {
a^{^xa}=
a \,+\,
lna.a.a_1 \,.\, x \,+\,
(lna.a.a_2 + \frac{lna^2.a}{2}.a_1^2) \,.\, x^2 \,+\,
(lna.a.a_3 + lna^2.a.a_1.a_2 +\frac{lna^3.a}{6}.a_1^3) \,.\, x^3 \,+\,
(lna.a.a_4 + lna^2.a.a_3.a_1 + \frac{lna^2.a}{2}.a_2^2 +\frac{lna^3.a}{2}.a_2.a_1^2+\frac{lna^4.a}{24}.a_1^4) \,.\,x^4 \\
+(lna.a.a_5 + lna^2.a.a_1.a_4 + {lna^2.a}.a_3.a_2 + \frac{lna^3.a}{2}.a_3.a_1^2 + \frac{lna^3.a}{2}.a_1.a_2^2+\frac{ lna^4.a}{6}.a_2.a_1^3+\frac{lna^5.a}{120}.a_1^5) \,.\,x^5 \\
+(lna.a.a_6 + {lna^2.a}.a_1.a_5 + lna^2.a.a_2.a_4 + \frac{lna^3.a}{2}.a_4.a_1^2 + \frac{lna^2.a}{2}.a_3^2 + { lna^3.a}.a_1.a_3.a_2 + \frac{lna^3.a}{6}.a_2^3 +\frac{lna^4.a}{4}.a_2^2 . a_1^2+\frac{lna^4.a}{6}.a_3.a_1^3+\frac{lna^5.a}{24}.a_2.a_1^4+\frac{lna^6.a}{ 720}.a_1^6) \,.\,x^6 \\
} \)
\( {\color{Red} {
+a.( lna.a_7 + lna^2. a_4.a_3 + lna^2.a_5.a_2 + \frac{lna^3}{2}a_3.a_2^2 + lna^2.a_6.a_1 + \frac{lna^3}{2}a_3^2.a_1 + lna^3a_4.a_2.a_1 + \frac{lna^4}{6}a_2^3.a_1 + \frac{lna^3}{2}a_5.a_1^2 + \frac{lna^4}{2} a_3.a_2.a_1^2 + \frac{lna^4}{6}a_4.a_1^3 + \frac{lna^5}{12} a_2^2.a_1^3 + \frac{lna^5}{24}a_3.a_1^4 + \frac{lna^6}{120}a_2.a_1^5 + \frac{lna^7}{5040}.a_1^7 )\,.\,x^7
} \)
EDIT: more terms:
\( {\color{Red} {
+a.(ln.a.a_8 + \frac{lna^2}{2} .a_4^2 + lna^2 .a_3.a_5 + lna^2 .a_2 .a_6 + \frac{lna^3}{2} .a_2 .a_3^2 + \frac{lna^3}{2} .a_2^2 .a_4 + \frac{lna^4}{24} .a_2^4 + lna^2 .a_1 .a_7 + lna^3 .a_1 .a_3 .a_4 + lna^3 .a_1 .a_2 .a_5 + \frac{lna^4}{2} .a_1 .a_2^2 .a_3 + \frac{lna^3}{2} .a_1^2 .a_6 + \\
\frac{lna^4}{4} .a_1^2 .a_3^2 + \frac{lna^4}{2} .a_1^2 .a_2 .a_4 + \frac{lna^5}{12} .a_1^2 .a_2^3 + \frac{lna^4}{6}.a_1^3 .a_5 + \frac{lna^5}{6} .a_1^3 .a_2 .a_3 + \frac{lna^5}{24} .a_1^4 .a_4 + \frac{lna^6}{48} .a_1^4 .a_2^2 + \frac{lna^6}{120} .a_1^5 .a_3 + \frac{lna^7}{720} .a_1^6 .a_2 + \frac{lna^8}{40320} .a_1^8
)\,.\,x^8
} \)
\( {\color{Red} {
+a.(lna.a_9 + lna^2.a_4 .a_5 + lna^2 .a_3 .a_6 + \frac{lna^3}{6} .a_3^3 + lna^2.a_2 .a_7 + lna^3.a_2 .a_3 .a_4 + \frac{lna^3}{2}.a_2^2 .a_5 + \frac{lna^4}{6} .a_2^3 .a_3 + lna^2.a_1 .a_8 + \frac{lna^3}{2} .a_1 .a_4^2 + \\
lna^3 .a_1 .a_3 .a_5 + lna^3 .a_1 .a_2 .a_6 + \frac{lna^4}{2} .a_1 .a_2 .a_3^2 + \frac{lna^4}{2} .a_1 .a_2^2 .a_4 + \frac{lna^5}{24}.a_1 .a_2^4 + \frac{lna^3}{2}.a_1^2 .a_7 + \frac{lna^4}{2}.a_1^2 .a_3 .a_4 + \frac{lna^4}{2} .a_1^2 .a_2 .a_5 + \frac{lna^5}{4}.a_1^2 .a_2^2 .a_3 + \frac{lna^4}{6} .a_1^3 .a_6 +\\
\frac{lna^5}{12} .a_1^3 .a_3^2 + \frac{lna^5}{6} .a_1^3 .a_2 .a_4 + \frac{lna^6}{36}.a_1^3 .a_2^3 + \frac{lna^5}{24}.a_1^4 .a_5 + \frac{lna^6}{24}.a_1^4 .a_2 .a_3 + \frac{lna^6}{120}.a_1^5 .a_4 + \frac{lna^7}{240} .a_1^5 .a_2^2 + \frac{lna^7}{720} .a_1^6 .a_3 + \frac{lna^8}{5040}.a_1^7 .a_2 + \frac{lna^9}{362880}.a_1^9 ).x^9
} \)
This last expression seems to be related to the Taylor expansion of \( \\[15pt]
{^xa} \) around x=1.
So, each summatory of coefficients, for each power of x, is in some loose sense, "the derivative" of the summatory of the precedent power of x (taking da/dx=lna.a.a₁, and daᵢ/dx=aᵢ₊₁)
I wrote "in some loose sense", because there are multiplicative constants missing.
For example, the "derivative" of lna.a.a₁ "should be" (lna.a.a₂+lna².a.a₁²), but instead is (lna.a.a₂+lna².a.a₁²/n!), where n is the power exponent of a₁ in that term.
again the "derivative" of (lna.a.a₂+lna².a.a₁²/2!) "should be" (lna.a.a₃ + 2 . lna².a.a₁.a₂ +lna³.a.a₁³/2), but instead is (lna.a.a₃ + lna².a.a₁.a₂ +lna³.a.a₁³/3!)
¿Can you find the right rule for the coefficients of each power of x?
Here is an excel worksheet, which takes the coefficients calculated with Sheldonison Kneser.gp program (for base a=Pi), and compares the blue and red equations:
http://s000.tinyupload.com/index.php?fil...7813385010
Unfortunately, this website does not allow to upload most format files, so I uploaded it o Tinyupload.com, but there is no guarantee of permanence, so it may be deleted in the future.
We want the coefficients aᵢ of this Taylor expansion:
\( {^xa}=\sum_{n=0}^{\infty}{a_n .x^n} \)
They should match this equation:
\( {\color{Blue} {^{x+1}a} \)\( \,=\, \)\( {\color{Red} a^{^xa}} \)
Let's start with the left side. We compose the Taylor expansion (on the first equation) with (x+1):
\( {\color{Blue} {{^{(x+1)}a}=\sum_{n=0}^{\infty}{a_n.(x+1)^n}} \)
We use the binomial theorem
\(
\left.\begin{matrix}
(1+x)^n = \sum_{k=0}^n {n \choose k}x^k
\\
{n \choose k} = \frac{n!}{k!\,(n-k)!}
\end{matrix}\right\}
{^{(x+1)}a}=\sum_{n=0}^{\infty}{a_n.\sum_{k=0}^n {n \choose k}x^k}
\)
Each term on the summation is
\(
a_0.(1+x)^0 = a_0. \left({0 \choose 0}x^0 \right)\\
a_1.(1+x)^1 = a_1. \left({1 \choose 0}x^0 + {1 \choose 1}x^1\right)\\
a_2.(1+x)^2 = a_2. \left({2 \choose 0}x^0 + {2 \choose 1}x^1 + {2 \choose 2}x^2\right)\\
\vdots \\
a_n.(1+x)^n = a_n. \left({n \choose 0}x^0 + {n \choose 1}x^1 + {n \choose 2}x^2 + \cdots + {n \choose {n-1}}x^{n-1} + {n \choose n}x^n\right) \)
In resume:
\( {\color{Blue} {{^{(x+1)}a}=\sum_{n=0}^{\infty}\left({\sum_{m=n}^{\infty}{{m \choose n}.a_m}\right).x^n \)
\( ^{x+1}a=( a_{0} + a_{1} + a_{2} + a_{3} + a_{4} + a_{5} + a_{6} + a_{7} + a_{8} +a_{9} + a_{10} + ...)+\\
(a_{1} + 2a_{2} + 3a_{3} + 4a_{4} + 5a_{5} + 6a_{6} + 7a_{7} + 8a_{8} + 9a_{9} + 10a_{10} + ...)x +\\
(a_{2} + 3a_{3} + 6a_{4} + 10a_{5} + 15a_{6} + 21a_{7} + 28a_{8} + 36a_{9} + 45a_{10} + ...)x^2 +\\
(a_{3} + 4a_{4} + 10a_{5} + 20a_{6} + 35a_{7} + 56a_{8} + 84a_{9} + 120a_{10} + ...)x^3... \)
It makes a super nice linear system, with determinant=1.
I bet that it is well studied and known.
His inverse matrix is the same matrix, but each element has alternate sign.
Here is the inverse of the coefficients matrix, generated with PariGP
Code:
gp > m=matrix(Size, Size,n,i,binomial(i-1,n-1))
(20:22) gp > m^-1
[1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1 1 -1
[0 1 -2 3 -4 5 -6 7 -8 9 -10 11 -12 13 -14 15 -16 17 -18 19
[0 0 1 -3 6 -10 15 -21 28 -36 45 -55 66 -78 91 -105 120 -136 153 -171
[0 0 0 1 -4 10 -20 35 -56 84 -120 165 -220 286 -364 455 -560 680 -816 969
[0 0 0 0 1 -5 15 -35 70 -126 210 -330 495 -715 1001 -1365 1820 -2380 3060 -3876
[0 0 0 0 0 1 -6 21 -56 126 -252 462 -792 1287 -2002 3003 -4368 6188 -8568 11628
[0 0 0 0 0 0 1 -7 28 -84 210 -462 924 -1716 3003 -5005 8008 -12376 18564 -27132
[0 0 0 0 0 0 0 1 -8 36 -120 330 -792 1716 -3432 6435 -11440 19448 -31824 50388
[0 0 0 0 0 0 0 0 1 -9 45 -165 495 -1287 3003 -6435 12870 -24310 43758 -75582
[0 0 0 0 0 0 0 0 0 1 -10 55 -220 715 -2002 5005 -11440 24310 -48620 92378
[0 0 0 0 0 0 0 0 0 0 1 -11 66 -286 1001 -3003 8008 -19448 43758 -92378
[0 0 0 0 0 0 0 0 0 0 0 1 -12 78 -364 1365 -4368 12376 -31824 75582
Now, let's deal with the right side of the first equation:
This is the Taylor expansion for a^x:
\( a^x\,=\, \sum_{n=0}^{\infty}\frac{ln(a)^n}{n!}.x^n \)
Let's compose with the Taylor polynomial at the start of the post:
\( {\color{Red} {a^{^xa}\,=\, \sum_{n=0}^{\infty}\frac{ln(a)^n}{n!}.{\left(\sum_{m=0}^{\infty}{a_m .x^m}\right)}^n \)
(We should apply the multinomial theorem)
I made an expansion of a^x with 12 terms, replaced each x with a polynomial with 12 terms:
Code:
(01:40) gp > exp(lna*x)
%1 = 1 + lna*x + 1/2*lna^2*x^2 + 1/6*lna^3*x^3 + 1/24*lna^4*x^4 + 1/120*lna^5*x^5 + 1/720*lna^6*x^6 + 1/5040*lna^7*x^7 + 1/40320*lna^8*x^8 + 1/362880*lna^9*x^9 + 1/3628800*lna^10*x^10 + 1/39916800*lna^11*x^11
Code:
1 + lna*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)
+1/2*lna^2*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^2
+1/6*lna^3*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^3
+1/24*lna^4*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^4
+1/120*lna^5*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^5
+1/720*lna^6*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^6
+1/5040*lna^7*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^7
+1/40320*lna^8*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^8
+1/362880*lna^9*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^9
+1/3628800*lna^10*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^10
+1/39916800*lna^11*(1+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+a8*x^8+a9*x^9+a10*x^10+a11*x^11)^11
...and threw it to PariGp to expand. I got this 24Mb text file
http://s000.tinyupload.com/index.php?fil...4013545746
Since the series are truncated, I made some guess about what each term means:
\( a^{^xa}=
\left(\sum_{n=0}^{...}{\frac{lna^n}{n!}} \right) \,+\,
\left(a_1.\sum_{n=0}^{...}{\frac{lna^{n+1}}{n!}} \right).x \,+\,
\left(a_2 .\sum_{n=0}^{...}{\frac{lna^{n+1}}{n!}} \,+\,
a_1^2.\sum_{n=0}^{...}{\frac{lna^{2+n}}{2.n!}} \right).x^2\,+\, \\
\left(a_3 .\sum_{n=0}^{...}{\frac{lna^{n+1}}{n!}} \,+\,
a_1.a_2.\sum_{n=0}^{...}{\frac{lna^{2+n}}{n!}} \,+\,
a_1^3 .\sum_{n=0}^{...}{\frac{lna^{3+n}}{6.n!}} \right).x^3
\,+\,
\left(a_4 .\sum_{n=0}^{...}{\frac{lna^{n+1}}{n!}} \,+\,
a_2^2 .\sum_{n=0}^{...}{\frac{lna^{2+n}}{2.n!}} \,+\,
a_3.a_1 .\sum_{n=0}^{...}{\frac{lna^{2+n}}{n!}} \,+\,
a_2.a_1^2.\sum_{n=0}^{...}{\frac{lna^{3+n}}{2.n!}} \,+\,
a_1^4 .\sum_{n=0}^{...}{\frac{lna^{4+n}}{24.n!}}
\right).x^4 \,+\, ... \)
\( {\color{Red} {
a^{^xa}=
a \,+\,
lna.a.a_1 \,.\, x \,+\,
(lna.a.a_2 + \frac{lna^2.a}{2}.a_1^2) \,.\, x^2 \,+\,
(lna.a.a_3 + lna^2.a.a_1.a_2 +\frac{lna^3.a}{6}.a_1^3) \,.\, x^3 \,+\,
(lna.a.a_4 + lna^2.a.a_3.a_1 + \frac{lna^2.a}{2}.a_2^2 +\frac{lna^3.a}{2}.a_2.a_1^2+\frac{lna^4.a}{24}.a_1^4) \,.\,x^4 \\
+(lna.a.a_5 + lna^2.a.a_1.a_4 + {lna^2.a}.a_3.a_2 + \frac{lna^3.a}{2}.a_3.a_1^2 + \frac{lna^3.a}{2}.a_1.a_2^2+\frac{ lna^4.a}{6}.a_2.a_1^3+\frac{lna^5.a}{120}.a_1^5) \,.\,x^5 \\
+(lna.a.a_6 + {lna^2.a}.a_1.a_5 + lna^2.a.a_2.a_4 + \frac{lna^3.a}{2}.a_4.a_1^2 + \frac{lna^2.a}{2}.a_3^2 + { lna^3.a}.a_1.a_3.a_2 + \frac{lna^3.a}{6}.a_2^3 +\frac{lna^4.a}{4}.a_2^2 . a_1^2+\frac{lna^4.a}{6}.a_3.a_1^3+\frac{lna^5.a}{24}.a_2.a_1^4+\frac{lna^6.a}{ 720}.a_1^6) \,.\,x^6 \\
} \)
\( {\color{Red} {
+a.( lna.a_7 + lna^2. a_4.a_3 + lna^2.a_5.a_2 + \frac{lna^3}{2}a_3.a_2^2 + lna^2.a_6.a_1 + \frac{lna^3}{2}a_3^2.a_1 + lna^3a_4.a_2.a_1 + \frac{lna^4}{6}a_2^3.a_1 + \frac{lna^3}{2}a_5.a_1^2 + \frac{lna^4}{2} a_3.a_2.a_1^2 + \frac{lna^4}{6}a_4.a_1^3 + \frac{lna^5}{12} a_2^2.a_1^3 + \frac{lna^5}{24}a_3.a_1^4 + \frac{lna^6}{120}a_2.a_1^5 + \frac{lna^7}{5040}.a_1^7 )\,.\,x^7
} \)
EDIT: more terms:
\( {\color{Red} {
+a.(ln.a.a_8 + \frac{lna^2}{2} .a_4^2 + lna^2 .a_3.a_5 + lna^2 .a_2 .a_6 + \frac{lna^3}{2} .a_2 .a_3^2 + \frac{lna^3}{2} .a_2^2 .a_4 + \frac{lna^4}{24} .a_2^4 + lna^2 .a_1 .a_7 + lna^3 .a_1 .a_3 .a_4 + lna^3 .a_1 .a_2 .a_5 + \frac{lna^4}{2} .a_1 .a_2^2 .a_3 + \frac{lna^3}{2} .a_1^2 .a_6 + \\
\frac{lna^4}{4} .a_1^2 .a_3^2 + \frac{lna^4}{2} .a_1^2 .a_2 .a_4 + \frac{lna^5}{12} .a_1^2 .a_2^3 + \frac{lna^4}{6}.a_1^3 .a_5 + \frac{lna^5}{6} .a_1^3 .a_2 .a_3 + \frac{lna^5}{24} .a_1^4 .a_4 + \frac{lna^6}{48} .a_1^4 .a_2^2 + \frac{lna^6}{120} .a_1^5 .a_3 + \frac{lna^7}{720} .a_1^6 .a_2 + \frac{lna^8}{40320} .a_1^8
)\,.\,x^8
} \)
\( {\color{Red} {
+a.(lna.a_9 + lna^2.a_4 .a_5 + lna^2 .a_3 .a_6 + \frac{lna^3}{6} .a_3^3 + lna^2.a_2 .a_7 + lna^3.a_2 .a_3 .a_4 + \frac{lna^3}{2}.a_2^2 .a_5 + \frac{lna^4}{6} .a_2^3 .a_3 + lna^2.a_1 .a_8 + \frac{lna^3}{2} .a_1 .a_4^2 + \\
lna^3 .a_1 .a_3 .a_5 + lna^3 .a_1 .a_2 .a_6 + \frac{lna^4}{2} .a_1 .a_2 .a_3^2 + \frac{lna^4}{2} .a_1 .a_2^2 .a_4 + \frac{lna^5}{24}.a_1 .a_2^4 + \frac{lna^3}{2}.a_1^2 .a_7 + \frac{lna^4}{2}.a_1^2 .a_3 .a_4 + \frac{lna^4}{2} .a_1^2 .a_2 .a_5 + \frac{lna^5}{4}.a_1^2 .a_2^2 .a_3 + \frac{lna^4}{6} .a_1^3 .a_6 +\\
\frac{lna^5}{12} .a_1^3 .a_3^2 + \frac{lna^5}{6} .a_1^3 .a_2 .a_4 + \frac{lna^6}{36}.a_1^3 .a_2^3 + \frac{lna^5}{24}.a_1^4 .a_5 + \frac{lna^6}{24}.a_1^4 .a_2 .a_3 + \frac{lna^6}{120}.a_1^5 .a_4 + \frac{lna^7}{240} .a_1^5 .a_2^2 + \frac{lna^7}{720} .a_1^6 .a_3 + \frac{lna^8}{5040}.a_1^7 .a_2 + \frac{lna^9}{362880}.a_1^9 ).x^9
} \)
This last expression seems to be related to the Taylor expansion of \( \\[15pt]
{^xa} \) around x=1.
So, each summatory of coefficients, for each power of x, is in some loose sense, "the derivative" of the summatory of the precedent power of x (taking da/dx=lna.a.a₁, and daᵢ/dx=aᵢ₊₁)
I wrote "in some loose sense", because there are multiplicative constants missing.
For example, the "derivative" of lna.a.a₁ "should be" (lna.a.a₂+lna².a.a₁²), but instead is (lna.a.a₂+lna².a.a₁²/n!), where n is the power exponent of a₁ in that term.
again the "derivative" of (lna.a.a₂+lna².a.a₁²/2!) "should be" (lna.a.a₃ + 2 . lna².a.a₁.a₂ +lna³.a.a₁³/2), but instead is (lna.a.a₃ + lna².a.a₁.a₂ +lna³.a.a₁³/3!)
¿Can you find the right rule for the coefficients of each power of x?
Here is an excel worksheet, which takes the coefficients calculated with Sheldonison Kneser.gp program (for base a=Pi), and compares the blue and red equations:
http://s000.tinyupload.com/index.php?fil...7813385010
Unfortunately, this website does not allow to upload most format files, so I uploaded it o Tinyupload.com, but there is no guarantee of permanence, so it may be deleted in the future.
I have the result, but I do not yet know how to get it.