The Unapologetic Mathematician

Mathematics for the interested outsider

The Integral as a Function of the Interval

Let’s say we take an integrator \alpha of bounded variation on an interval \left[a,b\right] and a function f that’s Riemann-Stieltjes integrable with respect to \alpha over that interval. Then we know that f is also integrable with respect to \alpha over the subinterval \left[a,x\right]\subseteq\left[a,b\right]. Let’s use this to define a function F on \left[a,b\right] by

\displaystyle F(x)=\int\limits_{\left[a,x\right]}fd\alpha

We can immediately say some interesting things about this function. First of all, F is, like \alpha, of bounded variation. Next, wherever \alpha is continuous, so is F. Finally, if \alpha is increasing, then F is differentiable wherever \alpha is differentiable and f is continuous. At such points, we have F'(x)=f(x)\alpha'(x). Notice that, as usual, the first two results will hold if we show that they hold for increasing integrators.

These results are similar to those we get from the Fundamental Theorem of Calculus, and we can use some of the same techniques to prove them. In particular, we call on the Integral Mean-Value Theorem for Riemann-Stieltjes integrals. If we take points x and y in \left[a,b\right], this tells us that

\displaystyle F(y)-F(x)=\int\limits_{\left[a,y\right]}fd\alpha-\int\limits_{\left[a,x\right]}fd\alpha=\int\limits_{\left[x,y\right]}fd\alpha=f(x_0)\left(\alpha(y)-\alpha(x)\right)

where x_0 is some point between x and y.

Now we let M be the supremum of |F| on \left[a,b\right]. For any partition of \left[a,b\right] we calculate the variation

\displaystyle\sum\limits_{i=1}^n|F(x_i)-F(x_{i-1})|\leq M\sum\limits_{i=1}^n|\alpha(x_i)-\alpha(x_{i-1})|\leq MV_\alpha(a,b)

thus giving an upper bound to the variation of F.

Similarly, let x be a point where \alpha is continuous. Then given an \epsilon>0 we can find a \delta so that |y-x|<\delta implies that |\alpha(y)-\alpha(x)|<\frac{\epsilon}{M}. Thus |F(y)-F(x)|<\epsilon, and we see that F is continuous at x as well.

Finally, we can divide by y-x to find

\displaystyle\frac{F(y)-F(x)}{y-x}=f(x_0)\frac{\alpha(y)-\alpha(x)}{y-x}

Then as y gets closer to x, x_0 gets squeezed in towards x as well. If f is continuous at x and \alpha is differentiable there, then the limit of this difference quotient exists, and has the value stated.

March 31, 2008 Posted by John Armstrong | Analysis, Calculus | | 2 Comments

Graded Rings in Grade School

Over at the new, improved Mathematics Under the Microscope (off of Blogger, onto WordPress), Alexandre makes an interesting point. word problems take place in a graded ring.

March 31, 2008 Posted by John Armstrong | Uncategorized | | 1 Comment

Sunday Samples 62

Sting is right: he used to be kind of cool once.

From their 1980 album Zenyattà Mondatta, The Police with “Don’t Stand So Close to Me”.
Read more »

March 31, 2008 Posted by John Armstrong | Sunday Samples | | No Comments

Two Mean Value Theorems

We’ve got two different analogues of the integral mean value theorem for the Riemann-Stieltjes integral.

The first one says that if \alpha is increasing on \left[a,b\right] and f is integrable with respect to \alpha, with supremum M and infimum m in the interval, then there is some “average value” c between m and M. This satisfies

\displaystyle\int\limits_{\left[a,b\right]}fd\alpha=c\int\limits_{\left[a,b\right]}d\alpha=c\left(\alpha(b)-\alpha(a)\right)

In particular, we should note that if f is continuous then the intermediate value theorem tells us that there is some x_0 with f(x_0)=c. That is, there is some x_0 such that

\displaystyle f(x_0)=\frac{1}{\alpha(b)-\alpha(a)}\int\limits_{\left[a,b\right]}fd\alpha

When \alpha(x)=x this gives us the old integral mean value theorem back again.

So why does this work? Well, if \alpha(a)=\alpha(b) then both sides are zero and the theorem is trivially true. Now, the lowest lower sum is L_{\alpha,\{a,b\}}(f)=m\left(\alpha(b)-\alpha(a)\right), while the highest upper sum is U_{\alpha,\{a,b\}}(f)=M\left(\alpha(b)-\alpha(a)\right). The integral itself, which we’re assuming to exist, lies between these bounds:

