dsolve(diff(diff(y(x),x),x)*y(x)-diff(y(x),x)^2-(-1+a*y(x))*diff(y(x),x)+2*a^2*y(x)^2-2*b^2*y(x)^3+a*y(x)=0,y(x))