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