The stress-energy tensor for the massless spin 1/2 field is numerically computed outside and on the event horizons of both charged and uncharged static nonrotating black holes, corresponding to the Schwarzschild, Reissner-Nordström, and extreme Reissner-Nordström solutions of Einstein's equations. The field is assumed to be in a thermal state at the black hole temperature. Comparison is made between the numerical results and previous analytic approximations for the stress-energy tensor in these spacetimes. For the Schwarzschild (charge zero) solution, it is shown that the stress energy differs even in sign from the analytic approximation. For the Reissner-Nordström and extreme Reissner-Nordström solutions, divergences predicted by the analytic approximations are shown not to exist.