裏 RjpWiki

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

ダメ出し:簡明直截なプログラミングをすれば,誤りもなくなる

2012年08月13日 | ブログラミング

円を描く」にて

t <- seq(0, 2, pi/1000) * pi
x <- cos(t)
y <- sin(t)
plot(x, y, type = "l", asp = 1)

などと書いているが,

なんで,t <- seq(0, 2, pi/1000) * pi なのだろうか?意味不明である。
引数か何か誤解があるのでは?

それに,描かれた円は,閉じていない!!!

座標(1,0) のあたりを拡大すると,以下のようになっている!!

こんなことがないように,

t <- seq(0, 2*pi, length=1000)

でよいだろ

t <- seq(0, 2*pi, by=0.001)

は端点の処理という点で(t に 2*pi が含まれない),元のプログラムと同じで,今回の目的には合わない。

> t <- seq(0, 2, pi/1000) * pi
> sin(t[length(t)])
[1] -0.00611687              # 0に近くないとマズイ

> t <- seq(0, 2*pi, length=1000)
> sin(t[length(t)])
[1] -2.449294e-16            # ちゃんとやればこうなる。こうならないとだめ。

コメント    この記事についてブログを書く
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする
« ダメ出し:ちゃんとしたグラ... | トップ | ダメ出し:初心者の手本にな... »
最新の画像もっと見る

コメントを投稿

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