Articulo de referencia

ecuación de Monge

En la teoría matemática de las ecuaciones diferenciales parciales , una ecuación de Monge , que recibe su nombre de Gaspard Monge , es un tipo de ecuación diferencial parcial de...

En la teoría matemática de las ecuaciones diferenciales parciales , una ecuación de Monge , que recibe su nombre de Gaspard Monge , es un tipo de ecuación diferencial parcial de primer orden .

Una ecuación de Monge es una función de tipoF(,q1:norte,pag1:norte):R2norte+1R{\displaystyle F(u,q^{1:n},p_{1:n}):\mathbb {R} ^{2n+1}\to \mathbb {R} }. El problema consiste en encontrar soluciones de tipo(q1:norte):RnorteR{\displaystyle u(q^{1:n}):\mathbb {R} ^{n}\to \mathbb {R} }, de tal manera queF(,q,q)=0.{\displaystyle F\left(u,q,\partial _{q}u\right)=0.}En notación moderna, es una ecuación en una variedad diferenciable.METRO{\displaystyle M}definido por una funciónF:R×TMETROR{\displaystyle F:\mathbb {R} \times T^{*}M\to \mathbb {R} }, dóndeTMETRO{\displaystyle T^{*}M}es el fibrado cotangente . El problema consiste en encontrar soluciones de tipo(q):METROR{\displaystyle u(q):M\to \mathbb {R} }, de tal manera queF(,q,d)=0.{\displaystyle F\left(u,q,du\right)=0.}La ecuación de Hamilton-Jacobi es un ejemplo particularmente importante.

Solución

La ecuación de Monge se suele resolver mediante el método de las características . Específicamente mediante el cono de Monge .

caso cuasilineal

Monge estudió primero el caso en el que F es lineal en las derivadas:A0(,q)+i=0norteAi(,q)qi=0{\displaystyle A^{0}(u,q)+\sum _{i=0}^{n}A^{i}(u,q)\partial _{q^{i}}u=0}dóndeA0:norte{\displaystyle A^{0:n}}son funciones de,q{\displaystyle u,q}Este caso se denomina ecuación diferencial parcial cuasilineal de primer orden . CuandoA1,,Anorte{\displaystyle A^{1},\dots ,A^{n}}no depender de{\displaystyle u}, es una EDP semilineal de primer orden .

Consideremos la gráfica de una superficie de solución enR×METRO{\displaystyle \mathbb {R} \times M}La tangente a la gráfica en un punto fijo.,q{\displaystyle u,q}tiene ecuacionesa0d+aidqi=0{\displaystyle a_{0}du+a_{i}dq^{i}=0}donde los coeficientes de valor real satisfacenAμaμ=0{\displaystyle A^{\mu }a_{\mu }=0}Esta es una ecuación lineal y un conjunto de todas las soluciones.{aμ:aμAμ=0}{\displaystyle \{a_{\mu }:a_{\mu }A^{\mu }=0\}}conforman un subespacio de codimensión 1 en el espacio cotangente.T(,q)(R×METRO){\displaystyle T_{(u,q)}^{*}(\mathbb {R} \times M)}.

Dualmente, el aniquilador de{aμ:aμAμ=0}{\displaystyle \{a_{\mu }:a_{\mu }A^{\mu }=0\}}es un subespacio de dimensión 1 en el espacio tangenteT(,q)(R×METRO){\displaystyle T_{(u,q)}(\mathbb {R} \times M)}, llamado eje de Monge . Está generado por el vector característico.A0+Aiqi{\displaystyle A^{0}\partial _{u}+A^{i}\partial _{q^{i}}}. Cualquier gráfica de solución que pase por el punto,q{\displaystyle u,q}necesariamente debe contener el vector característico.

De esta forma, obtenemos el campo vectorial característico de la ecuación enR×METRO{\displaystyle \mathbb {R} \times M}. Integrando el campo, obtenemos una fibración deR×METRO{\displaystyle \mathbb {R} \times M}en unRnorte{\displaystyle \mathbb {R} ^{n}}Familia parametrizada de curvas unidimensionales. Estas son las curvas características de la ecuación. Dado que resolver estas curvas características solo requiere integrar una EDO , se dice que se ha "reducido a cuadraturas" o que "es integrable ".

