[Exercise] A deal of Uniqueness-critrion:Gamma-functionas iteration
#1
I considered the iteration of the multiplication in the sense of the interpolation of the factorial. I adapted the idea of a truncation of a matrixoperator in the same way as Andy has described it with his slog-solution. The occuring powerseries are "automatically" that of the gamma/lngamma-function - it seems, that method provides a proper "uniqueness-criterion" implicitely (that's then also a backup for the possible justifications of Andy's method).

When I stuck at a certain problem I discussed this in the math.SE-forum, see the link here: math.stackexchange.com

I'll also put my results in my text on Uncompleting the Gamma (PDF) which I'll upload here if I'm finished as update of the provisorical chapter 4.

Gottfried
Gottfried Helms, Kassel
Reply
#2
(05/19/2011, 05:49 AM)Gottfried Wrote: I considered the iteration of the multiplication in the sense of the interpolation of the factorial. I adapted the idea of a truncation of a matrixoperator in the same way as Andy has described it with his slog-solution. The occuring powerseries are "automatically" that of the gamma/lngamma-function - it seems, that method provides a proper "uniqueness-criterion" implicitely (that's then also a backup for the possible justifications of Andy's method).    

When I stuck at a certain problem I discussed this in the math.SE-forum, see the link here: math.stackexchange.com

I'll also put my results in my text on Uncompleting the Gamma (PDF) which I'll upload here if I'm finished as update of the provisorical chapter 4.

Gottfried

The incomplete gamma function is one of the most interesting ones !!

I had the idea of removing the singularities of the sexp in a similar way as you removed the poles and zero's of the gamma.
Not sure if that will work. What do you think ??

If G(z) is the uncomplete ( or incomplete ) gamma function.

G(z) =/= 0 implies that there is an entire logG(z) analogue to loggamma !

Also logG(exp(z)) must be interesting.

Notice that all derivatives of G(z) has positive derivatives.
Im referring to fake function theory again yes.

I wonder how close G(z) , fakegamma(z) and fakeG(z) are.
And how many zero's the fake will have. Only finite zero's ??

Forgive me, maybe I said all that before.

the pxp analogue of G(z) is also interesting.

continu sums and products relate too.

A kind of superfactorial is also worth considering.

For the superfunction of G(z) we could perhaps try using the recent auxilliary function strategy :

F(s) = G( F(s-1) - exp(-s) )

So clearly it connects to many things here at the tetration forum.

regards

tommy1729
Reply
#3
(03/07/2021, 11:26 PM)tommy1729 Wrote:
(05/19/2011, 05:49 AM)Gottfried Wrote: I considered the iteration of the multiplication in the sense of the interpolation of the factorial. I adapted the idea of a truncation of a matrixoperator in the same way as Andy has described it with his slog-solution. The occuring powerseries are "automatically" that of the gamma/lngamma-function - it seems, that method provides a proper "uniqueness-criterion" implicitely (that's then also a backup for the possible justifications of Andy's method).    

When I stuck at a certain problem I discussed this in the math.SE-forum, see the link here: math.stackexchange.com

I'll also put my results in my text on Uncompleting the Gamma (PDF) which I'll upload here if I'm finished as update of the provisorical chapter 4.

Gottfried

The incomplete gamma function is one of the most interesting ones !!

I had the idea of removing the singularities of the sexp in a similar way as you removed the poles and zero's of the gamma.
Not sure if that will work. What do you think ??

If G(z) is the uncomplete ( or incomplete ) gamma function.

G(z) =/= 0 implies that there is an entire logG(z) analogue to loggamma !

Also logG(exp(z)) must be interesting.

Notice that all derivatives of G(z) has positive derivatives.
Im referring to fake function theory again yes.

I wonder how close G(z) , fakegamma(z) and fakeG(z) are.
And how many zero's the fake will have. Only finite zero's ??

Forgive me, maybe I said all that before.

the pxp analogue of G(z) is also interesting.

continu sums and products relate too.

A kind of superfactorial is also worth considering.

For the superfunction of G(z) we could perhaps try using the recent auxilliary function strategy :

F(s) = G( F(s-1) - exp(-s) )

So clearly it connects to many things here at the tetration forum.

regards

tommy1729

Another question : 

Is G(z) a superfunction ? of what ?

Now we can define F(z) := G( invG(z) + 1 ) and then plug in approximations of G and invG.
But I guess more can be done. 
Is F(z) entire and similar questions would be imho interesting.

Also finding tetration by infinite compositions such as f(z+1) = exp(f(z) - G(-z)) have crossed my mind.

