我觉得你这个问题没有解。因为-1.95和-2.95都小于0，那么如果d的绝对值比较大，对应的图像的值就会比较小，这样积分不会到1。这是因为半个图像的积分大致在1.25左右。

也可以在两边取导数，按照变上限积分计算，这样计算得到的d=0.411，但是0.411不能满足积分值为1

我觉得是这样，欢迎讨论。

我的原始问题是：
R是一个k×k阶的矩阵，其中每个元素Pij可以表示为：
Pij=Ф[((j*d)-(1-r)*(i-0.5)*d)/(r)]- Ф[((j-1)*d-(1-r)*(i-0.5)*d)/(r)] ；Ф为标准正态概率分布函数，其中，i,j=1,2,…,k;r=0.1
Y=P(I-R)-11
其中，P=[0,…,1,…,0]1×k,P(1,k/2)=1;I为k阶单位阵，1为所有元素均为1的k×1的列向量.
现在已知y的值为500，
对Y进行矩阵计算得到Rk×k=I-1/ARL*Ek×1*P1×k
这样我得到了R矩阵，已知里面任意一个元素，再根据Pij的计算表达式返回去求解d的值，是不是就转化成 这样一个表达式（m1,m2和a分别为常数），我想用通过二分法等求解非线性方程的方法求出d的值，但是首先得找到关于d的目标函数，麻烦各位帮忙看看，我对于matlab了解不多，不知道这个思路对不对？恳请各位指点！！！

你的问题我没看懂，你说你有一个k阶方阵R，后来又写Pij = ... 你是不是写错了啊。

我的问题连接在下面附件的word文档里面了，那个Pij就是一个代号，因为和马尔科夫链有关，这个Pij是转移概率的表示方法，代表R方阵中的每一个元素是怎么得来的，麻烦大侠再给我看看， 指点一下！！！

先 通过变换t=x/d把积分限的d去掉，然后对被积分函数做幂级数展开，最后解关于d的高次方程方程，。。。满意吗？