Each result cited below is universally quantified over the data in its own statement.
Conventions. For a natural number l≥1 let Bl and λl be the σ-algebra and measure on Rl built in Finite Products of Lebesgue Measure and Coordinate Integration on Rl. By claim 5 of The Borel Sigma-Algebra of a Euclidean Space as a Product, and Measurability of Projections, Sequentially Continuous Maps, and Open and Closed Sets, Bl=B(Rl), and by Lebesgue Measure on Rn this λl is Lebesgue measure on B(Rl); for l=1 it is Lebesgue measure λ on B1=B(R). By claim 1 of Finite Products of Lebesgue Measure and Coordinate Integration on Rl each λl is σ-finite. For l≥2 we use the identification of Rl with Rl−1×R fixed there, writing x=(x′,xl) with x′=(x1,…,xl−1); under it Bl=Bl−1⊗B(R) and λl=λl−1⊗λ. A Borel f:Rl→[0,∞] is one with {f>a}∈B(Rl) for every real a (Measure Spaces and the Lebesgue Integral: Standing Notation §measurable).
Claim A (Composition). If S:Rp→Rm is measurable with respect to B(Rp) and B(Rm) and h:Rm→[0,∞] is Borel, then h∘S is Borel: for real a, {h∘S>a}=S−1({h>a})∈B(Rp) by Measurable Function and Real-Valued Measurable Function.
Claim B (Arithmetic in [0,∞]). For positive reals a,b and I∈[0,∞] we have a(bI)=(ab)I and a(a−1I)=I: for finite I this is field arithmetic, and for I=∞ both sides are ∞ by the convention a⋅∞=∞ for 0<a of Measure Spaces and the Lebesgue Integral: Standing Notation §extended (with 0<ab and 0<a−1 by claims 5 and 7 of Elementary Order Arithmetic in an Ordered Field).
Step 1 (Proof of Lebesgue Measure under a Triangular or Symmetric Positive Definite Linear Map §borel). Let T be a real d×d matrix. By claim 3 of Linearity, Compatibility with the Matrix Product, and a Norm Bound for the Matrix-Vector Product there is a real C≥0 with ∥Tv∥≤C∥v∥ for all v∈Rd, and Tx−Ty=T(x−y) by claim 1 there, so ∥Tx−Ty∥≤C∥x−y∥. Since dE(x,y)=∥x−y∥ (claim 2 of Elementary Properties of the Euclidean Norm on Rn), x↦Tx is Lipschitz with constant C, hence continuous by A Lipschitz Map is Uniformly Continuous. For i∈[d] the component x↦(Tx)i satisfies ∣(Tx)i−(Ty)i∣≤∥Tx−Ty∥≤C∥x−y∥ by claim 4 of Elementary Properties of the Euclidean Norm on Rn. On R=R1 the Euclidean distance is ∣s−t∣: for u∈R1, ∥u∥2=u2=∣u∣2 by claim 1 of Elementary Properties of the Euclidean Norm on Rn and claim 1 of Properties of the Absolute Value in an Ordered Field, both ∥u∥ and ∣u∣ are nonnegative, so ∥u∥=∣u∣ by claim 3 of Monotonicity of Squaring on the Nonnegative Elements of an Ordered Field. Hence each component is Lipschitz, so continuous (A Lipschitz Map is Uniformly Continuous), from (Rd,dE) to (R,dE), and therefore measurable with respect to Bd and B(R) by claim 3(a) of The Borel Sigma-Algebra of a Euclidean Space as a Product, and Measurability of Projections, Sequentially Continuous Maps, and Open and Closed Sets. By claim 2 there (with (X,F)=(Rd,Bd)), x↦Tx is measurable with respect to Bd and Bd, that is, Borel, as Bd=B(Rd).
Step 2 (Positivity of triangular determinants). Let T be a lower or upper triangular real l×l matrix with 0<Tkk for all k∈[l]. By claim 1 or claim 2 of The Determinant of a Triangular Matrix is the Product of its Diagonal Entries, whose determinant is that of Determinant of a Real Square Matrix by Row Properties of the Determinant, detT=∏k=1lTkk. By the recursion in claim 1 of Properties of Finite Products and claim 5 of Elementary Order Arithmetic in an Ordered Field, induction on the number of factors gives 0<detT. This proves the first assertion of Lebesgue Measure under a Triangular or Symmetric Positive Definite Linear Map §triangular. If l≥2, let T′ be the (l−1)×(l−1) matrix with Tij′=Tij for i,j∈[l−1]; it is triangular of the same kind with the same first l−1 diagonal entries, and by The Determinant of a Triangular Matrix is the Product of its Diagonal Entries and claim 1 of Properties of Finite Products (restriction and recursion),
detT=detT′⋅Tll.
Step 3 (Splitting off the last coordinate). Let l≥2, T and T′ as in Step 2, and x=(x′,xl)∈Rl. By Matrix-Vector Product and the recursion in claim 1 of Properties of Finite Sums, (Tx)i=∑j=1l−1Tijxj+Tilxl for every i∈[l], and ∑j=1l−1Tijxj=(T′x′)i for i≤l−1.
If T is lower triangular, then Til=0 for i<l, so, with β(x′)=∑j=1l−1Tljxj,
Tx=(T′x′,Tllxl+β(x′)).
If T is upper triangular, then Tlj=0 for j<l, so ∑j=1l−1Tljxj=0 (claim 3 of Properties of Finite Sums with λ=0), and, with b=(T1l,…,Tl−1,l)∈Rl−1 and the sum and scalar multiple of Rl−1 taken coordinatewise,
Tx=(T′x′+xlb,Tllxl).
Step 4 (Proof of the integral identity in Lebesgue Measure under a Triangular or Symmetric Positive Definite Linear Map §triangular). We prove by induction on l≥1 the statement P(l): for every lower or upper triangular real l×l matrix T with positive diagonal entries and every Borel f:Rl→[0,∞], detT∫Rlf(Tx)dλl(x)=∫Rlfdλl. Here f∘T is Borel by Claim A and Step 1, so the left integral is defined. Then P(d) is the claim.
Base case l=1. Here T=(t) with t=T11>0, Tx=tx by Matrix-Vector Product and claim 1 of Properties of Finite Sums, and detT=t by Step 2 and claim 1 of Properties of Finite Products. By claim 2 of Scaling of Lebesgue Measure and the Lebesgue Integral on Rn with n=1 and c=t, where ∣t∣1=∣t∣=t by Natural Number Power of an Element of a Field, claim 1 of Properties of Finite Products and Absolute Value in an Ordered Field,
∫Rf(tx)dλ(x)=t−1∫Rfdλ,
and multiplying by t gives P(1) by Claim B.
Inductive step. Let l≥2 and assume P(l−1). Let T and f be as in P(l), let T′ be as in Step 2, and put H=f∘T, Borel on Rl=Rl−1×R. We apply Tonelli and Fubini Theorems with (X,F,μ)=(Rl−1,Bl−1,λl−1) and (Y,G,ν)=(R,B(R),λ), both σ-finite, to H and to f, which are Bl−1⊗B(R)-measurable.
Lower triangular T. By the Tonelli part of Tonelli and Fubini Theorems and Step 3,
∫RlHdλl=∫Rl−1(∫Rf(T′x′,Tllt+β(x′))dλ(t))dλl−1(x′).
Fix x′. The section k(s)=f(T′x′,s) is B(R)-measurable by the sections part of Tonelli and Fubini Theorems. By claim 2 of Translation and Reflection Invariance of Lebesgue Measure on Rn with n=1 and a=β(x′), k1(s)=k(s+β(x′)) is measurable with ∫Rk1dλ=∫Rkdλ, and by claim 2 of Scaling of Lebesgue Measure and the Lebesgue Integral on Rn with n=1 and c=Tll (as in the base case), ∫Rk1(Tllt)dλ(t)=Tll−1∫Rk1dλ. As k1(Tllt)=f(T′x′,Tllt+β(x′)), the inner integral equals Tll−1g(T′x′), where g(z′)=∫Rf(z′,s)dλ(s) defines a Bl−1-measurable g:Rl−1→[0,∞] by the Tonelli part of Tonelli and Fubini Theorems. The map g∘T′ is Borel (Claim A and Step 1), so claim 1 of Linearity and Monotonicity of the Lebesgue Integral with the constant Tll−1 gives ∫RlHdλl=Tll−1∫Rl−1g(T′x′)dλl−1(x′). By P(l−1) applied to T′ and g, and the Tonelli part of Tonelli and Fubini Theorems applied to f,
detT′∫Rl−1g(T′x′)dλl−1(x′)=∫Rl−1gdλl−1=∫Rlfdλl.
Using detT=detT′⋅Tll (Step 2) and Claim B,
detT∫Hdλl=detT′⋅Tll(Tll−1∫g∘T′dλl−1)=detT′∫g∘T′dλl−1=∫fdλl.
Upper triangular T. By the Tonelli part of Tonelli and Fubini Theorems, in the other order, and Step 3,
∫RlHdλl=∫R(∫Rl−1f(T′x′+tb,Tllt)dλl−1(x′))dλ(t).
Fix t. The section k(z′)=f(z′,Tllt) is Bl−1-measurable by the sections part of Tonelli and Fubini Theorems. By claim 2 of Translation and Reflection Invariance of Lebesgue Measure on Rn with n=l−1 and a=tb, k2(z′)=k(z′+tb) is measurable with ∫k2dλl−1=∫kdλl−1. The inner integrand is k2(T′x′), and k2∘T′ is Borel (Claim A, Step 1); by P(l−1) applied to T′ and k2,
detT′∫Rl−1k2(T′x′)dλl−1(x′)=∫Rl−1kdλl−1=g(Tllt),g(s)=∫Rl−1f(z′,s)dλl−1(z′),
where g:R→[0,∞] is B(R)-measurable by the Tonelli part of Tonelli and Fubini Theorems. By Claim B the inner integral equals (detT′)−1g(Tllt). The map t↦g(Tllt) is measurable by claim 2 of Scaling of Lebesgue Measure and the Lebesgue Integral on Rn, so by claim 1 of Linearity and Monotonicity of the Lebesgue Integral, claim 2 of Scaling of Lebesgue Measure and the Lebesgue Integral on Rn (with n=1, c=Tll) and the Tonelli part of Tonelli and Fubini Theorems,
∫RlHdλl=(detT′)−1∫Rg(Tllt)dλ(t)=(detT′)−1(Tll−1∫Rgdλ)=(detT′)−1(Tll−1∫Rlfdλl).
Multiplying by detT=detT′⋅Tll and using Claim B gives detT∫Hdλl=∫fdλl. This proves P(l) and completes the induction.
Step 5 (Proof of Lebesgue Measure under a Triangular or Symmetric Positive Definite Linear Map §spd). Let A∈S(d) be positive definite. By Cholesky Factorisation of a Symmetric Positive Definite Real Matrix there is a lower triangular real d×d matrix L with 0<Lii for all i and A=LL⊤. Put U=L⊤; by Transpose of a Real Matrix, Uij=Lji, which is 0 when j<i, and Uii=Lii>0, so U is upper triangular with positive diagonal. By The Determinant is Multiplicative, Step 2 and claim 5 of Elementary Order Arithmetic in an Ordered Field, detA=detL⋅detU>0.
By claim 2 of Linearity, Compatibility with the Matrix Product, and a Norm Bound for the Matrix-Vector Product, Ax=(LU)x=L(Ux). Let f:Rd→[0,∞] be Borel; f∘L is Borel by Claim A and Step 1. By Lebesgue Measure under a Triangular or Symmetric Positive Definite Linear Map §triangular (proved in Step 4) applied to U and f∘L, and then to L and f,
detU∫Rdf(Ax)dλd(x)=∫Rdf∘Ldλd,detL∫Rdf∘Ldλd=∫Rdfdλd,
so by Claim B, detA∫f(Ax)dλd(x)=detL(detU∫f(Ax)dλd(x))=∫fdλd.
By Invertibility of Symmetric Positive Definite Matrices, A is invertible. For y∈Rd, claim 2 of Linearity, Compatibility with the Matrix Product, and a Norm Bound for the Matrix-Vector Product gives A(A−1y)=(AA−1)y=Idy=y and likewise A−1(Ay)=y, where (Idy)i=∑jδijyj=yi by Identity Matrix, Matrix-Vector Product and claim 7 of Properties of Finite Sums. Hence x↦Ax is a bijection of Rd onto Rd with inverse y↦A−1y.
Let B∈B(Rd). Then y∈A(B) if and only if A−1y∈B, so A(B) is the preimage of B under y↦A−1y, which is Borel by Step 1; thus A(B)∈B(Rd). The indicator f=1A(B) is Borel (each set {f>a} is Rd, A(B) or ∅), and f(Ax)=1B(x) for all x, by injectivity. The identity just proved and The Integral of an Indicator Function is the Measure of the Set give detAλd(B)=λd(A(B)).