This paper presents a numerical algorithm for BIBO stability testing of a certain class of the so-called fractional-delay systems. The characteristic function of the systems under consideration is a multi-valued function of the Laplace variable s which is defined on a Riemann surface with finite number of Riemann sheets where the origin is a branch point. The stability analysis of such systems is not straightforward because there is no universally applicable analytical method to find the roots of the characteristic equation on the right half-plane of the first Riemann sheet. The proposed method is based on the Rouche’s theorem which provides the number of the zeros of a given function in a given simple closed contour. One advantage of the proposed method over previous works is that it gives the number and the location of the unstable poles. The algorithm has a reliable result which is illustrated by several examples.