Cualquier gráfica de solución debe ser fibrada en curvas características. Por el contrario, cualquierRnorte1{\displaystyle \mathbb {R} ^{n-1}}Una subfamilia parametrizada de curvas características conforma un gráfico de solución y puede obtenerse tomando cualquier superficie suave de codimensión 2.SR×METRO{\displaystyle S\subset \mathbb {R} \times M}y tomando la unión de todas las curvas características que pasan por esta superficie. Esta superficie representa la condición de contorno de Cauchy para la ecuación. Esto tiene algunas excepciones:

  • Cáustica : Si no se cumplen las condiciones del teorema de Picard-Lindelöf , distintas curvas características pueden converger o divergir del mismo punto; es decir, las curvas son soluciones singulares . Un gráfico de soluciones que contiene una curva de solución singular se vuelve indeterminado más allá del punto singular. Varias características pueden colapsar en una sola ( onda de choque ), y una característica puede divergir en varias ( rarefacción ).
  • Si la superficieS{\displaystyle S}que codifica la condición de contorno de Cauchy no es transversal al campo vectorial característico, entonces la superficie barrida porS{\displaystyle S}contendría curvas de solución duplicadas y, por lo tanto, perdería al menos una dimensión. La condición de que el límite de Cauchy sea transversal al campo vectorial característico es la condición de rango no característico .
  • Si(A1,,Anorte)=0{\displaystyle (A^{1},\dots ,A^{n})=0}En algún punto de la curva característica, entonces la curva se movería directamente verticalmente, lo que significa que una superficie que contiene dicha curva, incluso si es suave y está bien definida, ya no será la gráfica de una función diferenciable de tipoMETROR{\displaystyle M\to \mathbb {R} }. Esta es una singularidad de tiempo finito . Esto corresponde a la degeneración de la ecuación diferencial original enA0(,q)=0{\displaystyle A^{0}(u,q)=0}.

Banda característica

Propagación de bicaracterísticas en una franja característica. Las líneas negras representan las curvas características de la función, separadas infinitesimalmente. Las líneas rojas representan las curvas de contorno de la función, también separadas infinitesimalmente. Los gradientes rojos representan los covectores que se propagan a lo largo de una curva característica. La bicaracterística es la propagación combinada del covector y la curva característica.

Una franja característica es un engrosamiento infinitesimal de una curva característica, construida como la sección de la superficie de la solución obtenida a partir de una superficie de condición inicial infinitesimal de codimensión 2.δST(,q)(R×METRO){\displaystyle \delta S\subset T_{(u,q)}(\mathbb {R} \times M)}que es transversal al campo vectorial característico en ese punto. De forma equivalente, puede entenderse como la trayectoria en el espacio de frentes de onda infinitesimales (ondículas).R×TMETRO{\displaystyle \mathbb {R} \times T^{*}M}Desde esta perspectiva, el conjunto de ecuacionesA0(,q)+Ai(,q)pagi=0,dpagidqi=0{\displaystyle A^{0}(u,q)+A^{i}(u,q)p_{i}=0,\quad du-p_{i}dq^{i}=0}es una ecuación de onda. Una superficie de solución{\displaystyle u}se descompone en una familia (n-1) dimensional de frentes de onda infinitesimales que se propagan independientemente según la ecuación de onda. [ 1 ]

Caso general

Cuando el caso no es cuasilineal, el problema sigue siendo encontrar franjas características deF(,q,pag)=0,dpagidqi=0{\displaystyle F(u,q,p)=0,\quad du-p_{i}dq^{i}=0}Defina la 1-formaα=dpagidqi{\displaystyle \alpha =du-p_{i}dq^{i}}. La condiciónα=0{\displaystyle \alpha =0}define una estructura de contacto en la variedad de ondículas de (2n+1) dimensiones.R×TMETRO{\displaystyle \mathbb {R} \times T^{*}M}, escrito como(R×TMETRO,α){\displaystyle (\mathbb {R} \times T^{*}M,\alpha )}y llamado colector de 1 chorro , mientras que la condiciónF(,q,pag)=0{\displaystyle F(u,q,p)=0}define una hipersuperficie . El problema consiste en encontrar una subvariedad n -dimensional que sea tangente a la estructura de contacto mientras se permanece dentro de la subvariedad 2n -dimensional.S:=F1(0){\displaystyle S:=F^{-1}(0)}En otras palabras, se trata de encontrar subvariedades legendrianas enS{\displaystyle S}Cualquier función diferenciable:METROR{\displaystyle u:M\to \mathbb {R} }eleva de forma única a una subvariedad legendriana y, a la inversa, cualquier subvariedad legendriana que no sea vertical en ningún punto es la elevación de una única:METROR{\displaystyle u:M\to \mathbb {R} }obtenido al tomar el mapa de proyecciónR×TMETROR×METRO{\displaystyle \mathbb {R} \times T^{*}M\to \mathbb {R} \times M}(es decir, al olvidarpag{\displaystyle p}).

