pub fn cayley_transform(z: Complex) -> Complex
Cayley transform: maps the upper half-plane to the unit disk.
w = (z - i) / (z + i)