Not only because G(-z) grows to 0 faster than exp(-z) but also because of its behavior on the complex plane.
Afterall G does not have a period and does not have a fetish for making copies.

A riemann surface for invG and similar would be nice to see.

Btw the (normal) inverse gamma function is often approximated by using lambertW.

I think I posted a formula somewhere here ...

regards

tommy1729
Reply
#4
Heh, this is fairly interesting. I'd like to give my two cents on what you are doing briefly by taking a detour through fractional calculus. This relates very largely to the indefinite sum, fractional iteration, differsums, and recurrence relations in general. It's largely why I got so into fractional calculus in the first place--which seems to me, you're stumbling upon one of the fundamental analytic continuations I use.

To begin, observe, one I consider, the simplest way to analytically continue the Gamma function.

\(
\Gamma(z) = \int_0^\infty e^{-x}x^{z-1}\,dx\\
\)

Break the integral into \( \int_0^1 + \int_1^\infty \). When we look at,

\(
\int_0^1 e^{-x}x^{z-1}\,dx = \sum_{n=0}^\infty \frac{(-1)^n}{n!} \int_0^1 x^{n+z-1}\,dx\\
= \sum_{n=0}^\infty \frac{(-1)^n}{n!(n+z)}\\
\)

And this converges everywhere in \( \mathbb{C}/\{-\mathbb{N}\} \). The second part of the integral,

\(
\int_1^\infty e^{-x}x^{z-1}\,dx\\
\)

Converges everywhere in the complex plane, giving us the analytic continuation of \( \Gamma(z) \) into \( \mathbb{C}/\{-\mathbb{N}\} \),

\(
\Gamma(z) = \sum_{n=0}^\infty \frac{(-1)^n}{n!(n+z)} + \int_1^\infty e^{-x}x^{z-1}\,dx\\
\)

It's difficult who to attribute this analytic continuation to, I have heard it attributed to Euler, so it goes back a fair way.

The second player to the game is Ramanujan, who although didn't provide an absolutely rigorous description of the following phenomena, but still displayed the correct answer in a lot of cases (we can now rigorously say where and when it works). This is called Ramanujan's master theorem. The manner in way Ramanujan would write it is very indicative of the general theory.

Let us consider \( E \) a type of operator (really a function, but this was the original language),

\(
e^{-Ex} = \sum_{n=0}^\infty E^n \frac{(-x)^n}{n!}\\
\)

And then note that,

\(
\int_0^\infty e^{-Ex}x^{z-1}\,dx = E^{-z} \Gamma(z)\\
\)

Which of course, is nothing but a formalism, which does display some truth. The correct way, as would write it now, is to write,

\(
e^{-Ex} = f(x) = \sum_{n=0}^\infty E(n)\frac{(-x)^n}{n!}\\
\)

And, Ramanujan wrote, sort of unjustified, that,

\(
\int_0^\infty f(x) x^{z-1}\,dx = E(-z)\Gamma(z)\\
\)

Now, this clearly won't work for arbitrary functions \( E \). It will only work for certain functions, and to determine which functions is surprisingly not that hard. Let us take a function \( E \) holomorphic for \( \Re(z) > -1 \), let us assume that \( |E(z)| \le C e^{\tau|\Im(z)|} \) for \( \tau < \pi/2 \) as \( |\Im(z)| \to \infty \). And to keep things simple, lets assume that \( |E(z)| \le C e^{\rho|\Re(z)|} \) as \( \Re(z) \to \infty \) for arbitrary \( \rho \).

Now write the inverse mellin transform as,

\(
f(x) = \int_{\frac{1}{2} - i\infty}^{\frac{1}{2} + i\infty} \Gamma(z)E(-z)x^{-z}\,dz\\
\)

Now, the poles of this integrand in the left half plane, are at the negative integers, and their residues are,

\(
\text{Res} (z=-n;\Gamma(z)E(-z)x^{-z}) = E(n)\frac{(-x)^n}{n!}\\
\)

Doing a standard exercise in contour integration, and thanks to the nice asymptotics of the Gamma function (think Sterling), and the exponential bounds on \( E \), we get,

\(
f(x) = \sum_{n=0}^\infty E(n) \frac{(-x)^n}{n!}\\
\)

Now, by way of the mellin transform (thanks to the fourier transform), we know instantly that \( f(x) = \mathcal{O}(1/x) \) and therefore,

\(
\int_0^\infty |f(x)|x^{\sigma - 1}\,dx < \infty\,\,\text{for}\,\,0 < \sigma < 1\\
\)

