Abstract: The Robin problem for the Helmholtz equation in normal-polar annuli is addressed by using a suitable Fourier-Hankel series technique. Attention is in particular focused on the wide class of domains whose boundaries are defined by the so-called superformula introduced by Gielis. A dedicated numerical procedure based on the computer algebra system Mathematica© is developed in order to validate the proposed methodology. In this way, highly accurate approximations of the solution, featuring properties similar to the classical ones, are obtained.