Javascriptで数値計算その1(ルンゲクッタ法でローレンツ方程式)
Javascriptの実験を繰り返してきたけど、ようやく本題に入れる。数値計算をやってみよう。このブログではExcel VBAで8次とか14次のルンゲクッタ法を計算したり、複素数の構造体作ってゼータ関数を計算したりとあんまりよそでやってないことをやってきたけど、それはExcelがたいていどのPCにも入っていて超気軽に試せるのと、あまり誰もやってないからプログラムで文句をつけられることもない。。。ということ。
で、さらに気軽にやろうとするとExcelすらインストールされてないPCでもブラウザくらいはあるだろうと。でそこにはJavascriptは使えるはず、ということでしばらくはJavascriptで数値計算やってみよう。
まずは、ローレンツ方程式くらいだろうなあ。
dx/dt = -σ*(x-y)
dy/dt = -y-x*z+r*x
dz/dt =x*y-b*z
をルンゲクッタ4次(普通。。。まずはこれで試させて)でやってみよう。OKを押すと計算して図示する。
※Canvasが現れるまでOK押さないで。
« ココログでJavascriptでCanvas機能(excanvas.js)を使う。 | トップページ | Javascriptで数値計算その2(Runge-Kutta8次のDormand-Prince法) »
「パソコン・インターネット」カテゴリの記事
- Interface2025年8月号Pythonで体験!はじめての暗号を買った。上杉暗号からRSA、AES、DHなど、特に楕円曲線暗号についてはコードも実際に動かすところまで詳しくかかれていた。耐量子暗号や聞いたことなかったY-00暗号や関数型暗号も記載。(2025.07.10)
- Gemini CLIが使えるようになっていたので早速VSCodeのターミナルから使って、JavaScriptで連立一次方程式を計算するコードを書いてもらった。普通にガウスの消去法で計算するhtmlを作ってくれた。(2025.06.27)
- Google ColabのJulia言語で1次元のGray-Scottモデル(∂u/∂t=u²v-(F+k)u+Du∂²u/∂x²,∂v/∂t=-u²v+F(1-v)+Dv∂²v/∂x²)を計算してパルスが次々分裂する様子を見る。空間6次の差分、時間8次のルンゲクッタ法で計算。(2025.07.08)
「学問・資格」カテゴリの記事
- Interface2025年8月号Pythonで体験!はじめての暗号を買った。上杉暗号からRSA、AES、DHなど、特に楕円曲線暗号についてはコードも実際に動かすところまで詳しくかかれていた。耐量子暗号や聞いたことなかったY-00暗号や関数型暗号も記載。(2025.07.10)
- 高周波・RFニュース 2025年7月8日 NordicとSercommのセルラーIoTモジュール、iFixitがFairphone 6を分解、スコアは10/10、RCR wireless newsのウェビナー2件(6GとIndustry4.0)、SEMCOが高耐圧C0G MLCCを車載急速充電に提案(2025.07.09)
- 高周波・RFニュース 2025年7月2日 5G Americasが6Gに向けセンシングと通信ホワイトペーパー発行、KYOCERA AVXが3dBハイブリッドカプラ発表、TDKが車載薄膜インダクタ発表、Nordicが1次電池向けPMIC発表、ローデ・シュワルツの6GとAI/ML解説記事(2025.07.02)
- 高周波・RFニュース 2025年6月30日 QualcommがAIを用いた6Rxアンテナ解説、Next G Allianceと日本のXGMFが5G,6Gで協力、5G Americasが25Q1で5G加入者増加と発表、TechInsigtsがHuawei Pura 80 Pro+分解、Qorvoが5-7GHzのWi-Fi 7 FEM発表(2025.06.30)
- Gemini CLIが使えるようになっていたので早速VSCodeのターミナルから使って、JavaScriptで連立一次方程式を計算するコードを書いてもらった。普通にガウスの消去法で計算するhtmlを作ってくれた。(2025.06.27)
« ココログでJavascriptでCanvas機能(excanvas.js)を使う。 | トップページ | Javascriptで数値計算その2(Runge-Kutta8次のDormand-Prince法) »
コメント