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