So my fellow numerical computational peeps it may be that I am suffering from sleep deprivation but I'm struggling to numerically compute a function x-> u(x) defined implicitly as follows:
x^h(x) + x^-h(x) = a x
To be clear, that's just saying x raised to h(x) etc, and this should hold for all x. We can just consider x>0 and h(x)>0 actually.