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