Skip to content

Instantly share code, notes, and snippets.

@colomon
Created January 3, 2014 14:58
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save colomon/8239209 to your computer and use it in GitHub Desktop.
Save colomon/8239209 to your computer and use it in GitHub Desktop.
sub postfix:<!>($n) {
[*] 1..$n;
}
sub infix:<choose>($n, $k) {
$n! div ($k! * ($n - $k)!);
}
sub Lupas-ln2-term($n) {
((2 * $n) choose $n).FatRat * ((-1) ** ($n - 1) * (5 * $n + 1)) / (16 ** $n * $n * ($n + 1/2));
}
sub Lupas-ln2($k) {
3/4 - ([+] (1..$k).map(&Lupas-ln2-term)) / 8;
}
for (1..30) {
my $lupas = Lupas-ln2($_);
say $lupas ~ " " ~ 2.log - $lupas;
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment