s1rout3.m

Industrial Mathematics

% Routine 1.3
q = 0;
m=20;
lT = 15;
p = 1;
for k=1:m
  p = p*lT/k;
  q = q + p;
  end;
q = q*exp(-lT);
1-q