Sunday Function

Yes, it's Tuesday. Busy weekend, including a Relay for Life (you should find your local one and donate!). Today we'll get back on schedule a bit and make this one a rather more utilitarian than usual Sunday Function, playing into the last step we'll need in order to complete the Bose-Einstein condensation discussion.

In the course of our journey we will meet this somewhat alarming creature:


It's a function that takes two parameters: a more traditional variable z, and a number (lowercase Greek letter nu) denoting the particular subspecies of function g. The capital Greek letter gamma is itself a function.

We're only going to have to deal with this for nu = 3/2, so we may as well plot the stuff under the integral sign. Taking z = 1/2 just for the random heck of it:


As z changes that curve will change shape and correspondingly the area under the curve will change as well. Is there some z that will maximize that area? Sometimes you can save yourself trouble if you just plot the graph, and you can save the rigorous proof work for later. We'll do that. This is the actual graph of g(z):


It attains a maximum at z = 1, and we shall see that in our case the physical meaning of z requires that it not be able to go above 1 anyway.

Mystery upon mystery at the moment, I know. Unless you've done this before, of course. But we'll put it all together shortly and I think it'll be worth the wait.

More like this

Let's start with a pretty simple function. It's not this week's official Sunday Function, but we'll use it to get there. Take the number 1, divide it by x. Pretty easy. Now imagine putting a dot at the point x = 1 on the x axis, and draw a line straight up from there to the curve formed by…
We've done a lot of discussion of the concept of integrals of a function here on this blog. Their definitions and applications are so broad as to defy any one-sentence description, but one of the most basic is the idea of the area under a curve. More precisely it's the signed area under the curve…
Let's do two functions today. As sometimes happens, in this case we're not so interested in the functions themselves as the fact that these functions happen to be part of a general class of functions. Just as we can classify the real numbers as even, odd, or neither (numbers like pi, 1/2, and the…
Here's a very simple function: You give it a positive real number, it gives you the square root. It maps 81 to 9, 100 to 10, 2 to 1.414..., and so on and so forth. It's pretty much the only one-argument function that's built into most pocket calculators, which says something about its utility…

The apparent reason for the restriction on z is obvious: for z > 1 you have a simple pole in that integrand. Does the function have an analytic continuation, or is it completely undefined?

The z = 1 case also has a restriction to ν > 1, in order to make the singularity integrable. Without seeing the physical setting in which this function arises, I can't say what that physically means, but you say that we will get the specific version of this function with ν = 3/2.

By Eric Lund (not verified) on 28 Apr 2009 #permalink

I plugged in z= 1/2 and nu = 3/2 and did the integral and found g(1/2) = .624. After doing the integral the dummy variable x disappears. SO you write that you are using z = 1/2 to plot this but isn't z the variable here?
For z = 1/2, g(z) = a constant.
What exactly are you plotting g vs x? or g vs z?
How command do i use to plot this mathematica?

By Just being dumb here (not verified) on 28 Apr 2009 #permalink

Just to make it clear. I understand how you get the second graph where g->1 as z->infinity.
But I just don't understand the first plot.
The integral of g(z) when nu = 3/2 is
PolyLog[3/2,z] and plotting this from 0 to 5 gives the second graph. But what about the first? when z = 1/2 g(z) should be a constant and you get a lognormal-type graph.

By just being dumb here (not verified) on 28 Apr 2009 #permalink

We integrate a function of x and evaluate that function with the appropriate z values to get g(z). So what I've plotted in the first graph is just the function we're integrating to get g(z). This, in other words:

1/Gamma[3/2] x^((3/2) - 1)/((1/2)^-1 (Exp[x] - 1))

You have rediscovered social engineering! The first 40% is interpolated linear, the next 30% is interpolated linear with a slightly larger slope. The last 30% extrapolated, legislated, and enforced by jackbooted State compassion is observed as disaster. Economics wll espy your exposition and weep tears of self-realization, as will Baby Boomers expecting Social Security and Medicare.

Like infinity, centrally administered social engineering gets real exciting toward the end.

g(z) looks a lot like the riemann zeta function (in one of its forms). Is that just a coincidence, or is there something more significant in that it looks like zeta(nu) (with a 1/z in the denominator of the integrand)? And is there a name for g(z)?

By Anonymous (not verified) on 28 Apr 2009 #permalink