where Fi and Pi are the feed and product flow rates of size i, Hi is the
total mass of size i held in the mill, ri is the specific breakage rate, aij
is called appearance function, representing the fraction of size j
material which reports to the i size fraction after breakage, and di
is the specific discharge rate