mode(-1);

deff('y=f1(x)','y=(exp(x)-1)./x')     // Implémentation "naive" 
deff('y=f2(x)','z=exp(x);if (z==1) then y=1; else y=(z-1)./log(z);end') // Implémentation plus stable

x=logspace(-6,-15,10);      // 10 valeurs entre 1e-6 et 1e-15

r1=[x', f1(x)' (f1(x)-1)'];
r2=[x', f2(x)' (f2(x)-1)'];

mprintf('Version naturelle\n     x \t            f1(x)  \t         f1(x)-1\n')
mprintf(' %.2e \t   %.15f \t  % .4e\n', r1)

mprintf('\n\nVersion stable\n     x \t            f2(x)  \t         f2(x)-1\n')
mprintf(' %.2e \t   %.15f \t  %.4e\n', r2)
printf('\n');