diff(diff(y(x),x),x)*y(x)-diff(y(x),x)^2+(-1+a*y(x))*diff(y(x),x)-y(x)*(1+y(x))*(b^2*y(x)^2-a^2)