def METHOD_NAME(n, charge, *args, **kwargs):
    '''Becke, JCP 88, 2547 (1988); DOI:10.1063/1.454033'''
    if charge == 1:
        rm = BRAGG_RADII[charge]
    else:
        rm = BRAGG_RADII[charge] * .5
    t, w = numpy.polynomial.chebyshev.chebgauss(n)
    r = (1+t)/(1-t) * rm
    w *= 2/(1-t)**2 * rm
    return r, w