We fit axisymmetric 3-integral dynamical models to NGC 3379 using the line-of-sight velocity distribution obtained from HST/FOS spectra of the galaxy center and ground-based long-slit spectroscopy along four position angles, with the light distribution constrained by WFPC2 and ground-based images. We have fitted models with inclinations from 29 • (intrinsic galaxy type E5) to 90 • (intrinsic E1) and black hole masses from 0 to 10 9 M ⊙ .The best-fit black hole masses range from 6 × 10 7 − 2 × 10 8 M ⊙ , depending on inclination. The preferred inclination is 90 • (edge-on); however, the constraints on allowed inclination are not very strong due to our assumption of constant M/L V . The velocity ellipsoid of the best model is not consistent with either isotropy or a two-integral distribution function. Along the major axis, the velocity ellipsoid becomes tangential at the innermost bin, radial in the mid-range radii, and tangential again at the outermost bins. The rotation rises quickly at small radii due to the presence of the black hole. For the acceptable models, the radial to tangential ((σ 2 θ + σ 2 φ )/2) dispersion in the mid-range radii ranges from 1.1 < σ r /σ t < 1.7, with the smaller black holes requiring larger radial anisotropy. Compared with these 3-integral models, 2-integral isotropic models overestimate the black hole mass since they cannot provide adequate radial motion. However, the models presented in this paper still contain restrictive assumptions-namely assumptions of constant M/L V and spheroidal symmetry-requiring yet more models to study black hole properties in complete generality.