
%zgamma.m

function  g=zgamma(n,m,s,c,a)

g0=0; g1=m+2*c/a;

for i=1:n
    g2=(g1+g0)/2;
    k0=gamma(g0,m,s,c,a);
    k2=gamma(g2,m,s,c,a);

    if k0*k2<0 g1=g2;
    else g0=g2;
    end
end g=g2; 
