// calcul de l'epsilon machine
d = 1.0;
while ( 1.0 + d <> 1.0)
    eps = d;
    d = d/2;
end
printf(" epsilon calculé : %20.10e", eps);
printf("epsilon founi par la machine : %20.10e", %eps);