  A new Monte Carlo model is proposed for radiative corrections to Bhabha
scattering by extending QEDPS developed for multi-photon emission in muon pair
production in $e^+e^-$ annihilation. This is the QED version of the model known
as parton shower in QCD. The main difference between muon pair production and
Bhabha scattering is that the latter cross section shows the singularity of
$1/t^2$. A shower algorithm is constructed on the radiator formalism modified
in a suitable form for this singularity. Some results of the model are
presented and compared with $O(\alpha)$ corrections.
