公告:1)网站程序升级:Q2A升级到1.8.6,Wordpress升级到5.7.2
2)修复了头像加载慢与提交问题反应慢等问题
2021-06-16

欢迎来到 Mathematica 问答社区

提问时请贴上文本代码

语法高亮:在编辑器中点击

被禁止的话题:广告破解

请阅读:《提问的智慧》

备用域名:mma.ooo

支持LaTex数学公式
行内公式标识符:\$ 或“$\backslash ($”+“$\backslash )$”,
行间公式标识符:\$\$ 或 “$\backslash [$”+“$\backslash ]$”

社区建议QQ群:365716997

分类

0 投票
478 浏览

根据一篇论文里面的代码来求解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,...}}, <>]}}

 

分类:方程 | 用户: wssunyu123 (11 分)
修改于 用户:野鹤
……楼主来先把这个看了:http://tieba.baidu.com/p/2964416898 。顺便提示一下你这段代码有约7个错误。

1个回答

0 投票
sqrt 换成 Sqrt , 再试试看
用户: imma (81 分)
不是啊,说什么东西让保护着呢~
...