>

restart;

>	eq := exp(-x)*x=eps;

>	L1 := ln(1/eps);

>	L2 := ln(L1);

>	eps := 0.01;

>	exact := fsolve(eq, x=L1);

>	x1 := L1; x2 := L1+L2; x3 := L1+L2+L2/L1;

>	rel_err := (%-exact)/exact;

>