Each result cited is universally quantified over the data in its own statement. Integrals against measures on (Rn,B(Rn)) are those of Measure Spaces and the Lebesgue Integral: Standing Notation §integral, as fixed for probability measures in Probability Measures on Euclidean Space and Random Vectors: Standing Notation §measures, and a Borel map into [0,∞] is measurable in the sense of Measure Spaces and the Lebesgue Integral: Standing Notation §measurable (Probability Measures on Euclidean Space and Random Vectors: Standing Notation §borel-maps); so Nonnegative Combinations of Two Finite Measures, and the Average of Two Couplings §combination applies with E=Rn and E=B(Rn). For ρ′∈P(Rn) and B∈B(Rn), ρ′(B)≤ρ′(Rn)=1 by claim 2 of Basic Properties of a Measure, so ρ′(B) is a real number in [0,1]; in particular all sums below are finite sums of real numbers.
Step 1 (Induction for clause 1). For r∈[M] let γr:B(Rn)→R be given by γr(B)=∑i=1rciρi(B). We prove by induction on r∈[M] the statement H(r): (a) γr is a measure on (Rn,B(Rn)) with γr(Rn)<∞; (b) for every Borel f:Rn→[0,∞], if ∫fdρi<∞ for every i∈[r] then ∫fdγr=∑i=1rci∫fdρi, and if ∫fdρi=∞ for some i∈[r] with 0<ci then ∫fdγr=∞; (c) every Borel f:Rn→R that is integrable with respect to every ρi, i∈[r], is integrable with respect to γr, and ∫fdγr=∑i=1rci∫fdρi.
Base case r=1. By claim 1 of Properties of Finite Sums, γ1(B)=c1ρ1(B)=c1ρ1(B)+0⋅ρ1(B). Apply Nonnegative Combinations of Two Finite Measures, and the Average of Two Couplings §combination with α=β=ρ1 (finite, as ρ1(Rn)=1), s=c1 and t=0: γ1 is a measure with γ1(Rn)<∞, giving (a). For Borel f≥0 the same clause gives ∫fdγ1=c1∫fdρ1+0⋅∫fdρ1 in [0,∞], where the second term is 0 even when ∫fdρ1=∞, by the convention 0⋅∞=0 of Measure, Measure Space, and Probability Measure; so ∫fdγ1=c1∫fdρ1. If ∫fdρ1<∞ this is ∑i=11ci∫fdρi by claim 1 of Properties of Finite Sums; if ∫fdρ1=∞ and 0<c1 it is c1⋅∞=∞ by Measure Spaces and the Lebesgue Integral: Standing Notation §extended. This is (b). For (c), the last assertion of Nonnegative Combinations of Two Finite Measures, and the Average of Two Couplings §combination (with the same data) shows that f is integrable with respect to γ1 with ∫fdγ1=c1∫fdρ1+0=∑i=11ci∫fdρi.
Inductive step. Let r∈[M] with r+1∈[M], and assume H(r). By claim 1 of Properties of Finite Sums, γr+1(B)=γr(B)+cr+1ρr+1(B)=1⋅γr(B)+cr+1ρr+1(B). Apply Nonnegative Combinations of Two Finite Measures, and the Average of Two Couplings §combination with α=γr (a measure with γr(Rn)<∞ by H(r)(a)), β=ρr+1, s=1 and t=cr+1. It gives (a) for r+1, and for Borel f≥0
∫fdγr+1=1⋅∫fdγr+cr+1∫fdρr+1in [0,∞].
If ∫fdρi<∞ for every i∈[r+1], then ∫fdγr=∑i=1rci∫fdρi is real by H(r)(b), and the right side equals ∑i=1r+1ci∫fdρi by claim 1 of Properties of Finite Sums. If ∫fdρi=∞ for some i∈[r+1] with 0<ci, then either i∈[r], in which case ∫fdγr=∞ by H(r)(b) and the first term is 1⋅∞=∞, or i=r+1, in which case the second term is cr+1⋅∞=∞; in both cases the sum is ∞ by the conventions of Measure Spaces and the Lebesgue Integral: Standing Notation §extended. This is (b) for r+1. For (c), if f is integrable with respect to every ρi with i∈[r+1], then it is integrable with respect to γr by H(r)(c) and with respect to ρr+1, so the last assertion of Nonnegative Combinations of Two Finite Measures, and the Average of Two Couplings §combination shows that f is integrable with respect to γr+1 with ∫fdγr+1=∫fdγr+cr+1∫fdρr+1=∑i=1r+1ci∫fdρi, by H(r)(c) and claim 1 of Properties of Finite Sums.
Step 2 (Conclusion of clause 1). By Step 1, H(M) holds, and γM is the function ∑i=1Mciρi of the statement. It is a measure by H(M)(a), and γM(Rn)=∑i=1Mciρi(Rn)=∑i=1Mci=1, so it is a probability measure. The two assertions about Borel f:Rn→[0,∞] are H(M)(b), and the assertion about real-valued f is H(M)(c).
Step 3 (Clause 2: the weights). The set F is finite, so by Pairs of Euclidean Points: Coordinate Projections, Pairings, the Product Measure on a Euclidean Space, Borel Norm Functions and Finite Sets §finite-sets the sets F, Rn∖F and {zi} for i∈[M] belong to B(Rn). The points zi being pairwise distinct, the sets {z1},…,{zM} are pairwise disjoint with union F, so claim 1 of Basic Properties of a Measure gives ρ(F)=∑i=1Mρ({zi}), a sum of real numbers. By claim 3 of the same lemma (differences, ρ being finite), ρ(F)=ρ(Rn∖(Rn∖F))=ρ(Rn)−ρ(Rn∖F)=1−0=1. Hence ∑i=1Mρ({zi})=1. Put wi=ρ({zi}), a nonnegative real number. Each δzi belongs to P(Rn) by Dirac Measures on Euclidean Space: Probability Measure, Integrals, Push-Forwards and the Coupling of Two Dirac Measures §measure, so by clause 1 the function σ=∑i=1Mwiδzi is a probability measure on Rn.
Step 4 (Clause 2: the representation). Let B∈B(Rn). The sets B∩F and B∖F are Borel, disjoint, with union B, and ρ(B∖F)≤ρ(Rn∖F)=0 by claim 2 of Basic Properties of a Measure, so ρ(B∖F)=0 and, by claim 1 of that lemma, ρ(B)=ρ(B∩F). For i∈[M] let Ci={zi}∩B, a Borel set; the sets C1,…,CM are pairwise disjoint with union B∩F, so claim 1 of Basic Properties of a Measure gives ρ(B∩F)=∑i=1Mρ(Ci). If zi∈B, then Ci={zi} and ρ(Ci)=wi=wiδzi(B); if zi∈/B, then Ci=∅ and ρ(Ci)=0=wiδzi(B), by the definition of δzi in The Dirac Measure at a Point of Euclidean Space §dirac. Therefore
ρ(B)=i=1∑Mwiδzi(B)=σ(B)
for every B∈B(Rn), that is, ρ=∑i=1Mρ({zi})δzi.