stirling2.m

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

%stirling2.m
clear
for n=1:150
   f(n)=gamma(n+1);
   s(n)=exp(-n)*(n^50)*(n^(n-50));
end
subplot(1,2,1)
semilogy(f-s)
grid
xlabel('n')
ylabel('Absolute Error')
title('n! - [n^ne^-^n]')
subplot(1,2,2)
plot(f./s)
grid
xlabel('n')
ylabel('Relative Error')
title('n!/[n^ne^-^n]')
figure(1)