# Liouville's Formula - Proof of Liouville's Formula

Proof of Liouville's Formula

We omit the argument x for brevity. By the Leibniz formula for determinants, the derivative of the determinant of Φ = (Φi, j)i, j ∈ {0,...,n} can be calculated by differentiating one row at a time and taking the sum, i.e.

$(detPhi)'=sum_{i=1}^ndetbegin{pmatrix} Phi_{1,1}&Phi_{1,2}&cdots&Phi_{1,n}\ vdots&vdots&&vdots\ Phi'_{i,1}&Phi'_{i,2}&cdots&Phi'_{i,n}\ vdots&vdots&&vdots\ Phi_{n,1}&Phi_{n,2}&cdots&Phi_{n,n} end{pmatrix}.$

(2)

Since the matrix-valued solution Φ satisfies the equation, we have for every entry of the matrix Φ'

or for the entire row

$(Phi'_{i,1},dots,Phi'_{i,n}) =sum_{j=1}^n a_{i,j}(Phi_{j,1},ldots,Phi_{j,n}), qquad iin{1,ldots,n}.$

When we subtract from the i th row the linear combination

of all the other rows, then the value of the determinant remains unchanged, hence

$detbegin{pmatrix} Phi_{1,1}&Phi_{1,2}&cdots&Phi_{1,n}\ vdots&vdots&&vdots\ Phi'_{i,1}&Phi'_{i,2}&cdots&Phi'_{i,n}\ vdots&vdots&&vdots\ Phi_{n,1}&Phi_{n,2}&cdots&Phi_{n,n} end{pmatrix} =detbegin{pmatrix} Phi_{1,1}&Phi_{1,2}&cdots&Phi_{1,n}\ vdots&vdots&&vdots\ a_{i,i}Phi_{i,1}&a_{i,i}Phi_{i,2}&cdots&a_{i,i}Phi_{i,n}\ vdots&vdots&&vdots\ Phi_{n,1}&Phi_{n,2}&cdots&Phi_{n,n} end{pmatrix} =a_{i,i}detPhi$

for every i ∈ {1, . . ., n} by the linearity of the determinant with respect to every row. Hence

(3)

by (2) and the definition of the trace. It remains to show that this representation of the derivative implies Liouville's formula.

Fix x0I. Since the trace of A is assumed to be continuous function on I, it is bounded on every closed and bounded subinterval of I and therefore integrable, hence

is a well defined function. Differentiating both sides, using the product rule, the chain rule, the derivative of the exponential function and the fundamental theorem of calculus, we obtain

due to the derivative in (3). Therefore, g has to be constant on I, because otherwise we would obtain a contradiction to the mean value theorem (applied separately to the real and imaginary part in the complex-valued case). Since g(x0) = det Φ(x0), Liouville's formula follows by solving the definition of g for det Φ(x).