Esto tiene una solución estándar en geometría de contacto . Comencemos con la variedad de contacto de (2n+1) dimensiones.(R×TMETRO,α){\displaystyle (\mathbb {R} \times T^{*}M,\alpha )}y definir la 2-formaω=dα{\displaystyle \omega =d\alpha }y el campo vectorial de ReebR:={\displaystyle R:=\partial _{u}}Las franjas características son las curvas integrales del campo vectorial bicaracterístico .incógnitaF{\displaystyle X_{F}}que satisfacen las ecuaciones del hamiltoniano de contacto : [ 2 ]α(incógnitaF)=F,ω(incógnitaF,)=dF(RF)α{\displaystyle \alpha (X_{F})=F,\quad \omega (X_{F},\cdot )=dF-(RF)\alpha }Las curvas integrales del campo vectorial bicaracterístico son las curvas bicaracterísticas , que son trayectorias trazadas por las franjas características. En coordenadas,incógnitaF=(incógnitaF,,incógnitaF,q,incógnitaF,pag)=(pagpagF,pagF,qF(F)pag){\displaystyle X_{F}=(X_{F,u},X_{F,q},X_{F,p})=(p\cdot \partial _{p}F,\partial _{p}F,-\partial _{q}F-(\partial _{u}F)p)}Cualquier subvariedad legendriana enS{\displaystyle S}se convierte en fibra en unRnorte1{\displaystyle \mathbb {R} ^{n-1}}-familia de curvas bicaracterísticas. Por el contrario, dado cualquierRnorte1{\displaystyle \mathbb {R} ^{n-1}}-subvariedad parametrizada enS{\displaystyle S}que es transversal aincógnitaF{\displaystyle X_{F}}, se extiende a una subvariedad legendriana enS{\displaystyle S}tomando la unión de todas las curvas bicaracterísticas que pasan por ella.

Una condición de contorno de Cauchy es unaRnorte1{\displaystyle \mathbb {R} ^{n-1}}-subvariedad parametrizada enR×METRO{\displaystyle \mathbb {R} \times M}:(s1,,snorte1)(0(s),q01(s),,q0norte(s)){\displaystyle (s^{1},\dots ,s^{n-1})\mapsto (u_{0}(s),q_{0}^{1}(s),\dots ,q_{0}^{n}(s))}La solución se obtiene elevándola a unaRnorte1{\displaystyle \mathbb {R} ^{n-1}}-subvariedad parametrizada enR×METRO{\displaystyle \mathbb {R} \times M}:(s1,,snorte1)(0(s),q01(s),,q0norte(s),pag0,1(s),,pag0,norte1(s)){\displaystyle (s^{1},\dots ,s^{n-1})\mapsto (u_{0}(s),q_{0}^{1}(s),\dots ,q_{0}^{n}(s),p_{0,1}(s),\dots ,p_{0,n-1}(s))}que satisfaceF(0(s),q0(s),pag0(s))=0,sk0i=1nortepagiskq0i=0,k=1,,norte1{\displaystyle F(u_{0}(s),q_{0}(s),p_{0}(s))=0,\quad \partial _{s^{k}}u_{0}-\sum _{i=1}^{n}p_{i}\partial _{s^{k}}q_{0}^{i}=0,\quad k=1,\dots ,n-1}La solución existe y es única, con algunas excepciones:

  • Si la condición de contorno no es transversal aincógnitaF{\displaystyle X_{F}}, entonces se rompe la condición de rango no característico. En este caso, la ecuación no logra determinar cómo se debe levantar la condición de contorno.
  • Si en algún momento(,q,pag)S{\displaystyle (u,q,p)\in S}, la superficie es tangente akerα{\displaystyle \ker \alpha }, entoncesincógnitaF(,q,pag)=0{\displaystyle X_{F}(u,q,p)=0}, lo que indica que la ecuación no logra determinar cómo se propaga el frente de onda infinitesimal más allá de este punto.

cono Monge

