a=0.8;
n=100;
s=0;
for k=0:n
    for j=0:k
        Fs=(-1)^k*3^(-j)*gamma(a)*gamma(a);
        Js=gamma(j+1)*gamma(a-j)*gamma(k-j+1)*gamma(a-k+j);
        s=s+Fs/Js;
    end
    s=s*(3/2)^a;
end
s