And so we can take the mellin transform, and by the invertible nature,

\(
E(-z)\Gamma(z) = \int_0^\infty f(x)x^{z-1}\,dx\,\,\text{for}\,\,0 < \Re(z) < 1\\
\)

Now, we enter in the above analytic continuation we did for the Gamma function (which works in the exact same manner as above), to arrive at,

\(
E(-z)\Gamma(z) = \sum_{n=0}^\infty E(n)\frac{(-1)^n}{n!(n+z)} + \int_1^\infty f(x)x^{z-1}\,dx\\
\)

Which is valid for \( \Re(z) <1 \). And upon here, we enter the third player, the Riemann-Liouville differintegral. Now there are a bunch of differintegrals, but the one we want is the Riemann-Liouville because it is really the mellin transform in disguise.  We write this, in a slightly nonstandard form, but it's really nothing but a change of variables. Here, we change our function of interest to remove the pesky alternation. Let us write,

\(
\vartheta(w) = \sum_{n=0}^\infty E(n) \frac{w^n}{n!}\\
\)

\(
\frac{d^{-z}}{dw^{-z}}\vartheta(w) = \frac{1}{\Gamma(z)} \int_0^\infty \vartheta(w-x)x^{z-1}\,dx\\
\)

Now, we have an equivalence here. I'm going to skip a bit, but bear with me. Assume that \( |\vartheta(w)| \) has decay of the form \( \mathcal{O}(|w|^{-b}) \) for \( |\arg(-w)| < \theta < \pi/2 \) as \( |w|\to\infty \) in this sector. Then there exists a function \( E(z) \) which is holomorphic for \( \Re(z) > -b \) and satisfies \( |E(z)| \le Ce^{(\pi/2 - \theta)|\Im(z)|} \) and \( E(z)\Gamma(-z) \to 0 \) in a fast enough manner to be integrated, and for the sum to be entire.

With this we have an isomorphism between differintegration, and very well behaved functions in the right half plane \( E \). Now, it's important to remember the rules of differintegration. This is that,

\(
\frac{d^{z}}{dw^z}\frac{d^{z'}}{dw^{z'}} = \frac{d^{z+z'}}{dw^{z+z'}}\\
\frac{d^n}{dw^n}\,\,\text{is just differentiation}\\
\frac{d^{-1}}{dw^{-1}} = \int_{-\infty}^w \vartheta(x)\,dx\\
\)

So that integration is inherently integration with the lower limit at negative infinity (this is important to ensure consistency). And it is hereupon where we can return to Ramanujan's notation a bit more rigorously. I'll take the example of fractional iteration, as it's a nice one to start with (though the same heuristic works for indefinite sum, indefinite product, and the sort). Let us write,

\(
\mathcal{C}_\phi g(\xi) = g(\phi(\xi))\\
\)

Which is a linear operator. Now, we write,

\(
e^{\mathcal{C}_\phi w} g(\xi) = \sum_{n=0}^\infty g(\phi^{\circ n}(\xi)) \frac{w^n}{n!}\\
\)

Let's assume that this is differintegrable (which is the real challenge, in showing this happens). For the case that \( \phi \) has a fixed point with multiplier \( 0 < \lambda < 1 \) this is doable, thanks to Schroder's equation (and is doable in other scenarios as well). And using Ramanujan's abuse of notation, we write,

\(
\frac{d^z}{dw^z} e^{\mathcal{C}_\phi w} = \mathcal{C}_\phi^z e^{\mathcal{C}_\phi w}\\
\)

And, let us combine this with the semi-group property of the differintegral,

\(
\mathcal{C}_\phi^{z+z'} e^{\mathcal{C}_\phi w}=\frac{d^{z+z'}}{dw^{z+z'}}e^{\mathcal{C}_\phi w} = \frac{d^z}{dw^z}\frac{d^{z'}}{dw^{z'}} e^{\mathcal{C}_\phi w} = \mathcal{C}_\phi^z \mathcal{C}_\phi^{z'} e^{\mathcal{C}_\phi w}\\
\)

Which tells us the semi-group property of the differintegral preserves when iterating operators! Thus, all we really need is to prove the auxiliary function is differintegrable to express it in a manner that only depends on its natural iterates. We use \( E^n \) to create \( E^z \); so long as we have a differintegrable condition (which proves very hard to guarantee from first principles unfortunately).  Also, next to this, we are guaranteed a uniqueness condition. I refer to this as Ramanujan's identity theorem, a much stronger version was proven by Carlson (earning him a fields medal)--but we don't need the strength of Carlson.

