paradox.m

一些常被用于教学或者参考的概率论的实例的源代码

%paradox.m/created by PJNahin for "Duelling Idiots"(3/1/99)
%
%
clear s
clear r
a=input('What is average lifetime of single system: ')
a=1/a;
t=1:200;
for n=1:1:200
    s(n)=1-exp(-a*n);
    r(n)=s(n)*s(n)*(3-2*s(n));
end
grid
plot(t,s,'-',t,r,'.')
ylabel('probability of failure by time t')
xlabel('time t, in arbitrary units')
grid
title('Fig.20.1 - Average Lifetime of a system is 100')
legend('system','SYSTEM')
figure(1)