計算誤差については人一倍気を遣っていたつもりなのだけど,別の問題をチェックしていてキモを冷やした。
以下を見て,何とも思わない人は幸せな人なんだろうか。
> c(sin=sin(pi), cos=cos(pi), tan=tan(pi))
sin cos tan
1.224647e-16 -1.000000e+00 -1.224647e-16
> c(sin=sin(pi/2), cos=cos(pi/2), tan=tan(pi/2))
sin cos tan
1.000000e+00 6.123234e-17 1.633124e+16
> c(sin=sin(pi*3/2), cos=cos(pi*3/2), tan=tan(pi*3/2))
sin cos tan
-1.000000e+00 -1.836970e-16 5.443746e+15
> c(sin=sin(pi*4/2), cos=cos(pi*4/2), tan=tan(pi*4/2))
sin cos tan
-2.449294e-16 1.000000e+00 -2.449294e-16
> c(sin=sin(pi*(4/2)), cos=cos(pi*(4/2)), tan=tan(pi*(4/2)))
sin cos tan
-2.449294e-16 1.000000e+00 -2.449294e-16
> c(sin=sin(pi*2), cos=cos(pi*2), tan=tan(pi*2))
sin cos tan
-2.449294e-16 1.000000e+00 -2.449294e-16
※コメント投稿者のブログIDはブログ作成者のみに通知されます