为什么要用fsolve那么复杂？这个积分可以自己手算的诶……

把所有p和r带下标的玩意统统视为常数，也就是被积分函数是
\[\frac{2b - ar + r + \sqrt{ar - r - C + 4abrD}}{2br}\]
然后，将分子拆开，从而变成
\[\frac{1}{r} + \frac{1 - a}{2b} + \frac{\sqrt{(a - 1 + 4abD) r - C}}{2br}\]
前两项都可以求出积分来，下面考虑第三项

那么记常数\[\\k = a - 1 + 4 abD\\\]，并做变量代换得\[\\t = \sqrt{(a - 1 + 4abD)r - C}\\\]
因此
\[r = \frac{t^2 + C}{k}\]
\[dr = 2ktdt\]
上面被积函数第三项变成\[\int 2k^2 \frac{t^3}{t^2 + C} dt\]
这个积分你要是算不出我也没办法了……

BTW，有时候多算一点，matlab毕竟是数值运算，直接解误差一般比较大~

非常感谢！你这个处理方法对我的文章非常有借鉴之处！我准备使用！！
程序昨晚做出来了
function ff=djf(p1)
pe=10;pw=4;la=0.001;lc=0.2;a=(la-lc)/la;b=1/la;re=100;rw=0.1;

  g=@(r,p1)((2*b-a*r+r+sqrt((a*r-r-b*(p1-pw)/log(re/rw)).^2+4*a*b*r*(p1-pw)/log(re/rw)))./(2*b*r));
   ff=pe-pw-quadl(@(r)g(r,p1),0.1,100);
end

   
   

%主程序调用文件
% [x,feval]=fsolve(@djf,6,optimset('Display','off','TolX',1e-8,'TolFun',1e-8))