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