我试图在数值上求解偏微分方程,其中非齐次项是另一个函数的积分。像这样:
NDSolve [{D [f [x,y],x] == NIntegrate [h [x,y + y2],{y2,x,y}],f [0,y] == 0}, f,{x,0,1},{y,0,1}]
其中h [x,y]是先前定义的众所周知的函数。 但似乎Mathematica不知道如何评估积分。
我不经常使用Mathematica,所以我确信有一个简单的解决方案。 有人能告诉我我做错了吗?
感谢。
答案 0 :(得分:0)
问题并不是很清楚,但我遇到了类似的问题,我的解决方案就是遵循wolfram forums的建议,将积分放在一个额外的函数中并强制实际输入。 / p>
所以在你的情况下,这将是
integral[x_Real,y_Real] := NIntegrate[h[x,y+y2],{y2, x, y}];
NDSolve[{D[f[x, y], x] == integral[x,y], f[0,y] == 0}, f, {x, 0, 1}, {y,0,1}]