dn y/dxn =-ky


By Julian Pulman on Wednesday, May 08, 2002 - 07:49 pm:

dy/dx = -ky , when solved gives an Exponential Decay Relationship.
d2 y/dx2 =-ky , when solved gives a trigonometric relationship.

Both of these are revelent in modelling certain events, such as Radioactive Decay, and Harmonic Oscillators (respectively).
Does the solution to d3 y/dx3 = -ky exist? and if so, is the solution of use in modelling some aspect of physics?
What about dn y/dxn =-ky as n surpasses 3?
What methods can be used to solve such equations? I'm only aware of what I learnt in Further Maths, for 2nd order solutions.

Thanks,

Julian


By Dan Goodman on Thursday, May 09, 2002 - 01:09 am:
Do you know about complex numbers? If so, you'll know that the trigonometric solutions and exponential decay solutions are, in some sense, the some sort of thing.

For y(n)+k y=0 (here y(n) means the nth derivative dn y/d xn), we try a solution y=elx. Substitute this in and you get ln elx+k elx=0. Divide by elx to get ln +k=0.

There are n solutions to this equation (if you use complex numbers), and each gives rise to a different solution. Linear combinations of these give you the general solution.

For example, if n=3 and k=-2 then we want to solve l3=2. One solution is obvious, l = 21/3. So, if we check it we do get y=e21/3x is a solution, because then y ' =21/3y, y ' ' =22/3 y and y ' ' ' =2y, so y ' ' ' -2y=0.

If you know about complex numbers, you'll see the same thing works with the other two solutions l = 21/3 e2ip/3 and l = 21/3eip/3. This will give us things like e21/3cos(2p/e)xcos(21/3sin(2p/3)x), yuch.

If you don't know about complex numbers, the thing to do is to try solutions like eA xcos(B x), eA xsin(B x) and simply eA x and see where you can get.

The same method works for general nth order homogeneous linear ordinary differential equations, i.e. yn+an-1yn-1+¼+a0 y=0. You end up needing to solve the equation ln +an-1ln-1+¼+ a1l+a0=0, which is difficult in general but easy if all the ai are 0 except a0. In case you wanted to know the level of difficulty, at Cambridge this is taught in the first term of the first year of maths. So, pretty difficult but not too difficult (even if I've made it so by posting too late).

I don't know of any applications in physics I'm afraid, but no doubt there are some :-)

Sorry if that was a bit confused (or a bit too difficult), I'm up too late. Let me know if you didn't follow.


By Yatir Halevi on Friday, May 10, 2002 - 10:51 am:

for even n, the sinh and cosh functions fit.

Yatir


By Julian Pulman on Saturday, May 11, 2002 - 11:54 am:

Thanks everyone, I'm confident now I can use the same techniques outlined by Dan to solve the nth order ones. I'll get back if I have any trouble, thanks :-)