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