We deal with the discretization of a 2D radiating panel. After dimensioning the panel by a properly developed approach, a non-uniform array can be a feasible solution for its practical implementation. In the proposed scheme, the determination of the number, the positions and the excitations of the array elements is accomplished by exploiting a Gaussian quadrature rule. Considering that the most significant singular functions supported on the region of interest are related to those supported on the equivalent panel by a radiation integral, the quadrature allows to replace the radiation integrals by weighted summations. The quadrature also provides a unique set of nodes, namely the elements positions, and a unique set of weights, which are involved in the definition of the excitation coefficients. Numerical results assess the effectiveness of the technique.