Let w(x)=e−xβxα, w¯(x)=xw(x) and let {pm(w)}m, {pm(w¯)}m be the corresponding sequences of orthonormal polynomials. Since the zeros of pm+1(w) interlace those of pm(w¯), it makes sense to construct an interpolation process essentially based on the zeros of Q2m+1:=pm+1(w)pm(w¯), which is called “Extended Lagrange Interpolation”. In this paper the convergence of this interpolation process is studied in suitable weighted L1 spaces, in a general framework which completes the results given by the same authors in weighted Lup((0,+∞)), 1≤p≤∞ (see [31], [28]). As an application of the theoretical results, an extended product integration rule, based on the aforesaid Lagrange process, is proposed in order to compute integrals of the type∫0+∞f(x)k(x,y)u(x)dx,u(x)=e−xβxγ(1+x)λ,γ>−1,λ∈R+, where the kernel k(x,y) can be of different kinds. The rule, which is stable and fast convergent, is used in order to construct a computational scheme involving the single product integration rule studied in [23]. It is shown that the “compound quadrature sequence” represents an efficient proposal for saving 1/3 of the evaluations of the function f, under unchanged speed of convergence.