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.000000ArcTan(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