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 du logarithme, outre le projet GNU (HaypoCALC). Bien qu'il existe de nombreuses fonctions de logarithme du langage Pascal, il peut être intéressant d'effectuer les calculs par nous même :

  1. WITH TEXT_IO;
  2.  
  3. PROCEDURE Logarithme IS 
  4.    
  5.    USE TEXT_IO;
  6.    
  7.    I:Float;  
  8.    
  9.    FUNCTION Sqrt(X:in Float) RETURN Float IS 
  10.       A,B,M,XN:FLOAT;
  11.    BEGIN
  12.       IF X=0.0 THEN 
  13.          RETURN 0.0;
  14.       ELSE
  15.          M:=1.0;
  16.          XN:=X;
  17.          WHILE XN>=2.0 LOOP
  18.             XN:=0.25*XN;
  19.             M:=2.0*M;
  20.          END LOOP;
  21.          WHILE XN<0.5 LOOP
  22.             XN:=4.0*XN;
  23.             M:=0.5*M;
  24.          END LOOP;
  25.          A:=XN;
  26.          B:=1.0-XN;
  27.          LOOP
  28.             A:=A*(1.0+0.5*B);
  29.             B:=0.25*(3.0+B)*B*B;
  30.             EXIT WHEN B<1.0E-15;
  31.          END LOOP;
  32.          RETURN A*M;
  33.       END IF;
  34.    END Sqrt;
  35.    
  36.    FUNCTION Ln(X:IN Float) RETURN Float IS
  37.       Negatif:Boolean;      
  38.       Fois,I:Integer;
  39.       Xt,Ajout,Savx,Xp,Quotient,Dl:Float;
  40.    BEGIN      
  41.       Negatif := False;      
  42.       Fois := 1;      
  43.       Ajout := 0.0;      
  44.       Xt := X;
  45.       IF Xt <= 0.0 THEN         
  46.          RETURN 0.0;         
  47.       END IF;      
  48.       IF Xt < 1.0 THEN         
  49.          Negatif := True;         
  50.          Xt := 1.0 / Xt;         
  51.       END IF;      
  52.       WHILE Xt >= 10.0 LOOP         
  53.          Xt := Xt / 10.0;         
  54.          Ajout := Ajout + 2.302585092994046;         
  55.       END LOOP;      
  56.       WHILE Xt >= 1.1 LOOP         
  57.          Xt := Sqrt(Xt);         
  58.          Fois := Fois * 2;         
  59.       END LOOP;
  60.       Xt := Xt - 1.0;      
  61.       Savx := Xt;      
  62.       I := 2;      
  63.       Xp := Xt * Xt;
  64.       Quotient := (Xp / Float(I));      
  65.       Dl := Xt - Quotient;      
  66.       WHILE 1.0E-15 < Quotient LOOP
  67.          I := I + 1;         
  68.          Xp := Xp * Xt;         
  69.          Dl := Dl + (Xp / Float(I));         
  70.          I := I + 1;         
  71.          Xp := Xp * Xt;         
  72.          Quotient := (Xp / Float(I));         
  73.          Dl := Dl - Quotient;         
  74.       END LOOP;     
  75.       Dl := Dl * Float(Fois);      
  76.       Dl := Dl + Ajout;      
  77.       IF(Negatif)THEN 
  78.          Dl := - Dl;
  79.       END IF;
  80.       RETURN Dl;      
  81.    END Ln;
  82.    
  83. BEGIN        
  84.    I:=0.0;   
  85.    WHILE I <= 2.0 LOOP      
  86.       PUT_LINE("Ln(" & FLOAT'IMAGE(I) & ")=" & FLOAT'IMAGE(Ln(I)));      
  87.       I := I + 0.1;      
  88.    END LOOP;            
  89. END Logarithme;

on obtiendra le résultat suivant :

Ln(0.0)=0.0
Ln(0.1)=-2.302585092994046
Ln(0.2)=-1.6094379124341056
Ln(0.3)=-1.2039728043259357
Ln(0.4)=-0.916290731874156
Ln(0.5)=-0.6931471805599471
Ln(0.6)=-0.5108256237659916
Ln(0.7)=-0.3566749439387316
Ln(0.8)=-0.22314355131420963
Ln(0.9)=-0.10536051565782642
Ln(1.0)=0.09531017980432469
Ln(1.2)=0.18232155679395437
Ln(1.3)=0.2623642644674894
Ln(1.4)=0.3364722366212136
Ln(1.5)=0.40546510810816594
Ln(1.6)=0.47000362924573785
Ln(1.7)=0.5306282510621684
Ln(1.8)=0.5877866649021186
Ln(1.9)=0.6418538861723971

Voir également

Science - Mathématique

Dernière mise à jour : Samedi, le 25 août 2012