A field-first method is developed for locating equilibrium points in three-body gravitational systems by treating gravity as a continuous vector field on R^3. Each mass contributes a spherically symmetric acceleration field, and the net field is obtained by superposition. Equilibria are defined as zeros of the net acceleration. A Newton-type refinement based on the 3×3 Jacobian of the acceleration field (equivalently, the Hessian of the potential) yields fast local convergence. Robust initialization is obtained by a spherical-shell exploration stage: on each shell, directions with small field magnitude are identified by minimizing ||g|| on the sphere, and the corresponding Cartesian points seed the Jacobian-based refinement. The Jacobian eigenstructure provides a local classification of the equilibria. The resulting pipeline is naturally parallel and compatible with high-throughput computation. Quantum acceleration is potentially relevant to the global search component, while the 3×3 refinement solves are already inexpensive classically.