Solution (source code)

= Solution

Multiplying $(z+\sqrt{z^2-c^2})(z-\sqrt{z^2-c^2})=c^2$ gives the reciprocal identity and hence $z=c(\zeta+\zeta^{-1})/2$.