Assume that \( E \) is one of these well bounded functions. And let's assume that \( E(n) = 0 \) for all \( n \in \mathbb{N} \). Well by the isomorphism, and the interchange, our function,

\(
e^{Ew} = \sum_{n=0}^\infty E(n) \frac{w^n}{n!} = 0\\
\)

And therefore, when taking the differintegral, we prove that \( E= 0 \). Which implies that if \( E(n) = 0 \), is bounded in this nice manner, it must be zero everywhere. From this, we can derive, if we have two functions \( E_1 \) and \( E_2 \) bounded in the appropriate manner, and satisfy \( E_1(n) = E_2(n) \) then \( E_1 - E_2 \) is bounded the same and is zero on the naturals, therefore \( E_1 - E_2 = 0 \) everywhere and \( E_1 = E_2 \). Which lets us conclude we have a uniqueness condition on this space.

This is infinitely helpful when talking about iterations. For the indefinite sum, for instance, all we need is that it's bounded in this manner--it's unique. For the fractional iteration, for instance, all we need is that it's bounded in this manner--it's unique. Etc etc...

And here is where we come to what I think you are beginning to uncover. Though, you are looking at the Gamma function, and rather than the Gamma function as multiplier. Let us take a function \( f \) which exists in this bounded space, and let us assume that \( f(n) =1 \) for the natural numbers. Then we are guaranteed that \( f=1 \) everywhere, because it is in this space. Now, this is actually a uniqueness condition for the Gamma function in disguise. Where we are instead declaring the uniqueness of the equation,

\(
1 = \frac{1}{\Gamma(z)} \int_0^\infty e^{-x}x^{z-1}\,dx\\
\)

Which appears trivial, but it really isn't. Because there can be no function \( f \) and \( \vartheta \) such that,

\(
f(-z) = \frac{1}{\Gamma(z)} \int_0^\infty \vartheta(-x)x^{z-1}\,dx
\)


where \( f(n) =1 \) and \( f\neq 1 \).


Now how this all ties to Andy's slog, seems to be that we are going to work in the same manner. Now I can't speak for Andy's slog (I do not know enough), but if it's holomorphic for \( \Re(z) > -1 \) and it satisfies these bounds (which makes perfect sense because of how slow the slog function grows). We would declare a uniqueness condition as,

\(
\text{slog}(z) = \frac{d^z}{dw^z}|_{w=0} \sum_{n=0}^\infty \text{slog}(n)\frac{w^n}{n!}\\
\)

Which means there is only one function in this space which interpolates the values of Andy's slog at \( n \in \mathbb{N} \)--and this solution happens to solve the slog equation. Now I do not know if this works, because I don't know where Andy's slog is holomorphic (or if it even is). I do not know if any slog will work tbh, because I've never studied their domains of holomorphy. But the bounds shouldn't be a problem, just whether they're holomorphic in a right half plane. This doesn't help constructively much, because finding the values \( \text{slog}(n) \) are rather impossible without solving \( \text{slog} \) first.

This means, the manner you are using, truncating matrices and approximating solutions--may be equivalent to the Ramanujan method using the differintegral. It certainly is using the gamma function. Which may be helpful for me to re-explain. Assume that,

\(
g(z)
\)

interpolates the factorial and satisfies \( zg(z) = g(z+1) \) Then, additionally assume that \( \frac{g(z)}{\Gamma(z)} \) is in this nice space, then \( \frac{g(z)}{\Gamma(z)} = 1 \). Now, this may seem trivial, again. But the manner in which you are breaking up the Gamma function in your paper, is actually pretty much all of the above in disguise.

I hope this might be helpful. Anyway, had fun writing it. Least of all, it should explain a lot of the work I've been showing you lately.

Regards, James.


EDIT:

I also thought it may be helpful to describe how the functional relationship of Indefinite Summation can be interpreted in the fractional calculus method. Let us assume that \( f \) is in this bounded space; and define,

\(
\vartheta(w) = \sum_{n=0}^\infty f(n) \frac{w^n}{n!}\\
\Theta(w) = \sum_{n=0}^\infty \sum_{j=0}^n f(j) \frac{w^n}{n!}\\
\)

Now, these functions satisfy the relationship,

\(
\frac{d}{dw} \Theta(w) = \Theta(w) + \frac{d}{dw}\vartheta(w)\\
\)

Now, we know that,

\(
f(z) = \frac{d^z}{dw^z}\vartheta(w)|_{w=0}\\
\)

