ex0110.m

《精通MATLAB科学计算》一书的实例源代码

function ex0110()
r = 80; 
M = 160;
p = exp(-r);
S1 = 0;
for k = 1:M
    p=p*r/k;% 
    S1=S1+p;
end
S1
S2= 0;
for k = 1:M
    p = r^k/factorial(k);%
    S2 = S2 + p;
end
S2*exp(-r);
S2