dsolve(diff(y(x),x) = 1/2*(x^(3/2)+2*F(y(x)-1/6*x^3))*x^(1/2),y(x))