Helium and heavy-element diffusion are both included in precise calculations of solar models. In addition, improvements in the input data for solar interior models are described for nuclear reaction rates, the solar luminosity, the solar age, heavy-element abundances, radiative opacities, helium and metal diffusion rates, and neutrino interaction cross sections. The effects on the neutrino fluxes of each change in the input physics are evaluated separately by constructing a series of solar models with one additional improvement added at each stage. The effective 1 cr uncertainties in the individual input quantities are estimated and used to evaluate the uncertainties in the calculated neutrino fluxes and the calculated event rates for solar neutrino experiments. The calculated neutrino event rates, including all of the improvements, are 9.3* J'\ SNU for the 37 C1 experiment and 137*7 SNU for the 71 Ga experiments. The calculated flux of 7 Be neutrinos is 5.1 (1.00ig;Sf)X10 9 cm -2 s _1 and the flux of 8 B neutrinos is 6.6(1.00^) X 10 6 cm -2 s _1 . The primordial' helium abundance found for this model is Y=0.278. The present-day surface abundance of the model is Y s = 0.247, in agreement with the helioseismological measurement of 7^ = 0.242±0.003 determined by Hernandez and Christensen-Dalsgaard (1994). The computed depth of the convective zone is R -0.712R 0 , in agreement with the observed value determined from /?-mode oscillation data of R = 0.713±0.003/?© found by Christensen-Dalsgaard et al. (1991). Although the present results increase the predicted event rate in the four operating solar neutrino experiments by almost l