Thus we include the backreaction of the negative energy Hawking flux in the interior geometry of the collapsing star and solve the full 4-dimensional Einstein and hydrodynamical equations numerically.