円周率を計算するPerlプログラムを github で公開している人がいた。
https://gist.github.com/djsas/2083243
===ここから===
use strict;
use warnings;
use bigint;
my($t1,$t2);
my $n=2000;
my $a;
my $b;
my $i;
$t1=times;
$a = $b = 10 ** $n;
for($i = $n * 8 + 1; $i >= 3; $i -= 2) {
$a = int(int($i / 2) * ($a + $b * 2) / $i);
}
$a -= $b;
print("3.$a\n");
$t2=times;
print "<".($t2-$t1).">";
円周率の公式と計算法
http://www.kurims.kyoto-u.ac.jp/~kenkyubu/kokai-koza/H16-ooura.pdf
円周率を計算してみよう
http://w3e.kanazawa-it.ac.jp/e-scimath/contents/t16/textbook_t16_all.pdf
https://gist.github.com/djsas/2083243
===ここから===
use strict;
use warnings;
use bigint;
my($t1,$t2);
my $n=2000;
my $a;
my $b;
my $i;
$t1=times;
$a = $b = 10 ** $n;
for($i = $n * 8 + 1; $i >= 3; $i -= 2) {
$a = int(int($i / 2) * ($a + $b * 2) / $i);
}
$a -= $b;
print("3.$a\n");
$t2=times;
print "<".($t2-$t1).">";
===ここまで===
円周率の公式と計算法
http://www.kurims.kyoto-u.ac.jp/~kenkyubu/kokai-koza/H16-ooura.pdf
円周率を計算してみよう
http://w3e.kanazawa-it.ac.jp/e-scimath/contents/t16/textbook_t16_all.pdf
コメント
コメントを投稿