def METHOD_NAME():
    # TODO: first implementation to generate wake overlap lookup table
    # (currently supplied by turbopark_lookup_table.mat.)
    # However, the result of this function doesn't generate the same
    # interpolant as the .mat file, so if used, needs to be corrected.
    dist = np.arange(0, 10, 1.0)
    radius_down = np.arange(0, 20, 1.0)
    overlap_gauss = np.zeros((len(dist), len(radius_down)))
    for i in range(len(dist)):
        for j in range(len(radius_down)):
            if radius_down[j] > 0:
                def fun(r, theta):
                    return r * np.exp(
                        -1 * (r ** 2 + dist[i] ** 2 - 2 * dist[i] * r * np.cos(theta)) / 2
                    )
                out = integrate.dblquad(fun, 0, radius_down[j], lambda x: 0, lambda x: 2 * np.pi)[0]
                out = out / (np.pi * radius_down[j] ** 2)
            else:
                out = np.exp(-(dist[i] ** 2) / 2)
            overlap_gauss[i, j] = out
    return dist, radius_down, overlap_gauss