# -*- coding: utf-8 -*-
"""
Created on Mon Dec 08 11:13:26 2014

@author: grivet

Listing 2.3 : calcul de la fonction exponentielle
par son développement en série

"""
from pylab import *

eps = 10**(-8); kmax = 40
terme = 1.0; somme = 1.0; k = 1
x = float(input('valeur de x: '))
while k < kmax and abs(terme) > eps:
    terme = x*terme/k
    somme = somme + terme
    print ("%d\t%12.6e\t%12.6f" %(k, terme, somme))
    k = k+1
print ("exp(",x,") = %14.8f" %somme)
print ("valeur fournie par numpy: %14.8f" % (exp(x)))