高精度計算サイト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)

ですよ。