Conos de Monge a lo largo de una superficie de solución. El cono en cada punto es tangente al plano tangente que pasa por el punto. La línea tangente es la dirección bicaracterística. Cada cono es barrido por todos los planos tangentes posibles que resuelvenF=0{\displaystyle F=0}en ese momento.

Monge resolvió la ecuación mediante una construcción geométrica llamada cono de Monge . Consideremos la totalidad de todas las superficies de solución que pasan por un punto fijo.(,q)R×METRO{\displaystyle (u,q)\in \mathbb {R} \times M}Su superficie envolvente es un cono con vértice(,q){\displaystyle (u,q)}. Este es el cono de Monge en este punto. Dada una condición de contorno de Cauchy especificada por unRnorte1{\displaystyle \mathbb {R} ^{n-1}}Subvariedad de dimensión enR×METRO{\displaystyle \mathbb {R} \times M}La superficie de la solución es la extensión tangente a todos los conos en todas las ubicaciones. En otras palabras, la superficie de la solución se construye tomando una rama de la envolvente de los conos.

Los conos de Monge a lo largo de una condición de contorno de Cauchy barren las ramas de una solución. Las curvas tangentes son las características bifásicas.

Las intersecciones entre los conos y la superficie de la solución son tangentes a las curvas características y se elevan a curvas bicaracterísticas enR×TMETRO{\displaystyle \mathbb {R} \times T^{*}M}Cuando se cumple la condición de rango no característico, la envoltura de los conos tiene solo un número finito de ramas, y por lo tanto, solo hay un número finito de superficies de solución.

En el caso especial de una ecuación cuasilineal, el cono degenera en una sola línea, el eje de Monge.

En el lenguaje de la geometría de contacto, cada punto(,q)R×METRO{\displaystyle (u,q)\in \mathbb {R} \times M}se eleva a un frente de onda infinitesimal(,q,pag)R×TMETRO{\displaystyle (u,q,p)\in \mathbb {R} \times T^{*}M}. La condiciónF(,q,pag)=0{\displaystyle F(u,q,p)=0}restringe la elección depag{\displaystyle p}a un subconjunto de (n-1) dimensiones. Para cualquier talpag{\displaystyle p}La curva bicaracterística produce una posible trayectoria del frente de onda.s((s),q(s),pag(s)){\displaystyle s\mapsto (u(s),q(s),p(s))}, que se proyecta a un vector tangente(d/ds,dq/ds)R×TqMETRO{\displaystyle (du/ds,dq/ds)\in \mathbb {R} \times T_{q}M}. El conjunto generado por todos esos vectores tangentes es el cono de Monge en este punto.

Teoría de la transformación

Sophus Lie estudió la ecuación de Monge en el marco de la teoría de grupos de Lie y la geometría de contacto . [ 3 ]

Una transformación de contacto es una transformación del espacio de 1-jet.R×TMETRO{\displaystyle \mathbb {R} \times T^{*}M}que preserve la estructura de contacto. Por ejemplo, cualquier difeomorfismo de la subyacente.R×METRO{\displaystyle \mathbb {R} \times M}Se eleva a una transformación de contacto mediante prolongación. Estas transformaciones de contacto se denominan transformaciones puntuales . La transformación de Legendre es una transformación de contacto que no es una transformación puntual.

Dada una transformación de contactogramo{\displaystyle g}enR×TMETRO{\displaystyle \mathbb {R} \times T^{*}M}transforma una ecuación de MongeF(,q,pag)=0{\displaystyle F(u,q,p)=0}aF(gramo(,q,pag))=0{\displaystyle F(g\cdot (u,q,p))=0}. Las transformaciones de contacto mapean elementos de superficie a elementos de superficie (es decir, hiperplanos infinitesimales enR×METRO{\displaystyle \mathbb {R} \times M}), superficies de solución a superficies de solución, conos de Monge a conos de Monge y bicaracterísticas a bicaracterísticas.

El grupo de simetría de una ecuación de Monge es el grupo de transformaciones de contacto tales queF(,q,pag)=0{\displaystyle F(u,q,p)=0}yF(gramo(,q,pag))=0{\displaystyle F(g\cdot (u,q,p))=0}son la misma superficie.

