"Testing tgamma"
clear
for(k,1,3,check(tgamma(k)==factorial(k-1)))
