Sorry for the two-day delay. The personal business to which I alluded kept me out until yesterday morning, delaying this post until today. I hope it’s decent for all that!

This week, the same function as last week. There’s nothing so special about the function itself, but it does serve as a nice Rorschach blot to project our own wandering mathematical thoughts upon. As with last time, the specific “f(x) = stuff” form of the function is not particularly enlightening. It’s just a few lines tacked together. The graph is all we need.

Last time we set ourselves the task of writing this function as a trigonometric series, offhandedly mentioning the power series as a more usual way of expanding a function. A power series, as you may know, is when we figure out the constants “a_{n}” that make the following relationship true, if indeed such a thing can be done.

This can usually be done with most functions of physical interest. Sometimes it will work only for some x within a particular interval, and more rarely it may in fact work for no x (or at least almost no x). Fortunately telling which is the case is not just a shot in the dark. If a function is complex-differentiable in the neighborhood of x, the power series converges appropriately at that point. What about our function?

Well we can forget about complex differentiability right off the bat. It’s not even regular differentiable everywhere. For those who are not aficionados of the calculus, remember that a derivative represents the slope of the graph at a particular point. At the sharp points our function is not differentiable. It just doesn’t have a slope at all. And so if you try to find a power series that works everywhere you won’t succeed. There just isn’t one.

The Fourier series worked because its conditions for convergence are less onerous. To ruin the nice technical conditions with simplification, I can say that if the function is integrable and not deliberately stupid the Fourier series will converge. But come on, there’s got to be *some* polynomial of degree n that approximates the function better than the rest of the degree n polynomials. Can we find it? Of course!

It turns out that Fourier series aren’t the only game in town if you’re looking for complete and orthogonal sets of functions. There’s tons of ’em. An infinite number, in fact. One of the best known and most useful – we have in fact used it for a very similar purpose before – is the Legendre polynomials.

They’re just ordinary looking polynomials. To list a few,

P1 = 1

P2 = x

P3 = 0.5*(x^2 – 1)

A series in the Legendre polynomials would look like this:

Over the interval [-1,1], they’re orthogonal in the sense we talked about last time (feel free to take a gander at the last Sunday Function to refresh your memory). They work in the same way, and we find the coefficients in the same way. We do need to shift to the interval [0,2], but this isn’t hard. In the style of TV chefs everywhere I’ll pull our necessary formula for the coefficients pre-baked out of the oven. We have to pre-bake the individual polynomials themselves too. This isn’t hard either, just replace x with x – 1. For the coefficients, we get:

Don’t forget that the Pn are themselves functions of x. We then can write a series in the Legendre polynomials. But this will just be a sum of powers of x – how is this not just a power series in disguise? One important reason is that the coefficients in front of each individual power of x don’t stay constant as you add in higher-order Legendre polynomials. This has its price. Unlike the standard power series we don’t have beautiful results like Taylor’s theorem, nor can we say nearly as much in general about what the Legendre series is doing about any given point. But frequently the trade-off is worth it.

Shall we see how this series looks? The convergence turns out not to be so pretty or fast as the Fourier series. But functional and ugly works too.

Only the first non-zero term:

The first two non-zero terms:

The first six non-zero terms:

I’m including a

in case you want to play around with this a little bit without having to write it yourself. Not that it matters immensely, but consider the code licensed so that you can do whatever you want with it so long as you attribute and share alike.

legendre.nb by Matthew Springer is licensed under a Creative Commons Attribution-Share Alike 3.0 United States License.

Based on a work at scienceblogs.com.