#!/usr/bin/env python # # taylor series for log(1+x) # import math def myexp(x): if x < 0 : return 1./myexp(-x) sum=0.0 term=1.0 i=1 while (sum + term) != sum: sum+=term term*=x/i i+=1 return sum for x in [0.0,3.0,10.0,30.0,100.0,300.0]: print myexp(x) , math.exp(x) print myexp(-x) , math.exp(-x)