dobrodeev1970.py 文件源码

def __init__(self, n):
        self.name = 'Dobrodeev'
        self.degree = 7
        self.dim = n

        A = fr(1, 8)
        B = fr(19-5*n, 20)
        alpha = 35*n * (5*n - 33)
        C = fr((alpha + 2114)**3, 700 * (alpha+1790.0) * (alpha+2600.0))
        D = fr(729, 1750) * fr(alpha + 2114, alpha + 2600)
        E = fr(n * (n-1) * (n - 4.7), 3) - 2*n * (C + D) + fr(729, 125)

        a = sqrt(fr(3, 5))
        b = a
        c = sqrt(fr(3, 5) * fr(alpha+1790, alpha+2114))
        data = [
            (A, fsd(n, (a, 3))),
            (B, fsd(n, (b, 2))),
            (C, fsd(n, (c, 1))),
            (D, fsd(n, (1.0, 1))),
            (E, z(n)),
            ]

        self.points, self.weights = untangle(data)
        self.weights /= fr(729, 125 * 2**n)
        return