It is my impression that if we find a function f(z) that satisfies
$$f(f(z)) = e^z $$
there is only one point z that satisfies the relation.
This dawned on me when I noticed that the pesky z that kept popping up in my attempts to look at the problem was the one my book proposed I start with, to wit: $z_o = 0.318 + 1.337i.$ So the joke was on me.
Now I would like to prove this. I would instinctively begin by assuming there was a $z \neq z_o$ and deriving a contradiction. Hopefully I will make some progress before an answer is posted but I am sure I will miss nuances. Maybe it's as simple as showing that $\log^nz$ has a fixed point, which I don't know to be true.
Thanks for any insights.
Daniel,
One solution would be the half iterate generated from real valued tetration, but one can also start with the $z_0\approx0.318+1.337i$ fixed point, and develop the half iterate directly from there. Of course, such a solution is not real valued at the real axis. Then, $z_0$ is defined such that $\exp(z_0)=z_0$. Then there is a Schroeder function $s(z)$, and its inverse, $s^{-1}(z)$. The Schroeder function for $\exp(z)$ has $\lambda=z_0\approx0.318+1.337i$; where $s(z_0)=0$, and s(z) is a taylor series developed in the neighborhood of $z_0$. And the inverse of the Schroder function, $s^{-1}(0)=z_0$.
$s(\exp(z))=\lambda s(z)\;\;\;\;\;\;\;\;\;\;\;s(z) = (z-z_0) + a_2(z-z_0)^2 + a_3(z-z_0)^3 + ... $
$s^{-1}(\lambda z)=\exp(s^{-1}(z))\;\;\;\;s^{-1}(z) = z_0 + z + b_2 z^2 + b_3 z^3 + ...$
Finally, the half iterate of f(z) that you might be looking for would be $h(z) = s^{-1}(s(z) \times \lambda^{0.5} )$ $h(h(z)) = \exp(z)$
The solution for h(z) isn't real valued at the real axis, and h(z) has a singularity at z=0. Nonetheless, the half iterate of say a number like 1/2 is defined. $h(0.5)\approx 0.99303919280011+0.1311428457124i$
A completely different solution involves real valued tetration, which actually involves both $z_0$ and $\overline{z_0}$. This is Kneser's solution, which involves both fixed points, and a Riemann mapping. I wrote a pari-gp program that purports to calculates Kneser's solution. You can download the pari-gp code from http://math.eretrandre.org/tetrationforum/showthread.php?tid=486 Call the abel function for the slog, or inverse of tetration, $\alpha(z)$, and the inverse abel function $\alpha^{-1}(z)$ is tetration itself.
Then $h(z)=\alpha^{-1}(\alpha(z)+0.5)$. With this solution, $h(0)$ is defined, and the half iterate of 0 is: $h(0)\approx 0.498563287941114434679619$
$h(h(z))=\exp(z)$