 function [f,df] = fgamma(dummy,g);

    f  =  1./gamma(g);
    df =  -psi(g)./gamma(g);
    
return;
