En Ada, il n'existe pas toujours une fonction «Atn», «ATan» ou «ArcTan» permettant de calculer le ArcTangente d'un cercle. Pour remédier à se problème, il suffit de créer une fonction ressemblant à ceci avec un préalable la fonction de racine carrée :
- WITH TEXT_IO;
-
- PROCEDURE ArcTangente IS
-
- USE TEXT_IO;
-
- R:Float;
-
- FUNCTION Sqrt(X:in Float) RETURN Float IS
- A,B,M,XN:FLOAT;
- BEGIN
- IF X=0.0 THEN
- RETURN 0.0;
- ELSE
- M:=1.0;
- XN:=X;
- WHILE XN>=2.0 LOOP
- XN:=0.25*XN;
- M:=2.0*M;
- END LOOP;
- WHILE XN<0.5 LOOP
- XN:=4.0*XN;
- M:=0.5*M;
- END LOOP;
- A:=XN;
- B:=1.0-XN;
- LOOP
- A:=A*(1.0+0.5*B);
- B:=0.25*(3.0+B)*B*B;
- EXIT WHEN B<1.0E-15;
- END LOOP;
- RETURN A*M;
- END IF;
- END;
-
- FUNCTION ArcTan(X:IN Float) RETURN Float IS
- A,B:Float;
- N:Integer;
- BEGIN
- A := 1.0 / Sqrt(1.0 + (X * X));
- B := 1.0;
-
- N:=1;
- WHILE N<=11 LOOP
- A := (A + B) / 2.0;
- B := Sqrt(A * B);
- N:=N+1;
- END LOOP;
-
- RETURN X / (Sqrt(1.0 + (X * X)) * A);
- END;
-
- BEGIN
- R:=0.0;
- WHILE R<=1.1 LOOP
- PUT_LINE("ArcTan(" & FLOAT'IMAGE(R) & ")=" & FLOAT'IMAGE(ArcTan(R)));
- R:=R+0.1;
- END LOOP;
- END ArcTangente;
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 25 août 2012