A numerical study of lattice QED with dynamical fermions