zi = np.linspace(0, 3, 47)