1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにルンゲクッタ8次のDOP853を移植したもので計算してみた。なかなか面白い動きをする。
だいぶ前に、
という記事を書いた。
このときも感心していたのだが、先日こんな論文があるのを知った。
なんじゃこれ！力技で1223個も見つけたのか！
とりあえず検証してみよう。ルンゲクッタ8次のDormand&Princeルーチン、DOP853を使っていると論文には書いてあるが、私もちょうどそのフォートランルーチンをExcelのVBAに移植しているのだ（この記事の後にリンクつけます）
では計算しよう。どれがいいかな、、、1223個もあるし、、、
とりあえずなんとなくみかんのあみあみ袋みたいなI.A68(0.5)というのにしよう。
論文の結果と同じになった。GIFアニメにもしてみよう。
1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにDOP853移植ルーチンを使って計算、GIFアニメにしてみた。今回はI.A68(0.5)というもの。えぐい動きをしてるな、、、計算の詳細は続くツイートにて。 pic.twitter.com/BjCu7uJFIQ
— tomo (@tonagai) 2017年9月23日
うーん、なんちゅう動きをしてるんだ（ちょっと表示の刻みが粗すぎるか、、、）
次はIID2(2)ってやつ。
1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAへのDOP853移植ルーチンを使って計算、GIFアニメにしてみた。今回はII.D2(2)というもの。 pic.twitter.com/npOPHjfBSe
— tomo (@tonagai) 2017年9月23日
いろいろ計算しようかと思ったけど1223個か、、、萎える、、、
---
関連リンク：
« #せやねん メチャ売れのコーナー、9/23の特集はエースコックの新食感焼そば モッチッチ。 | トップページ | 味仙監修の名古屋名物台湾ラーメンのカップ麺を食べた。なかなかの辛さでした。 »
「学問・資格」カテゴリの記事
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(3/16) IEEE Microwave Magazineは女性マイクロ波研究者特集、Microwave Journalで285GHz帯で30GHz帯域のOTA測定！QorvoがUWB室内ナビのデモ、STMとonsemiのBluetooth新商品、IDTechExの6Gレポート、など。(2023.03.16)
- KeysightのADSで位相を±180°に限らずに連続にする関数をよく聞かれるがいつも忘れる…unwrap()だ。PythonのNumPyでもあるので（matlabにもある）いい加減に覚えたい。とりあえずPythonでやってみて記憶する。(2023.03.02)
- 50万人が毎年受ける試験で採用、“謎”のプログラミング言語「DNCL」という記事を見た。大学入試センターで2022年1月に出している仕様（日本語プログラムぽい）と、令和７年度試験の問題作成の方向性，試作問題等で出している仕様（Pythonを日本語にしたっぽい）違うのか…(2023.03.01)
- LC共振回路のモンテカルロシミュレーションでL,Cを一様乱数で振って共振周波数を見る、、、と全然一様じゃない。そりゃそうだ。何を勘違いしてたんだろう…(2023.02.27)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(2/21) IEEE Microwave Magazineで学生の設計コンテストの結果（がれきに埋もれた被災者を見つける技術）、Microwave Journalは衛星通信（もうStarlinkとOneWebで4000基も飛んでる！）MWC2023あれこれ（Mobandi, Qualcomm,アナデバ）など。(2023.02.22)
「日記・コラム・つぶやき」カテゴリの記事
- BingのAIチャットに最大の素数はいくつ？と聞くと、存在しないが今知られている最大の素数は2^89589933-1だと答えた。が、これは引用している記事が誤植していて、実際は2^82589933-1だ...間違った記事を引用すると間違った答えになる。英語で聞くと１つ前のメルセンヌ素数に…(2023.03.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース(3/16) IEEE Microwave Magazineは女性マイクロ波研究者特集、Microwave Journalで285GHz帯で30GHz帯域のOTA測定！QorvoがUWB室内ナビのデモ、STMとonsemiのBluetooth新商品、IDTechExの6Gレポート、など。(2023.03.16)
- 今日は円周率の日（3.14)ということでChatGPTとBingのAIチャットに円周率を求める面白い方法を教えてくださいと聞いてみた。ChatGPTが訳の分からない方法をRegenerateするたびに提案してくる...ピザ押し。Bingは検索していることもあってまあこんな感じか。(2023.03.15)
- Nokiaの新ロゴ、シンプルというか書きかけというかアイフルの新ロゴに似てる気がする。そういやKIAとかミューズのロゴもこんな感じになってる。流行ってるのかな。(2023.03.06)
- KeysightのADSで位相を±180°に限らずに連続にする関数をよく聞かれるがいつも忘れる…unwrap()だ。PythonのNumPyでもあるので（matlabにもある）いい加減に覚えたい。とりあえずPythonでやってみて記憶する。(2023.03.02)
トラックバック
この記事へのトラックバック一覧です： 1223個の新しい三体問題の周期解が見つかったということで、私もExcel VBAにルンゲクッタ8次のDOP853を移植したもので計算してみた。なかなか面白い動きをする。:
« #せやねん メチャ売れのコーナー、9/23の特集はエースコックの新食感焼そば モッチッチ。 | トップページ | 味仙監修の名古屋名物台湾ラーメンのカップ麺を食べた。なかなかの辛さでした。 »
コメント