I need to code up the correct way my computer (in python) views the function $\sqrt(z^2+1)$. I want the branch cuts to extend from $i \to i\infty$ and from $-i \to -i\infty$.
I was hoping someone could suggest a way to write some pseudo code on how to do this, or give me a hint. It is more the implementation part I am having troubles understanding.
Edit: I do have some code I, but not sure this is the best place to post such a question. If someone could point me in the right direction that would also be helpful- i.e which stackexchange site.