import  numpy as np
def f(z): return (np.exp(z) - 1 - z)/z**2
print(f(1e-8))




def g(z):
	N = 16
	ws = z + np.exp(2j*np.pi*np.arange(N)/N)
	return sum(f(ws))/N
    
    
print(g(1e-8))