用MMA求解一个KDV方程

Question

根据一篇论文里面的代码来求解KDV方程，但是没有成功。附上论文图，和代码，希望高手给看看。

代码：

kdv[{M_}, t_, x_] := D[u[x, t], t] + u*D[u[x, t], x] + M*D[u[x, t,], {x, 3}] == 0
S[{M_, D_, x0_}, t_, x_] := (3 D) sech[1/2*sqrt[D/M] (x - x0 - Dt)]^2

jie=S[{1,3,5},0,x]+S[{1,1,20},0,x]+
S[{1,3,5+50},0,x]+S[{1,1,20+50},0,x]
+S[{1,3,5-50},0,x]+S[{1,1,20-50},0,x];

S1=NDSolve[kdv[{1},t,x]&&u[0,x]=jie&&u[t,0]=u[t,50],u,{t,0,12,0.5},Method
→StiffnessSwitching]
{{u→InterpolatingFunction[{{0,12},{...,0,50,...}}, <>]}}

 

Comments

……楼主来先把这个看了：http://tieba.baidu.com/p/2964416898 。顺便提示一下你这段代码有约7个错误。

Answer 1

sqrt 换成 Sqrt , 再试试看

Comments

不是啊，说什么东西让保护着呢~