solve y'(t) = ay(t) - b/(1+e^(2t)) * (1-y(t)^2)