Section courante

A propos

Section administrative du site

Après avoir effectué des recherches dans de nombreux livres comme Scientific Pascal, Dictionnaire mathématique,..., je n'ai jamais trouvé aucun livre fournissant une réponse correct du calcul de l'exponentiel (soit l'inverse de Ln ou la méthode exp), outre le projet GNU (HaypoCALC). Je me suis donc basé sur sa formule pour arriver au résultat contenu dans cette page. Voici un code source Turbo C effectuant un calcul correct de cette formule en se basant sur les savants calculs de ce projet :

  1. #include <stdio.h>
  2. #include <stdlib.h>
  3. #include <math.h>
  4.  
  5. double Exponentiel(double x) {
  6.  double q;
  7.  int inverse = 0;
  8.  double n = 0, dl = 1;
  9.  int i = 1;
  10.  if(x<0) {
  11.   inverse = 1;
  12.   x = -x;
  13.  }
  14.  while(x >= 2) {
  15.      x /= 2;
  16.      n++;
  17.  }
  18.  x /= 16;
  19.  n += 4;
  20.  q = x;
  21.  while (q > 1.0E-15) {
  22.     dl += q;
  23.     i++;
  24.     q = q*x/i;
  25.  }
  26.  for(i=1;i<= n; i++) dl=dl*dl;
  27.  if(inverse) dl = 1/dl;
  28.  return dl;
  29. }
  30.  
  31. int main()
  32. {
  33.     double I;
  34.     for(I=0;I<=2.0;I+=0.1) {
  35.           printf("Exp(%f)=%f %f\n",I,Exponentiel(I),exp(I));
  36.     }
  37.     return 0;
  38. }

on obtiendra le résultat suivant :

Exp(0.000000)=1.000000 1.000000
Exp(0.100000)=1.105171 1.105171
Exp(0.200000)=1.221403 1.221403
Exp(0.300000)=1.349859 1.349859
Exp(0.400000)=1.491825 1.491825
Exp(0.500000)=1.648721 1.648721
Exp(0.600000)=1.822119 1.822119
Exp(0.700000)=2.013753 2.013753
Exp(0.800000)=2.225541 2.225541
Exp(0.900000)=2.459603 2.459603
Exp(1.000000)=2.718282 2.718282
Exp(1.100000)=3.004166 3.004166
Exp(1.200000)=3.320117 3.320117
Exp(1.300000)=3.669297 3.669297
Exp(1.400000)=4.055200 4.055200
Exp(1.500000)=4.481689 4.481689
Exp(1.600000)=4.953032 4.953032
Exp(1.700000)=5.473947 5.473947
Exp(1.800000)=6.049647 6.049647
Exp(1.900000)=6.685894 6.685894

Voir également

Science - Mathématique

Dernière mise à jour : Dimanche, le 26 juillet 2015