function F=myfun(t)     
  ubound=10000;     
  x=0:0.1:ubound;     
  y=sqrt(x)./(1+exp(x-t));     
  F=trapz(x,y);%梯形积分

clear,clc
F=1.4893 ;   %给定的F值
ta=0;tb=100;
fa=F-myfun(ta);
fb=F-myfun(tb); %fa*fb<0即可，不满足的话可以适当放大区间
error=inf;
while error>0.001   %自定义精度     
    tc=(ta+tb)/2;     
    fc=F-myfun(tc);     
    if fa*fc<0         
        tb=tc;     
    else         
        ta=tc;     
    end     
    error=tb-ta;
end
t=(ta+tb)/2;
t

11楼: Originally posted by yz457694 at 2016-08-20 00:37:51
代码没缩进了，重发一下
函数代码：

function F=myfun(t)     
  ubound=10000;     
  x=0:0.1:ubound;     
  y=sqrt(x)./(1+exp(x-t));     
  F=trapz(x,y);%梯形积分

求解代码：

clear,clc
F= ...

12楼: Originally posted by 星——星 at 2016-08-20 12:30:18
嗯嗯，非常感谢~我还有一个问题想请教下您，那要是正常的从0到无穷这个范围内对x进行积分，是不是也要用这种梯形积分的方法？...