<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 def _resample_residuals(self, stars, epsf):
 """ Compute normalized residual images for all the input stars. Parameters stars : `EPSFStars` object The stars used to build the ePSF. epsf : `EPSFModel` object The ePSF model. Returns ------- star_imgs : 3D `~numpy.ndarray` A 3D cube containing the resampled residual images. """


        shape = (stars.n_good_stars, epsf.shape[0], epsf.shape[1])
        star_imgs = np.zeros(shape)
        for i, star in enumerate(stars.all_good_stars):
            star_imgs[i, :, :] = self._resample_residual(star, epsf)

        return star_imgs