dsolve(diff(y(x),x)+a*diff(phi(x),x)*y(x)^3+6*a*phi(x)*y(x)^2+(2*a+1)*y(x)*diff(diff(phi(x),x),x)/diff(phi(x),x)+2*a+2 = 0,y(x))