function f=factorgamma(n)
    f = gamma(n+1)
endfunction