# -*- coding: utf-8 -*-
"""
Created on Mon Dec 08 11:29:10 2014

@author: grivet

Listing 2.4 : Approximation de Padé pour l'exponentielle
"""
import matplotlib.pyplot as plt
import numpy as np

x = np.linspace(0.0,1.0,100)
poly = 1.0 + x + x*x/2 + x*x*x/6 + x*x*x*x/24
pade = ((x+6.0)*x+12.0)/((x-6.0)*x+12) 
plt.plot(x,poly-np.exp(x),'-b',label="erreur pour le polynôme")
plt.plot(x,pade-np.exp(x),'-r',label="erreur pour l'approximant de Padé")
plt.legend()