An algorithm to get an optimal choice for the number of symmetric quadrature points is given to find symmetric quadrature points is given to find symmetric quadrature for-mulas over a unit disk with a minimal number of points even when a high degree of polynomial precision is required. The symmetric quad-rature formulas for numerical integration over a unit disk of complete polynomial functions up to degree 19 are presented.