dv(t)/dt = g(T - mu) + mu v(t)^2 /r , v(0)=0