We describe a numerical method to measure the location and strength of
complex singularities of a Hele-Shaw flow. We are able to detect the
singularities with exponent (close to) -4/3, predicted by perturbative
calculations. We also give numerical evidence for a new mechanism of
singularity generation.