\displaystyle m\left(\alpha(b)-\alpha(a)\right)\leq\int\limits_{\left[a,b\right]}fd\alpha\leq M\left(\alpha(b)-\alpha(a)\right)

So we can divide through by \int_{\left[a,b\right]}d\alpha=\alpha(b)-\alpha(a) to get the result we seek.

We can get a similar result which focuses on the integrator by using integration by parts. Let’s assume \alpha is continuous and f is increasing on \left[a,b\right]. Our sufficient conditions tell us that the integral of f with respect to \alpha exists, and the integration by parts formula says

\displaystyle\int\limits_{\left[a,b\right]}fd\alpha=f(b)\alpha(b)-f(a)\alpha(a)-\int\limits_{\left[a,b\right]}\alpha df

But the first integral mean value theorem tells us that the integral on the right is equal to \alpha(x_0)\left(f(b)-f(a)\right) for some x_0. Then we can rearrange the above formula to read

\displaystyle\int\limits_{\left[a,b\right]}fd\alpha=f(b)\alpha(b)-f(a)\alpha(a)-\alpha(x_0)\left(f(b)-f(a)\right)
\displaystyle=f(a)\left(\alpha(x_0)-\alpha(a)\right)+f(b)\left(\alpha(b)-\alpha(x_0)\right)
\displaystyle=f(a)\int_{\left[a,x_0\right]}d\alpha+f(b)\int_{\left[x_0,b\right]}d\alpha

So there is some point x_0 so that the integral of f is the same as the integral of the step function taking the value f(a) until x_0 and the value f(b) after it.

March 28, 2008 Posted by John Armstrong | Analysis, Calculus | | No Comments

On a Lighter Note

March 28, 2008 Posted by John Armstrong | Uncategorized | | No Comments

Hoax!

Isabel at God Plays Dice cites (uncharacteristically credulously) this story of a 13th-century monk who supposedly discovered the Mandelbrot set. How surprising, since we didn’t see it until we had computer graphics. But there’s a lot more anachronism just beneath the surface. There are so many topics that feed into the standard representation of the Mandelbrot set that we have so completely internalized we hardly know how to think without them.

First, this story asks us to believe that Brother Udo worked with complex numbers 300 years before Tartaglia and Cardano hesitantly advanced them. They were motivated by problems that could not be solved without them, while here they’re just being used for multiplications that can be perfectly well defined over the reals. Why would he have even thought to use complex numbers? Sure, maybe he was using a geometric description that happened to exactly correspond, but he wouldn’t have thought of it as complex multiplication. Still, more on that possibility later.

Next the story asks us to believe that Brother Udo graphed algebraic relations 400 years before Descartes introduced analytic geometry. So he somehow had the idea that algebra and geometry were interlinked so far in advance of everybody else and didn’t tell anyone? Again, maybe the construction was wholly geometrical.

Next the story asks us to believe that Brother Udo was not only graphing any old algebraic relations, but using the graphical representation of complex numbers 500 years before Argand described that interpretation. Again, this is something we internalize in high school (or earlier) and so we forget how late in the game it actually came along.

Okay, so could there have been a totally geometrical construction that these historians are simplifying into the function z\mapsto z^2+c? I doubt it. Mostly because you run it 70 times and then.. what? In rendering of the Mandelbrot set you know whether to keep or toss a point by whether you eventually go outside the circle of radius 2. So now we need to posit an unspoken proof that once a point leaves that circle it can never return. No reference is made to that result either.

Of course, the page cites sources. But the only source directly on point for the main content of the story is this one:

[5] Schipke, R.J. and Eberhardt, A. “The forgotten genius of Udo von Aachen”, Harvard Journal of Historical Mathematics, 32, 3 (March 1999), pp 34-77.

but this “Harvard Journal of Historical Mathematics” doesn’t seem to exist either. In fact, a Google search on that journal’s name gives… a bunch of references back to this very story! Let me be clear about this:

Just because someone includes a bibliography doesn’t mean they’re not forging the bibliography too!

And, of course, the page itself says it was published April 1, 1999. It’s not even a new hoax. But as I show above, you don’t need to see that cue to realize it’s a joke. You just have to stop, think, and really understand what someone is trying to sell you. You have to be skeptical about any assertion of fact. And it seems a lot of people aren’t.

This story got included in course materials for a course on fractal geometry taught at Yale by my teaching mentor Michael Frame, with help by Benoît Mandelbrot himself. John Allen Paulos wrote it up for ABC news. Sure, maybe they knew it was a hoax and passed it along as such. Paulos wrote his article on April 1, 2001, after all. But then it gets swept up by woo like “Life Technology”, “The All-Oneness Hadron Materia”, and badly-researched psychology texts. And nobody tags it as the hoax it is.

