> | 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; |
> |