solve y(t)' = a - b*y(t)^2