円周率は3ではだめなのか？

概要
Rでさまざまな円周率で円を書いて比較してみた。
元ネタは飲み会での先輩の一言。

コード
Rを選んだのは

  • プロットを簡単にしてくれる
  • 既にインストールされていた

Octaveでも良かったけど二つ目の条件に合致しなかったため。

x<-1:1000
p<-3.1415926535
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3.141592653
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3.14159265
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3.1415926
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3.141592
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3.14159
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3.1415
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3.141
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3.14
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3.1
plot(cos(p*x/500),sin(p*x/500),type='l')
p<-3
plot(cos(p*x/500),sin(p*x/500),type='l')

実行結果
円周率3.1415926535　以降 3.1415まで同様の結果なので省略
f:id:fakeyakuza:20110724221843p:image
円周率3.141　(1,0)あたりで微妙に縁が切れているのに注目
f:id:fakeyakuza:20110728000217p:image
円周率3.14
f:id:fakeyakuza:20110728000624p:image
円周率3.1
f:id:fakeyakuza:20110728000219p:image
円周率3
f:id:fakeyakuza:20110728000220p:image

結論
円周率は3桁未満にすると大変。3.14は妥当なライン。