Program to subtract stars from an image by modelling the PSF