无法在mathematica中找到二维非线性热方程的解

时间:2013-06-03 13:24:24

标签: wolfram-mathematica 2d equation heat

请帮助我在Wolfram Mathematica计划中找到解决方案。 我有几次检查输入数据的准确性。他们是真的。解决方案必须存在。但Wolfram提供空图(任何时间点)或错误“NDSolve :: eerr”。这是我的代码:

    solution = 
 NDSolve[{D[fun[t, x, y], t] == 
    Exp[-t]*Cos[Pi*y] + D[fun[t, x, y], {x, 2}] + 
     D[fun[t, x, y], {y, 2}], fun[t, 0, y] == 0, fun[t, 1, y] == 0, 
   fun[0, x, y] == 0, (D[fun[t, x, y], y] /. y -> 0) == 
    0, (D[fun[t, x, y], y] /. y -> 1) == 0}, 
  fun[t, x, y], {t, 0, 5}, {x, 0, 1}, {y, 0, 1}]


Plot3D[Evaluate[First[fun[5, x, y] /. solution]], {x, 0, 1}, {y, 0, 
  1}, PlotRange -> All, Mesh -> None, PlotPoints -> 40]

这是错误代码

  

NDSolve :: eerr:警告:缩放的局部空间误差估计值   140.65851971330582 at t = 5.在自变量x的方向上远大于规定的误差容差。格   与15个点的间距可能太大而无法达到预期的效果   准确性或精确度。奇点可能已形成或较小的网格   可以使用MaxStepSize或MinPoints方法指定间距   选项。

请告知在这种情况下可以做些什么。非常感谢提前!

0 个答案:

没有答案