Throughout, the values of a finite measure are real numbers ≥0, on which the sum and the product of [0,∞] are those of R (routine facts and background of the setting). Let μ=a1μ1+a2μ2.
Step 1: a finite measure. The values of μ are real numbers ≥0, μ(∅)=0, and μ(X)<∞. For a pairwise disjoint sequence (Am) in F, countable additivity of μ1 and μ2 and the additivity and homogeneity of sums of sequences in [0,∞] give
μ(m∈N⋃Am)=a1m∈N∑μ1(Am)+a2m∈N∑μ2(Am)=m∈N∑(a1μ1(Am)+a2μ2(Am))=m∈N∑μ(Am).
This proves (i).
Step 2: nonnegative simple functions. For s∈S+, the simple integral is a finite sum over the image s[X] whose terms are real, and
Iμ(s)=y∈s[X]∑yμ(s−1({y}))=a1y∈s[X]∑yμ1(s−1({y}))+a2y∈s[X]∑yμ2(s−1({y}))=a1Iμ1(s)+a2Iμ2(s).
Step 3: nonnegative measurable functions. Let g∈M+, and let (sm) be its sequence of dyadic approximations: sm∈S+, sm≤sm+1 pointwise, and g(x)=supmsm(x) for every x∈X. For each of the measures ρ=μ,μ1,μ2, the integral of a nonnegative simple function is its simple integral, and the monotone convergence theorem gives ∫gdρ=supmIρ(sm), the sequence (Iρ(sm))m being nondecreasing. Since the supremum of the sum of two nondecreasing sequences in [0,∞] is the sum of their suprema, and suprema commute with multiplication by ai∈[0,∞], Step 2 gives
∫gdμ=msup(a1Iμ1(sm)+a2Iμ2(sm))=a1msupIμ1(sm)+a2msupIμ2(sm)=a1∫gdμ1+a2∫gdμ2.
Step 4: integrable functions. Let f∈L1(μ1)∩L1(μ2). Then f is measurable, f+,f−∈M+, and their integrals against μ1 and μ2 are real numbers. By Step 3 with g=f+ and g=f−,
∫f±dμ=a1∫f±dμ1+a2∫f±dμ2<∞,
so f∈L1(μ) by Integrable Functions and Their Lebesgue Integral (i), and subtracting the identity for f− from that for f+ gives, by Integrable Functions and Their Lebesgue Integral (iii),
∫fdμ=a1(∫f+dμ1−∫f−dμ1)+a2(∫f+dμ2−∫f−dμ2)=a1∫fdμ1+a2∫fdμ2.
With Step 3 this proves (ii).