The evaporation of the collapsing star can be equivalently described as a negative energy flux of radiation travelling radially inwards towards the center of the star. Here, we are interested in the evolution of the star during its collapse. 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.