y'(t)=-sqrt(a)*sqrt(1/y(t)-1/b)