So here I am. This story is false. Unquestionably false, like faking the moon landing or a face on Mars. False like psychic powers and divine cameos on corn chips. False like that girl who got so high on LSD she microwaved the baby she was taking care of. Those of you who love linking to big-name skepti-bloggers might want to point them here so we can squelch this thing once and for all.

March 28, 2008 Posted by John Armstrong | rants | | 11 Comments

2008 Abel Prize

As Chris Hillman just pointed out in a comment, the 2008 Abel Prize went to John Griggs Thompson and Jacques Tits “for their profound achievements in algebra and in particular for shaping modern group theory”. The comment went on my recent throwaway post about the 7×7x7 Rubik’s Cube, but a more appropriate one might have been this one from over a year ago, in which I discuss the Feit-Thompson theorem in passing.

Incidentally, I think I’ve met both of the winners. Tits I’m sure of, since I tried and failed horribly to take a short course he gave at Yale on “buildings”. Thompson I believe showed up for Walter Feit’s memorial, but I could be wrong about that. I wish I could say I was particularly close to one or the other, but I suppose that will have to wait until Adams and Vogan win the prize.

March 27, 2008 Posted by John Armstrong | Uncategorized | | 19 Comments

Step Function Integrators

Now that we know how a Riemann-Stieltjes integral behaves where the integrand has a jump, we can put jumps together into more complicated functions. The ones we’re interested in are called “step functions” because their graphs look like steps: flat stretches between jumps up and down.

More specifically, let’s say we have a sequence of points

\displaystyle a\leq c_1<c_2<...<c_n\leq b

and define a function \alpha to be constant in each open interval \left(c_{i-1},c_i\right). We can have any constant values on these intervals, and any values at the jump points. The difference \alpha(c_k^+)-\alpha(c_k^-) we call the “jump” of \alpha at c_k. We have to be careful here about the endpoints, though: if c_1=a then the jump at c_1 is \alpha(c_1^+)-\alpha(c_1), and if c_n=b then the jump at c_n is \alpha(c_n)-\alpha(c_n^-). We’ll designate the jump of \alpha at c_k by \alpha_k.

So, as before, the function \alpha may or may not be continuous from one side or the other at a jump point c_k. And if we have a function f discontinuous on the same side of the same point, then the integral can’t exist. So let’s consider any function f so that at each c_k, at least one of \alpha or f is continuous from the left, and at least one is continuous from the right. We can chop up the interval \left[a,b\right] into chunks so that each one contains only one jump, and then the result from last time (along with the “linearity” of the integral in the interval) tells us that

\displaystyle\int\limits_{\left[a,b\right]}fd\alpha=\sum\limits_{i=1}^nf(c_k)\alpha_k

That is, each jump gives us the function at the jump point times the jump at that point, and we just add them all together. So finite weighted sums of function evaluations are just special cases of Riemann-Stieltjes integrals.

Here’s a particularly nice family of examples. Let’s start with any interval \left[a,b\right] and some natural number n. Define a step function \alpha_n by starting with \alpha_n(a)=a and jumping up by \frac{b-a}{n} at a+\frac{1}{2}\frac{b-a}{n}, a+\frac{3}{2}\frac{b-a}{n}, a+\frac{5}{2}\frac{b-a}{n}, and so on. Then the integral of any continuous function on \left[a,b\right] gives

\displaystyle\int\limits_{\left[a,b\right]}fd\alpha_n=\sum\limits_{i=1}^nf\left(a+(i-\frac{1}{2})\frac{b-a}{n}\right)\frac{b-a}{n}

But notice that this is just a Riemann sum for the function f. Since f is continuous, we know that it’s Riemann integrable, and so as n gets larger and larger, these Riemann sums must converge to the Riemann integral. That is

\displaystyle\lim\limits_{n\to\infty}\int\limits_{\left[a,b\right]}fd\alpha_n=\int\limits_{\left[a,b\right]}fdx

But at the same time we see that \alpha_n(x) converges to x. Clearly there is some connection between convergence and integration to be explored here.

March 27, 2008 Posted by John Armstrong | Analysis, Calculus | | 1 Comment

Integrating Across a Jump

In the discussion of necessary conditions for Riemann-Stieltjes integrability we saw that when the integrand and integrator are discontinuous from the same side of the same point, the integral can’t exist. But how close can we come to that situation? It turns out that as long as one of the two functions is continuous from each side, then things generally work out.

