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