Section courante

A propos

Section administrative du site

En C, il existe déjà une fonction permettant de connaitre l'ArcTangente. Cependant, il peut être intéressant d'en reproduire une pour notre plaisir personnel. Pour remédier à se problème, il suffit de créer une fonction ressemblant à ceci avec un préalable la fonction de racine carré :

  1. #include <stdio.h>
  2.  
  3. double SquareRoot(double X) {
  4.   double A,B,M,XN;
  5.   if(X==0.0) {
  6.      return 0.0;
  7.   } else {
  8.      M=1.0;
  9.      XN=X;
  10.      while(XN>=2.0) {
  11.         XN=0.25*XN;
  12.         M=2.0*M;
  13.      }
  14.      while(XN<0.5) {
  15.         XN=4.0*XN;
  16.         M=0.5*M;
  17.      }
  18.      A=XN;
  19.      B=1.0-XN;
  20.      do {
  21.         A=A*(1.0+0.5*B);
  22.         B=0.25*(3.0+B)*B*B;
  23.      } while(B>=1.0E-15);
  24.      return A*M;
  25.   }
  26. }
  27.  
  28. double ArcTan(double X) {
  29.   double A = 1.0 / SquareRoot(1.0 + (X * X)),B = 1.0;
  30.   int N = 1;
  31.   while(N<=11) {
  32.      A = (A + B) / 2.0;
  33.      B = SquareRoot(A * B);
  34.      N++;
  35.   }
  36.   return X / (SquareRoot(1.0 + (X * X)) * A);
  37. }
  38.  
  39. void main() {
  40.     double R=0.0;
  41.     while(R<1.1) {
  42.        printf("ArcTan(%.6f)=%.6f\n",R,ArcTan(R));
  43.        R+=0.1;
  44.     }
  45. }

on obtiendra le résultat suivant :

ArcTan(0.00000)= 0.000000
ArcTan(0.10000)= 0.099668
ArcTan(0.20000)= 0.197396
ArcTan(0.30000)= 0.291457
ArcTan(0.40000)= 0.380506
ArcTan(0.50000)= 0.463648
ArcTan(0.60000)= 0.540420
ArcTan(0.70000)= 0.610726
ArcTan(0.80000)= 0.674741
ArcTan(0.90000)= 0.732815
ArcTan(1.00000)= 0.785398

Dernière mise à jour : Samedi, le 22 août 2015