dsolve(diff(y(x),x) = (b^3+y(x)^2*b^3+2*y(x)*b^2*a*x+x^2*b*a^2+y(x)^3*b^3+3*y(x)^2*b^2*a*x+3*y(x)*b*a^2*x^2+a^3*x^3)/b^3,y(x))