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 y1 and y2 is y1y2'-y1'y2.

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)=y1y2'' + y1'y2'-y1'y2'-y1''y2=y1y2''-y1''y2
then W'(t) + p W(t) = y1(y2''+py2') - y2(y1''+py1')
But since y1 and y2 are solutions of (*) then y2''+py2' = -qy2 and y1''+py1' = -qy1
thus W'(t) + pW(t) = -y1y2q + y2y1q=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 y1 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:
y1y2-y1'y2=W(t)
y2/y1-y1'y2/y12 = W(t)/y12
Noticing that the left hand side is the derivative of y2/y1 yields y2=y1at W(s)/y1(s)2ds
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)y12

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 y1 and y2. If h=0, then the solutions are the functions y = λy1+μy2 , λ and μ constants. The idea behind this method is to try and find a particular solution that is a nonlinear combination of y1 and y2, ie y = f y1+g y2, 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' y1+g' y2=0 we then have :
y=f y1+g y2
y'=f y1'+g y2' (using our condition on f and g)
y''=f' y1' +g' y2'+f y1''+g y2''

this gives y''+py'+qy = f (y1''+py1'+qy1)+g (y2''+py2'+qy2) +f' y1'+g' y2'=h
y1 and y2 both solve the homogenous equation, so the first 2 terms disappear. we are left with :
f'y1'+g'y2'=h Multiplying through by W/h gives:
(f'W/h)y1'+(g'W/h)y2'=W Identifying this with the expression of the Wronksian in terms of y1,y1',y2 and y2' gives :
f=-∫hy2/W
h=∫hy1/W
Which gives you a particular solution. Use it with care.

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