convergenceOfNfactorial.py

import numpy as numpy
"""This program takes a series
and checks for convergence value
of 1/n!"""
currSum = 1.0
t = 1e-3 # tolerance
prevSum = 0
perDiff = 1.0
n = 0
seriesTerm = 1.0

while perDiff > t:
	n += 1
	prevSum = currSum
	seriesTerm *= 1/n # This saves the individual
	#factorial value terms.
	currSum += seriesTerm # This adds them up.
	perDiff = abs(currSum-prevSum)/currSum
	print(n, seriesTerm, currSum, perDiff)

print("The series converges to ", currSum, " at the ", n, "th iteration.")