リッカチ方程式をドルマン・プリンスDOP853で計算(Excel VBA)
だいぶ前にリッカチ/リカッチ論争にけりがついたリッカチ(Riccati)方程式ですが、
具体的に解こうと思ってすっかり忘れていた。
まずは、狭義のリッカチ方程式、
y' + a * y^2 = b * x^m
を計算しよう。mが特別な値、
m=-2, 4k/(1-2k) (ただしkは整数)
のときは厳密に初等的に計算できる。(たぶん、こういう話が「リッカチのひ・み・つ」に載っているんだろうと思うが、何せ読んでない。。。ベッセル関数を使えば一般解が求められるが、それもめんどくさい。。。)
でめんどくさいので、まずは数値的に計算。
« マラゴニー効果って何?(トイレスタンプクリーナー)、というかバブル君効果って何? | トップページ | 神足神社でお参り »
「学問・資格」カテゴリの記事
- TensorFlow(Keras)でモデルをsaveで保存してload_modelで読み込むときに突然エラー(utf8で読めないとかディレクトリが存在しないとか)が出始めた。なんで?といろいろやると、単に日本語がフォルダ名に使われているときにエラーになるだけだった…(Windowsネイティブの場合)(2023.05.28)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(4) 10進31桁まで計算できることを生かして4段4次のルンゲクッタ法でローレンツ方程式を計算して図示してみる。テキストベース(アスキーアート)で!(2023.05.23)
- 時代に逆行してCOBOL(GnuCOBOL)を学んでみる(3) ロジスティック写像の分岐図を31桁まで10進計算ができることを生かして描いてみる。テキストベースで!今回は2次元配列の練習。(2023.05.19)
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース(5/18) IEEE Microwave Magazineはグラフェンのマイクロ波応用、Microwave Journalはテラヘルツ波の発生、Next G allianceの新6Gレポート、Hexa-Xの新動画、ESAが軌道上に6Gラボ、De-embeddingって何?など。(2023.05.18)
- カシオの高精度計算サイトkeisan.casio.jpが復旧してログインや自作式の作成ができるようになった。ということで作った自作式のアクセスベスト10とお勧めを紹介。(2023.05.17)
« マラゴニー効果って何?(トイレスタンプクリーナー)、というかバブル君効果って何? | トップページ | 神足神社でお参り »
コメント