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)