And, as I showed in my paper that I posted in the other thread, the function \( \Theta \) is also differintegrable. Therefore define,

\(
F(z) = \frac{d^z}{dw^z}\Theta(w) |_{w=0}\\
\)

Then,

\(
F(z+1) = \frac{d^{z+1}}{dw^{z+1}}\Theta(w) |_{w=0}= \frac{d^z}{dw^z} \frac{d}{dw}\Theta(w) |_{w=0}\\
= \frac{d^z}{dw^z} \Theta(w) + \frac{d}{dw}\vartheta(w)|_{w=0}\\
= F(z) + \frac{d^{z+1}}{dw^{z+1}}\vartheta(w)|_{w=0}\\
=F(z) + f(z+1)\\
\)

Which is the functional equation of indefinite summation. We can do this with many functional equations; and it seems pretty universal in its effectiveness. The same procedure works with indefinite products (but it's a tad less effective). It works with Matrix operators as well, if we consider \( A \) a diagonalizable square nxn matrix; and assuming we have a certain kind of distribution of the eigenvalues (I don't know the word for it, I can write it out if you want), then the iteration,

\(
A^z e^{Aw}= \frac{d^z}{dw^z} e^{Aw}\\
\)

Is just as convergent in the matrix space. Giving a rather quick way of iterating matrices... This may help in your truncating matrix method and iterating them. But, don't ask me how--it just occurred to me it might help.
Reply
#5
Hi James -

 this seems really great what, and how, you're explaining/discussing that things. Even it takes me in derivations, which I'd always avoided due to lack of original/basic insight (my formal math-edication has been limited to linear algebra, while I completely missed any progress in the analysis-course, for students of "general informatik" in the early 1970ies). But how you expose consideration & arguing is really a delight, and perhaps I'll find the energy to step into it even deeper. Unfortunately since my retiring two years ago my mental energy for those things degrades already and I think I'll have to accept, that not many new things shall get entry in my mind... But really, your style of discourse is attractive again.     