Specifically, let’s consider jump discontinuities. These are especially useful to understand, since they’re the only sort a function of bounded variation can have. So let’s say we take an interior point c\in\left(a,b\right) and define as simple a function \alpha as we can with a jump there. We let it be constant on either side, with \alpha(x)=\alpha(a) for x\in\left[a,c\right) and \alpha(x)=\alpha(b) for x\in\left(c,b\right]. We’ll let \alpha(c) be anything at all. Generally it will be discontinuous from both sides at c, but if \alpha(c)=\alpha(a) then we’ll have continuity from the left, and similarly on the right. Of course, we could have \alpha(a)=\alpha(b), with c being the only point with a different value for \alpha.

Now let’s let f be any other function on \left[a,b\right]. We know that we can’t let it be discontinuous from the left at c if \alpha is, or from the right either, so let’s assume it’s continuous from the left or the right at c to satisfy these conditions, but put no other assumptions on it. I assert that f is then integrable with respect to \alpha on \left[a,b\right], and has the value

\displaystyle\int\limits_{\left[a,b\right]}fd\alpha=f(c)\left(\alpha(c^+)-\alpha(c^-)\right)

Where \alpha(c^+) is the limit of \alpha as we approach c from the right. In our situation this will be \alpha(b), but we’ll telegraph a bit by writing it like this. Similarly, \alpha(c^-) is the limit of \alpha as we approach c from the left, which here is \alpha(a).

To see that this is the case, take a tagged partition x, and if it doesn’t already contain c just refine it by throwing the new point in. Now every term in the Riemann-Stieltjes sum

\displaystyle f_{\alpha,x}=\sum\limits_{i=1}^nf(t_i)\left(\alpha(x_i)-\alpha(x_{i-1})\right)

is zero except for the one on either side of c=x_k. We can then write the difference f_{\alpha,x}-f(c)\left(\alpha(c^+)-\alpha(c^-)\right) as

\displaystyle\left(f(t_{k-1})-f(c)\right)\left(\alpha(c)-\alpha(x_{k-1})\right)+\left(f(t_k)-f(c)\right)\left(\alpha(x_{k+1})-\alpha(c)\right)

Since either f or \alpha is continuous from the left, the first term above must go to zero. Similarly, because at least one is continuous from the right, the second term must also go to zero. Thus the sums f_{\alpha,x} converge to the value asserted.

March 26, 2008 Posted by John Armstrong | Analysis, Calculus | | 2 Comments

Necessary Conditions for Integrability

We’ve talked about sufficient conditions for integrability, which will tell us that a given integral does exist. Now we’ll consider the situation where we know that a given integral exists, and see what we can tell about its integrand or integrator.

First, let’s consider an increasing integrator \alpha over an interval \left[a,b\right], and an interior point c\in\left(a,b\right). Let’s assume that both \alpha and the integrand f are discontinuous from the right at c. Then the integral \int_{\left[a,b\right]} f\,d\alpha cannot exist.

Consider what the assumption of discontinuity means: there is some \epsilon so that for every \delta>0, no matter how small, there are points x,y\in\left(c,c+\delta\right) so that |f(x)-f(c)|\geq\epsilon and |\alpha(y)-\alpha(c)|\geq\epsilon. That is, as we approach c from the right we can always find a sample point that differs from the value we want by at least \epsilon, both in the integrand and in the integrator. So let’s see what happens when we throw c in as a partition point.

Suddenly we’re stuck! Riemann’s condition tells us to look at the sum

\displaystyle U_{\alpha,x}(f)-L_{\alpha,x}(f)=\sum\limits_{i=1}^n\left(M_i(f)-m_i(f)\right)\left(\alpha(x_i)-\alpha(x_{i-1})\right)

where each term is positive. So to make this sum go down to zero in Riemann’s condition we need to make each term go down to zero. But for the subinterval whose left endpoint is c, we can’t make this happen! The change in \alpha will always be at least \epsilon for some smaller subinterval, and the difference in sample values of f must always exceed \epsilon as well. Thus the difference between the upper and lower sum will always be at least \epsilon^2>0, violating Riemann’s condition, and thus integrability.

Of course, we could run through the same argument approaching c from the left. We could also use an integrator \alpha of bounded variation, as usual.

What this tells us is that while we may allow integrands or integrators to be discontinuous, they can’t be discontinuous from the same side at the same point. For example, any discontinuity in an integrator of bounded variation is fine as long as we’re pairing it off with a continuous integrand, as we did last time.

March 25, 2008 Posted by John Armstrong | Analysis, Calculus | | 1 Comment