real scalar comb1(n, k) {
	return(exp(lnfactorial(n)-lnfactorial(k)-lnfactorial(n-k)))
}

real scalar perm(n, k) {
	return(exp(lnfactorial(n)-lnfactorial(n-k)))
}