mathematica-8 - 在 NDSolve 中使用 NIntegrate

Question

我试图数值求解一个偏微分方程，其中非齐次项是另一个函数的积分。像这样的东西:

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，所以我相信有一个简单的解决方案。
有人能告诉我我做错了什么吗？

谢谢。

Answer 1

从问题中并不是很清楚，但我遇到了类似的问题，在我的情况下，解决方案是遵循 wolfram forums 的建议。将积分放在一个额外的函数中并强制实际输入。

所以在你的情况下，这将是

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}]