next up previous contents
Next: . Up: L'agorithme de CORDIC sur Previous: .   Contents

$\exp$.

On écrira cette fois:

\begin{displaymath} \exp(x)= (1+1)^{n_0} \times (1+1/10)^{n_1}\times \cdots \times (1+1/10^6)^{n_6} \times (1+\epsilon) \end{displaymath}

k= 0; y= 1;        // y est le produit partiel précédent
TantQue (k <= 6)
 TantQue (x <= L[k])
  x= x-L[k] ;
  y= y+ytimes10^(-k);
 FinTant;
 k= k+1;           // à la fin on a:
FinTant; // exp(x)=(1+ epsilon) donc epsilon= exp(x)-1

Rendre y+ytimes x; //le terme suivant est + y * x^2/2!
                   //réponse exacte: y * exp(x)

Faq de fr.sci.maths 2003-12-14