Scratch(プログラム言語)でルンゲクッタ法(4段4次)を実装した。まずはローレンツ方程式を計算。
Scratchで数値計算シリーズもようやくルンゲクッタまで来たよ。
うちのブログではおなじみ、8次(Dormand&Prince DOP853)や14次のルンゲクッタ法を実装しようかと思ったが、あまりにめんどくさいので断念。かわりにおなじみの4段4次のふつーーーのRunge-Kuttaを実装した。
こちら。
http://scratch.mit.edu/projects/11698142/
今回、初めてMake a blockという自分でブロックを作れる機能で、呼び出す関数を実装してみた。
最初に解くのは、、、もちろん、ローレンツ方程式。
dx/dt = -σ*(x-y)
dy/dt=-y-x*z+r*x
dz/dt=x*y-b*z
これだけいろんな人がプロジェクト作っていて、誰もルンゲクッタを実装していないというのは微分方程式ファンの私には驚き、、、
さて、これを応用していろいろやってみるか(続く)。
« 天神社でお参り。 | トップページ | 仮面ライダーのプロデューサー、平山亨さんがお亡くなりに、、、ご冥福をお祈りします。 »
「パソコン・インターネット」カテゴリの記事
- RF Weekly Digest (Gemini 3.1 Pro・Google AI Studio BuildによるAIで高周波・RF情報の週刊まとめアプリ)2026/4/12-4/19(2026.04.19)
- Qwen3.6-35B-A3Bが発表され、Ollamaでも使える。そこで電子レンジの動作原理(2.45GHzは水分子の共振周波数でない)と隕石が大気圏突入で燃える原理(摩擦熱ではない)を聞くと、誘電緩和と断熱圧縮について正しく答えられた。今までのローカルLLMで一番賢い回答と思う。(2026.04.17)
- ExcelのOfficeスクリプト(TypeScript)で数値計算ライブラリmath.jsを使う(1) Officeスクリプトは外部API呼び出せるし、math.jsは RESTful APIで呼び出せることがわかった。まずは選択したセルのデータを読み、行列演算。LU分解で一次方程式を解き、逆行列と行列式を求める。(2026.04.17)
- RF Weekly Digest (Gemini 3.1 Pro・Google AI Studio BuildによるAIで高周波・RF情報の週刊まとめアプリ)2026/4/5-4/12(2026.04.12)
- GLM-5.1(Ollamaから利用)でPythonのscikit-rfを使ってTouchstoneフォーマットのSパラメータファイルを読んでdB, 位相, スミスチャート, TDRを表示するGUIアプリを作ってもらった。5分など長く考えた後、Gemma 4:31bよりさらに出来が良く、思った通りのものができた。(2026.04.09)
« 天神社でお参り。 | トップページ | 仮面ライダーのプロデューサー、平山亨さんがお亡くなりに、、、ご冥福をお祈りします。 »



コメント