« 9/24 20時過ぎからJRの三ノ宮駅に人が立ち入り、JR神戸線と京都線が遅れてる! | トップページ | 35段14次のルンゲクッタ法をPARI/GPに実装、ローレンツ方程式を計算し、通常の4次、そして8次、オイラー法と比較してみる。(再掲) »

2017年9月25日 (月)

Scratchで新しい三体問題の周期解(1223個)を計算してみる、、、計算できたが遅すぎる。8の字解を計算したものを代わりに。。。

この前、三体問題の周期解が1223個も見つかったという論文を見て、、、

ExcelのVBAにルンゲクッタ8次のDOP853 (Dormand&Prince)を移植したもので計算したりしていたが!

今度はScratchで計算してみよう。II.D2(2)というやつだが、、、
Scratch1223
計算はできたんですよ、計算は。
でもDOP853はさすがに移植できず、シンプレクティックオイラー法を使ったら
刻み幅をとんでもなく小さくしないとすぐに発散する。
なので、ターボモードですら遅すぎる、、、
そのためちょっとアップロードは保留。
代わりに8の字解を計算したものを再掲。こっちは速い!
Scratch8

« 9/24 20時過ぎからJRの三ノ宮駅に人が立ち入り、JR神戸線と京都線が遅れてる! | トップページ | 35段14次のルンゲクッタ法をPARI/GPに実装、ローレンツ方程式を計算し、通常の4次、そして8次、オイラー法と比較してみる。(再掲) »

パソコン・インターネット」カテゴリの記事

学問・資格」カテゴリの記事

日記・コラム・つぶやき」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

« 9/24 20時過ぎからJRの三ノ宮駅に人が立ち入り、JR神戸線と京都線が遅れてる! | トップページ | 35段14次のルンゲクッタ法をPARI/GPに実装、ローレンツ方程式を計算し、通常の4次、そして8次、オイラー法と比較してみる。(再掲) »

最近の記事

最近のコメント

2025年2月
            1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28  
フォト
無料ブログはココログ