Abstract: The interior and exterior Robin problems for the Helmholtz equation in starlike planar domains are addressed by using a suitable Fourier-like technique. Attention is in particular focused on normal-polar domains whose boundaries are defined by the so-called superformula introduced by J. Gielis. A dedicated numerical procedure based on the computer algebra system Mathematica© is developed in order to validate the proposed approach. In this way, highly accurate approximations of the solution, featuring properties similar to the classical ones, are obtained. The computed results are found to be in good agreement with the theoretical findings on Fourier series expansion presented by L. Carleson.