stirling1.m

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

%stirling1.m
clear
for n=1:170
   f(n)=log(gamma(n+1));
   s(n)=n*log(n)-n;
end
subplot(1,2,1)
plot(f-s)
grid
xlabel('n')
ylabel('Absolute Error')
title('ln(n!) - [nln(n)-n]')
subplot(1,2,2)
semilogx(f./s)
grid
xlabel('n')
ylabel('Relative Error')
title('ln(n!)/[nln(n)-n]')
figure(1)