function Ga=Gamma_q(q,alpha,k)
s=1;
m=(1-q^(alpha));
for i=1:k
    s=s*(1-q^i);
    m=m*(1-q^(alpha+i));
end
Ga=(s/m)*(1-q)^(1-alpha);
end