Lie demostró que en el caso deMETRO=R2{\displaystyle M=\mathbb {R} ^{2}}:

  • Si el grupo de simetría de la ecuación contiene un grupo de Lie abeliano con 3 dimensiones, entonces se puede transformar en una ecuación de la formaF(q1,q2)=0{\displaystyle F(\partial _{q^{1}}u,\partial _{q^{2}}u)=0}.
  • Si el grupo de simetría de la ecuación contiene un grupo de Lie abeliano con 2 dimensiones, entonces se puede transformar en una ecuación de la formaF(,q1,q2)=0{\displaystyle F(u,\partial _{q^{1}}u,\partial _{q^{2}}u)=0}.
  • Si el grupo de simetría de la ecuación contiene un grupo de Lie abeliano de 1 dimensión, entonces se puede transformar en una ecuación de la formaF(q1,q2,q1,q2)=0{\displaystyle F(q^{1},q^{2},\partial _{q^{1}}u,\partial _{q^{2}}u)=0}.

Para ver cómo se obtienen estos resultados, consideremos el caso del grupo de Lie abeliano tridimensional.

Sea el grupo generado por los 3 campos vectoriales.incógnita0,incógnita1,incógnita2{\displaystyle X_{0},X_{1},X_{2}}Dado que el grupo es abeliano, los campos vectoriales tienen corchetes de Lie nulos.[incógnitai,incógnitaj]=0{\displaystyle [X_{i},X_{j}]=0}, por lo tanto, son integrables. Es decir, existe un sistema de coordenadas.(Q1,Q2,U,PAG1,PAG2){\displaystyle (Q^{1},Q^{2},U,P_{1},P_{2})}en el cualincógnita0=U,incógnita1=Q1,incógnita2=Q2{\displaystyle X_{0}=\partial _{U},X_{1}=\partial _{Q^{1}},X_{2}=\partial _{Q^{2}}}Según el teorema de Darboux , las otras dos variablesPAG1,PAG2{\displaystyle P_{1},P_{2}}se puede elegir de modo que el formulario de contacto seaθ=dUPAGidQi{\displaystyle \theta =dU-P_{i}dQ^{i}}Por lo tanto, la transformación es una transformación de contacto.

DesdeF=0{\displaystyle F=0}es invariante bajoincógnita0,incógnita1,incógnita2{\displaystyle X_{0},X_{1},X_{2}}, en las nuevas coordenadas debe tener la formaF(PAG1,PAG2)=0{\displaystyle F(P_{1},P_{2})=0}.

En el caso general, si el grupo de simetría contiene un grupo de Lie abeliano connorte{\displaystyle n}dimensiones, entonces se puede transformar en una ecuación de la formaF(q1,,qnorte)=0{\displaystyle F(\partial _{q^{1}}u,\dots ,\partial _{q^{n}}u)=0}. CuandoF=0{\displaystyle F=0}La ecuación de Hamilton-Jacobi corresponde a las coordenadas de acción-ángulo de un sistema totalmente integrable . El sistema, visto como el movimiento de una partícula, se mueve indefinidamente con momento constante, ya que las bicaracterísticas son líneas rectas. El sistema, visto como el movimiento de una onda, tiene soluciones de onda plana que nunca se dispersan , puesto que los planos planos son superficies de solución.

También son posibles los grupos de Mentiras No Abelianas.

Reducción de simetría

Cuando la dinámica de un sistema no se ve alterada por la acción de un grupo , esa simetría puede eliminarse calculando el cociente de la acción del grupo. Esto se conoce como reducción de simetría.

Si el grupo de simetría de la ecuación contiene un grupo de Lie abeliano con 3 dimensiones, se puede transformar a la forma deF(q1,q2)=0{\displaystyle F(\partial _{q^{1}}u,\partial _{q^{2}}u)=0}En este caso, los conos de Monge son los mismos en toda la(,q1,q2){\displaystyle (u,q^{1},q^{2})}espacio. Por lo tanto, el movimiento de una partícula a lo largo de una curva característica es el mismo en todas partes y puede reducirse fácilmente. En la dinámica reducida, la partícula permanece en el mismo lugar, y la única dinámica no trivial es su momento.(q1,q2){\displaystyle (\partial _{q^{1}}u,\partial _{q^{2}}u)}, que permanece constante en el tiempo. La dinámica se reduce a una posición estática en un círculo.

Si el grupo de simetría de la ecuación contiene un grupo de Lie abeliano con 2 dimensiones, entonces se puede transformar en una ecuación de la formaF(,q1,q2)=0{\displaystyle F(u,\partial _{q^{1}}u,\partial _{q^{2}}u)=0}En este caso, los conos de Monge son los mismos en toda la(q1,q2){\displaystyle (q^{1},q^{2})}plano, y varía sobre{\displaystyle u}Por lo tanto, la única dinámica no trivial de una partícula a lo largo de una curva característica es la interacción entre su altura{\displaystyle u}y el impulso(q1,q2){\displaystyle (\partial _{q^{1}}u,\partial _{q^{2}}u)}La dinámica se reduce a un flujo sobre un cilindro.

