File: stirling.py

package info (click to toggle)
pycode-browser 1%3A1.03-1
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 9,860 kB
  • sloc: python: 2,919; xml: 152; makefile: 76
file content (28 lines) | stat: -rw-r--r-- 953 bytes parent folder | download | duplicates (4)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
# Validity of stirling's approximation up to N=2000
import math

def stirling(N):  #approximation function. log(n!)  = nlog(n)-n
    return N*math.log(N) - N
errors=[]
factN=1
MAXNUM = 2001
print('N\tActual\tApprox\t%Error')
for N in range(1,MAXNUM): # loop from N=1 to N=99
	factN *= N  # factorial is automatically calculated and revised per iteration.	
	actual = math.log(factN)
	approx = stirling(N)
	
	if(N>2): #Actual is 0 at N=1, and error is quite high for N=2 => neglect
		err = 100*(actual-approx)/actual
		errors.append(err)
		#Print values to the screen every 5 iterations
		if N<5: print ('%d\t%.2f\t%.2f\t%.4f'%(N,actual, approx, err))
		elif N<100 and N%20==0: print ('%d\t%.2f\t%.2f\t%.4f'%(N,actual, approx, err))
		elif N%100==0:print ('%d\t%.2f\t%.2f\t%.4f'%(N,actual, approx, err))

#Plot N vs %error 
from pylab import *
# x = [3 ... MAXNUM-1] , y = errors
plot(range(3,MAXNUM),errors)
xlabel('N');ylabel('percentage error')
show()