I had tried to step into this fractional calculus some years ago, did some exercises, but when it came to the "half-derivative of zeta at zero" I gave up (even had tried to get help on this in MSE; that has been 7 years ago, I've been much fitter...). I'll surely read your answer again with hope to be able to catch up some more insight, let's see. And nice to see, that one gets up this explorative adventure with the "uncompleting the gamma" with a qualified access at all. <heart-warming> ... :-)    

At the moment I'm absorbed with the other problem (discussing x --> x-1/x) and properties of periodic points im comparision with what I've learned about periodic points of the exp(z) function, and writing up an essay on my explorations and discoveries. The occuring problem to make a vision of a continuous iteration for instance through a 3-periodic set in the exp(z) compared with a 3-periodic set in the x-->x-1/x problem is absorbing me most massively at the moment.        

Just wanted -at least- respond to your engagement timely enough to show that I'm valuating it very much!

Gottfried
Gottfried Helms, Kassel
Reply
#6
Hey, so I read your post trying to find the half-derivative of zeta. It's very important to remember in fractional calculus that the differintegral always depends on a base point.

The standard form of any differintegral is written as,

\(
I_a^{-z} f(x) = \frac{1}{\Gamma(z)}\int_a^x f(t)(x-t)^{z-1}\,dt
\)

Now each, \( I_a \) provides different solutions for each value \( a \). Particularly, the Gamma function derivative of monomials works as,

\(
I_a^{-z} (x-a)^n = \frac{\Gamma(n+1)}{\Gamma(n+z + 1)}(x-a)^{n+z}\\
\)

But, for \( b \neq a \),

\(
I_a^{-z} (x-b)^n \neq \frac{\Gamma(n+1)}{\Gamma(n+z + 1)}(x-b)^{n+z}\\
\)

Additionally, we get a lot of problems when considering transcendental functions for \( a \) finite. Which is epitomized by,

\(
I_a^{-z} e^{x} \neq e^x\\
\)

To get the above identity, we use Riemann-Liouville's differintegral, which is when we set \( a = -\infty \) (or more accurately, we set \( a = \infty \) the point on the Riemann Sphere). This differintegral will not converge on monomials without A LOT of massaging and discussion of distribution theory (something I think I wasted way too much time on). But, if we write,

\(
\frac{d^z}{dw^z} \vartheta(w) = \frac{1}{\Gamma(-z)}( \sum_{n=0}^\infty \vartheta^{(n)}(w)\frac{(-1)^n}{n!(n-z)} + \int_1^\infty \vartheta(w-x)x^{-z-1}\,dx)\\
\)

This is probably the best analytic continuation of the Riemann Liouville I came up with. And if you note,

\(
\frac{d^z}{dw^z} e^{\lambda w} = \lambda^z e^{\lambda w}\\
\)

Which means it's the differintegral which preserves the exponential. Even if \( \lambda \) is complex, we can choose different branches of \( \lambda^z \) which corresponds to different branches of the Riemann-Liouville integral (which is sort of like integrating in different directions and rotating the direction collecting powers of \( e^{2\pi i k z} \)).

Now, as to the zeta function, you have two choices. Expand the zeta function into a Taylor series, and differintegrate term wise. Or consider the zeta function as a transcendental function (which is a much better choice in my opinion).

The trouble is, the zeta function has decay as \( \Re(s) \to \infty \) rather than \( \Re(s) \to -\infty \). This isn't much of a problem really, as we just rotate the integration and add a exponential factor.

\(
\frac{d^{-z}}{dw^{-z}}\zeta(w) = \frac{e^{(2k+1)\pi i z}}{\Gamma(z)}\int_0^\infty \zeta(w+x)x^{z-1}\,dx\\
\)

For some \( k \in \mathbb{Z} \). This will almost equal the expected result if it weren't for the pesky first term of the zeta function, the constant \( 1 \). Let's ignore it briefly and write,

\(
\frac{d^{z}}{dw^{z}}\zeta(w) -1= \sum_{n=2}^\infty \frac{e^{(2k+1)\pi i z}\log(n)^z}{n^{w}}\\
\)

To uncover the value at \( 1 \) we have to use the above analytic continuation, this is,

\(
\frac{d^z}{dw^z} 1 = \frac{1}{\Gamma(-z)}(\sum_{n=0}^\infty 1^{(n)} \frac{(-1)^n}{n!(n-z)} + \int_1^\infty x^{-z-1}\,dx)\\
= \frac{1}{\Gamma(-z)}( \frac{-1}{z} + \frac{1}{z}) = 0\,\,\text{for}\,\,z\neq 0\\
\)

Which is where we have to use pesky distribution theory. If \( \delta(0) = 1 \) and \( \delta(z) = 0 \) for all other \( z \in \mathbb{C} \), then we can conclude,

\(
\frac{d^z}{dw^z} \zeta(w) = \delta + \sum_{n=2}^\infty \frac{e^{(2k+1)\pi i z}\log(n)^z}{n^{w}}\\
\)

Which tells us that \( \frac{d^z}{dw^z} \zeta(w) \) is holomorphic for \( \Re(w) > 1 \) and \( \Re(z) > 0 \) (which avoids all that pesky \( \delta \) nonsense).

Now, you're safest bet to uncover the fractional derivative at zero using the Riemann Liouville differintegral is to analytically continue this Dirichlet series (which isn't hard at all considering the complexity of the dirichlet coefficients). This will give you a function,

\(
\frac{d^z}{dw^z}\zeta(w) \,\,\text{is holomorphic for}\,\,w \neq 1\,\Re(z) > 0\\
\)

Remembering there is a choice of which branch you are using, and there is a different differintegral for each \( k \in \mathbb{Z} \).

Now the other way, is fairly simpler, but I absolutely hate it for power series. Largely because it will not behave well on the prototypical power series \( e^w \). This is much plainer and it is just differintegrating termwise (which is boring).

Anywho, that's my take on the fractional derivative of \( \zeta \).

Regards, James

EDIT: It's important to remember that the Riemann-Liouville integral is largely meant for when our function has decay to zero in some direction, as it goes to infinity. If the function decays to a constant, like how \( \zeta \to 1 \) then the differintegration can't be holomorphic at \( z=0 \), we will have a discontinuity here. This also happens if we have \( \vartheta(w) = \mathcal{O}(w^k) \) for monomials, which would mean the differintegration is holomorphic for \( \Re(z) > k \). I thought it may be helpful to write out the formula for the differintegrations of monomials using Riemann-Liouville.

\(
\Gamma(-z)\frac{d^z}{dw^z}w^k = \sum_{n=0}^\infty \frac{k!}{(k-n)!} w^{k-n} \frac{(-1)^n}{n!(n-z)} + \int_1^\infty (w-x)^{k}x^{-z-1}\,dx\\
= \sum_{n=0}^k \binom{k}{n} (-1)^n \frac{w^{k-n}}{n-z} + \sum_{n=0}^k \binom{k}{n} (-1)^{n} \frac{w^{k-n}}{z-n}\\
\)

