# Fórmula de Liouville

En matemáticas, la fórmula de Liouville, también conocida como fórmula de Abel-Liouville, es una identidad que expresa el determinante de una matriz cuadrada que soluciona un sistema de ecuaciones diferenciales lineales homogéneas de primer orden en función de la suma de los coeficientes de la diagonal del sistema. La fórmula debe su nombre al matemático francés Joseph Liouville.

## Enunciado

Sea $\Phi (t)$ una matriz cuadrada de dimensión $n$ x $n$ que verifica la siguiente ecuación diferencial homogénea de primer orden:
    $\Phi '(t)=A(t)\Phi (t),\qquad t\in I$,
donde $I$ es un intervalo de la recta real y $A(t)$ es una matriz cuadrada de dimensión $n$ x $n$ con coeficientes reales o complejos.
Entonces, si la traza de $A(t)$ es integrable en $I$, se cumple la siguiente relación con el determinante de $\Phi (t)$:
    $\det \Phi (t)=\det \Phi (t_{0})\,\exp \left(\int _{t_{0}}^{t}trA(s)ds\right),\qquad t,t_{0}\in I.$

## Demostración

Denotamos $a_{ij}$ y $\phi _{ij}$ como los elementos individuales de las matrices $A(t)$ y $\Phi (t)$ respectivamente. Por brevedad se omite la variable $t$ en estas matrices y sus coeficientes.
Por la fórmula de Leibniz para el cálculo de determinantes se cumple que
| $(\det \Phi )'=\sum _{i=1}^{n}\det {\begin{pmatrix}\phi _{11}&\phi _{12}&\cdots &\phi _{1n}\\\vdots &\vdots &&\vdots \\\phi '_{i1}&\phi '_{i2}&\cdots &\phi '_{in}\\\vdots &\vdots &&\vdots \\\phi _{n1}&\phi _{n2}&\cdots &\phi _{nn}\end{pmatrix}}=\det {\begin{pmatrix}\phi '_{11}&\phi '_{12}&\cdots &\phi '_{1n}\\\vdots &\vdots &&\vdots \\\phi _{i1}&\phi _{i2}&\cdots &\phi _{in}\\\vdots &\vdots &&\vdots \\\phi _{n1}&\phi _{n2}&\cdots &\phi _{nn}\end{pmatrix}}+\cdots +\det {\begin{pmatrix}\phi _{11}&\phi _{12}&\cdots &\phi _{1n}\\\vdots &\vdots &&\vdots \\\phi '_{i1}&\phi '_{i2}&\cdots &\phi '_{in}\\\vdots &\vdots &&\vdots \\\phi _{n1}&\phi _{n2}&\cdots &\phi _{nn}\end{pmatrix}}+\cdots +\det {\begin{pmatrix}\phi _{11}&\phi _{12}&\cdots &\phi _{1n}\\\vdots &\vdots &&\vdots \\\phi _{i1}&\phi _{i2}&\cdots &\phi _{in}\\\vdots &\vdots &&\vdots \\\phi '_{n1}&\phi '_{n2}&\cdots &\phi '_{nn}\end{pmatrix}}.$ | (1) |
En el $i$-ésimo sumando se aplica una combinación lineal sobre su $i$-ésima fila del resto de sus filas, lo que no altera su valor. Usando la ecuación diferencial de la hipótesis, que en términos de $a_{ij}$ y $\phi _{ij}$ se escribe
    $\phi '_{ik}=\sum _{j=1}^{n}a_{ij}\phi _{jk}\,,\qquad i,k\in \{1,\ldots ,n\},$
se obtiene la expresión del $i$-ésimo término de la suma anterior en función del determinante de $\Phi$:
$\det {\begin{pmatrix}\phi _{11}&\phi _{12}&\cdots &\phi _{1n}\\\vdots &\vdots &&\vdots \\\phi '_{i1}&\phi '_{i2}&\cdots &\phi '_{in}\\\vdots &\vdots &&\vdots \\\phi _{n1}&\phi _{n2}&\cdots &\phi _{nn}\end{pmatrix}}=\det {\begin{pmatrix}\phi _{11}&\phi _{12}&\cdots &\phi _{1n}\\\vdots &\vdots &&\vdots \\\phi '_{i1}-\sum _{\scriptstyle j=1 \atop \scriptstyle j\neq i}^{n}a_{ij}\phi _{j1}&\phi '_{i2}-\sum _{\scriptstyle j=1 \atop \scriptstyle j\neq i}^{n}a_{ij}\phi _{j2}&\cdots &\phi '_{in}-\sum _{\scriptstyle j=1 \atop \scriptstyle j\neq i}^{n}a_{ij}\phi _{jn}\\\vdots &\vdots &&\vdots \\\phi _{n1}&\phi _{n2}&\cdots &\phi _{nn}\end{pmatrix}}=\det {\begin{pmatrix}\phi _{11}&\phi _{12}&\cdots &\phi _{1n}\\\vdots &\vdots &&\vdots \\a_{ii}\phi _{i1}&a_{ii}\phi _{i2}&\cdots &a_{ii}\phi _{in}\\\vdots &\vdots &&\vdots \\\phi _{n1}&\phi _{n2}&\cdots &\phi _{nn}\end{pmatrix}}=a_{ii}\det \Phi .$ Usando esto en la fórmula (1), se obtiene la siguiente ecuación diferencial para el determinante de $\Phi$:
    $(\det \Phi )'=\sum _{i=1}^{n}a_{ii}\det \Phi =\mathrm {tr} \,A\,\det \Phi$.
Se trata de una ecuación diferencial ordinaria de primer orden separable, cuya solución es
    $\det \Phi (t)=\det \Phi (t_{0})\,\exp \left(\int _{t_{0}}^{t}trA(s)ds\right),\qquad t,t_{0}\in I.$