.On utilise une représentation de type7.11
k= 0; y= 1; r= 0;
TantQue (k leq 4)
TantQue (x < y * 10^(-k)
k= k+1;
FinTant;
xp= x-y*10^(-k);
y = y+x*10^(-k);
x= xp;
r= r+A[k]; //où A[k] = arctan(10^(-k))
FinTant; //à ce stade: arctan(x/y)= epsilon.
Rendre r+(x/y); //le terme suivant est -(x/y)^3/3
//réponse exacte: r+arctan(x/y)