裏 RjpWiki

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

フィボナッチ数列の一般項

2019年06月20日 | ブログラミング

「フィボナッチ数列の一般項の実数演算結果がExcelとLibreOffice Calcで違う」ということだ。

まあ,そうなんだろうけど

議論の展開とは別だが,実際問題として,結果は整数になって欲しいんだから,演算誤差があろうと,整数に丸めてやれば何の問題もない。というか,整数に丸めない方が異常だ。

ちゃんとやっていれば,実数演算結果が違うってことに気づかなかったはずだワン。

n = 70
A = 1:n
B = round((((1+sqrt(5))/2)^A-(((1-sqrt(5))/2)^A))/sqrt(5))
C = numeric(n)
C[1] = C[2] = 1
for (i in 3:n) {
    C[i] = C[i-2] + C[i-1]
}
B==C
B[n]

n = 70 まで正しい。

コメント    この記事についてブログを書く
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする
« 防衛省の調査ミス(その3) | トップ | フィボナッチ数列の一般項(... »
最新の画像もっと見る

コメントを投稿

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