10 lines
413 B
Text
10 lines
413 B
Text
sample(f,a,b)=for(i=1,4, my(n1=random(b-a)+a,n2=random(b-a)+a); [n1,n2]=[max(n1,n2),min(n1,n2)]; print(n1", "n2": "f(n1,n2)))
|
|
permExact(m,n)=factorback([m-n+1..m]);
|
|
combExact=binomial;
|
|
permApprox(m,n)=exp(lngamma(m+1)-lngamma(m-n+1));
|
|
combApprox(m,n)=exp(lngamma(m+1)-lngamma(n+1)-lngamma(m-n+1));
|
|
|
|
sample(permExact, 1, 12);
|
|
sample(combExact, 10, 60);
|
|
sample(permApprox, 5, 15000);
|
|
sample(combApprox, 100, 1000);
|