高精度計算サイトkeisan.casio.jpに協力してみる(その5)。 - 2個のソリトンの衝突(KdV方程式)
よく本には書いてあるんだけど、実際に計算するのが超めんどくさいKdV方程式の解(2個のソリトンがある場合)をUPしてみた。
超めんどくさいので当然Maximaで計算。答えは...
2*((((k1-k2)^2*(2*k2+2*k1)^2*%e^(2*(k2*x-4*k2^3*t+d2)+2*(k1*x-4*k1^3*t+d1)))/(k2+k1)^2+4*k2^2*%e^(2*(k2*x-4*k2^3*t+d2))+4*k1^2*%e^(2*(k1*x-4*k1^3*t+d1)))/(((k1-k2)^2*%e^(2*(k2*x-4*k2^3*t+d2)+2*(k1*x-4*k1^3*t+d1)))/(k2+k1)^2+%e^(2*(k2*x-4*k2^3*t+d2))+%e^(2*(k1*x-4*k1^3*t+d1))+1)-
(((k1-k2)^2*(2*k2+2*k1)*%e^(2*(k2*x-4*k2^3*t+d2)+2*(k1*x-4*k1^3*t+d1)))/(k2+k1)^2+2*k2*%e^(2*(k2*x-4*k2^3*t+d2))+2*k1*%e^(2*(k1*x-4*k1^3*t+d1)))^2/(((k1-k2)^2*%e^(2*(k2*x-4*k2^3*t+d2)+2*(k1*x-4*k1^3*t+d1)))/(k2+k1)^2+%e^(2*(k2*x-4*k2^3*t+d2))+%e^(2*(k1*x-4*k1^3*t+d1))+1)^2)
ですよ。
« Excelで拡散方程式を解いてみる(その1)-陽解法 | トップページ | ロジスティック写像の計算比較(Keisan.casio.jpとExcel) »
「学問・資格」カテゴリの記事
- Wordleを100回やってみて正解率の分布を見た。1回目で当たる確率をpとして、2回目からはk倍ずつ確率が上がっていく(kp,k²p,k³p…)とすると、私の結果はp=0.0089, k=2.4になった。(2022.05.27)
- 高周波回路シミュレータQucsStudioを使ってみる(その5)マルチポート(4ポートだとs4pみたいな)のSパラメータデータをTouchstoneファイルに出力しようとするとExport to SnPでは2ポートになってしまう。Project→Import Dataからポート数を手で入れないと。(2022.05.24)
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2022年5月19日 IEEE Microwave Magazineで磁性体を使わない非可逆デバイス、Microwave Journalは5Gスモールセル、THz、Movanoがミリ波血圧&血糖値計発表、三菱の3Dプリンタ衛星アンテナ、など。(2022.05.19)
- シン・ウルトラマンをIMAXで観てきた。面白かった!物理学者の有岡君がラグランジアン,AdS,プランクブレーンとかしゃべってる!Randall–Sundrumも。マグカップは標準理論+重力で、物理監修は予想通り橋本幸士さんでした!ブラックホールの構造もタイムリー!(2022.05.14)
- 個人の正解確率pとして多数決をとったときの正解確率が、pが0.5に近づくとどのくらいの人数で0.9を超えるか計算してみた。p=0.50005だと日本の人口くらいの多数決。(神とさざなみの密室読んで正当性確率が気になったので)(2022.05.10)
« Excelで拡散方程式を解いてみる(その1)-陽解法 | トップページ | ロジスティック写像の計算比較(Keisan.casio.jpとExcel) »
コメント