代码没骗你，exp(5+para*x)这个函数的定积分是
(e^(5+para*2)-e^(5-para*2))/para，把para看作自变量，这个函数的曲线类似一个开口向上的抛物线，没有最大值。

你把第一行的函数换一个试试，maxLik()是可以出结果的，说明maxLik(), integrate()这些函数都没有问题。

换句话说，你的函数找不到最优参数

谢谢你的回复，确实是我简化的函数设定的问题，我把简化函数改了，但还是会有问题，可以麻烦你再看一下吗？我的目的是估计“para”的值，使U_Cur/（U_Cur+U_Tgt）最大，其中U_Cur和U_Tgt是exp(A)/(exp(A)+exp(B))的形式。这是我重新设定的函数：我通过调试，发现报错的原因是U_Cur或U_Tgt有时可以是inf或0。我把U_Cur/（U_Cur+U_Tgt）加个log再变换，或者变成把U_Cur=exp(A)/(exp(A)+exp(B))变成1-(exp(B)/(exp(A)+exp(B)))等方法都试了。还是会报错。请问还有别的什么办法吗？

我试了新的函数，也用你的初始参数做了函数图像和积分图像，然后把你的函数又作了简化，用一个参数时exp(x)/(exp(x)+exp(a*x))，似乎有结果。改成用两个参数时，报错。我也不知道为什么。抱歉，能力有限。
类似exp(x)/(exp(x)+exp(a*x))的函数，似乎只在某个区间内可以做定积分的，结果基本上是区间上限的值。
猜，可能还是参数对函数积分结果的影响无法确定。

嗯嗯，好的，还是谢谢你，麻烦了。我后来直接给exp(A)/(exp(A)+exp(B))限定了范围，如果出现(-)inf，我就赋值为(-)2^1000，不过运行结果还是会报错"The integral is probably divergent"。我目前打算直接用GA算法求参数值，但就是没办法得到参数的t检验结果了。。再次谢谢你，耽误你时间和精力了，祝好~