dsolve(diff(y(x),x)-a*(x^n-x)*y(x)^3-y(x)^2 = 0,y(x))