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