Numerical studies of electron-phonon mediated superconductors