Nchimy
Solving y=x^x for x does have a solution in terms of Lambert W Function. Although this reply may have come late, one could still make use of it later.
To make the problem easy, let's assume that x and y are real. Then we can proceed as follows:
If y = x^x, ...(1)
then since we also have
x = exp[ln[x]], ...(2)
we may conviniently express (1) as
y=x^exp[ln[x]]. ...(3)
Now, taking logs on both sides of (3) (and noting that ln[d^c]=cln[d]) gives
ln[y] = ln[x] * exp[ln[x]] ... (4).
Imediately, we see that (4) can be solved for ln[x], using Lambert W Function, as
ln[x] = W[ln[y]] ... (5)
so that
x = exp[W[ln[y]]] ... (6) .
The solution in (6) is valid for y > 0 because ln[y] for real values of y makes sense only within this range. From (6) (and taking note that W[0]=1while W[e] =1) we see that
(a) x = 0 when y = 1,
(b) x = 1 when y = e
(c) x is only real when ln[y] >= -(1/e) (or equivalently x is real for y >= exp[-(1/e)) but complex and multivalued otherwise.
(d) x increases monotonically with increasing value of y.
(e) The exists a taylor series expansion of x, about ln[y]=0, with a radius of converges equal to -(1/e).
For more details on th Lambert W function, please refer to Corless R M et al, “On the Lambert W function”, Adv. Comput. Math, Vol. 5, pp.329-359.
Cheers...