在NDSolve中调用含参变量的NIntegrate

自我回复：
我没有找到直接的解决方案，但是我用一个间接的方法实现了我的目标。
方案是使用插值函数：
List1 = Table[f[i], {i, 0, 30, 1]；
f1 = ListInterpolation[List1, {{0, 30}}];
sol=NDSolve[{a'[t]==f1[t],a[0]==1},{a},{t,0,1}];
其中的f1是使用插值函数得到的数值函数，因此可以在ndsolve中使用。

报错是NIntegrate得到非数值结果，出现原因是R作为符号而非数值被传入，解决方法很简单，给f添加一个NumericQ判断即可
ClearAll[f];
f[R_?NumericQ] := NIntegrate[R, {r, 0, 10}]
sol = NDSolve[{a'[t] == f[t], a[0] == 1}, {a}, {t, 0, 1}]
另：在NDSolve中调用NIntegrate的时间开销可能会非常大，所以你的间接解决方案在某种意义上其实是更好的方案，只要注意取点够密让插值函数误差不要太大就可以