Wronskians are handy for dealing with some second order linear ordinary differential equations.
For example consider the equation y'' + p(t)y' +q(t)y=0 (*)

where p and q are non constant functions that are everywhere non zero.

**Calculating the Wronskian**
The wronskian of 2 linearly independant solutions y_{1} and y_{2} is y_{1}y_{2}'-y_{1}'y_{2}.

If we want to be able to do anything useful with the wronskian, we're going to want to be able to calculate it without knowing 2 linearly independant solutions.

Consider W'(t)=y_{1}y_{2}'' + y_{1}'y_{2}'-y_{1}'y_{2}'-y_{1}''y_{2}=y_{1}y_{2}''-y_{1}''y_{2}

then W'(t) + p W(t) = y_{1}(y_{2}''+py_{2}') - y_{2}(y_{1}''+py_{1}')

But since y_{1} and y_{2} are solutions of (*) then y_{2}''+py_{2}' = -qy_{2} and y_{1}''+py_{1}' = -qy_{1}

thus W'(t) + pW(t) = -y_{1}y_{2}q + y_{2}y_{1}q=0

We now have a first order linear equation for W, which causes us to rub our hands with glee and shout "Huzzah", as we know how to solve these:

W(t)=k exp(-∫p(s)ds), with k constant (for our purposes non zero)
This also shows that the wronskian is everywhere non null.

**Uses**

We now have a nice shiny expression for the Wronskian and I can tell some of you just can't wait to taste its awesome power. One use is if, by luck, the inspiration of the Holy Spirit or some other method, you have managed to guess what y_{1} is, and you would like to find a second linearly independant solution, as you will then have all the solutions of the equation (the set of solutions of an nth order linear differential equation is a vector space of dimension n). To do this, you just need to write down what the wronskian is again:

y_{1}y_{2}-y_{1}'y_{2}=W(t)

y_{2}/y_{1}-y_{1}'y_{2}/y_{1}^{2} = W(t)/y_{1}^{2}

Noticing that the left hand side is the derivative of y_{2}/y_{1} yields y_{2}=y_{1}∫_{a}^{t} W(s)/y_{1}(s)^{2}ds

Of course some times this integral will be so horrible that you've more or less wasted your time, but you will probably be able to obtain a series expansion by expanding the integrand. You could also define a new function which is defined to be the antiderivative of W(t)y_{1}^{2}

Another application of wronskians is guessing those pesky particular solutions. Our equation is now y'' + p(t)y' +q(t)y=h(t). As before, we'll have p and q everywhere non zero, and to make things interesting, non constant. I'll also assume h everywhere non-zero(if p and q constant you should try functions "like" h). I'll assume that, by hook or by crook you've solved the homogenous equation and so you have 2 linearly independant solutions y_{1} and y_{2}.
If h=0, then the solutions are the functions y = λy_{1}+μy_{2} , λ and μ constants. The idea behind this method is to try and find a particular solution that is a nonlinear combination of y_{1} and y_{2}, ie y = f y_{1}+g y_{2}, for some functions f and g of t. To make this useful, we are going to add an extra condition on f and g : f' y_{1}+g' y_{2}=0
we then have :

y=f y_{1}+g y_{2}

y'=f y_{1}'+g y_{2}' (using our condition on f and g)

y''=f' y_{1}' +g' y_{2}'+f y_{1}''+g y_{2}''

this gives y''+py'+qy = f (y_{1}''+py_{1}'+qy_{1})+g (y_{2}''+py_{2}'+qy_{2}) +f' y_{1}'+g' y_{2}'=h

y_{1} and y_{2} both solve the homogenous equation, so the first 2 terms disappear. we are left with :

f'y_{1}'+g'y_{2}'=h Multiplying through by W/h gives:

(f'W/h)y_{1}'+(g'W/h)y_{2}'=W
Identifying this with the expression of the Wronksian in terms of y_{1},y_{1}',y_{2} and y_{2}' gives :

f=-∫hy_{2}/W

h=∫hy_{1}/W

Which gives you a particular solution. Use it with care.

In case you were wondering, the wronskian is named after Mr. Wronski.