The expression is $(a - j \mu_1 b) \div (z - \sqrt{z^2 - a^2 - \mu_1^2 b^2})$
$z = a$, the above expression reduces to $(a - j \mu_1 b) \div (a - \sqrt{- \mu_1^2 b^2})$
One way to refactor the above expression is using $j^2 = -1$, the above expression becomes $(a - j \mu_1 b)\div(a - j \mu_1 b)$, which should be 1 for any value of $\mu_1$...
But when I evaluate numerically using Python (cmath.sqrt()), I get 1 for $\mu_1 = 1, 2, 3....$ but not for $\mu_1 = 2+2j$ or any other complex number.