M832 did not really fit in the "big picture" of my study plan this year, not only because I am notoriously bad at numerical calculations. Having invested a lot of effort and sleepless nights in developing intuition for the behaviour of analytic functions I was suddenly confronted by the cubic splines which seemed to have all those properties, that well-mannered functions would be never allowed to possess.

Nicely, but artificially glued together of several pieces of cubics, smooth only up to the second derivative, vanishing on the entire intervals, now this is what seems really counter-intuitive. Be that as it may, the TMA deadlines had to be met.

The following is the kind of problem I got stuck with for a while. Obviously I am not replicating a TMA question here, giving an extended solution to one of the problems from the Course Notes. So the task is to express a function, say in terms of cubic B-splines on the entire real axis. I am omitting a lot of background material focusing on one particular idea that arises in the solution.

Since has a supporting interval of length 4 outside which it vanishes, we can start by expressing the function in terms of on and then try to extend the result. Calculating the expressions for the splines on :

multiplying by the respective coefficients, summing and equating powers of on each side we arrive at the following system of equations:

Having the solution (guaranteed by the Schoenberg-Whitney theorem):

Now we want to find all coefficients on each of the intervals for the points . From the general expression for the B-spline it can be deduced that

which for leads to the following recurrence relation:

or, after changing the index

Now here is the trick that I came up with and that was not (at least not explicitly) described in the Course Notes or the set book. The last expression can be thought of as a "second-order linear inhomogeneous recurrence relation". The advantage of this approach is that the structure of the solution instantly becomes clear.

The general solution of the corresponding homogeneous relation

is derived in the course notes, using the standard method of solving this type of recurrencies and is given by the following expression:

It can also be found using generating functions. Not surprisingly it depends on 2 arbitary constants, as it takes 2 initial terms, and to reconstruct the whole sequence from the three-term recurrency. Applying the general ideas from the linear systems we deduce that in order to obtain the general solution of the inhomogeneous recurrency we have to add a particular solution to the expression above.

Since the RHS is the quadratic polynomial it makes sence to look for the particular solution in the form:

Substituting this into the original recurrency and gatherig together the powers of we obtain:

which after equating powers gives the solution

Thus the general solution of the inhomogeneous equation is given by the following formula:

Now we can use the values of and to determine the constants (bearing in mind that ):

which gives . Thus finally:

which is the solution of the original problem.