  We present and apply rigorous dynamical modeling with which we infer
unprecedented constraints on the stellar and dark matter mass distribution
within our Milky Way (MW), based on large sets of phase-space data on
individual stars. Specifically, we model the dynamics of 16,269 G-type dwarfs
from SEGUE, which sample 5 < R_GC/kpc < 12 and 0.3 < |Z|/kpc < 3. We
independently fit a parameterized MW potential and a three-integral,
action-based distribution function (DF) to the phase-space data of 43 separate
abundance-selected sub-populations (MAPs), accounting for the complex selection
effects affecting the data. We robustly measure the total surface density
within 1.1 kpc of the mid-plane to 5% over 4.5 < R_GC/kpc < 9. Using metal-poor
MAPs with small radial scale lengths as dynamical tracers probes 4.5 < R_GC/kpc
< 7, while MAPs with longer radial scale lengths sample 7 < R_GC/kpc < 9. We
measure the mass-weighted Galactic disk scale length to be R_d = 2.15+/-0.14
kpc, in agreement with the photometrically inferred spatial distribution of
stellar mass. We thereby measure dynamically the mass of the Galactic stellar
disk to unprecedented accuracy: M_* = 4.6+/-0.3+3.0x(R_0/kpc-8)x10^{10}Msun and
a total local surface density of \Sigma_{R_0}(Z=1.1 kpc) = 68+/-4 Msun/pc^2 of
which 38+/-4 Msun/pc^2 is contributed by stars and stellar remnants. By
combining our surface density measurements with the terminal velocity curve, we
find that the MW's disk is maximal in that V_{c,disk} / V_{c,total} =
0.83+/-0.04 at R=2.2 R_d. We also constrain for the first time the radial
profile of the dark halo at such small Galactocentric radii, finding that
\rho_{DM} (r;near R_0) \propto 1 / r^\alpha with \alpha < 1.53 at 95%
confidence. Our results show that action-based distribution-function modeling
of complex stellar data sets is now a feasible approach that will be fruitful
for interpreting Gaia data.