And again we have to use a type of distribution where \( \frac{d^z}{dw^z}w^k = \frac{k!}{(k-n)!}w^{k-n} \) for \( z = n\,\,0 \le n \le k \) and it equals zero everywhere else. Which is a perfectly natural interpolation of \( \frac{k!}{(n-k)!}w^{k-n} \) but looks a lot uglier than the typical formula for a monomial. But, when you enter the distribution field of fractional calculus this formula can be very helpful, especially when correcting errors, and interpreting things like \( w^k + e^w \)--it essentially gives us justification for throwing away \( w^k \), especially when we let \( \Re(z) > k \). It preserves the semi-group property and the relationship \( \frac{d}{dw}\frac{d^z}{dw^z} = \frac{d^{z+1}}{dw^{z+1}} \); however it really only does this trivially. This is rather unconventional--and is because we are analyticially continuing the operator as a semi-group and sort of ignoring the differential relationship--I haven't investigated it much but it's nice. So it disagrees, somewhat, with what we expect a differintegral to do--especially for negative real values (when we fractionally iterate the integral); because, well, it doesn't integrate at all--it just nulls it.

The alternative way to do this, is to no longer insist that \( \frac{d^z}{dw^z} \) interpolates differentiation necessarily, but is instead just a semigroup. And part of the kernel of this semi-group is the space of all polynomials,

\(
\frac{d^z}{dw^z} : \mathbb{C}[w] \to \{0\}\\
\)

I like this interpretation because it essentially restricts the Riemann-Liouville differintegral to non-polynomials; and only works with functions expressible by power-series with an infinite number of terms. Both interpretations are equally valid though; largely because in the distribution sense we are only worried about a measure zero set of values (so we can throw it away almost everywhere); they're equivalent almost everywhere.

Non-integer monomials behave nice though, thought I'd add this too.

\(
\Gamma(-z)\frac{d^z}{dw^z}\sqrt{w} = \sum_{n=0}^\infty \frac{\Gamma(\frac{1}{2}+1)}{\Gamma(\frac{1}{2}+1-n)} w^{\frac{1}{2}-n} \frac{(-1)^n}{n!(n-z)} + \int_1^\infty \sqrt{w-x}x^{-z-1}\,dx\\
= \sum_{n=0}^\infty \binom{1/2}{n} w^{\frac{1}{2}-n}\frac{(-1)^n}{n-z} + \sum_{n=0}^\infty \binom{1/2}{n} w^{\frac{1}{2}-n}\frac{(-1)^n}{z-n}\\
=0\\
\)

Which means even more functions are in the kernel. This produces some trouble when talking about the logarithm, but only if we think about it dramatically. Remembering these are semigroups, we still have,

\(
\Gamma(-z)\frac{d^z}{dw^z} \log(w) = \sum_{n=0}^\infty \log^{(n)}(w) \frac{(-1)^n}{n!(n-z)} + \int_1^\infty \log(w)w^{-z-1}\,dw\\
\)

Which is holomorphic for \( \Re(z) > 0 \) and is a semi-group. However we lose right side associativity, \( \frac{d^z}{dw^z} \frac{d}{dw} = \frac{d^{z+1}}{dw^{z+1}} \). This isn't a problem, we just have to remember to be very careful once we analytically continue the Riemann-Liouville integral. The correct manner of interpreting this, is that MOST OF THE TIME \( \frac{d^z}{dw^z} = \frac{d^n}{dw^n} \) when \( z=n \) but NOT ALWAYS. Especially once the differintegral is extended to a larger space of functions. However, the semi-group property is still conserved. And the reason it works for the logarithm, is because we are NOT performing/allowing the differintegral identity \( \frac{d^{z-1}}{dw^{z-1}} \frac{1}{w} = \frac{d^z}{dw^z}\log(w) \) (recall that \( \frac{d^{-1}}{dw^{-1}} \) is integration with initial point at infinity, which is infinite for \( \frac{1}{w} \), this is meaningless no matter how we slice it).

So try to remember that \( \frac{d^z}{dw^z} \) is a semi-group acting on a space of integrable functions at infinity (with certain decay conditions)--and in this space interpolates the derivative. When we extend the operator (analytically continue the operator), it is still a semi-group but it no longer necessarily interpolates the derivative. For all of our purposes this is great news. We don't really care that it interpolates the derivative as much as we care that the semi-group property holds.

