DSolve[Derivative[1][y][x] == y[x]/x - F[x]*(-(a*x^2) + y[x]^2),y[x],x]