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