Section courante

A propos

Section administrative du site

En MinGW, il existe déjà une fonction permettant de connaître 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é:

#include <iostream.h>

double SquareRoot(double X) {
  double A,B,M,XN;
  if(X==0.0) {
     return 0.0;
  } else {
     M=1.0;
     XN=X;
     while(XN>=2.0) {
        XN=0.25*XN;
        M=2.0*M;
     }
     while(XN<0.5) {
        XN=4.0*XN;
        M=0.5*M;
     }
     A=XN;
     B=1.0-XN;
     do {
        A=A*(1.0+0.5*B);
        B=0.25*(3.0+B)*B*B;
     } while(B>=1.0E-15);
     return A*M;
  }
}

double ArcTan(double X) {
  double A = 1.0 / SquareRoot(1.0 + (X * X)),B = 1.0;
  int N = 1;
  while(N<=11) {
     A = (A + B) / 2.0;
     B = SquareRoot(A * B);
     N++;
  }
  return X / (SquareRoot(1.0 + (X * X)) * A);
}

int main() {
    double R=0.0;
    while(R<1.1) {
       std::cout << "ArcTan(" << R << ")=" << ArcTan(R) << std::endl;
       R+=0.1;
    }
    return 0;
}

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 7 avril 2018