def pink(N, depth=80):
    """
    N-length vector with (approximate) pink noise
    pink noise has 1/f PSD
    """
    a = []
    s = iterpink(depth)
    for n in range(N):
        a.append(next(s))
    return a