« 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) »

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

コメント

コメントを書く

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

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

最近の記事

最近のコメント

2019年10月
    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 31    
フォト
無料ブログはココログ