8楼: Originally posted by dbb627 at 2012-07-31 08:52:14
s的取值是多少呢？，如果alpha依次取0~97%对应计算出t的值是：t<20,t=21,t=23;t=25;t=27;t=29;t=31;t=33;t=35;t=37;t=39，那么计算结果应该不会解不出来...

function yan
global s phi lam alpha
s=3.2;
phi=1;
lam=3;
alpha=0.6;
[X,FVAL,EXITFLAG,OUTPUT]=fsolve(@sumt,10)

function y=sumt(t)
global s phi lam alpha
y=0;
for N=0:floor(s*t-phi)
    y=y+exp(-lam*t)*(lam*t)^N/gamma(N+1);
end
y=y-alpha;