function G=gamma(i)

% gamma.m --> Returns the value of the sectionalized breakup
% kernerl, gamma, for the particle population balance
% design problem

Go=1e-8*100;			% constant, s/m

global d0 mfr T Po MW mu smf rhos kb R rhog Vo SECT rhomix nu startnum D L k m 
global uz A eps

if (i==1)
   G=0;
else
   G=Go*(eps(m)/nu)^(3/2)*(3*Vo/2)^(1/3)*2^((i-1)/3);
end