Section courante

A propos

Section administrative du site

En Visual C++, 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é :

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

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