« Excelで拡散方程式を解いてみる(その1)-陽解法 | トップページ | ロジスティック写像の計算比較(Keisan.casio.jpとExcel) »

2008年10月15日 (水)

高精度計算サイトkeisan.casio.jpに協力してみる(その5)。 - 2個のソリトンの衝突(KdV方程式)

よく本には書いてあるんだけど、実際に計算するのが超めんどくさいKdV方程式の解(2個のソリトンがある場合)をUPしてみた。

2個のソリトンの衝突(KdV方程式の解)

超めんどくさいので当然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) »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/42799358

この記事へのトラックバック一覧です: 高精度計算サイトkeisan.casio.jpに協力してみる(その5)。 - 2個のソリトンの衝突(KdV方程式):

« Excelで拡散方程式を解いてみる(その1)-陽解法 | トップページ | ロジスティック写像の計算比較(Keisan.casio.jpとExcel) »

最近のコメント

2018年11月
        1 2 3
4 5 6 7 8 9 10
11 12 13 14 15 16 17
18 19 20 21 22 23 24
25 26 27 28 29 30  
フォト
無料ブログはココログ