In this paper, triangular orthogonal functions (TF’s) method is applied as a basis in collocation method, to solve the population balance integro-differential equation, which arise in the description of particle-size distribution of a continuous, mixed-suspension, mixed product removal crystallizer with taking account of the effect of particle breakage. To demonstrate the validity of this method obtained results have been compared with other methods.