def linfun(i):
        return np.linspace(0, np.pi / 2., 2 ** i + 1)