Abstract
Cubature schemes, in this case for uniformly weighted integration over the unit disk, enable exact evaluation of numerical integrals of polynomials but have been explicitly constructed for only low or moderate degrees. In this paper, cubature formulae are discovered for a wider range of degrees by leveraging numerical optimization. These results include a degree-17 cubature scheme with fewer points than existing solutions and up to a degree-77 solution with 1021 cubature points. Optimization heuristics and patterns in the distributions of cubature points are discussed, which serve as vital guides in this work. For example, these heuristics leverage a connection to circle-packing configurations to facilitate the discovery of fully symmetric cubature schemes. (C)& nbsp;2022 The Authors. Published by Elsevier B.V.& nbsp; & nbsp;