Les Chevaliers de Malte/Développeur - Perl - Atn/ATan/ArcTan/ArcTangente


En Perl, 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é:
#!/usr/bin/perl

sub SquareRoot($X) {
  my($A,$B,$M,$XN);
  if(@_[0]==0.0) {
     return 0.0;
  } else {
     $M=1.0;
     $XN=@_[0];
     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;
  }
}

sub ArcTan($X) {
   my($A,$B,$N);
   $A = 1.0 / &SquareRoot(1.0 + (@_[0] * @_[0]));
   $B = 1.0;
   $N = 1;
   while($N<=11) {
     $A = ($A + $B) / 2.0;
     $B = &SquareRoot($A * $B);
     $N++;
   }
   return @_[0] / (&SquareRoot(1.0 + (@_[0] * @_[0])) * $A);
}

$R=0.0;
while($R<=1.1) {
   print "ArcTan(".$R.")=".&ArcTan($R)."\n";
   $R+=0.1;
}

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: Mardi, le 7 février 2006