This whole edit (which I must've done ten or twelve times so far, lol) is an argument that The Riemann Liouville integral can be analytically continued and has a very robust kernel. Also, it's meant as an argument, that this function is infinitely more important/effective for transcendental/infinite power series functions. As a nice end to end, not all algebraic functions suffer this problem.

\(
\Gamma(-z) \frac{d^z}{dw^z} \frac{1}{1-w} = \int_0^\infty \frac{x^{-z-1}}{1-w+x}\,dx\\
\)

Equals a transcendental function and is not \( 0 \) everywhere. It equals something with the \( \sin \) function I can't remember now, but,

\(
\int_0^\infty \frac{x^{-z-1}}{1+x}\,dx =\frac{-\pi}{\sin\pi z}\\
\)

And we can use a power series expansion to derive the general formula.

\(
\frac{d^z}{dw^z} \frac{1}{1-w} = \sum_{n=0}^\infty \lim_{q\to n} (\frac{-\pi}{\Gamma(-z-q)\sin\pi (z+q)}) \frac{w^n}{n!}\\
\)

This has a closed form but I'm too lazy to go through it right now, lol. But it does equal \( \frac{n!}{(1-w)^{n+1}} \) for \( z = n \).

PS:

I'm sorry to hear you're having problems doing math like you used to. That must be very hard, and I can't even imagine. I appreciate you reading what I've written. I'd say just keep at the math, because even if it's hard, it's the type of thing that can only help your brain Smile.

Ironically enough, I'm rather horrible with matrices and linear algebra (at best, I have a working rudimentary knowledge), analysis and integrals and the sort we're always my strong suit. It's always lovely to see matrix people coming to similar results as analysis people. I like to think of the history of Schrodinger's mathematical formalism and Heisenberg's mathematical formalism; one with complex analysis, the other with matrices; and the eventual equivalence proven by Neumann. When I read what you wrote in the initial post I found it strikingly similar to the uniqueness condition imposed by Ramanujan's master theorem; especially when you pulled out the incomplete Gamma function. It took me a very long time to make sense of the equations, and the exact language needed. But if,

\(
f(z) = \frac{1}{\Gamma(-z)}(\sum_{n=0}^\infty f(n) \frac{(-1)^n}{n!(n-z)} + h(s))\\
\)

For a holomorphic function \( h \). There is only ONE function \( h \) which makes our function \( f \) belong to the space we want. Which, if we change \( h \) to \( \tilde{h} \) a holomorphic function, to get a function \( \tilde{f} \); we still get \( \tilde{f}(n) = f(n) \), but \( \tilde{f} \neq f \). This implies that \( \tilde{f} \) is not in our nice space. And that there is only one function \( h \) which gives us our uniqueness condition. And this function happens to be the incomplete Mellin transform/(incomplete Gamma function). And then, the added benefit is that in this space it must preserve the semi-group property! And the semi-group property of the differintegral is easily translatable into functional equations!
Reply
#7
see also : 

https://math.eretrandre.org/tetrationfor...p?tid=1309

where I SERIOUSLY consider this uncompleted gamma for tetration.
Reply


Possibly Related Threads…
Thread Author Replies Views Last Post
  A related discussion on interpolation: factorial and gamma-function Gottfried 9 20,364 07/10/2022, 06:23 AM
Last Post: Gottfried
  Uniqueness of fractionally iterated functions Daniel 7 2,117 07/05/2022, 01:21 AM
Last Post: JmsNxn
  Universal uniqueness criterion? bo198214 57 121,330 06/28/2022, 12:00 AM
Last Post: JmsNxn
  Fractional iteration of x^2+1 at infinity and fractional iteration of exp bo198214 17 34,020 06/11/2022, 12:24 PM
Last Post: tommy1729
  A question concerning uniqueness JmsNxn 4 11,010 06/10/2022, 08:45 AM
Last Post: Catullus
  [exercise] fractional iteration of f(z)= 2*sinh (log(z)) ? Gottfried 4 3,952 03/14/2021, 05:32 PM
Last Post: tommy1729
  A conjectured uniqueness criteria for analytic tetration Vladimir Reshetnikov 13 27,841 02/17/2017, 05:21 AM
Last Post: JmsNxn
  Uniqueness of half-iterate of exp(x) ? tommy1729 14 34,529 01/09/2017, 02:41 AM
Last Post: Gottfried
  Removing the branch points in the base: a uniqueness condition? fivexthethird 0 3,812 03/19/2016, 10:44 AM
Last Post: fivexthethird
  Tommy's Gamma trick ? tommy1729 7 15,270 11/07/2015, 01:02 PM
Last Post: tommy1729



Users browsing this thread: 1 Guest(s)