Ejemplos

Casos de baja dimensionalidad

Cuandonorte=1{\displaystyle n=1}, la ecuaciónF(,q1,q1)=0{\displaystyle F(u,q^{1},\partial _{q^{1}}u)=0}se reduce a una EDO para una función de tipo:RR{\displaystyle u:\mathbb {R} \to \mathbb {R} }y en cada punto(,q1){\displaystyle (u,q^{1})}, el cono de Monge consiste en líneas rectas que pasan por el punto con pendientesk{\displaystyle k}satisfactorioF(,q1,k)=0{\displaystyle F(u,q^{1},k)=0}Una curva de solución permanece tangente a cada cono de Monge.

Cuandonorte=2{\displaystyle n=2}, en cada punto(,q1,q2){\displaystyle (u,q^{1},q^{2})}, la ecuaciónF(,q1,q2,q1,q2)=0{\displaystyle F(u,q^{1},q^{2},\partial _{q^{1}}u,\partial _{q^{2}}u)=0}produce una variedad unidimensional de posibles elecciones de(q1,q2){\displaystyle (\partial _{q^{1}}u,\partial _{q^{2}}u)}En general, la variedad se asemeja a la unión de varias curvas diferenciables cerradas. Cada curva diferenciable genera una familia de planos, que a su vez generan un cono como envoltura. La unión de todos estos conos es el cono de Monge.

ecuación de Eikonal

La ecuación totalmente no lineal más simple es la ecuación eikonal . Esta tiene la forma

||2=1,{\displaystyle |\nabla u|^{2}=1,}

de modo que la función F viene dada por

F(incógnita,y,,incógnita,y)=incógnita2+y21.{\displaystyle F(x,y,u,u_{x},u_{y})=u_{x}^{2}+u_{y}^{2}-1.}

El cono dual consta de 1-formasadincógnita+bdy+dodz{\displaystyle adx+bdy+cdz}satisfactorio

a2+b2do2=0.{\displaystyle a^{2}+b^{2}-c^{2}=0.}

Esto define un cono circular. El dual convexo de un círculo es un círculo, por lo que el cono de Monge en cada punto es un cono circular.

Esto se generaliza a la ecuación eikonal para geodésicas en una variedad riemanniana y, más generalmente, en una variedad de Finsler .

Véase también

Referencias

  1. Gray, Jeremy (2021), "Monge y soluciones a ecuaciones diferenciales parciales" , Cambio y variaciones , Springer Undergraduate Mathematics Series, Cham: Springer International Publishing, pp. 95–110 , doi : 10.1007/978-3-030-70575-6_8 , ISBN  978-3-030-70574-9
  2. Bravetti, Alessandro (11 de octubre de 2017). "Dinámica hamiltoniana de contacto: el concepto y su uso" . Entropía . 19 (10): 535. Bibcode : 2017Entrp..19..535B . doi : 10.3390/e19100535 . ISSN 1099-4300 . 
  3. Hawkins, Thomas (2000), "Los orígenes geométricos de la teoría de Lie" , Surgimiento de la teoría de los grupos de Lie , Nueva York, NY: Springer New York, pp. 1–42 , doi : 10.1007/978-1-4612-1202-7_1 , ISBN  978-1-4612-7042-3
  • Rozov, N. Kh. (1990). "Ecuación de Monge". En Hazewinkel, Michiel (ed.). Enciclopedia de Matemáticas . vol.  6. Editores académicos de Kluwer . ISBN 1-55608-005-0.
  • David Hilbert y Richard Courant (1989). Métodos de física matemática, Volumen 2. Wiley Interscience.
  • Ivanov, AB (2001) [1994], "Monge cone" , Encyclopedia of Mathematics , EMS Press
  • Monge, Gaspard (1787). "Mémoire sur le Calcul Intégral des Équations aux Différentialles Partielles". Histoire de l'Académie royale des sciences avec les mémoires de mathématique et de physique tirés des registres de cette Académie (en francés). 1784 : 118-192 .
  • Monge, Gaspard (1850). Aplicaciones de l'Analyse à la Géométrie (en francés) (5ª  ed.). París: Bachiller.