A surface integral method developed for one-dimensional multivalued surfaces has been extended to include multiple profiles. This paper is concerned with the numerical implementation of the parametrization of the surface profile, and the development of an algorithm in order to determine whether a point is located inside or outside the material. This formulation is especially suitable for arbitrary-shaped particles and for particles placed on flat surfaces. The result is a useful and versatile tool for studying the scattering by irregular particles in both the near and far fields. This method can be used to investigate optical effects associated with morphologydependent resonances and multiple scattering by particles of arbitrary shape.