Scratchで新しい三体問題の周期解（1223個）を計算してみる、、、計算できたが遅すぎる。8の字解を計算したものを代わりに。。。
この前、三体問題の周期解が1223個も見つかったという論文を見て、、、
ExcelのVBAにルンゲクッタ8次のDOP853 (Dormand＆Prince)を移植したもので計算したりしていたが！
今度はScratchで計算してみよう。II.D2（2）というやつだが、、、
計算はできたんですよ、計算は。
でもDOP853はさすがに移植できず、シンプレクティックオイラー法を使ったら
刻み幅をとんでもなく小さくしないとすぐに発散する。
なので、ターボモードですら遅すぎる、、、
そのためちょっとアップロードは保留。
代わりに8の字解を計算したものを再掲。こっちは速い！
« 9/24 20時過ぎからJRの三ノ宮駅に人が立ち入り、JR神戸線と京都線が遅れてる！ | トップページ | 35段14次のルンゲクッタ法をPARI/GPに実装、ローレンツ方程式を計算し、通常の4次、そして8次、オイラー法と比較してみる。（再掲） »
「パソコン・インターネット」カテゴリの記事
- Visual Basic （VB.NET)でC#用の数値計算ライブラリMath.NET Numericsを使う(1)複素行列を定義して一次方程式や逆行列、行列式などを計算する。(2023.03.24)
- 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)
- BingのAIチャット、2変数の連立方程式も解けないし、指摘したら切れられた。WolframAlpha先生を見習ってほしい。いやWolframさんはChatGPTとの融合に可能性を感じているということなので、そのうち計算できるようになるんだろうな。(2023.02.18)
「学問・資格」カテゴリの記事
- Visual Basic （VB.NET)でC#用の数値計算ライブラリMath.NET Numericsを使う(1)複素行列を定義して一次方程式や逆行列、行列式などを計算する。(2023.03.24)
- 高周波（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)
「日記・コラム・つぶやき」カテゴリの記事
- Visual Basic （VB.NET)でC#用の数値計算ライブラリMath.NET Numericsを使う(1)複素行列を定義して一次方程式や逆行列、行列式などを計算する。(2023.03.24)
- 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)
« 9/24 20時過ぎからJRの三ノ宮駅に人が立ち入り、JR神戸線と京都線が遅れてる！ | トップページ | 35段14次のルンゲクッタ法をPARI/GPに実装、ローレンツ方程式を計算し、通常の4次、そして8次、オイラー法と比較してみる。（再掲） »
コメント