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