One approach is to use Newton's method on a function known to have a zero at $e$—for instance, $f(x) = \ln x - 1$. Then $f'(x) = 1/x$, and we can iterate as
$$
x_{n+1} = x_n (2-\ln x_n)
$$
Of course, this isn't very satisfying, since the natural log is sitting right there in the expression. (I also haven't looked at the interval of convergence.) But it is available on most scientific calculators.
I'll have to give this more thought to see if something clever can be done to avoid the natural log (or exp, for that matter).
Also, as I mentioned in the comments, this question has been looked at before on Math.SE, with the nice result
$$
x_{n+1} = x_n - \frac1n x_{n-1}
$$
with $x_1 = x_2 = 1$. This is not exactly in the form desired by the OP, but it's pretty nice all the same.