A method for constructing the generalized Gaussian quadrature rules for Müntz polynomials on (0, 1) is given. Such quadratures possess several properties of the classical Gaussian formulae (for polynomial systems), such as positivity of the weights, rapid convergence, etc. They can be applied to the wide class of functions, including smooth functions, as well as functions with end-point singularities, such as those in boundary-contact value problems, integral equations with singular kernels, complex analysis, potential theory, etc. The constructive method is based on an application of orthogonal Müntz polynomials introduced by Badalyan [Akad. Nauk Armyan. SSR. Armenian summary)] and studied intensively by Borwein, Erdélyi, and Zhang [Trans. Amer. Math. Soc., 342 (1994), pp. 523-542], as well as by Milovanović on a numerical procedure for evaluation of such polynomials with high precision [Müntz orthogonal polynomials and their numerical evaluation, in