f <- function(h){f15.3.1(1,2,3, z=2, h=h)}
# Winding number [around 1/z] matters:
f(0.5)
f(c(1-1i, 1+1i, -2i))
# Accuracy isn't too bad; compare numerical to analytical result :
f(0.5) - (-1+1i*pi/2)Run the code above in your browser using DataLab