裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

Julia の小ネタ--030 円周率 π の高精度表示

2021年06月28日 | ブログラミング

Julia では π は特殊な型を持つ

typeof(π)Irrational{:π}

なので,

π # π = 3.1415926535897...

という表示になる。 ... に意味がある。

big(π) あるいは BigFloat(π)

3.141592653589793238462643383279502884197169399375105820974944592307816406286198

になるが,もっと桁数が欲しい(?)ときは,

setprecision(BigFloat, 4096) do
    BigFloat(pi)
end

などとすると,

3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679821480865132823066470938446095505822317253594081284811174502841027019385211055596446229489549303819644288109756659334461284756482337867831652712019091456485669234603486104543266482133936072602491412737245870066063155881748815209209628292540917153643678925903600113305305488204665213841469519415116094330572703657595919530921861173819326117931051185480744623799627495673518857527248912279381830119491298336733624406566430860213949463952247371907021798609437027705392171762931767523846748184676694051320005681271452635608277857713427577896091736371787214684409012249534301465495853710507922796892589235420199561121290219608640344181598136297747713099605187072113499999983729780499510597317328160963185950244594553469083026425223082533446850352619311881710100031378387528865875332083814206171776691473035982534904287554687311595628638823537875937519577818577805321712268066130019278766111959092164201989380952572010654858632788659361533818279682303019520353018529689957736225994138912497217752834791315155748572424541506959508295331168617278558890750983817546374649393192550604009277016711390098488240128583616035637076601047101819429564

が得られる。4096 が計算に使うビット数ということで,これを大きくすればいくらでも桁数を増やせる。

指定するビット数を n とすれば,ほぼ n*log10(2) 桁が表示される。

コメント    この記事についてブログを書く
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする
« Julia の小ネタ--029 実数演... | トップ | Julia で五数要約値(R の fi... »
最新の画像もっと見る

コメントを投稿

ブログラミング」カテゴリの最新記事