It is well known that the fragmentation equation admits self-similar solutions for evolving particle-size distributions (PSD); i.e., if the shape of PSD is independent of time after an initial transient period. Although an analytical derivations of the self-similar PSD cases have been studied extensively, results for cases requiring numerical solutions are rare. The aim of the present work is to fill this gap for the case of homogeneous breakage functions. The known analytical and approximate solutions for the self-similar PSD are reviewed and a general algorithm for the numerical solution is proposed. Results for a broad range of breakage functions (kernel and rate) are presented. Further, the work is focused on the sensitivity of the relation between self-similar PSD and breakage kernel and its influence on the inverse breakage problem, i.e., that of estimating the breakage kernel from experimental self-similar PSDs. Useful suggestions are made for tackling the inverse problem. (c) 2004 Elsevier Inc. All rights reserved.