Calculo Variacional PDF

Cálculo Variacional.
Condiciones necesarias y suficientes de extremo. La

ecuación de Euler. Condiciones de Jacobi, Legendre y
Weiertrass.
01/01/2010
Universidad de Almería.
Carmen Gádor Garzón Escamilla y Melina Gorini.
Contenido.
1. Introducción…………………………………………………………………………………………4
1.1. Introducción histórica…………………………………………………………….…..4
1.2. Definición de Cálculo Variacional………………………………………………...5
1.3. Problemas importantes………………………………………………………….…...7
1.3.1. Problema isoperimétrico………………………………………….………..7
1.3.2. Principio de acción mínima……………………………………….……….7
1.3.3. Problema de la Braquistocrona………………………………………….8
1.3.4. El problema de las Geodésicas……………………………………….….10
1.3.5. La Catenaria……………………………………………………………….……11
1.3.6. Conclusión…………………………………………………………...……….…11
1.4. Métodos de resolución de los problemas variacionales…………………..12
1.4.1. Métodos indirectos……………………………………………….………….12
1.4.2. Métodos directos……………………………………………….…………….12
2. Condiciones necesarias de extremo. Ecuación de Euler…………………………15
2.1. Funciones que dependen de funciones de varias
Variables independientes…………………….……………..…………………….28
3. Condiciones suficientes de Extremo…………………………………………..………..33
4. Condición de Legendre…………………………………………………………..…………..39
4.1. Condición necesaria de Legendre para la realización de un
mínimo de un funcional del tipo variacional………..………..……………39
4.2. Condición suficiente de Legendre de inclusión de un
extremal de un funcional en un campo de extremales…….……….....40
5. Condición de Jacobi……………………………………………………………………………42
5.1. Condición necesaria de Jacobi…………………….……………….……………42
5.2. Condición suficiente de Jacobi bajo las cuales un funcional
tiene un mínimo………………………………………………………………………49
5.3. Relación entre la condición de Jacobi y la teoría de
formas cuadráticas………………….………………………………………………..51
5.4. Condición suficiente de Jacobi de inclusión de una extremal
en un campo de extremales centrales………………………………………..56
6.1. Función (, , , ′). Condición de Weiertrass………….……….………..59

6. Condición de Weiertrass……………………………………………………………………..59
7. Anexo: Métodos directos en el cálculo variacional. …………………….………..72

7.1. Método de Euler de diferencias finitas……………………..………..………72
7.2. Método de Ritz……………………………………………………………..………….74
7.3. Método de Kantoróvich………………………………………………….…………79
7.4. Métodos variacionales de búsqueda de valores propios y
funciones propias…………………………………………………………….…….….81
7.4.1. El problema de Sturm-Liouville………………..………………...…..81
7.4.2. Principio de Rayleigh……………………………...……………….….….85
8. Biografías…………………………………………………………………………………………..88
8.1. Leonhard Euler……………….………………………………………………………88
8.2. Karl Gustav Jacobi…………………….………………………………………….….90
Universidad de Almería. | Cálculo Variacional. 2

8.3. Adrien Marie Legendre…………………………………………………….…….…91
8.4. Karl Weiertrass……………………………….………………………….………….…91
9. Bibliografía……………………………………………………………………………………..…93

1. Introducción
1.1. Introducción histórica

El cálculo de variaciones o cálculo variacional es una rama clásica y
fundamental de las matemáticas. No es una exageración afirmar que el
desarrollo de esta rama de las matemáticas ha ido a la par con el desarrollo de
los conceptos centrales del análisis matemático y sus aplicaciones. En lo que
respecta a las aplicaciones, muchos de los conceptos centrales de la física teórica
están en estrecha relación con el cálculo variacional.
Las raíces del cálculo variacional se extienden a tiempos anteriores a la

Grecia clásica. Uno de los problemas más antiguos del cálculo variacional, y de
las matemáticas en general, es el problema isoperimétrico. Este problema
está relacionado con la legendaria Dido fundadora de la ciudad fenicia de
Cartago (buena parte de la leyenda de Dido se encuentra en la Eneida de
Virgilio, aunque por otras fuentes se sabe que fue un personaje histórico).
Cuenta la leyenda que Dido y un grupo de seguidores llegaron a las costas de lo
que ahora es Túnez y solicitaron un pedazo de tierra a los habitantes locales.
Dido pidió la tierra que pueda ser encerrada por la piel de un toro. Desde luego
la petición no parecía muy ambiciosa así que le fue esto concedido. Dido corto la
piel en tiras muy delgadas formando así una cuerda muy larga. Utilizó entonces
esta cuerda para rodear una extensión de tierra en la costa que pasó a
convertirse en la ciudad de Cartago. Independientemente de la veracidad de la
leyenda no es difícil aceptar que el problema de abarcar la mayor área posible
dada una cuerda de longitud fija apareció hace mucho tiempo en la historia. El
filósofo Zenodoros (200 a.n.e.) planteó de manera precisa éste y otros
problemas matemáticos relacionados con encontrar figuras “óptimas", que hoy
podemos considerar problemas clásicos del cálculo variacional. Hay otros
problemas clásicos que son parte del cálculo de variaciones que fueron
planteados y estudiados por Aristóteles y Pappus.
Hasta aquí hemos hablado de problemas de cálculo variacional, pero no

hemos definido esta rama de las matemáticas. De hecho no lo haremos ahora
sino que postergaremos la definición del cálculo de variaciones. Esto no nos
impide notar que en los problemas de cálculo de variaciones siempre se requiere
encontrar curvas, figuras, procesos, “óptimos".
Se le atribuye a Pierre de Fermat, matemático francés del siglo XVII, el

principio físico de tiempo mínimo, el cual establece que la trayectoria que toma
la luz entre dos puntos es la trayectoria que puede ser recorrida en el menor
tiempo. Este principio está relacionado con el principio de distancia
mínima de Herón de Alejandría, filósofo griego del siglo I (la luz sigue la
trayectoria entre dos puntos que resulta ser la más corta). En 1662, Fermat
utilizó su principio de tiempo mínimo para deducir la ya entonces conocida ley
de Snell que describe la refracción de la luz al pasar de un medio a otro. Es a
partir de este momento que se empiezan a utilizar métodos analíticos para la
resolución de problemas de “optimización" (anteriormente estos problemas se
habían abordado por métodos puramente geométricos).

El tratamiento de Fermat de este problema es considerado por varios
historiadores del cálculo variacional como el comienzo del mismo, precisamente
por el uso de técnicas analíticas similares a las que se usarían más tarde en el
análisis matemático. Es interesante mencionar que éstas técnicas jugaron un
papel importante en el desarrollo del cálculo unos años después.
En 1696 Johann Bernoulli publica un desafío para los matemáticos de su

tiempo: el llamado problema de la Braquistocrona. Varios matemáticos
dieron respuesta al desafío. Entre las respuestas destacadas se encuentran la de
su hermano Jakob, la de Newton (publicada de manera anónima) y la de
Leibniz.
Jakob Bernoulli utiliza un método similar al de Fermat, pero más refinado,
para dar respuesta al problema. El desarrollo y generalización de estos métodos
por Euler, y después Lagrange, llevan a un método sistemático para estudiar
este tipo de problemas y éste al cálculo de variaciones. Fue precisamente Euler
quién acuño el término.
Johann Bernoulli también estudio geodésicas en varias superficies. Este es

otro problema clásico del cálculo variacional. La geodésica es la curva más corta
sobre cierta superficie que une a dos puntos de esa superficie.
Durante el siglo XIX los trabajos de Euler y Lagrange son formalizados y

generalizados para conformar lo que es el cálculo de variaciones hoy en día. Es
de destacarse las contribuciones de Weierstrass en la formalización de la teoría.
El desarrollo del cálculo variacional está relacionado con el desarrollo de la

física. Esto es así por el marco conceptual en el que se han desarrollado las ideas
sobre el “comportamiento" de la realidad. Con innegable influencia religiosa el
pensamiento físico ha considerado que los procesos naturales se desarrollan de
manera “óptima". Durante la evolución de los procesos algo se minimiza o
maximiza (Dios o la naturaleza deben ser perfectos). Así las leyes de la física
deben ser el producto de principios variacionales. Es así como surgen la
mecánica analítica y la mecánica hamiltoniana y de ahí la formalización de la
mecánica cuántica. Es también notable que la Teoría General de la Relatividad
también esté relacionada con el cálculo de variaciones.
1.2. Definición de Cálculo Variacional

En una serie de problemas de la física y de la matemática nos encontramos
con funciones definidas sobre un conjunto cuyos elementos también son
funciones de una o varias variables.
Las funciones definidas sobre un conjunto cuyos elementos son funciones, se

llaman funcionales.
Por ejemplo, la longitud l del arco de una curva plana que une dos puntos
dados A( x0 , y 0 ) y B( x1 , y1 ) , es un funcional. La magnitud l puede calcularse si se
da la ecuación de la curva y = y(x ) . Entonces

x1
l [ y ( x )] = ∫ 1 + ( y ') dx
2
x0
El área S de cierta superficie es también un funcional, puesto que se

determina escogiendo la superficie, es decir, escogiendo la función z = z ( x, y ) de
la superficie. Como es sabido,
2 2
 dz   dz 
S [z ( x, y )] = ∫ ∫ 1 +   +   dxdy
D
 dz   dy 
donde D es la proyección de la superficie en el plano OXY .
Los momentos de inercia, los momentos estáticos, las coordenadas del

centro de gravedad de cierta curva o superficie homogénea, son también
funcionales, puesto que sus valores se determinan eligiendo la curva o la
superficie, es decir, las funciones contenidas en la ecuación de dicha curva o
superficie.
En todos estos ejemplos se tiene una dependencia que es característica para

los funcionales: a una función (escalar o vectorial) le corresponde un número,
mientras que al dar una función z = f ( x ) a un número le correspondería otro
número.
En el Cálculo Variacional se consideran los métodos para hallar el valor

máximo o el mínimo de un funcional. Los problemas en que se exige investigar
el máximo o el mínimo de un funcional, se denominan problemas
variacionales.
Durante más de dos siglos, el Cálculo de Variaciones ha sido una de las

principales ramas del Análisis. Es un instrumento de gran utilidad que se puede
aplicar en muy diversos problemas, como ya he dicho antes, tanto en
Matemáticas, como en Física.
Es fácil captar el interés del tema si se toman en cuenta algunos de sus

problemas típicos. A continuación, introducimos algunos ejemplos clásicos del
Cálculo de Variaciones, en los que se muestran los elementos fundamentales del
problema de optimización. Éstos son:
1. Un espacio de funciones V , tal que u : Ω → ℜ q , donde Ω es un

abierto, normalmente acotado, de ℜ n , de frontera, Γ , regular.
2. Restricciones sobre el conjunto de soluciones, que pueden imponerse

bien sobre la frontera Γ , bien sobre el dominio Ω . Por ejemplo u = 0
en Γ , u ≥ Ψ en Ω , etc. El conjunto de funciones que satisfacen estas
restricciones es, en general, un subconjunto, U de V .
3. Un funcional J : V → ℜ de la forma siguiente:

J (u ) := ∫ L( x, u ( x ), u ' ( x ))dx (1)
Ω
La hipótesis sobre V y L deben asegurar la existencia de J sobre V , o al

menos sobre U .
El problema de optimización consiste en hallar el mínimo, u ∈ U , del

funcional J .
1.3. Problemas importantes
1.3.1. Problema isoperimétrico
De entre todas las curvas de longitud λ dada, que unen el punto (0,0 ) con un
punto variable (ξ ,0 ) , encontrar aquella que, junto con el eje OX , encierra una
superficie máxima. El problema es, pues, el de hallar una función, u y un
número, ξ tales que u (0 ) = 0 , u (ξ ) = 0 , u ≥ 0 y que minimicen el funcional
ξ
J (u , ξ ) := − ∫ u
0
y satisfagan la restricción
ξ
∫
2
1 + u' = λ
0
1.3.2. Principio de acción mínima
En Mecánica Clásica, cuando una partícula se mueve bajo la acción de un

potencial V ( x ) , el movimiento real es el dado por las ecuaciones de Newton, que
expresan la aceleración de la partícula en términos de las fuerzas. Cuando las
fuerzas derivan de un potencial V ( x ) , el movimiento real t → x (t ) satisface la
ecuación diferencial:
d 2 x(t ) dV (x (t ))
m 2
=−
dt dx
cuya solución determina el movimiento real que sigue una partícula que en un
instante inicial t1 sale del punto x1 , se mueve bajo la acción del potencial, y
llega en un instante final t 2 al punto x 2 .
Una pregunta interesante es: ¿Podemos singularizar el movimiento real

dado por las soluciones de esta ecuación, entre todos los movimientos que la

partícula podría seguir, para ir desde el punto inicial x1 en el instante t1 al
punto final x 2 en el instante t 2 ?
La respuesta a esta pregunta es un principio básico en Física, que en

Mecánica se denomina principio de Hamilton, o principio de mínima acción.
Este principio caracteriza a los movimientos reales entre todos los movimientos
imaginables que llevarían a la partícula del estado inicial (posición x1 en el
instante t1 ) al estado final (posición x 2 en el instante t 2 ), ambos dados.
La caracterización dada por el principio de Hamilton asocia una cantidad,

denominada acción a cada movimiento imaginable. La acción es una cantidad
de naturaleza bastante diferente a las cantidades que usualmente describen el
estado de la partícula, como posición y/o velocidad. A diferencia de ellas, la
acción no se asocia al estado, sino a la historia completa de la partícula entre
dos instantes inicial y final. Para cada movimiento imaginable, descrito por
t → x (t ) con las condiciones x(t1 ) = x1 , x (t 2 ) = x 2 , la acción de ese movimiento se
define como:
t2  1  dx (t )  2 
S [x (t )] = ∫  m  − V ( x (t ))dt
t1
 2  dt  
El principio de la mínima acción dice: entre todos los movimientos

imaginables, la propiedad que distingue al movimiento real es que el valor de
la acción S [x (t )] es menor para el movimiento real que para cualquier otro.
¿Cuál es la relación entre este principio y la forma newtoniana de planear las

ecuaciones del movimiento? Resulta que ambas maneras de describir el
movimiento son equivalentes. Para verlo, necesitamos abordar el problema de la
búsqueda de la función x (t ) con las condiciones requeridas, que minimice el
valor de la acción. No se trata de un problema ordinario de mínimo, ya que la
acción depende del movimiento como un todo, esto es, depende de la función
x (t ) .
1.3.3. Problema de la Braquistocrona
El problema de la braquistocrona, o curva de descenso más rápido, es uno

de los problemas más antiguos del cálculo de variaciones. La primera solución
fue dada por Johann Bernoulli en 1696, aunque también dieron soluciones
algunos contemporáneos como Jacob Bernoulli, Leibniz y Newton.
Entre todas las curvas que unen los puntos A y B , se desea hallar aquella a
lo largo de la cual un punto material, moviéndose bajo la fuerza de la gravedad
desde A llega al punto B en el menor tiempo.
Para resolver este problema debemos considerar todas las posibles curvas
que unen A y B . A una determinada curva, γ , le corresponderá un valor
determinado, T , del tiempo invertido para el descenso del punto material a lo
largo de ella.. El tiempo, T , dependerá de la elección de γ . De todas las curvas
que unen A con B debemos hallar aquella a la que corresponda el menor valor
de T . El problema puede plantearse de la siguiente forma.
Tracemos un plano vertical que pase por los puntos A y B . La curva de más
rápido descenso debe evidentemente estar en él, así que podemos restringirnos
a curvas sobre dicho plano. Tomemos el punto A como el origen de
coordenadas, el eje OX apuntando en la dirección de la gravedad y sea
B = ( x1 , y1 ) , con x1 > 0 y y1 ≥ 0 . Consideremos una curva arbitraria descrita por
la ecuación
y = y(x ) 0 ≤ x ≤ x1 (2)
donde y es una función regular. Como la curva pasa por A y B , la función y

debe verificar
0 = y (0) , y (x1 ) = y1 (3)
El movimiento de la masa puntual puede describirse por medio de la ley de

la conservación de la energía, E c + E p = cte. , del siguiente modo: en el punto A ,
en el que asumimos que la velocidad inicial es nula, se tiene
E c + E p = E p = mgh A = E
donde E > 0 es una constante y h A es la altura a la que se encuentra el punto A .

En cualquier punto por debajo será
1 2
mv + mgh = E
2
luego
v 2 = 2 g (h A − h )
y tomando la coordenada vertical como x = h A − h , deducimos que la velocidad

del movimiento del punto material es
ds
v≡ = 2 gx
dt
siendo s una parametrización de la trayectoria del punto material. Deducimos

que
ds
dt =
2 gx
y como la longitud de arco de la curva viene dada por
ds = 1 + y ' ( x ) dx
2

tenemos que el tiempo empleado a lo largo de la curva y viene dado por
1
x1  1 + y ' ( x ) 
2 2
J ( y ) = ∫   dx
 (4)
 2 gx
0

Hallar la braquistocrona es equivalente a resolver el siguiente problema de

mínimos: entre todas las posibles funciones (2) que verifican las condiciones
(3), hallar la que corresponda al menor valor de la integral (4).
1.3.4. El problema de las Geodésicas
Las geodésicas son aquellas curvas contenidas en una superficie regular que
minimizan la distancia entre dos puntos de la misma. Enunciaremos este
problema de dos formas:
1. Consideremos una superficie regular S ⊂ ℜ3 definida por la

parametrización:
x = x(u, v ) , y = y (u, v ) , z = z (u, v ) ,
con (u , v ) ∈ [u 0 , u1 ]× [v 0 , v1 ] . Cualquier curva contenida en S puede

parametrizarse en la forma t : [t1 , t 2 ] → (u (t ), v(t )) .
El elemento de arco de las curvas contenidas en S está determinado

por la primera forma fundamental:
ds 2 := Eu ' 2 +2 Fu ' v'+Gv ' 2

con
E := xu2 + y u2 + z u2 , F := xu xv + y u y v + z u z v , G := xv2 + y v2 + z v2
De modo que la longitud del arco entre los puntos correspondientes a

los valores t1 y t 2 es
t2
J (u , v ) = ∫ Eu ' 2 +2 Fu ' v'+Gv ' 2 dt
t1
que es el funcional a minimizar.
2. Si la superficie viene dada de la forma implícita por ϕ ( x, y , z ) = 0 y

representamos una curva sobre ella de forma paramétrica,
(x(t ), y (t ), z (t )) , debemos minimizar el funcional

1
J ( x, y , z ) = ∫
t1
t 0
(x' (t ) 2
+ y ' (t ) + z ' (t ) )2 dt
2 2
Además, las funciones x , y y z deben someterse a la condición

ϕ ( x(t ), y (t ), z (t )) = 0 para t ∈ [t 0 t1 ] . Es lo que se llama un problema
variacional con restricciones de igualdad.
1.3.5. La Catenaria
¿Cómo cuelga un hilo inextensible y flexible, de longitud total L, suspendido

entre dos torres con separación horizontal d , y alturas dadas, A y B ?
Claramente, el principio que determina la forma de equilibrio del hilo es que su
energía potencial sea la menor posible. Cada forma posible del hilo está
descrita por una función x → z ( x ) que debe satisfacer las condiciones z (a ) = A ,
z (b ) = B (donde a y b son las coordenadas horizontales de las torres), d = b − a ,
y además otra condición importante, a saber, la longitud total del hilo debe ser
L ; esta condición se traduce en:
b
L=∫ 1 + z ' ( x ) dx
2
a
Veamos ahora cómo se expresa la energía potencial del hilo cuando su forma
es la función z (x ) . Suponiendo el hilo de densidad lineal ρ constante, la masa
del elemento entre las coordenadas x y x + dx es ρ 1 + z ' ( x ) dx y la energía
2
potencial de ese elemento es z ( x )gρ 1 + z ' (x ) dx . Así pues, la energía potencial

2
total es:
b
E = ρg ∫ z (x ) 1 + z ' (x ) dx
2
a
La forma real será aquella curva que, satisfaciendo la condición adicional de

tener longitud total L , haga mínima la energía potencial. Conviene notar que
este problema es más complicado que los anteriores, ya que interviene en él una
ligadura, o condición auxiliar.
1.3.6. Conclusión
En todos los casos, el problema propuesto se reduce a buscar, entre todas las
funciones f : x → f ( x ) definidas en un intervalo [a, b] , y con condiciones del
tipo f (a ) = A , f (b ) = B (además de otras condiciones de continuidad,
regularidad, etc. Que se precisarán a su tiempo), aquellas que minimizan o
maximizan una expresión del tipo

b
∫ Φ(x, f (x ), f ' (x ))dx
a
En algunos casos, la función f : x → f ( x ) debe satisfacer ciertas condiciones

adicionales, que pueden imaginarse como ligaduras; en todos los casos que
hemos discutido los ligaduras están expresadas también por condiciones del
tipo
b
∫ Ξ(x, f (x ), f ' (x ))dx = cte
a
Por otro lado veremos algunos métodos de resolución de problemas

variacionales.
1.4. Métodos de resolución de los problemas

variacionales
Existen dos aproximaciones fundamentales a la resolución de los problemas

variacionales.
1.4.1. Métodos indirectos
La primera de estas aproximaciones es la heredada de los métodos de

minimización de funciones (dimensión finita) vía el cálculo diferencial. Este
método proporciona condiciones necesarias y condiciones suficientes que dan
lugar a una base metodológica para la resolución de problemas variacionales, la
cual está íntimamente ligada a la teoría de ecuaciones diferenciales.
1.4.2. Métodos directos
La idea fundamental es la extensión del Teorema de Weierstrass a funciones

definidas de dimensión infinita, que tendrá un enunciado del tipo:
Teorema 1. Sea J : V → ℜ un funcional definido en un espacio de funciones

V dotado de cierta noción de convergencia para la que V es compacto y J es
semicontinuo inferiormente. Entonces existe un mínimo de J en V .
A partir de este teorema, se produce del siguiente modo:
1. Se elige la clase de funciones V junto con una noción adecuada de

convergencia para la que V sea completo.

2. Hay que mostrar que J está bien definido en V y que está acotado
inferiormente, de modo que inf u∈V J (u ) sea finito. Esto implica que
se puede construir una sucesión minimizante, u k ∈ V , tal que
J (u k ) → inf u∈V J (u ) .
3. Debemos probar que J es semicontinuo inferiormente

(secuencialmente), es decir, que u k → u implica
J (u ) ≤ lím J (u k )
k→∞
4. Finalmente, debemos demostrar que V es compacto

(secuencialmente) con respecto a la convergencia considerada en 1.
La hipótesis del Teorema de Weierstrass atañen a la función que se desea

minimizar (semicontinuidad inferior) y al conjunto en el cual se busca el
mínimo (compacto). En espacios de dimensión finita estas hipótesis son
relativamente fáciles de comprobar dado que la compacidad de un conjunto es
equivalente a que el mismo sea cerrado y acotado. La continuidad suele
deducirse de un análisis directo de la función a minimizar.
Sin embargo, el Teorema de Riesz establece que la bola unidad cerrada de un

espacio de Banach es compacta si y solo si la dimensión del espacio es finita.
Puesto que este criterio de compacidad falla en el caso de dimensión infinita, se
impone la investigación de nuevas condiciones sobre los subconjuntos de
espacios de dimensión infinita y sobre los funcionales definidos en estos
espacios que nos permitan usar una generalización del Teorema de Weierstrass.
Puesto que los conjuntos cerrados y acotados, en el sentido de la topología

fuerte, de un espacio de Banach no son compactos, puede esperarse que si se
reduce la cantidad de abiertos mediante la introducción de la nueva topología, la
cantidad de cerrados y, por tanto, de compactos, aumente. Esto resulta ser así.
En particular, cualquier subconjunto cerrado y acotado de un espacio de Banach
es relativamente compacto respecto a la topología débil (es la topología menos
fina que hace continuas a las aplicaciones lineales).
El problema que surge a continuación es el de la continuidad (respecto a la

topología débil) del funcional a minimizar. Claramente, al introducir una
topología con menos abiertos, la cantidad de funciones continuas también
disminuye y así, por ejemplo, la norma asociada a la topología fuerte no es una
función continua respecto a la topología débil. Cobra especial importancia en
este contexto la noción de semicontinuidad inferior.
Finalmente, observamos que aunque la introducción de la topología débil y

de los funcionales semicontinuos inferiormente respecto dicha topología nos
permiten asegurar la existencia de un mínimo sobre cualquier conjunto cerrado
y acotado respecto la topología débil, la verificación práctica de estas
propiedades dista de ser sencilla. Por ello, una de las cuestiones centrales es la

búsqueda de condiciones expresadas respecto a la topología fuerte que
impliquen las correspondientes respecto a la topología débil. En este contexto la
convexidad de conjuntos y funciones juega un papel fundamental.

2. Condiciones necesarias de extremo.
Ecuación de Euler.
Analicemos el extremo del funcional:

() = , (), ′ ()

si los puntos frontera de las curvas admisibles están fijos ( ) = e ( ) =
, tal y como podemos apreciar en la siguiente imagen:
Además (, , ′ ) se considera derivable tres veces. Sabemos que la
el funcional. Supongamos que en la curva = (), derivando dos veces, se

condición necesaria para que haya un extremo es la anulación de la variación de
tiene un extremo (exigiendo sólo la existencia de derivadas de primer orden de
también segunda derivada). Tomemos cierta curva admisible = () cercana a

las curvas admisibles, se puede demostrar por otro método que la curva posee
la curva = () e incluyamos ambas curvas en la familia mono-paramétricas

de curvas:
(, ) = () + (() − ())
cuando = 0, se obtiene la curva = (); para = 1, se tiene = (). Lo

podemos ver en la siguiente imagen:

Recordemos que () − () se llama variación de la función () y se
designa por . Esta función se puede derivar una o varias veces, siendo
()′ = () − ()= ′, es decir, la derivada de la variación es igual a la
de la variación. De este modo, consideremos la familia = (, ), donde

variación de la derivada. Análogamente podemos verlo en derivadas sucesivas
(, ) = () + , que contiene para = 0 la curva en la cual se alcanza el

extremo, y para = 1 cierta curva admisible cercana llamada curva de
comparación. Si consideramos los valores del funcional:

() = , (), ′ ()

sólo en las curvas de la familia = (, ), el funcional se transforma en una

función de :

(, ) = ()
ya que el valor del parámetro determina una curva de la familia = (, )

determinando también con esto el valor del funcional
(, ). Esta función
() tiene un extremo en = 0, ya que para dicho valor se obtiene = ()
teniendo el funcional, por hipótesis, un extremo con respecto a cualquier curva
= (, ) . La condición necesaria para que la función () tenga un extremo

cerca admisible y, en particular, con respecto a las curvas cercanas de la familia
en = 0, es la anulación de su derivada para = 0, es decir, ′ (0) = 0. Como

() = , (, ), ′ (, )

Derivando:

" " ′
′() = ! (, ) + ! ′ (, )#
" "
"
donde
! = , (, ), ′ (, ),
"
"
!′ = , (, ), ′ (, ),
"′
o, puesto que,
" "
(, ) =
() + =
" "
" " ′()
′(, ) = $ + ′ % = ′
" "
se obtiene

′ ()
= $! &, (, ), ′(, )' + !′ &, (, ), ′(, )'′ %

sustituyendo por = 0, obtenemos:

′ (0) = $! &, (), ′()' + !′ &, (), ′()'′ %

Como hemos visto ′ (0) se llama variación del funcional, y se designa por .
anulación de su variación: = 0. Para el funcional:

La condición necesaria para que la funcional v tenga un extremo consiste en la

() = &, , ′ '

esta condición tiene la forma:

$! + !′ ′% = 0

′ = ()′ obtenemos que:

Integrando el segundo sumando por partes y tomando en cuenta que

= $!′ % + )! − *
!′
= ( ) − ( ) = 0 ,
Pero
+
=
+ = ( ) − ( ) = 0 ,
=
Considerando que todas las curvas admisibles en el problema simple pasan

por puntos frontera fijos, se tiene que:

= )! − *
!′
De este modo, la condición necesaria de extremo toma la forma:

= )! − * = 0 (-)
!′
donde el primer factor ! − . !′ es una función continua dada en la curva

.
= () que realiza el extremo, y el segundo factor es una función arbitraria

que satisface que se anula en los puntos frontera = y = es continua y
derivable una o varias veces, o bien y ′ son pequeños en valor absoluto.
Para simplificar la condición necesaria de extremo ( I ) obtenida, aplicaremos el
siguiente lema, conocido como Lema Fundamental del Cálculo de Variaciones:
Lema 2. (Fundamental del Cálculo de Variaciones). Sea / ∈ 1 (

, ) = {/ ∈
1 (
, ): /( ) = 0, /(4 ) = 0}, una función continua tal que:

(/() ()) = 0

para todo ∈ 1 (

, ). Entonces /() = 0 para todo ∈
, .
Demostración. Suponiendo que en el punto = ̅ contenido en el segmento

≤ ≤ sea /() ≠ 0, se llega a una contradicción. Así es, como / ∈
1 (
, ) entonces se deduce que si /(̅ ) ≠ 0, entonces /() conserva el signo
en cierto entorno ≤ ≤
del punto ̅ . Pero entonces tomando una función
() que también conserve su signo en este entorno y sea igual a cero fuera del
mismo tal y como aparece en la siguiente imagen:
se obtiene:

(/() ()) = (/() ()) ≠ 0

ya que el producto /() () conserva su signo en el segmento ̅ ≤ ≤ y se
tanto /() ≡ 0. La función () puede escogerse de la siguiente manera:

anula fuera del mismo. De este modo, hemos llegado a una contradicción; por lo
4; (
() ≡ 0 fuera del segmento ≤ ≤
;
() = :( − )
−
)

4;
≤ ≤
segmento ,
en el

Es evidente que la función () satisface las condiciones consideradas
donde n es un entero positivo, y k, un factor constante.
anteriormente: es continua, tiene derivadas continuas hasta de orden 2= − 1, se

anula en los puntos y y puede hacerse tan pequeña como se quiera en valor
absoluto, conjuntamente con sus derivadas, disminuyendo el módulo del factor
∎
k.
Apliquemos ahora el lema 2 para simplificar la condición necesaria ( I ),

obtenida anteriormente, de extremo del funcional inicial con el que
comenzamos esta sección,

)! − * = 0
!′
extremo, el factor ! − !′ es función continua, y la variación es una

.
Todas las condiciones del lema se cumplen: en la curva que realiza el
.
ya previstas en el lema fundamental. Por lo tanto, ! − !′ ≡ 0 en la curva

.
función arbitraria a la cual se ha impuesto sólo limitaciones de carácter general,
.
= () que realiza el extremo del funcional considerado, es decir, = () es
solución de la ecuación diferencial de segundo orden

! − =0
!′
o también:

! , (), ′ () − , (), ′ () = 0 , ∀ ∈ ( , )
!′
o bien de la forma desarrollada:
! − !′ − !!′ ′ − !′!′ ′′ = 0
Esta ecuación se denomina Ecuación de Euler, publicada por primera vez por
variaciones. Las curvas integrales de la ecuación de Euler = (, 1 , 14 ) se

Euler en 1744. La ecuación de Euler juega un papel fundamental en el cálculo de
llaman extremales. Sólo en las extremales puede alcanzarse un extremo del

funcional:

() = , (), ′ ()

Para hallar la curva que realiza un extremo del funcional, se integra la
la solución general de esta ecuación, de las condiciones de frontera ( ) =

ecuación de Euler y se determinan las dos constantes arbitrarias, que figuran en
, ( ) = . Sólo en las extremales que satisfacen estas condiciones se puede

realizar un extremo del funcional.
Recordemos que el problema de la frontera:

! − = 0, ( ) = , ( ) =
!′
no siempre tiene solución, y si existe, puede no ser única.

A continuación veremos un par de problemas variacionales donde la
existencia de la solución es evidente en el sentido físico o geométrico del
problema, y si la ecuación de Euler que satisface las condiciones de frontera es
única, esta única extremal será la solución del problema variacional
considerado.

Ejemplo 1. Veamos en qué curvas puede alcanzar su extremo el funcional:
C
4 B

() = @& ′ ' − 4 A , (0) = 0, = 1 ?
4
2
Solución. La ecuación de Euler tiene la forma: ′′ + 6 = 0 donde =

1 cos + 14 IJ= . Utilizando las condiciones de frontera, se obtiene que 1 = 0,
14 = 1, por lo tanto, el extremo puede alzarse sólo en la curva = IJ=.
∎
Ejemplo 2. Veamos en qué curvas puede alcanzar su extremo el funcional:

() = @& ′ ' − 12A , (0) = 0, (1) = 1 ?
4
Solución. La ecuación de Euler tiene la forma: ′′ + = 0 donde = K +

1 + 14 . Utilizando las condiciones de frontera, se obtiene que 1 = 0, 14 = 0,
por lo tanto, el extremo puede alzarse sólo en la curva = IJ=.
∎
En estos dos ejemplos la ecuación de Euler fue integrada fácilmente; pero

esto no siempre ocurre, puesto que las ecuaciones diferenciales de segundo
orden se integran de forma finita, solo en casos excepcionales.
A continuación veamos algunos casos simples de integración de la ecuación
de Euler:
1) F no depende de ′:
= (, )
La ecuación de Euler tiene la forma ! (, ) = 0 puesto que ! ≡ 0. La

solución de la ecuación finita ! (, ) = 0 obtenida no contiene elementos
( ) = J ( ) = .
arbitrarios y, por esto, en general no satisface las condiciones de frontera
no existe. Sólo en casos excepcionales, cuando la curva ! (, ) = 0 pasa

Por lo tanto, la solución del problema variacional considerado en general
por los puntos frontera ( , ), ( , ) existe una curva en la que se
puede alcanzar un extremo.
2) F depende de ′ en forma lineal:
&, , ′ ' = L(, ) + M(, )′

() = L(, ) + M(, ) #

La ecuación de Euler tiene la forma:
"L "M ′
+ − M(, ) = 0
" "
"L "M ′ "M "M ′

ó
+ − − =0
" " " "
o bien:
"L "M
− =0
" "
pero ésta, al igual que en el caso anterior, de nuevo es una ecuación finita
− N = 0 no satisface, en general, las
NO NP
N!
y no diferencial. La curva
condiciones de frontera. Por lo tanto, el problema variacional, por regla
general, no tiene solución en la clase de funciones continuas. Si, en
cambio,
"L "M
− ≡0
" "
la expresión L + M es una diferencial total, entonces:

() = L+M # =
L + M

no depende del camino de integración con lo que el valor del funcional

es constante en las curvas admisibles. Por tanto el problema variacional
pierde el sentido.
3) F depende sólo de ′:
= ( ′ )
La ecuación de Euler tiene la forma !′!′ ′′ = 0, puesto que ! = ! ′ =

!!′ = 0. De aquí se obtiene que ′′ = 0, o bien, !′!′ = 0. Si ′′ = 0,
entonces = 1 + 14 , que es una familia biparamétricas de líneas rectas.
Si la ecuación ! ′ ! ′ & ′ ' = 0 tiene una o varias raíces reales ′ = :,
entonces = :Q + 1, y obtenemos una familia monoparamétrica de
rectas contenida en la familia biparamétrica = 1 + 14 obtenida

anteriormente. De esta forma, en el caso = ( ′ ) todas las líneas rectas
posibles = 1 + 14 son extremales.
4) F depende sólo de x e ′:
= (, ′)
! ′ &, ′ ' = 0, y por lo tanto, tiene

.
.
La ecuación de Euler toma la forma
la primera integral ! ′ &, ′ ' = 1 . Además, como la ecuación de primer
orden obtenida ! ′ &, ′ ' = 1 no contiene a , ésta puede integrarse o
bien resolviéndola directamente respecto a ′ e integrando, o bien
introduciendo un parámetro escogido en forma adecuada.
5) F depende sólo de e ′:
= (, ′ )
La ecuación de Euler tiene la forma ! − !!′ ′ − !′!′ ′′ = 0, puesto que

! ′ = 0. Si se multiplica esta ecuación miembro a miembro por ′,
transforma en la derivada exacta . ( − ′ !′ ).

entonces, como no es difícil comprobar, el primer miembro se
.

En efecto,
& − ′ ! ′ ' = ! ′ + !′ ′′ − !!′ ′4 − !!′ ′ ′′

= ′ &! − !! ′ ′ − ! R ! R SS '.
Por consiguiente, la ecuación de Euler tiene la primera integral:
− S !S = 1
a x, puede ser integrada resolviéndola con respecto a ′ y separando

además, como esta ecuación de primer orden no contiene explícitamente
variables, o introduciendo un parámetro.
Para obtener las condiciones necesarias de extremo del funcional de tipo más
general:

, 4 , K , … , ; = (, , 4 , … , S , 4S , … , ;S )

con condiciones de frontera dadas para todas las funciones:

( ) = , 4 ( ) = 4 , … … … … … … … … … , ; ( ) = ;
( ) = , 4 ( ) = 4 , … … … … … … … … … , ; ( ) = ;
variaremos sólo una de las funciones:
U () (V = 1,2 … , =)
dejando las demás invariables.
Entonces el funcional
, 4 , K , … , ; se transforma en un funcional que
depende sólo de una función variable, por ejemplo, de Q (),

, 4 , K , … , ; = W
Q
del tipo considerado en el apartado anterior.
Por tanto, la función que realiza el extremo debe satisfacer la ecuación de

Euler:

!X − R=0
!X
Como este razonamiento es aplicable a cualquier función Q (Y = 1,2, … =), se

obtiene el sistema de ecuaciones diferenciales de segundo orden:

!X − R=0 (Y = 1,2, … , =)
!X
curvas integrales en el espacio , , 4 , … , ; , que es la familia de extremales del

que determinan, en general, una familia dependiente de 2n parámetros de
problema variacional dado.
Si, en particular, el funcional depende sólo de dos funciones () y Z():

(), Z() = (, , Z, S , Z S ) ;

( ) = , Z( ) = Z ( ) = Z( ) = Z
o sea, se determina eligiendo la curva alabeada = (), Z = Z(), tal y como

observamos en el siguiente gráfico:

Entonces variando sólo () y fijando Z(), cambiamos nuestra curva de tal
todo el tiempo en el cilindro de proyección Z = Z() como podemos observar:

modo que su proyección en el plano xOz no varía, es decir, la curva permanece
Análogamente, fijando () y variando Z(), variamos la curva de modo que

ésta permanezca todo el tiempo en el cilindro de proyección = (). Entonces
obtenemos un sistema de dos ecuaciones de Euler:

! − =0
!S

\ − =0
\S
Ejemplo 3. Hallar los extremales del funcional:

C
4 B B

() =
( S )4 + (Z S )4 + 2Z, (0) = 0, = 1, Z(0) = 0, Z = −1
2 2
Solución. El sistema de ecuaciones diferenciales de Euler tiene la forma:
SS − Z = 0
Z SS − = 0

segunda se obtiene que Q] − = 0 integrando esta ecuación lineal con
Despejando por ejemplo, z, de la primera ecuación y sustituyendo en la
coeficientes constantes tendremos que:
= 1 J + 14 J ^ + 1K cos + 1_ IJ=
Como Z = ′′ entonces Z = 1 J + 14 J ^ − 1K cos − 1_ IJ= . Utilizando las

condiciones de frontera, que se exponen en el enunciado, se halla que:
1 = 0 14 = 0 1K = 0 1_ = 1
Y por lo tanto = IJ= , Z = −IJ= .

■
Analicemos el extremo del funcional:

, 4 , K , … , ; = &, (), S (), … , (;) ()'

donde la función F se considera derivable = + 2 veces con respecto a todos los

argumentos y supondremos que las condiciones de frontera tienen la forma:
( ) = , ′( ) = S , … … … … … … … … … , (;^) ( ) =

(;^)
( ) = , ′( ) = S , … … … … … … … … … , (;^) ( ) =

(;^)
funciones, sino también de sus derivadas hasta de orden (= − 1) inclusive.

Es decir, en los puntos frontera están dados los valores no sólo de las
Supongamos que el extremo se alcanza en la curva = (), derivable 2n veces.

Sea = () la ecuación de cierta curva de comparación, también derivable 2n
veces.
Consideramos la función monoparamétrica de funciones (, ) = () +

() − (), o bien, (, ) = () + . Para = 0, (, ) = (); para
= 1, (, ) = (). Si consideramos el valor del funcional
() sólo en las
curvas de la familia = (, ), entonces ésta se transforma en una función del
parámetro , la cual alcanza su extremo para = 0. Por lo tanto,
(, )
.
.
restringido a = 0, vale 0. Esta derivada se llama variación del funcional y se
designa por :

= ` &, (, ), S (, ), … , (;) (, )'a
=0

= ( ! + ! R S + ! RR SS + ⋯ + ! (c) (;) )

Integramos por partes una vez el segundo sumando del segundo miembro:

!S = $! R % −
S

!S


el tercer sumando, dos veces:

4
!SS = $SS
′% − # +

! RR
!SS
4 !SS
y así sucesivamente. El último sumando n veces, tenemos que:

! (c) (;) = $! (c) (;^) % − !(;) (;^4) # + ⋯ + (−1);

;
=

; ! (c)
variaciones = S = SS = ⋯ = (;^) = 0 para = y para = ,

Tomando en cuenta las condiciones de frontera, en virtud de las cuales las
obtenemos por último:

4 ;
= d! − !S + 4 !SS + ⋯ + (−1); ; ! (c) e

Como en la curva que realiza el extremo se tiene que:

4 ;
= d! − + + ⋯ + (−1) ;
(c) e = 0
!S 4 !SS ; !
para funciones arbitrarias, y como el primer factor bajo el símbolo integral es

función continua de x en la misma curva = (), entonces, debido al lema
fundamental, el primer factor es idénticamente nulo:
4 ;
! − !S + 4 !SS + ⋯ + (−1); ; ! (c) ≡ 0

De este modo, la función = (), que realiza el extremo del funcional:

() = &, (), S (), … , (;) ()'

debe ser solución de la ecuación
4 ;
! − !S + 4 !SS + ⋯ + (−1); ; ! (c) = 0

Esta ecuación diferencial de orden 2n recibe el nombre de ecuación de Euler-

Poisson, y sus curvas integrales se denominan extremales del problema
variacional considerado. La solución general de esta ecuación contiene 2n
constantes arbitrarias, las cuales pueden ser, en general, determinadas a partir
de las 2n condiciones de frontera:

( ) = , ′( ) = S , … … … … … … … … … , (;^) ( ) =
(;^)
( ) = , ′( ) = S , … … … … … … … … … , (;^) ( ) =

(;^)
Ejemplo 4. Hallar la extremal del funcional:

() = (1 + SS4 ), (0) = 0, ′(0) = 1, (1) = 1, ′(1) = 1

Solución. La ecuación de Euler-Poisson tiene la forma f (2 SS ) = 0 ó Q] = 0

.f
.
su solución general es = 1 K + 14 4 + 1K + 1_ . Utilizando las condiciones de
frontera, obtenemos:
1 = 0 14 = 0 1K = 1 1_ = 0
De esta manera, el extremo puede alcanzarse sólo en la recta = .

■
Ejemplo 5. Determinar la extremal del funcional:

C
4 B B

() =
SS4 − 4 + 4 , (0) = 1, ′(0) = 0, = 0, Z = −1
2 2
Solución. La ecuación de Euler-Poisson tiene la forma Q] − = 0 su solución

general es = 1 J + 14 J ^ + 1K cos + 1_ IJ= . Utilizando las condiciones de
1 = 0 14 = 0 1K = 1 1_ = 0
frontera, obtenemos:
De esta manera, el extremo puede alcanzarse sólo en la recta = cos .

■
Veamos que ocurre cuando el funcional tiene la forma:

(), Z() = &, , S , … , (;) , Z, Z S , … , Z (h) ' ;

variando sólo () y considerando Z() fija, se halla que las funciones () y
Z() que realizan el extremo deben satisfacer la ecuación de Euler-Poisson:
4 ;
! − !S + 4 !SS + ⋯ + (−1); ; ! (c) = 0;

Variando Z() y considerando () fija, obtenemos que las mismas funciones
deben satisfacer la ecuación:

4 h
\S + 4 \SS + ⋯ + (−1)h h \ (i) = 0;
\ −

De esta manera las funciones Z() e () deben satisfacer el sistema de
ecuaciones:
4 ;
! − + + ⋯ + (−1) ;
(c) = 0;
!S 4 !SS ; !
4 h
\ − \S + 4 \SS + ⋯ + (−1)h h \ (i) = 0;

De forma completamente análoga se puede razonar también al analizar el

extremo de un funcional que depende de un número arbitrario de funciones:

, 4 , … , h = , , S , … , , 4 , 4S , … , 4 , … , h , h , … , h i ;
(; ) (;f ) S (; )

Variando sólo alguna Q () y considerando invariables las demás, obtenemos
;X
la condición necesaria fundamental de extremo en la forma:
!X − ! R + ⋯ + (−1);X ; (cX) = 0 (Y = 1,2, … , j)
X X !X
2.1 Funciones que dependen de funciones de

varias variables independientes.
Analicemos el extremo del funcional
"Z "Z

Z(, ) = k ), , Z, , * ;
" "
l
Además en la frontera C de la región D los valores de la función Z(, ) están

dados, es decir, está dado un contorno alabeado 1̅ , por el cual deben pasar todas
las superficies admisibles, tal y como vemos en la siguiente figura:
= , = m. La función F se
N\ N\
N N!
considerará derivable tres veces. La superficie Z = Z(, ), en la cual se realiza el
Para abreviar la escritura, designamos
extremo, se supondrá derivable dos veces.

Z = Z(, , ) = Z(, ) + Z, siendo Z = Z̅(, ) − Z(, ), que contiene, cuando
Consideramos nuevamente la familia monoparamétrica de superficies
= 0, la superficie Z = Z(, ), en la cual se realiza el extremo y, para = 1,

cierta superficie admisible Z = Z̅(, ). En las funciones de la familia Z(, , ), el
funcional se transforma en una función de , la cual debe tener un extremo
cuando = 0; por lo tanto
Z(, )
= 0. Llamando variación del
N\
Nn =0
funcional a la derivada de
Z(, , ) con respecto a cuando = 0,
designándola por , tendremos que:
"Z

Z(, ) = o k (, , Z(, , ), (, , ), m(, , ))p
" =0
l
= k
\ Z + q + r m
l
Donde
Z(, , ) = Z(, ) + Z

"Z(, , )
(, , ) = = (, ) +
"
"Z(, , )
m(, , ) = = m(, ) + m
"
" "
Como
sq Zt = s tZ + q
" " q
"
sr Zt = s tZ + r m
" " r
Entonces
k(q + r m)

l
" " " "

=k sq Zt + sr Zt# − k sq t + s t# Z
" " " " r
l l
sq t es la llamada derivada parcial completa o total con respecto a x.
N
N
donde
Al calcularla, se considera fija, pero la dependencia de Z, , m de x se toma en

cuenta:
" "Z " "m

sq t = q + q\ + qq + qr
" " " "
y análogamente:
" "Z " "m
sr t = r + r\ + rq + rr
" " " "
En virtud de la fórmula de Green:
"M "L
k + # = k(M − L)
" "
l u
se obtiene
" "
k sq Zt + s Zt# = &q − r 'Z = 0
" " r
l u
Z = 0, puesto que todas las superficies admisibles pasan por el mismo contorno
La última integral es igual a cero, debido a que en el contorno C, la variación
alabeado 1̅ . Por lo tanto:
" "
k
\ Z + r m = − k sq t + s t# Z
" " r
l l
y la condición necesaria de extremo:
k
\ Z + q + r m = 0
l
toma la forma:
" "
k )\ − sq t − { }* Z = 0
" " r
l
Como la variación Z es arbitraria (a Z se le imponen limitaciones sólo de

carácter general con respecto a la continuidad y a la derivabilidad, anulación en
enunciado anteriormente, en la superficie Z = Z(, ) que realiza el extremo

el contorno C, etc.) y el primer factor es continuo por el lema fundamental
será:
" "
\ − sq t − { } ≡ 0
" " r
Por consiguiente, Z = Z(, ) es solución de la ecuación:
" "
\ − sq t − s t = 0
" " r

debe satisfacer la función Z(, ) que realiza el extremo, lleva el nombre de
Esta ecuación diferencial de segundo orden en derivadas parciales, a la cual
ecuación de Ostrogradski, en honor al eminente matemático ruso

M.V.Ostrogradski, quien la obtuvo por primera vez en el año 1834; sin embargo,
para las regiones D rectangulares se encontraba ya en los trabajos de L. Euler.
Dado el funcional de la siguiente forma:

Z( , 4 , … , ; ) = k … … ( , 4 , … , ; , Z, , 4 , … , ; ) 4 … ; ,
v
donde Q = en base a la condición necesaria fundamental de extremo = 0,

N\
NX
se obtiene, en forma completamente análoga, la siguiente ecuación de
Ostragradski:
;
"
\ − w s t = 0
"Q qX
Qx
a la cual debe satisfacer la función:
Z = Z( , 4 , … , ; )
que realiza el extremo del funcional .
Por ejemplo, para el funcional:
"Z "Z " 4 Z " 4 Z " 4 Z

Z(, ) = k (, , Z, , , , , )
" " " 4 "" " 4
l
se obtiene la ecuación
" " "4 "4 "4

\ − sq t − sr t + 4 {y } + {z } + 4 {{ } = 0
" " " "" "
donde
"Z "Z " 4Z " 4Z " 4Z

= ;m = ;| = 4 ;I = ;} = 4
" " " "" "
La función que realiza el extremo del funcional debe satisfacer esta

ecuación de cuarto orden en derivadas parciales.
Por ejemplo, para el funcional:

" 4Z " 4Z " 4Z
4 4 4
= ~d 4 e + d 4 e + 2 d e
" " " "
l

La función z que realiza el extremo debe satisfacer la llamada ecuación
" _Z " _Z " _Z

biarmónica
+ 2 +
" _ " 4 " 4 " _
Que habitualmente se escribe en forma compacta así: ∆∆Z = 0.
Para el funcional
" 4Z " 4Z " 4Z

4 4 4
= ~d 4 e + d 4 e + 2 d e − 2Z(, )
" " " "
l
La función z(x,y) que realiza el extremo debe satisfacer la ecuación ∆∆Z =

(, ).
4
Los problemas sobre el extremo funcional
" 4Z " 4Z
= d 4 e + d 4 e
" "
l
o del funcional de forma más general
" 4Z " 4Z " 4Z " 4Z " 4Z
4 4
= d 4 + 4 e − 2(1 − ) ~ 4 4 − d e
" " " " ""
l
donde es un parámetro, se reducen también a la ecuación biarmónica.

3. Condiciones suficientes de extremo.
Si en el plano (, ) por cada punto de cierta región v pasa una y sólo una
curva de la familia = (, 1), se dice que esta familia de curvas forma un
campo en la región v, o, más exactamente, un campo propio. El coeficiente
angular de la tangente (, )a la curva de la familia = (, 1) que pasa por el
punto (, ) se llama inclinación (o declive) del campo en el punto (, ).
Por ejemplo, las rectas paralelas = + 1 forman un campo dentro del

círculo x 2 + y 2 ≤ 1 (Figura 1) y su inclinación es (, ) = 1. Por el contrario, la
familia de parábolas y = ( x − a ) 2 − 1 (Figura 2) no forma un campo dentro del
mismo círculo, debido a que en su interior las parábolas de la familia
considerada se cortan.
Figura 1
Si todas las curvas de la familia = (, 1) pasan por cierto punto ( x0 , y 0 ) , es
campo propio en la región v si el centro del haz pertenece a ésta. Sin embargo,
decir, forman un haz de curvas, entonces éstas con seguridad no forman un
si las curvas del haz cubren toda la región v y se cortan en su interior sólo en el
centro del haz, es decir, se cumplen las condiciones impuestas al campo en
= (, 1) forma también un campo, llamado en este caso central, a diferencia

todos los puntos distintos del centro del haz, se die que la familia
del campo propio (Figura 3).
Figura 2
Por ejemplo, el haz de sinusoides y = Csen (x ) forma un campo central para
0 ≤ x ≤ a , a ≤ π (Figura 4). El mismo haz de sinusoides forma un campo propio
en un entrono suficientemente pequeño del segmento δ ≤ x ≤ a del eje de las
abscisas, donde δ > 0 , a < π (Figura 4). El mismo haz de sinusoides no forma
un campo en un entorno del segmento 0 ≤ x ≤ a1 , a1 > π , del eje de las abscisas
(Figura 4).
Si un campo central o propio está formado por una familia de extremales de

cierto problema variacional, se llama campo de extremales.
El concepto de campo se generaliza casi sin modificaciones también al caso

de un espacio de cualquier dimensión. La familia yi = y i ( x, C1 ,..., C n )
(i = 1,2,..., n) forma un campo en la región D del espacio x, y1 ,..., y n .
Figura 3 Figura 4
si por cada punto de dicha región pasa una y sólo una curva de la familia
yi = y i ( x, C1 ,..., C n ) . Se llaman funciones de inclinación del campo
pi = ( x, y1 , y 2 ,..., y n ) (i = 1,2,..., n) a las derivadas parciales de las funciones
yi = y i ( x, C1 ,..., C n ) con respecto a x, calculadas en el punto ( x, y1 , y 2 ,..., y n ) ; por
∂
consiguiente, para obtener pi = ( x, y1 , y 2 ,..., y n ) hay que tomar y i ( x, C1 ,..., C n )
∂x
y sustituir (C1 ,..., C n ) por sus expresiones mediante las coordenadas
x, y1 , y 2 ,..., y n . En forma análoga se define también el campo central.
Supongamos que la curva y = y (x ) es extremal del problema variacional

sobre el extremo del funcional simple
x1
v[ y ( x)] = ∫ F ( x, y, y ' )dx
x0
con puntos frontera A( x o , y 0 ) y B( x1 , y1 ) fijos. Se dice que la extremal y = y (x )

está incluida en un campo de extremales, si ha sido dada una familia de
extremales y = y ( x, C ) que forma un campo,

Figura 5 Figura 6
contiene a la extremal y = y (x ) para cierto valor C = C 0 y dicha extremal no

pertenece a la frontera de la región D en la cual la familia y = y ( x, C ) forma un
campo (Figura 5). Si el haz de extremales con centro en el punto A( x o , y 0 ) forma
un campo en un entrono de la extremal y = y (x ) , que pasa por este punto,
entonces con esto se halla un campo central que contiene a la extremal dada
y = y (x ) . Como parámetro de la familia se puede tomar en este caso el
coeficiente angular de la tangente a las curvas del haz en el punto A( x o , y 0 )
(Figura 6).
Ejemplo 6. Se da el funcional
∫ ( y'
2
− y 2 )dx
0
se pide incluir el segmento de la extremal y = 0 que une los puntos (0, 0) y (a,
0), donde 0 < a < π , en un campo central de extremales.
Solución. La solución general de la ecuación de Euler y ' '+ y = 0 tiene la forma

y = C1 cos( x) + C 2 sen( x) . De la condición de que las extremales pasen por el
punto (0, 0) se obtiene C1 = 0 , y = C 2 sen( x) ; las curvas de este haz forman un
campo central en el segmento 0 ≤ x ≤ a , a < π que incluye la extremal y = 0
para C 2 = 0 . El parámetro C2 de la familia es igual a la derivada y x en el punto
(0, 0). Si en este mismo problema fuera a ≥ π , la familia y = C 2 sen( x) no
formaría campo.
Es sabido que dos curvas infinitamente cercanas de la familia F ( x, y , C ) = 0 se

cortan en los puntos de la curva C-discriminante, la cual se determina por las
ecuaciones
∂F
F ( x, y , C ) = 0 =0
∂C
Figura 7
Recordemos que en la curva C-discriminante se incluyen, en particular, la

envolvente de la familia y los lugares geométricos de los puntos múltiples de las
curvas de dicha familia. Si F ( x, y , C ) = 0 es la ecuación de un haz de curvas, se
centro pertenece también a la curva C-discriminante. Por esto, si se toma un haz
de extremales y = y ( x, C ) que pasan por el punto ( x0 , y 0 ) y se determina su
curva C-discriminante Φ ( x, y ) = 0 .
En particular, las curvas de este familia próximas a la extremal considerada

y = y (x ) , que pasa por los puntos A( xo , y 0 ) y B( x1 , y1 ) , se cortarán en puntos
cercanos a los puntos de tangencia (o de intersección) de la curva y = y (x ) con
la curva C-discriminante (véase la Figura 7, en donde la curva C-discriminante
está representada por una línea gruesa). Si el arco AB de la extremal y = y (x ) no
tiene puntos comunes diferentes de A con la curva C-discriminante del haz de
extremales que incluye a la extremal dada, entonces las extremales del haz
suficientemente próximas al arco AB no se cortan, es decir, forman un campo
central que incluye al arco AB en un entorno de este arco (Figura 8).
Figura 8

Si el arco de la extremal y = y (x ) tiene un punto común A* , diferente del
punto , con la curva C-discriminante del haz y = y ( x, C ) entonces las curvas
las proximidades del punto A* se llama conjugado del punto .

del haz próximas a y = y (x ) pueden cortarse entre sí y con la curva y = y (x ) en
de extremales con centro en el punto que contenga el arco de la extremal,

El resultado obtenido se puede enunciar así: para construir un campo central
es suficiente que el punto A* conjugado del no pertenezca al arco . Esta

condición de la posibilidad de la construcción de un campo de extremales que
incluya a la extremal dada se llama condición de Jacobi.
y = y ( x, C ) la ecuación del haz de extremales con centro en el punto ; el

No es difícil formular esta condición también en forma analítica. Sea
parámetro 1 se puede considerar, para fijar ideas que coincide con el coeficiente
angular y’ de las extremales del haz en el punto . La curva C-discriminante se
determina por las ecuaciones
∂y ( x, C )
y = y ( x, C ); =0
∂C
∂y ( x, C )
A lo largo de cada curva fija de la familia la derivada = 0 es una
función sólo de . Denotaremos abreviadamente por dicha función:
∂C
1
∂y ( x, C ) ∂ 2 y ( x, C )
u= , donde está dado; de aquí u ' x = . Las funciones
∂C ∂C∂x
y = y ( x, C ) son soluciones de la ecuación de Euler por lo tanto,
d
Fy ( x, y ( x, C ), y ' x ( x, C )) − Fy ( x, y ( x, C ), y ' x ( x, C )) ≡ 0 .
dx
Derivando esta identidad con respecto a 1 y haciendo

∂y ( x, C )
= u , se obtiene
∂C
d
Fyy u + Fyy ' u '− ( Fyy u + Fy ' y ' u ' ) = 0 ,
dx
o bien
d d
( Fyy u − Fyy ' )u − ( Fy ' y ' u ' ) = 0 .
dx dx
Aquí Fyy ( x, y , y ' ) , F yy ' ( x, y , y ' ) y Fy ' y ' ( x, y , y ' ) son funciones conocidas de ,
debido a que el segundo argumento y es igual a la solución y = y ( x, C ) de la
. Esta ecuación lineal homogénea de segundo orden con respecto a u se

ecuación de Euler, tomada para el valor C = C 0 que corresponde a la extremal
llama ecuación de Jacobi.

∂y ( x, C )
Si la solución u = de este ecuación que se anula en el centro del haz
∂C
para x = x0 (el centro del haz siempre pertenece a la curva C-discriminante) se
conjugado de , que se determina por las ecuaciones

anula también en algún otro punto del intervalo x0 < x < x1 , entonces el punto
∂y ( x, C )
y = y ( x, C 0 ) y = 0 o bien u = 0 ,
∂C
Pertenece al arco de la extremal. Si, en cambio, existe una solución de la
del segmento x0 ≤ x ≤ x1 , entonces no hay puntos conjugados de en el arco ;

ecuación de Jacobi que se anule para x = x0 y no se anule en ningún otro punto
la condición de Jacobi se cumple y se puede incluir el arco de la extremal en

un campo central de extremales con centro en el punto .
que se alcance un extremo, es decir, para la curva que realiza un extremo el

Observación. Se puede demostrar que la condición de Jacobi es necesaria para
punto conjugado de no puede estar en el intervalo x0 < x < x1 .

4. Condición de Legendre
4.1.Condición necesaria de Legendre para la

realización de un mínimo de un funcional del
tipo variacional.
La condición necesaria para la realización de un mínimo de un funcional del
tipo variacional

() = &, (), S ()'

con condición de frontera fija, y ( x0 ) = y 0 y y ( x1 ) = y1 , es la condición de

Legendre
F y ' y ' ( x , y ( x ), y ' ( x )) ≥ 0 para todo x ∈ ( x 0 , x1 )
Legendre, en analogía al caso de dimensión finita, intentó demostrar, sin

éxito, que una condición suficiente para que v tenga un mínimo en y es que
satisfaga la desigualdad estricta
F y ' y ' ( x , y ( x ), y ' ( x )) > 0 para todo x ∈ ( x 0 , x1 )
Ahora, hallamos la siguiente expresión para la segunda variación
1 
x1 d  
δ y2 v (u ) =
2 ∫x0
  F yy −
 dx
F yy ' u 2 + F y ' y 'u ' 2 dx
 
con u ∈ C0 ([x0 , x1 ]) . Por brevedad, la escribimos como

1
1
∫ (Pu ' )dx
x1
δ y2 F (u ) = 2
+ Qu 2
2 x0
La idea de Legendre fue escribir esta expresión en la forma
1 x1
δ y2 F (u ) =
2 ∫x 0
(
Pu ' 2 + 2 wuu '+ (Q + w ')u 2 dx ) (5)
siendo w una función derivable arbitraria, y donde hemos usado la relación

d
0= ∫
x0
x1
dx
(
wu 2 = ) ∫ (w ' u x1
x0
2
)
+ 2 wuu ' dx
puesto que u ( x 0 ) = u ( x1 ) = 0 . Seguidamente, observó que la condición F y ' y ' > 0

sería suficiente si se pudiera encontrar una función w para la que el integrando
de (5) fuera un cuadrado perfecto. Sin embargo, esto no es siempre posible,
como el mismo Legendre demostró, puesto que w debería satisfacer la ecuación
P (Q + w') = w 2
que no posee, necesariamente, una solución global en todo el intervalo ( x 0 , x1 )

(Por ejemplo, si P = −1 y Q = 1 , la solución viene dada por w( x ) = tan (c − x ) . Si
x1 − x 0 > π , no hay solución en todo el intervalo ( x 0 , x1 ) , puesto que tan (c − x )
se hace infinita en algún punto de dicho intervalo). Legendre concluyó que
condiciones de tipo local como la ecuación de Euler o la condición estricta de
Legendre no podían ser las únicas condiciones suficientes para la realización de
un mínimo.
4.2 Condición suficiente de Legendre de

inclusión de un extremal de un funcional en
un campo de extremales.
La condición suficiente de inclusión de un extremal de un funcional

() = &, (), S ()'

con condición de frontera fija, y ( x 0 ) = y 0 y y ( x1 ) = y1 en un campo de
extremales es el cumplimiento de la condición reforzada de Legendre
F y ' y ' ( x , y ( x ), y ' ( x )) > 0

x ∈ ( x0 , x1 )
en todos los puntos de la extremal considerada (es decir,para todo )
Ejemplo 7. Mostrar que la extremal del problema variacional

4

() = (′_ + ′4 )

con condición de frontera y (0 ) = 1, y (2 ) = 5 se puede incluir en un campo de

extremales.

Solución. Las extremales son las rectas y = C1 x + C 2 . La extremal que satisface
la condición de contorno dadas es la recta y = 2 x + 1. En este caso

F y ' y ' = 12 y '2 + 2 y = 2 x +1 tenemos
y en todos los puntos de la extremal
F y ' y ' = 50 > 0 .
La condición reforzada de Legendre se cumple y, por tanto, la
y = 2 x +1 puede ser incluida en un campo de extremales.
extremal
Esto se ve sin necesidad de realizar cálculos. La extremal y = 2 x +1 pertenece
a la familia uniparamétrica de extremales y = 2 x + α , donde α es un parámetro,
las cuales forman un campo propio.
■

5. Condición de Jacobi
5.1. Condición necesaria de Jacobi.
En esta sección estudiaremos las condiciones bajo las cuales el funcional
x1
(
G (u ) = ∫ Pu ' 2 +Qu 2 dx
x0
) (6)
definido para u ∈ C0 ([x0 , x1 ]) , es definido positivo. Se puede observar, que en

1
relación con la función lagrangiana del problema de optimización se obtiene

este mismo funcional, en concreto
1 1 d 
P= Fy ' y ' , Q =  Fyy − Fy ' y ' 
2 2 dx 
De momento estudiaremos la positividad del funcional (6) como un problema

independiente.
Una condición necesaria para la no negatividad del funcional (6) es que
P( x ) ≥ 0 para todo x ∈ ( x 0 , x1 ) .
En esta sección asumiremos la condición
P(x ) > 0 para todo x ∈ ( x 0 , x1 ) .
y buscaremos condiciones necesarias y suficientes para que dicho funcional sea

definido positivo. Comenzamos escribiendo la ecuación de Euler asociada al
funcional (6):
d
− (Pu') + Qu = 0 (7)
dx
Esta es una ecuación diferencial lineal de segundo orden. Las condiciones de

frontera son u ( x 0 ) = u ( x1 ) = 0 . Este problema tiene la solución trivial u ≡ 0 . Sin
embargo, también puede tener soluciones no triviales (Por ejemplo, para
P = Q = 1 , la función u ( x ) = Csex es solución del problema en el intervalo (0, π ) ,
para toda constante C ). En este contexto, se introduce la siguiente definición:
Definición 3. Diremos que el punto x es un punto conjugado de x 0 si la

ecuación (7) tiene una solución que se anula en x 0 y x pero que no es
idénticamente nula.

Observación. Si u es una solución no idénticamente nula de (7) entonces
también lo será Cu , para cualquier C ≠ 0 . Hay varios criterios de
normalización que se pueden imponer para forzar la unidad de solución de este
problema. Aquí asumiremos que u ' (x 0 ) = 1 , que es siempre posible eligiendo C
adecuadamente. En efecto, si u ( x0 ) = 0 y v no es idénticamente nula entonces
u ' ( x0 ) no puede ser nulo, por el teorema de unicidad para la ecuación lineal (7).
Teorema 4. Supongamos que P ( x ) > 0 en [x 0 , x1 ] y que no hay puntos

conjugados en dicho intervalo. Entonces el funcional cuadrático
x1
(
G (u ) = ∫ Pu ' 2 +Qu 2 dx
x0
)
es definida positiva para toda u ∈ C0 ([x0 , x1 ]) .
1
Demostración. Para demostrar que el funcional G es definido positivo lo

reduciremos a la forma
x1
∫
x0
Pϕ 2 dx
donde ϕ 2 es cierta expresión cuya anulación implica que u ≡ 0 . Para conseguir

esto, comenzamos sumando a G la función cero expresada en la forma
d
∫x0
x1
dx
(
wu 2 dx )
A continuación, seleccionamos la función diferenciable w de modo que la
expresión
d
Pu ' 2 +Qu 2 +
dx
( )
wu 2 = Pu ' 2 +2 wuu '+ (Q + w')u 2 (8)
sea un cuadrado perfecto. Esto será cierto si w es solución de
P (Q + w') = w 2 (9)
ya que en tal caso (8) puede ser escrito como
2
 w 
P u '+ u 
 P 
De este modo, si (9) tiene una solución definida en todo el intervalo [x0 , x1 ],
entonces G puede escribirse como

2
x1  w 
∫x0 P u '+ P u  dx (10)
y es, por tanto, definido positivo. De hecho, si (10) se anula, entonces debe ser
w
u '+ u≡0
P
puesto que, por hipótesis, es P ( x ) > 0 en [x 0 , x1 ] . Pero esta ecuación de primer

orden, con la condición u ( x0 ) = 0 tiene por única solución a u ≡ 0 .
La demostración se reduce, pues, a comprobar que si no hay puntos

conjugados de x 0 en el intervalo [x 0 , x1 ] entonces la ecuación (9) tiene una
solución global en [x 0 , x1 ] . Esta ecuación diferencial, es una ecuación de Ricatti,
que puede ser reducida a una ecuación lineal de segundo orden mediante el
cambio
z'
w=− P (11)
z
donde z es una nueva incógnita. Con este cambio (9) se transforma en
d
− (Pz ') + Qz = 0 (12)
dx
que es, justamente, la ecuación de Euler de G . Ahora, puesto que no hay puntos
conjugados en [x 0 , x1 ] se sigue, por definición, que (12) tiene una solución que
no se anula en [x 0 , x1 ] . Por tanto, la ecuación (9) tiene una solución, dada por
(11), definida en todo el intervalo [x 0 , x1 ] .
■
A continuación veremos que la condición de que no existan puntos

conjugados de x 0 en el intervalo [x 0 , x1 ] no es solo una condición suficiente sino
también necesaria. Comenzamos con un lema que usaremos en la demostración.
Lema 5. Si la función u satisface la ecuación
d
− (Pz ') + Qz = 0
dx
y las condiciones de frontera
u ( x0 ) = u ( x1 ) = 0
entonces

∫ (Pu ' )
x1
2
+Qu 2 dx = 0
x0
Demostración. El lema es una consecuencia inmediata de la fórmula
x1  d  x1
0 = ∫  − (Pu ') + Qu udx = ∫ Pu ' 2 +Qu 2 dx ( )
x0
 dx  x0
que se obtiene por integración por partes y el uso de las condiciones de frontera.
■
Teorema 6. Supongamos que P ( x ) > 0 para todo x 0 en el intervalo [x 0 , x1 ] . Si

el funcional cuadrático G es definido positivo para toda u ∈ C0 ([x0 , x1 ])
1
entonces el intervalo [x 0 , x1 ] no posee puntos conjugados de x 0 .
Demostración. La idea de la demostración es la de construir una familia de

funcionales definidos positivos dependientes de un parámetro, t , tales que para
t = 1 se reduce al funcional G , mientras que para t = 0 nos da el funcional
cuadrático
x1
∫
x0
u ' 2 dx
el cual, claramente, no posee puntos conjugados de x 0 en el intervalo [x 0 , x1 ] .

Entonces, demostraremos que cuando variamos el parámetro t en [0,1] , no
pueden aparecer puntos conjugados en [x 0 , x1 ] .
Consideremos el funcional
∫ [t (Pu ' ) ]
x1
2
+ Qu 2 + (1 − t )u ' 2 dx (13)
x0
es definido positivo para todo t ∈ [0,1] , puesto que G lo es por hipótesis. La

ecuación de Euler correspondiente a este funcional es
d
− [(tP + (1 − t ))u '] + tQz = 0 (14)
dx
Sea u ( x, t ) una solución de (14) tal que u ( x 0 , t ) = 0 y u x ( x0 , t ) = 1 para todo

t ∈ [0,1] . Esta solución depende con continuidad del parámetro t , que para t = 1
se reduce a la solución v de la ecuación (7) con las condiciones

u ( x0 ,1) = u ( x1 ,1) = 0 , y para t = 0 se reduce a la solución de u' ' = 0 con las
mismas condiciones de frontera, es decir, u ( x0 ,0 ) = x − x0 .
Supongamos ahora que el intervalo [x 0 , x1 ] contiene un punto x conjugado

de x 0 . Necesariamente será x < x1 porque, si x = x1 , el lema anterior implica
que existe una v no idénticamente nula tal que G (v ) = 0 , contradiciendo la
hipótesis de positividad de G . Por tanto, la demostración se reduce a
comprobar que no puede haber un punto conjugado, x , en el interior de [x 0 , x1 ]
Para demostrarlo, consideremos el conjunto de puntos
C = {(x, t ) ∈ [x0 x1 ]× [0,1] : u (x, t ) = 0}
Si mostramos que en los puntos en los que u ( x, t ) = 0 no puede tenerse

u x (x 0 , t ) = 0 entonces podremos deducir el teorema de la función implícita que
el conjunto C representa una curva diferenciable en el plano xt . Tenemos que,
si u ( x*, t *) = 0 en algún ( x*, t *) entonces debe ser u x ( x*, t *) ≠ 0 ya que para
cualquier t fijo u ( x, t ) satisface la ecuación (14), y si se tuviera
u (x*, t *) = u x (x*, t *) = 0 entonces debería ser u ( x, t *) = 0 para todo x ∈ [x 0 , x1 ]
debido al teorema de unicidad para ecuaciones diferenciales lineales. Pero esto
es imposible puesto que u x es una función continua en [x 0 , x1 ]× [0,1] y
u x (x0 , t ) = 1 para todo t ∈ [0,1] .
Tenemos entonces que el teorema de la función implícita nos asegura la

existencia de una curva, x(t ) , tal que u ( x(t ), t ) = 0 en un entorno de cada punto
( )
de C . Por hipótesis, el punto x,1 pertenece a dicha curva. Partiendo de este
punto tenemos que:
A. La curva no puede terminar en le interior de [x 0 , x1 ]× [0,1] , pues sería

una condición de la dependencia continua de u ( x, t ) respecto del
parámetro t .
B. La curva no puede cortar el segmento {x = x1 , t ∈ [0,1]} , puesto que, por
el mismo argumento que el del lema anterior, pero aplicado a la
ecuación (14), las condiciones de frontera u (x 0 , t ) = u (x1 , t ) = 0 y el
funcional (13), se tendría una contradicción de la positividad del
funcional para todo t .
C. La curva no puede cortar el segmento t = 1 , x ∈ [x 0 , x1 ] , puesto que
tendríamos u ( x, t ) = u x ( x, t ) = 0 para algún ( x, t ) .
D. La curva no puede cortar el segmento {t = 0, x ∈ [x 0 , x1 ]} , puesto que
para t = 0 la ecuación (14) se reduce a u ' ' = 0 cuya solución sólo se
anula en x = x0 .

E. La curva no puede aproximarse al segmento {x = x 0 , t ∈ [0,1]} puesto
que tendríamos u x ( x 0 , t ) = 0 para algún t , contrario a nuestras
hipótesis.
Podemos ver esto en la siguiente figura:
Se sigue que tal curva no puede existir, con lo que se concluye la

demostración.
■
Si reemplazamos la condición de que el funcional G sea definido positivo

por la de que sea no negativo obtenemos el siguiente resultado.
Corolario 7. Si G(u ) , con P ( x ) > 0 para todo x 0 en el intervalo [x 0 , x1 ] , es

definido no negativo entonces el intervalo [x0 , x1 ) no contiene puntos
conjugados de x 0 .
Demostración. La única diferencia con la demostración del teorema anterior

es que no podemos asegurar que el funcional auxiliar dado por (13) sea definido
positivo en t = 1 . Así, no se puede excluir la posibilidad de que x = x1 .
■
Los teoremas anteriores se combinan en el siguiente enunciado.
Teorema 8. El funcional cuadrático
∫ (Pu ' )
x1
2
+Qu 2 dx
x0
con P ( x ) > 0 para todo x 0 en el intervalo [x 0 , x1 ] , es definido positivo para toda

u ∈ C01 ([x0 , x1 ]) si y solo si el intervalo [x 0 , x1 ] no contiene puntos conjugados en
x0 .

A continuación aplicaremos estos resultados al problema de optimización
x1
J (u ) = ∫ L (x , u (x ), u ' (x ))dx
x0
(15)
con las condiciones de frontera u ( x0 ) = u 0 y u ( x1 ) = u1 . Recordemos que la

segunda variación de este funcional en un extremo viene dada por
∫ (Pv' )
x1
2
+Qv 2 dx (16)
x0
con
1 1 d 
P= Lu 'u ' , Q =  Luu − Lu 'u ' 
2 2 dx 
Definición 9. La ecuación de Euler

d
− (Pv ') + Qv = 0
dx
del funcional cuadrático (16) es llamada ecuación de Jacobi del funcional
original dado por (15).
Definición 10. Se dice que el punto x es conjugado de x 0 con respecto al

funcional (15) si es el conjugado de x 0 con respecto al funcional cuadrático (16).
A continuación enunciamos la condición necesaria de Jacobi:
Teorema 11. Si u es un mínimo del funcional

() = &, (), S ()'

(17)
para el cual F y ' y ' > 0 , entonces el intervalo (x0 , x1 ) no contiene puntos
conjugados de x 0 .
Demostración. Como ya se ha dicho anteriormente, una condición necesaria

para la realización de F es la no negatividad de la segunda variación evaluada
en dicho mínimo. Un colorarlo anterior, asegura que si el funcional cuadrático
(16) es no negativo entonces el intervalo ( x 0 , x1 ) no puede contener puntos
conjugados en x 0 .
■

5.2.Condiciones suficientes de Jacobi bajo las
cuales un funcional tiene un mínimo.
En esta sección formularemos las condiciones suficientes bajo las cuales un
funcional de la forma

() = &, (), S ()'

(18)
con condiciones de frontera fija, y ( x0 ) = y 0 y y ( x1 ) = y1 tiene un mínimo en la

curva y . Veremos que estas condiciones se parecen mucho a las condiciones
necesarias obtenidas en las secciones anteriores. Las condiciones necesarias
fueron consideradas separadamente, puesto que cada una de ellas es necesaria
en sí misma. Sin embargo, las condiciones suficientes deben considerarse en
conjunto puesto que la presencia de un mínimo está asegurada solo si se
satisfacen todas las condiciones simultáneamente. Demostraremos previamente
un lema que usaremos en la demostración del teorema de suficiencia.
Lema 12. Sea u una función diferenciable con y ( x 0 ) = y 0 , y ( x1 ) = y1 y

u ∈ C01 ([x0 , x1 ]) . Entonces si F es tres veces diferenciable con continuidad
respecto todos sus argumentos, se tiene
∫ ( ) ∫ (ξ u )
x1 x1
v(y + u ) − v(y ) = F Py ' 2 + Qy dx + 2
+ η u ' 2 dx
x0 x0
donde ξ ,η → 0 cuando u 1 → 0 .
Demostración. El desarrollo de Taylor de orden dos nos proporciona la

identidad
1
v(y + u ) − v(y ) = ∫ (
x1
x0
F u +yF u ' )dx +
2∫
(F
y'
x1
x0
yy )
u 2 + 2 F yy 'uu '+ F y '' y 'u ' 2 dx + ε
donde el resto, ε , puede escribirse como
∫ (ε u )
x1
2
ε = 1 + ε 2 uu '+ ε 3 u ' 2 dx (19)
x0
Debido a la continuidad de las derivadas F yy , F yy ' y F y ' y ' , se sigue que

ε i → 0 cuando u 1 → 0 , para i = 1,2,3 . Ahora integrando por partes y usando
las condiciones de frontera de h , podemos escribir (19) como
∫ (ξ u )
x1
2
+ η u ' 2 dx
x0
con ξ y η satisfaciendo ξ ,η → 0 cuando u 1 → 0 . ■

Teorema 13. Supongamos que el funcional (18) evaluado en u satisface las
siguientes condiciones:
1. La curva u es un extremo, es decir, satisface la ecuación de Euler
d
Fy (x, y ( x ), y ' (x )) − Fy ' ( x, y ( x ), y ' (x )) = 0 para todo x ∈ (x0 , x1 )
dx
2. A lo largo de la curva u se satisface la condición estricta de Legendre, es
decir
1
P(x ) = = Fy ' y ' (x, y ( x ), y ' ( x )) > 0 para todo x ∈ (x0 , x1 )
2
3. El intervalo [x 0 , x1 ] no contiene puntos conjugados de x 0 .
Entonces el funcional (18) tiene un mínimo en y .
Demostración. Si el intervalo [x 0 , x1 ] no contiene puntos conjugados de x 0 , y

si P ( x ) > 0 en dicho intervalo, entonces por la continuidad de la solución de la
Ecuación de Jacobi y de la función P , tenemos que tampoco habrá puntos
conjugados en un intervalo mayor [x 0 , x1 + ε ] , en el cual también podemos
asumir que P > 0 . Consideremos ahora el funcional cuadrático
∫ (Pu ' +Qu )dx − α ∫

x1 x1
2 2 2
u ' 2 dx (19)
x0 x0
que tiene por ecuación de Euler
d
−
dx
[( )]
P − α 2 u ' + Qu = 0 (20)
Puesto que P > 0 en [x 0 , x1 + ε ] y, por tanto, tiene una cota inferior positiva
en este intervalo, y como la solución de (20) que satisface las condiciones
iniciales u ( x 0 ) = 0 , u ' (x 0 ) = 1 depende con continuidad del parámetro α , se
sigue que:
1. P ( x ) − α 2 > 0 para todo x ∈ [x 0 , x1 ] .

2. La solución de (2.15) que satisface las condiciones de frontera u ( x 0 ) = 0 ,
u ' (x 0 ) = 1 no se anula en ( x 0 , x1 ] .
Por el segundo teorema de este punto se sabe que el funcional (19) es

definido positivo para todo α suficientemente pequeño. Es decir, existe una
constante c > 0 tal que

∫ (Pu ' )
x1 x1
2
+Qu 2 dx > c ∫ u ' 2 dx (21)
x0 x0
A partir de (21) deducimos que el mínimo es, efectivamente, y . En efecto,

sea h una curva tal que y + h está suficientemente próximo a y . Entonces
tenemos que (por un lema anterior):
∫ ( ) ∫ (ξ h )
x1 x1
v(y + h ) − v(y ) = F Pu ' 2 + Qu dx + 2
+ η h ' 2 dx
x0 x0
donde ξ ,η → 0 en [x 0 , x1 ] cuando h 1 → 0 . Además, usando la desigualdad de

Schwarz obtenemos
2
h 2 ( x ) =  ∫ h ' dx  ≤ ( x − x 0 )∫ h ' 2 dx ≤ ( x − x 0 )∫ h ' 2 dx
x x x1
 x0  x0 x0
es decir,
x1 (x1 − x 0 )2 x1
∫ ∫
2
h dx ≤ h ' 2 dx
x0 2 x0
que implica que
 ( x1 − x 0 )
2
 x1 2
∫ ( )
x1
ξ h + η h ' dx ≤ ε  1 +
2 2
 ∫ h ' dx
 x0 (22)
x0
 2 
si ξ ≤ ε y η ≤ ε . Puesto que ε > 0 puede tomarse arbitrariamente pequeño, se

sigue de (21) y (22) que v( y + h ) > v( y ) para toda h con h 1
suficientemente
pequeña.
■
5.3. Relación entre la condición de Jacobi y la

teoría de formas cuadráticas.
El funcional cuadrático
∫ (Pu ' )
b
2
+Qu 2 dx (23)
a
donde hemos cambiado la notación [x 0 , x1 ] por [a, b ] , y donde P ( x ) > 0 para

todo x ∈ [a, b ] , es definido positivo para todo u ∈ C 0 ([a, b]) si y solo si el
1
intervalo [a, b ] no contiene puntos conjugados de a . El funcional (23) es el

análogo a una forma cuadrática en dimensión finita. Por tanto, es natural

comenzar estudiando las condiciones de este tipo de formas en un espacio n-
dimensional y luego tomar el límite n → ∞ . Esto puede hacerse de la siguiente
manera: introducimos la partición
a = x0 , x1 ,..., xn , xn +1 = b
del intervalo [a, b ] , en la cual, por comodidad, suponemos los nodos

equiespaciados, ∆x = (b − a ) (n + 1) , a continuación consideramos la forma
cuadrática
n   u i +1 − u i  2 
∑ 
i =0 
Pi 
 ∆ x


+ Q u 2
i i  ∆x (24)
 
donde Pi , Qi y v i son los valores de las funciones P , Q y v en los nodos x i .

Esta forma cuadrática proporciona una aproximación finito dimensional
cuadrático (23), agrupando términos similares y teniendo en cuenta que
u 0 = u (a ) = 0 , u n +1 = u (b ) = 0 , podemos escribir (24) como
n
 Pi −1 − Pi  2 P 
∑  Q ∆x + i
∆x 
u i − 2 i −1 u i −1u i 
∆x
(25)
i =1  
En otras palabras, el funcional cuadrático (23) puede aproximarse por una

forma cuadrática de n variables cuya matriz viene dada por
 
 a1 b1 0 ... 0 0 0 
 
b a2 b2 ... 0 0 0 
 1 
 
0 b2 a3 ... 0 0 0 
 
  (26)
 .......... .......... ... 
 
0 0 0 ... bn − 2 a n −1 bn −1 
 
 0 0 0 ... 0 bn −1 a n 
 
donde
Pi −1 + Pi
ai = Qi ∆x + para i = 1,2,..., n
∆x
y
Pi
bi = − para i = 1,2,..., n − 1
∆x

Un matriz como (26) en la cual todos los electos excepto la diagonal
principal y sus dos diagonales adyacentes se anulan es llamada matriz de
Jacobi, y su forma cuadrática asociada es llamada forma de Jacobi. Para
cualquier matriz de Jacobi existe una fórmula de recurrencia para el cálculo de
los menores principales dados por:
para i = 1,2,..., n . En efecto, expandiendo Di con respecto a los elementos de la

última fila, obtenemos la fórmula
Di = ai Di −1 − bi2−1 Di − 2 (27)
que nos permite determinar D3 ,..., Dn en términos de los dos primeros

menores. De hecho, si tomamos D0 = 1 y D −1 = 0 , entonces esta fórmula es
válida para todo i = 1,2,..., n .
El criterio de Silvester asegura que una forma cuadrática simétrica es

definida positiva si y solo si todos los menores Di son positivos. Podemos así
obtener un criterio para que el funcional cuadrático (23) sea definido positivo
haciendo n → ∞ en la fórmula (27). Sustituyendo la expresión de los
coeficientes ai y bi en dicha fórmula, obtenemos
 Pi −1 + Pi  Pi −21
Di =  Q i ∆ x +  Di −1 − D (28)
 ∆x  (∆x )2 i − 2
para i = 1,2,..., n . Es imposible pasar directamente al límite n → ∞ en esta
expresión, puesto que los coeficientes de Di −1 y Di −2 se hacen infinitos. Para
evitar esta dificultad realizamos el cambio

P1 ...Pi Z i +1 Z1
Di = , D0 = = 1 , D−1 = Z 0 = 0
(∆x )i +1 ∆x
para i = 1,2,..., n . La fórmula (28) se escribe entonces, en términos de las

variables Z i como
P1 ...Pi Z i +1  Pi −1 + Pi  P1 ...Pi −1 Z i Pi −21 P1 ...Pi − 2 Z i −1

=  Q ∆x +  −
(∆x )i +1  i ∆x  (∆x )i (∆x )2 (∆x )i −1
es decir,
Qi Z i (∆x ) + Pi −1 Z i + Pi Z i − Pi Z i +1 − Pi −1 Z i −1 = 0
2
(29)
o
1  Z i +1 − Z i Z − Z i −1 
Qi Z i −  Pi − Pi −1 i =0
∆x  ∆x ∆x 
para i = 1,2,..., n . Pasando al límite n → ∞ obtenemos la ecuación diferencial
d
− (PZ ') + QZ = 0
dx
que es justamente la ecuación de Jacobi.
La condición de que las cantidades Di sean positivas es equivalente a la

condición de que las cantidades Z i que satisfacen la ecuación en diferencias
(29) sean positivas ya que el factor
Pi ...Pi
(∆x )i +1
es siempre positivo (ya que P ( x ) > 0 ). De modo que hemos probado que la
forma cuadrática (23) es definida positiva si y solo si todas, excepto las primeras
n + 2 cantidades Z 0 ,..., Z n +1 que satisfacen la ecuación en diferencias (29) son
positivas.
Ahora, si consideramos la línea poligonal Π n con vértices
(a0 , Z 0 ), (x1 , Z1 ),..., (b, Z n+1 )

la condición de que Z 0 = 0 y Z i > 0 para i = 1,2,..., n + 1 significa que Π n no
corta el intervalo [a, b ] excepto en el punto a . Así, cuando ∆x → 0 , la ecuación
en diferencias (29) se transforma en la ecuación de Jacobi, y la línea poligonal
Π n tienda a una solución no trivial de dicha ecuación, la cual satisface la
condición inicial
Z1 − Z 0 ∆x
Z (a ) = Z 0 = 0 , Z ' (a ) = lím = lím =1
x→0 ∆x x → 0 ∆x
y además, dicha solución no se anula en (a, b ] . En otras palabras, cuando

n → ∞ , la forma de Jacobi converge al funcional cuadrático (23), y la condición
de que (25) sea definida positiva se traduce en la condición de que (23) sea
definida positiva, que es equivalente a que el intervalo [a, b ] no contenga puntos
conjugados de a .
Ejemplo 8. ¿Se cumple la condición de Jacobi para calcular el extremo del

a
( )
funcional v = ∫ y ' − y dx que pasa por los puntos A(0,0 ) y B (a,0 ) ?
0
2 2
Solución. La ecuación de Jacobi tiene la forma
d
− 2y − (2 y ') = 0 , o bien y ' '+ y = 0
dx
de donde
y = C1 sen ( x − C 2 ) .
Como y (0 ) = 0 , entonces C 2 = 0 , y = C1 senx . La función y se anula en los

puntos x = kπ , donde k es un número entero.
Por tanto, si a ∈ (0, π ) , la función y se anula en el intervalo x ∈ [0, a ] solo en el

punto x = 0 , y la condición de Jacobi se cumple. Si, en cambio, a ∈ (π , ∞ ) , la
función y se anula en el intervalo x ∈ [0, a ] también por lo menos en el punto
x = π , y la condición de Jacobi no se cumple.
■
Ejemplo 9. ¿Se cumple la condición de Jacobi para calcular el extremo del

a
0
2
( 2 2
)
funcional v = ∫ y ' + y + x dx que pasa por los puntos A(0,0 ) y B (a,0 ) ?
Solución. La ecuación de Jacobi tiene la forma
y ' '+ y = 0
Tomemos la solución general en la forma
y = C1 shx + C 2 chx
De la condición y (0 ) = 0 se halla C 2 = 0 , por lo que y = C1 shx . Las curvas

del haz y = C1 shx cortan el eje OX solo en el punto x = 0 .
La condición de Jacobi se cumple para todo a .
5.4.Condición suficiente de Jacobi de inclusión

de una extremal en un campo de
extremales central.
Para que un arco AB de una extremal se pueda incluir en un campo de

extremales central con centro en un punto A( x 0, y 0) , es suficiente que el arco AB
no contenga puntos A * conjugados del punto A.
Ejemplo 10. Consideremos el funcional

a
v[ y( x)] = ∫ ( y '2 − 9 y 2 + e x −1)dx;
2
y (0) = 0, y (0) = 0, y (a ) = 0.
0
Determinar si la extremal y = 0 puede ser incluida en un campo central de

extremales con centro en el punto O (0, 0).
Solución. La ecuación de Euler del funcional dado es y ''+ 9 y = 0 y su solución

general, y ( x) = C1sen(3 x) + C2 cos(3 x).
Si a ≠ kπ / 3 (k es un número entero), entonces la extremal que satisface las
condiciones de contorno dadas es la recta y = 0. Considerando la familia
uniparamétrica de extremales y1 = C1sen(3 x), es fácil verificar que su C-

discriminante está formado por los puntos ( kπ / 3, 0). Por eso, si a < π / 3, en la
extremal y = 0 no existen puntos conjugados del punto (0, 0), por lo cual esta
extremal, evidentemente, puede ser incluida en el campo central de extremales
con centro en el punto O (0, 0). Pero si a ≥ π / 3, entonces en la extremal y = 0
hay al menos un punto conjugado del punto O (0, 0), es decir, la condición
suficiente de Jacobi no se cumple. En este caso las extremales y = C1sen(3 x) no
forman un campo.
■
Forma analítica de la condición de Jacobi.
Consideremos el problema variacional elemental

x1
v[ y ( x)] = ∫ F ( x, y, y ')dx; y ( x0 ) = y0 , y( x1 ) = y1.
x0
Si la solución u = u ( x ) de la ecuación deJacobi
d d
( Fyy − Fyy ' )u − ( Fy ' y ' , u ') = 0
dx dx
que satisface la condición u ( x0 ) = 0 se anula en algún punto del intervalo

x0 < x < x1 , entonces el arco AB de la extremal (en el punto B tiene coordenadas
( x1 , y1 ) ) contiene un punto A * conjugado de A( x0 , y0 ).
Si existe una solución u ( x) de la ecuación de Jacobi que satisface la
condición u ( x0 ) = 0 y que no se anula en ningún punto del semiintervalo
x0 < x ≤ x1 , entonces en el arco AB no hay puntos conjugados de A. En este caso,
el arco AB de la extremal puede ser incluido en el campo central de extremales
con centro en el punto A( x0 , y0 ).
Ejemplo 11. Determinar si se cumple la condición de Jacobi para la extremal

del funcional
a
v[ y( x)] = ∫ ( y '2 + x 2 )dx
0
que pasa por los puntos O (0, 0) y B ( a, 3).
Solución. La ecuación de Jacobi es u '' = 0 y su solución general es

u ( x) = C1 x + C2 . De la condición u (0) = 0 resulta C2 = 0, de donde u = C1 x. Estas
soluciones u = C1 x ( C1 ≠ 0 ) no se anulan para ningún valor de a > 0. Es decir, en
el arco OB de la extremal no existe un punto conjugado del punto O (0, 0). Por
consiguiente, esta extremal puede ser incluida en un campo central de
extremales con centro en el punto O (0, 0). No es difícil verificar que la extremal

buscada es la recta y = 3 x / a, la cual, evidentemente, está incluida en el campo
central de extremales y = C1 x.
■

6. Condición de Weiertrass.
6.1. Función E ( x, y , p , y ' ) . Condición de

Weiertrass.
Supongamos que en el problema simple sobre el extremo del funcional
x1
v = ∫ F ( x, y, y ' )dx; y ( x0 ) = y 0 ; y ( x1 ) = y1
x0
se cumple la condición de Jacobi y, por consiguiente, la extremal C que pasa por
inclinación igual a (, ) (Figura 5).

los puntos A( x0 , y 0 ) y B( x1 , y1 ) puede ser incluida en un campo central con
Figura 5
Transformemos el incremente ∆v = ∫ F ( x, y, y ' )dx − ∫ F ( x, y, y ' )dx a una forma

C C
funcional v al pasar de la extremal 1 a cierta curva próxima admisible C . (Los

más cómoda para su estudio, para determinar el signo del incremento ∆v del
símbolos:
∫ F ( x, y, y' )dx y ∫ F ( x, y, y ' )dx

C
C
x1
representan los valores del funcional v = ∫ F ( x, y, y ' )dx tomados

x0
respectivamente por los arcos de las curvas C y C ).
Consideremos el funcional auxiliar:

  dy  
∫  F ( x, y, p) +  dx − p  F p ( x, y, p )dx ,

C
dy
que se transforma en ∫ F ( x, y, y ' )dx
C
en la extremal C , en virtud de que
dx
= p
en las extremales del campo. Por otro lado, el mismo funcional auxiliar
  dy  
∫  F ( x, y, p) +  dx − p  F p ( x, y , p )dx ,

C
o bien
∫ [ F ( x, y, p) − pF p ( x, y , p )]dx + F p ( x, y , p )dy (30)

C
es la integral de una diferencial total. En efecto, la diferencial de la función

v ( x, y ) , en la cual se transforma el funcional v[ y ( x )] en las extremales del
campo, tiene la forma
d v = [ F ( x, y , y ' ) − y ' Fy ' ( x, y , y ' )]dx + Fy ' ( x, y , y ' )dy ,
Y se diferencia sólo por la notación del coeficiente angular de la tangente a las

extremales del campo de la expresión subintegral en la instalar auxiliar (30)
considerada.
De este modo, la integral ∫ [ F ( x, y, p) + ( y '− p) F p ]dx coincide con la

C
∫ F ( x, y, y ' )dx
C
en la extremal C, y como el funcional ∫ [ F ( x, y , p ) + ( y '− p ) F p ]dx es
C
la integral de una diferencial total, (y por lo tanto, no depende del camino de
integración) entonces
∫ F ( x, y, y' )dx = ∫ [ F ( x, y, p) + ( y'− p) F

C
p ( x, y , p )]dx
C
no sólo para C = C , sino también para cualquier C .
Por lo tanto, el incremento
∆v = ∫ F ( x, y , y ' )dx − ∫ F ( x, y , y ' )dx

C C
Puede ser reducido a la forma siguiente:
∆v = ∫ F ( x, y , y ' )dx − ∫ [ F ( x, y , p ) + ( y '− p ) F p ( x, y , p )]dx .

C C

La función subintegral se llama función de Weierstrass, y se denota por
E ( x, y , p , y ' ) :
E ( x, y , p, y ' ) = F ( x, y , y ' ) − F ( x, y , p ) − ( y '− p ) F p ( x, y, p ) .
En estas notaciones
x1
∆v = ∫ E ( x, y , p, y ' )dx
x0
Es evidente que una condición suficiente para que el funcional v tenga un

mínimo en la curva C es que la función E no sea negativa, puesto que si E ≥ 0 ,
entonces también ∆v ≥ 0 ; y una condición suficiente para que tenga un máximo
será E ≤ 0 , debido a que en este caso también ∆v ≤ 0 . Para que haya un mínimo
débil es suficiente que la desigualdad E ( x, y , p, y ' ) ≥ 0 (ó E ≤ 0 en el caso de
máximo) se cumpla para valores de x e y próximos al valor de x e y en la
extremal C investigada, y para valores de y’ cercanos a p(x,y) en la misma
extremal; para que haya un mínimo fuerte la misma desigualdad debe ser válida
para las mismas x e y, pero para y’ arbitrarias, puesto que en el caso de un
extremo fuerte las curvas cercanas tendrán direcciones arbitrarias de las
tangentes, y en el caso de extremo débil los valores de y’ en las curvas cercanas
son próximos a los valores y ' = p en la extremal C.
Por consiguiente, las siguientes condiciones serán suficientes para que el

funcional v tenga un extremo en la curva C.
Para un extremo débil:
1. La curva C es una extremal que satisface las condiciones de frontera.

2. La extremal C puede ser incluida en un campo de extremales. Esta
condición se puede sustituir de la de Jacobi.
3. La función E ( x, y , p, y ' ) no cambia su signo en todos los puntos (x, y)
próximos a la curva C, y para valores de y’ cercanos a p(x, y). En el
caso de mínimo es E ≥ 0 , en el caso de máximo, E ≤ 0 .
Para un extremo fuerte:
1. La curva C es una extremal que satisface las condiciones de frontera.

2. La extremal C puede ser incluida en un campo de extremales. Esta
condición se puede sustituir por la de Jacobi.
3. La función E ( x, y , p, y ' ) no cambia su signo en todos los puntos (x, y)
próximos a la curva C, y para valores de y’. En el caso de mínimo es
E ≥ 0 , en el caso de máximo, E ≤ 0 .
Observación. Se puede demostrar que la condición de Weierstrass es

necesaria. Más exactamente, si en un campo central que incluya a la extremal C
la función E tiene signos opuestos en los puntos de la extremal para ciertas y’,

entonces no hay extremo fuerte. Si esta propiedad tiene lugar para valores de y’
arbitrariamente próximos a p, tampoco hay extremo débil.
Ejemplo 12. Investigar el extremo del funcional
a
v = ∫ y '3 dx y (0) = 0 y(a) = b
0
a > 0 b > 0.
Solución. Las rectas y = C1 x + C 2 son extremales. El extremo puede alcanzarse

b
sólo en la recta y = x . El haz de rectas y = C1 x con centro en el punto (0, 0)
a
b
forma un campo central que incluye a la extremal y = x (Figura 6).
a
Figura 6
La función E es igual a
E ( x, y , p, y ' ) = y ' 3 − p 3 − 3 p 2 ( y '− p ) = ( y '− p ) 2 ( y '+2 p ) .
b b
En la extremal y = x la inclinación del campo es p = > 0 , y si y’ toma
a a
b
valores próximos a p = , entonces E ≥ 0 y, en consecuencia se cumplen todas
a
las condiciones suficientes para que haya un mínimo débil. De este modo, en la
b
extremal y = x hay un mínimo débil. Si, en cambio, y’ toma valores
a
arbitrarios, entonces ( y '−2 p ) puede tener cualquier signo y, por lo tanto, la
función E no tiene signo constante; no se cumplen las condiciones suficientes
para que haya un mínimo fuerte. Entonces, se puede afirmar que no hay
b
mínimo fuerte en la recta y = x .
a

Ejemplo 13. Analizar el extremo del funcional
∫ (6 y '
2
− y ' 4 + yy' )dx; y (0) = 0 y(a) = b
0
a > 0 b > 0.
en la clase de las funciones continuas con derivada primera continua.
b
Solución. Las rectas y = C1 x + C 2 son extremales. La recta y =
x satisface las
a
condiciones de frontera, y se incluye en el haz de extremales y = C1 x que
forman un campo central.
La función E es igual a
E ( x, y, p, y ' ) = 6 y ' 2 − y ' 4 + yy'−6 p 2 + p 4 − yp − ( y '− p )(12 p − 4 p 3 + y ) =

.
[
= −( y '− p) 2 y ' 2 +2 py '−(6 − 3 p 2 ) . ]
El signo de E es opuesto a signo del último factor
y ' 2 +2 py '−(6 − 3 p 2 ) .
Este factor se anula y puede cambiar su signo sólo cuando ′ pasa por el valor
y ' = − p ± 6 − 2 p 2 . Para 6 − 2p2 ≤ 0 ó p ≥ 3 , para toda y’ tendremos
[ y ' 2 +2 py '−(6 − 3 p 2 )] cambia su signo. Si en este caso ′ se diferencia

[ y ' 2 +2 py '−(6 − 3 p 2 )] ≥ 0 ; si, en cambio, 6 − 2 p 2 > 0 ó p < 3 , la expresión
suficientemente poco de , la última expresión mantiene su signo positivo >

1, y negativo para < 1.
b
Por tanto, para p = < 1 ó b < a se tiene un mínimo débil, puesto que E ≥ 0
a
b
para valores de y’ próximos a p; para p = > 1 ó b > a se tiene un máximo
a
b
débil. Para p = ≥ 3 se tiene un máximo fuerte, ya que E ≤ 0 para valores
a
b
cualesquiera de y’. Para p = < 3 , no hay ni mínimo fuerte ni máximo fuerte
a
(figura 7).

Figura 7
del signo de la función causó ciertas dificultades; por esto, sería conveniente
Incluso en los ejemplos citados más arriba, sumamente simples, el estudio
sustituir la condición de que la función tenga sigo constante por otra, de más
veces respecto al argumento ′. Según la fórmula de Taylor,

fácil comprobación. Supongamos que la función F ( x, y , y ' ) es derivable tres
( y '− p ) 2
F ( x, y , y ' ) = F ( x, y , p ) + ( y '− p ) F p ( x, y , p ) + Fy ' y ' ( x, y, q ) ,
2!
donde m está contenido entre e ′.
La función
E ( x, y , p, y ' ) = F ( x, y , y ' ) − F ( x, y , p ) − ( y '− p ) F p ( x, y, p )
luego de sustituir F ( x, y , y ' ) por su desarrollo mediante la fórmula de Taylor,

toma la forma:
( y '− p ) 2
E ( x , y , p, y ' ) = Fy ' y ' ( x, y , q ) .
2!
De aquí se ve que la función conserva su signo se F y ' y ' ( x, y , q ) lo conserva.
valores de e en los puntos cercanos a los puntos de la extremal que se

Al estudiar el extremo débil, la función F y ' y ' ( x, y , q ) debe conservar su signo para
investiga, y para valores de m próximos a . Si F y ' y ' ( x, y , y ' ) ≠ 0 en los puntos de

la extremal 1, entonces en virtud de la continuidad esta derivada segunda
conserva su signo en los puntos próximos a la curva 1, y para valores de ′
próximos a los valores de y’ en la curva 1. De este modo, al estudiar un mínimo
débil la condición E ≥ 0 puede ser sustituida por la F y ' y ' < 0 en la curva 1. La
condición F y ' y ' > 0 (ó F y ' y ' < 0 ) se llama condición de Legendre).

F y ' y ' ( x, y , q ) ≥ 0 en los puntos (, ) próximos a los puntos de la curva 1 para
Al estudiar un mínimo fuerte la condición E ≥ 0 puede sustituirse por la
valores arbitrarios de m. En este caso se supone, claro está, que el desarrollo

mediante la fórmula de Taylor:
( y '− p ) 2
F ( x, y , y ' ) = F ( x, y , p ) + ( y '− p ) F p ( x, y , p ) + Fy ' y ' ( x, y, q )
2!
tiene lugar para ′ cualesquiera. Al estudiar un máximo fuerte se obtiene la

condición F y ' y ' ( x, y , q ) ≤ 0 , bajo las mismas suposiciones respecto a la región de
variación de los argumentos y a la validez del desarrollo de la función F ( x, y , y ' )
por la fórmula de Taylor.
a
v[ y ( x)] = ∫ ( y ' 2 − y 2 )dx a>0 y (0) = 0 y(a) = 0
0
Solución. La ecuación de Euler tiene la forma y ' '+ y = 0; su solución general es
y C 2 = 0 , si a ≠ kπ , donde : es un entero.
y = C1 cos( x) + C 2 sen( x). Utilizando las condiciones de frontera, se obtiene C1 = 0
De este modo, si a ≠ kπ , es extremo se puede alcanzar sólo en la recta y = 0 .

Si a < π , el haz de extremales y = C1 sen( x) con centro en el punto (0, 0) forma
un campo central. Si a > π , la condición de Jacobi no se cumple.
Como la función subintegral es derivable tres veces con respecto a ′ para ′

cualesquiera y Fy ' y ' = 2 > 0 para valores cualesquiera de y’, entonces en la recta
y = 0 , para a < π , hay un mínimo fuerte. Entonces, se puede afirmar que para
a > π no hay mínimo en la recta y = 0 .
x1
1 + y'2
v[ y ( x)] = ∫ dx y (0) = 0 y ( x1 ) = y1
0 y
Solución. Las extremales son las cicloides
x = C1 (t − sen(t ) + C 2
y = C1 (1 − cos(t )

El haz de cicloides x = C1 (t − sen(t ) + C 2 , y = C1 (1 − cos(t ) con centro en el
punto (0, 0) forma un campo central que incluye la extremal
x = a (t − sen (t )) y = a (1 − cos(t )) ,
Donde a se determina de la condición de que la cicloide pase por el segundo

punto frontera B( x1 , y1 ) (Figura 8), si x1 < 2πa
Figura 8
se tiene
y' 1
Fy ' = Fy ' y ' = >0
y 1 + y' 2
y (1 + y ' 2 ) 3 / 2
para y’ cualesquiera. Por consiguiente, para x1 < 2πa hay un mínimo fuerte en la
cicloide
x = a (t − sen (t )) y = a (1 − cos(t )) ,
a
v[ y ( x)] = ∫ y '3 dx y (0) = 0 y(a) = b a>0 b > 0.
0
Solución. Los extremales son líneas rectas. El haz y = Cx forma un campo

b b
central que incluye la extremal y = x . En la extremal y = x la derivada
a a
b b
segunda es Fy ' y ' = 6 y ' = 6 > 0 . En consecuencia, en la recta y = x hay un
a a
mínimo débil.

La derivada segunda Fy ' y ' = 6 y ' no mantiene su signo constante para ′
arbitrarias; por lo tanto, las condiciones suficientes para que haya un mínimo
fuerte indicadas más arriba no se cumplen. Sin embargo, de aquí no se puede
aún concluir que no hay un extremo fuerte.
a
y
v[ y ( x)] = ∫ dx y ( 0) = 1 y(a) = b a>0 0 < b < 1.
0 y' 2
Solución. La primera integral de la ecuación de Euler tiene la forma
y 2y
2
+ y ' 3 = C , o bien y ' 2 = 4C1 y
y' y'
Extrayendo la raíz cuadrada, separando variables e integrando, se obtiene

y = (C1 x + C 2 ) 2 , que es una familia de parábolas. De la condición y (0) = 1 se
halla C 2 = 1 . El haz de parábolas y = (C1 x + 1) 2 con centro en el punto A(0,1)
tiene la curva C1 -discriminante y = 0 (Figura 9).
Figura 9
Por el punto (, ) pasan dos parábolas de este haz. En el arco de una de
ellas ( L1 ) se encuentra el punto ∗ conjugado del ; en el otro ( L2 ) no hay
puntos conjugados y, por lo tanto, la condición de Jacobi se cumple en el arco
L2 , y en este arco de parábola puede haber un extremo. En un entrono de la
extremal estudiada se tiene Fy ' y ' = 4 > 0 para ′ arbitrarias; sin embargo, en
6y
y'

base a esto no su puede aún afirmar que en el arco L2 hay un mínimo fuerte,
y
debido a que la función F ( x, y, y ' ) = 2 no puede ser representada en la forma
y'
( y '− p ) 2
F ( x, y, y ' ) = F ( x, y, p ) + ( y '− p ) F p ( x, y, p ) − F y ' y ' ( x, y , q )
2!
Para valores de ′ arbitrarios, ya que la función F ( x, y , y ' ) es discontinua para
para valores de ′ cercanos a la inclinación del campo en la curva L2 tiene lugar

y ' = 0 . Se puede afirmar solamente que en L2 hay un mínimo débil, puesto que
dicho desarrollo de la función F ( x, y , y ' ) por la fórmula de Taylor.
Para un estudio completo del extremo de este funcional es necesario

considerar la función E ( x, y , p, y ' ) :
y y 2y y ( y '− p ) 2 (2 y '+ p )
E ( x, y, p, y ' ) = 2 − 2 + 3 ( y '− p ) = .
y' p p y' 2 p 3
Como el factor ( 2 y '+ p ) no mantiene su signo constante para ′ arbitrarias, se

puede afirmar, que no hay un mínimo fuerte en el arco L2 .
La teoría expuesta se generaliza sin cambios sustanciales a los funcionales de

la forma
x1
v[ y1 , y 2 ,..., y n ] = ∫ F ( x, y1 , y 2 ,..., y n , y '1 , y ' 2 ,..., y ' n )dx;

x0
yi ( x0 ) = y io y i ( x1 ) = yi1 (i = 1,2,...n)
La función tiene la forma:
E = F ( x, y1 , y 2 ,..., y n , y '1 , y ' 2 ,..., y ' n ) − F ( x, y1 , y 2 ,..., y n , p1 , p 2 ,..., p n ) −

n
− ∑ ( y i − pi ) F pi ( x, y1 , y 2 ,..., y n , p1 , p 2 ,..., p n )
i =1
donde p i son las funciones de inclinación del campo, al cual se le imponen

ciertas limitaciones (bajo estas limitaciones en campo se llama especial).
La condición de Legendre Fy ' y ' ≥ 0 se sustituye por las condiciones siguientes:
Fy '1 y '1 Fy '1 y '2

F y '1 y ' 2 ≥ 0 , ≥0
Fy '2 y '1 Fy '2 y '2

Las condiciones suficientes de mínimo débil se pueden obtener tanto para el
problema simple como para otros más complejos por otro método, basado en el
estudio del signo de la segunda variación.
El incremento del funcional en el problema simple se puede reducir,

mediante la fórmula de Taylor, a la forma siguiente:
x1
∆v = ∫ [F ( x, y + δy, y'+δy' ) − F ( x, y, y' )]dx =

x0
x1 x
1 1
= ∫ ( Fy δy, Fy 'δy ' )dx + ∫ ( Fyy δy 2 + 2 Fyy 'δyδy '+ Fy ' y 'δy ' 2 )dx + R
x0
2 x0
donde R tiene orden mayor que 2 respecto a δy y δy ' . Al investigar un extremo

débil δy y δy ' son suficientemente pequeñas, y en este caso el signo del
incremento ∆v se determina por el signo del término que se halla en el segundo
miembro y contiene las potencias menores de δy y δy ' . En las extremales la
primera variación será
x1
∫ ( F δy, F δy ' )dx = 0

x0
y y'
y, por lo tanto, el signo del incremento ∆v coincide en general con el de la

segunda variación
x1
δ 2 v = ∫ ( Fyyδy 2 + 2 Fyy 'δyδy '+ Fy ' y 'δy ' 2 )dx

x0
La condición de Legendre conjuntamente con la de Jacobi son precisamente

condiciones que aseguran que el signo de la segunda variación sea constante y,
conjuntamente con esto, sea constante el signo del incremento ∆v en el
problema sobre un extremo débil.
En efecto, consideremos la integral
x1
∫ (ω ' ( x)δy
2
+ 2ω ( x)δyδy ' )dx (31)
x0
donde ω (x) es una función derivable arbitraria. Esta integral es igual a cero:
x1 x1
∫ (ω ' ( x)δy + 2ω ( x)δyδy ' )dx = ∫ d (ωδy ) = [ω ( x)δy 2 ] xx10 = 0

2
x0 x0
Sumando la integral (31) a la segunda variación, se obtiene

x1
δ v = ∫ (( Fyy + ω ' )δy 2 + 2( Fyy ' + ω )δyδy '+ Fy ' y 'δy ' 2 )dx .
2
x0
La función ω (x) se escoge de modo que la función subintegral se transforme,

salvo un factor, en un cuadrado perfecto. Para esto la función ω (x) debe
satisfacer la ecuación
Fy ' y ' ( Fyy + ω ' ) − ( Fyy ' + ω ) 2 = 0 .
Con la función ω elegida de este modo, la segunda variación toma la forma
2
x1
 Fyy ' + ω 
δ v = ∫ Fy ' y '  δy '+
2
δy  dx
x0  F y ' y ' 
y, en consecuencia, su signo coincide con el de F y ' y ' .
Sin embargo, esta transformación es posible sólo bajo la hipótesis de que la

ecuación diferencial
Fy ' y ' (ω '+ Fyy ) − ( Fyy ' + ω ) 2 = 0
tenga una solución derivable ω (x) en el segmento ( x0 , x1 ) .
Transformando esta ecuación a nuevas variables mediante la sustitución
u'
ω = − Fyy ' − Fy ' y '
u
donde es una nueva función incógnita, se obtiene
 d  d
 Fyy − Fyy ' u − ( Fy ' y ' u ' ) = 0 ,
 dx  dx
que es la ecuación de Jacobi.
Si existe una solución de esta ecuación que no se anule para x 0 < x ≤ x1 , es

decir, si se cumple la condición de Jacobi, existe, para dichos valores de x, una
solución continua y derivable:
u'
ω ( x) = − Fyy ' − Fy ' y '
u
de la ecuación:
Fy ' y ' ( Fyy + ω ' ) − ( Fyy ' + ω ) 2 = 0

De este modo, las condiciones de Legendre y de Jacobi aseguran la
constancia del signo de la segunda variación y, por consiguiente, son
condiciones suficientes de mínimo ( Fy ' y ' > 0) o máximo ( Fy ' y ' < 0) débil.

7. ANEXO: Métodos directos en el cálculo
variacional.
7.1. Método de Euler de diferencias finitas.

Consideremos el siguiente problema variacional elemental: hallar el
extremo del funcional
b
v[ y ( x)] = ∫ F ( x, y, y ')dx; y(a) = A, y(b) = B.
a
En el método de Euler los valores del funcional anterior no se consideran

en curvas arbitrarias admisibles para el problema variacional dado, sino sólo en
líneas quebradas formadas por un número dado n de segmentos rectilíneos
cuyos vértices tienen abscisas
b−a
a + ∆x, a + 2∆x,…, a + (n − 1)∆x, donde ∆x = .
n
En estas líneas quebradas el funcional v[ y ( x)] se transforma en una

función Φ ( y1, y2 ,..., yn −1 ) de las coordenadas y1 , y2 ,...., yn −1 de los vértices de la
línea quebrada. Las ordenadas y1 , y2 ,...., yn −1 se eligen de tal manera que la
función Φ ( y1, y2 ,..., yn −1 ) alcance su extremo, es decir, se determinan del
sistema de ecuaciones
∂Φ ∂Φ ∂Φ
= 0, = 0, ....., = 0.
∂y1 ∂y2 ∂yn−1
La línea quebrada obtenida es una solución aproximada del problema

variacional inicial.
Ejemplo 18. Hallar la solución aproximada del problema del mínimo del
funcional
1
v[ y ( x )] = ∫ ( y '2 + 2 y )dx; y (0) = y (1) = 0.
0
1− 0
Solución. Tomemos ∆x = = 0,2 y hagamos
5
= (0) = 0, = (0,2), 4 = (0,4), K = (0,6),

_ = (0,8), = (1) = 0
Sustituyamos aproximadamente los valores de las derivadas mediante la

fórmula
−
S = S ( ) ≈
∆
Entonces,
− 0 4 − K − 4 _ − K
S (0) = , S (0,2) = , S (0,4) = , S (0,6) = ,
0,2 0,2 0,2 0,2
0 − _
S (0,8) = .
0,2
La integral se calcula por el método de los rectángulos:

() ≈ (() + ( ) + (4 ) + ⋯ + (;^ ))∆.

Entonces,
4 4 − 4 K − 4 4 _ − K 4 ,
Φ( , 4 , K , _ ) = `) * + ) * + 2 + ) * + 24 + ) * +
0,2 0,2 0,2 0,2
,+2K + )− _ * + 2_ a ∙ 0.2

4
0,2
Escribimos el sistema de ecuaciones para determinar las ordenadas , 4 , K , _

de la línea quebrada:
1 "Φ 4 −
∙ = − + 2 = 0,
0,2 " 0,02 0,02
1 "Φ 4 − K − 4
∙ = − + 2 = 0,
0,2 "4 0,02 0,02 ,
1 "Φ − −
∙ = − + 2 = 0,
K 4 _ K
0,2 "K 0,02 0,02
1 "Φ − _
∙ =
_ K
+ + 2 = 0,
0,2 "_ 0,02 0,02
o bien
2 − 4 = −0,04
− + 24 − K = −0,04,

−4 + 2K − _ = −0.04
−K + 2_ = −0,04

La solución de este sistema es
= −0,08 ; 4 = −0,12 ; K = −0,12 ; _ = −0,08
Los valores de la solución exacta
4 −
=
2
en los puntos respectivos coinciden con los valores de la solución aproximada.
7.2. Método de Ritz.

La idea del método consiste en que al hallar el extremo del funcional v[ y ( x)]
se consideran, en lugar del espacio de las funciones admisibles, sólo las
funciones que se pueden representar como combinaciones lineales de las
funciones admisibles:
n
y n ( x) = ∑ α j ϕ j ( x) (32)
j =1
donde α j son unas constantes y el sistema {ϕ j (x)}, llamado sistema de

funciones coordenadas, está formado por funciones ϕ j (x) que son linealmente
independientes y que constituyen un sistema completo de funciones en el
espacio considerado.
Hablado en términos generales, cuando pedimos que las funciones y n (x)

sean admisibles, imponemos a las funciones coordenadas ϕ j (x) ciertas
condiciones complementarias como, por ejemplo, limitaciones en cuanto a la
derivabilidad o en cuanto a la verificación de las condiciones de frontera.
En estas combinaciones lineales el funcional v[ y ( x)] se convierte en una

función de los argumentos α 1 , α 2 ,..., α n :
v[ y ( x)] = Φ (α1 ,α 2 ,...,α n ) .
Determinamos los valores α1 , α 2 ,..., α n que ofrecen extremo a la función

Φ (α 1 , α 2 ,..., α n ) ; para ello resolvemos el sistema de ecuaciones

∂Φ
=0 (i = 1,2,..., n) ,
∂α j
no lineales, como regla, respecto a α1 , α 2 ,..., α n , e introducimos en (32) los

valores encontrados para α j . La sucesión {y n (x)} que así resulta es una
sucesión minimizante, o sea, la sucesión de los valores del funcional {v[ y ( x)]}
obtenida a partir de ella converge hacia el mínimo o hacia la cota inferior del
funcional J [ y ( x )] .
Sin embargo, de
limv[ yn ( x)] = min v[ y ( x )]

n→∞
no se deduce aún que lim yn ( x ) = y ( x ) . La sucesión minimizante puede no

n→∞
converger hacia la función que realiza el extremo en la clase de las funciones
admisibles.
Se pueden indicar las condiciones que garanticen que el mínimo absoluto del
funcional exista y se alcance en las funciones {y n (x)}.
En el caso en el que se trate del extremo del funcional
x2
v[ y ( x)] = ∫ F ( x, y, y ')dx ;
x1
y ( x1 ) = y1 y( x2 ) = y 2
estas condiciones son:
1. La función F ( x, y, z ) es continua respecto al conjunto de sus argumentos

para cualquier z y para ( x, y ) ∈ D , donde D es un recinto cerrado del
plano x0 y al que pertenecen las líneas y n (x) .
2. Existan unas constantes α > 0, p > 1 y β tales que
p
F ( x, y , z ) ≥ z + β
cualquiera que sea z y para cualquier punto ( x, y ) ∈ D .
3. La función F ( x, y, z ) tiene la derivada parcial continua Fz ( x, y, z ) y esta

derivada es una función no decreciente de z ( −∞ < z < +∞ ) cualquiera que
sea el punto ( x, y ) ∈ D .
En particular, las condiciones enunciadas se cumplen para los funcionales

x2
v[ y ( x)] = ∫  p ( x) y '2 + q ( x) y 2 + 2r ( x) y  dx
x1
y( x1 ) = a y ( x2 ) = b
donde p ( x), q ( x) y r (x ) son funciones dadas, continuas en [ x1 , x2 ] , con la

particularidad de que existe la derivada continua p ' ( x ) de p (x ) y de que
p( x) > 0 y q( x) ≥ 0 .
Si por este método se determina el extremo absoluto del funcional, el valor

aproximado de su mínimo se obtiene por exceso y el valor aproximado de su
máximo, por defecto. Al aplicar este método, el éxito depende en gran medida
de la elección adecuada del sistema {ϕ j (x)} de funciones coordenadas.
En muchos casos basta tomar la combinación lineal de dos o tres funciones

ϕ j (x) para obtener una aproximación bastante satisfactoria de la solución
exacta.
Si hay que determinar el extremo aproximado del funcional v[ z ( x1 , x2 ,..., xn )]

que dependen de las funciones de varias variables independientes, se escoge un
sistema de funciones coordenadas
ϕ1 ( x1 , x 2 ,..., x n ) , ϕ 2 ( x1 , x 2 ,..., x n ) ,… ϕ n ( x1 , x 2 ,..., x n ) ,…
y la solución aproximada del problema variacional se busca en la forma:
m
z m ( x1 , x 2 ,..., x n ) = ∑ α k ϕ k ( x1 , x 2 ,..., x n ) ,
k =1
donde los coeficientes α k son unos números constantes. Para determinarlos se

forma, por analogía con lo que hemos explicado, el sistema de ecuaciones
∂Φ
= 0 ( k = 1,2,..., n) donde Φ (α 1 , α 2 ,..., α n ) es el resultado de introducir z m en
∂α k
el funcional v[ z ( x1 , x2 ,..., xn )] .
Ejemplo 19. Hallar la solución aproximada del problema sobre el mínimo del
funcional
1
v[ y ( x)] = ∫  y '2 − y 2 + 2 xy  dx (33)
0
y (0) = y (1) = 0
y compararla con la solución exacta.

Solución. Como sistema de funciones coordenadas ϕ k (x) tomamos
ϕ k ( x ) = (1 − x ) x k ( k = 1,2,...)
Es evidente que las funciones ϕ k (x) satisfacen las condiciones de frontera

ϕ k (0) = ϕ k (1) = 0 , son linealmente independientes y forman un sistema completo
en el espacio C1[0,1] .
Para k = 1 tenemos y1 ( x) = α1 ( x − x 2 ) . Introduciendo esta expresión de y1 ( x)

en el funcional (33), obtenemos
1
v[ y ( x)] = ∫ α12 (1 − 2 x) 2 − α12 ( x − x 2 ) 2 + 2 xα1 ( x − x 2 )  dx =
0
1
= ∫ α12 (1 − 4 x + 4 x 2 − x 2 + 2 x3 − x 4 ) + 2α1 ( x 2 − x3 )  dx =
0
3 2 1
= α1 + α1
10 6
El coeficiente α1 se determina de la ecuación
∂φ 3 1
= α1 + = 0
∂α1 5 6
5
de donde resulta α1 = − . Por consiguiente,
18
5 5
y1 ( x ) = − x + x2 .
18 18
Solución exacta. La ecuación de Euler del funcional considerado es
y ' '+ y = x
Resolviendo esta ecuación lineal no homogénea, encontramos
y = C1 cos( x) − C2 sen( x) + x .
Empleando las condiciones de frontera y (0) = y (1) = 0 , obtenemos

definitivamente

sen( x )
y = x− .
sen(1)
Comparemos las soluciones exacta y aproximada:
X Solución exacta Solución

aproximada
0,00 0 0
0,25 -0,044 -0,052
0,50 -0,070 -0,069
0,75 -0,060 -0,052
1,00 0 0
Ejemplo 20. Hallar la solución aproximada de la ecuación no lineal
3 2
y' ' = y
2
que satisfaga las condiciones y (0) = 4 , y (1) = 1 .
Solución. A este problema de contorno le corresponde el problema variacional
1
v[ y ( x)] = ∫  y '2 − y 3  dx y (0) = 4 y (1) = 1
0
Buscaremos la solución en la forma
y1 ( x ) = 4 − 3 x + α1 ( x − x 2 )
Es evidente que y1 ( x) satisface las condiciones de frontera dadas cualquiera que

sea el valor de α1 .
Tenemos:
1
v[ y ( x)] = ∫ α12 (1 − 2 x) − 3]2 + [4 − 3x + α1 ( x − x 2 )]3  dx ,
0
de donde
1
∂v[ y1 ( x)]
= ∫ (1 − 2 x)2[α1 (1 − 2 x) − 3] + 3( x − x 2 )[4 − 3x + α1 ( x − x 2 )]2  dx
∂α1 0

∂v[ y1 ( x)]
La condición = 0 toma la forma
∂α1
9α12 + 490α1 + 1407 = 0
y para α1 = −3,0413 obtenemos la solución del problema
y1 ( x ) = 3,0413 x 2 − 6,0413 x + 4
positiva en todos los puntos.
7.3. Método de Kantoróvich.

Este método ocupa un lugar intermedio entre la resolución exacta del
problema y el método de Ritz, y se utiliza en la investigación de los extremos de
funcionales
v[ z ( x1 , x2 ,...xn )],
dependientes de funciones de varias variables independientes ( n ≥ 2 ). Al igual

que en el caso del método de Ritz, elegimos un sistema coordenado de funciones
{ϕ ( x , x ,...x )} y buscamos una solución aproximada de la forma
k 1 2 n
m
zm = ∑ α k ( x j )ϕk ( x1 , x2 ,...xn ) ,
k =1
pero ahora los coeficientes α k ( x j ) son funciones desconocidas de una de las

variables independientes.
m
En las funciones de la forma zm = ∑α
k =1
k ( x j )ϕ k ( x1 , x2 ,...xn ) el funcional
v[ z ( x1 , x2 ,...xn )] se transforma en el funcional
v (α1 ( x j ),α 2 ( x j ),...,α m ( x j )),
dependiente de j funciones α1 ( x j ),α 2 ( x j ),...,α m ( x j ). Estas funciones se eligen

de tal manera que el funcional v alcance su extremo, y se determinan a partir de
las condiciones necesarias de extremo del funcional v .

{ϕ ( x , x ,...x )} y con el mismo
número j de términos de la aproximación, utilizando el método de Kantoróvich
Para las mismas funciones coordenadas k 1 2 n
obtenemos una solución aproximada más exacta, en general, que por el método
de Ritz.
Ejemplo 21. Hallar la solución aproximada de la ecuación de Poisson ∆z = −1

en el rectángulo D = {−a ≤ x ≤ a, −b ≤ y ≤ b} , con la condición z = 0 en la
frontera de la región D.
Solución: La ecuación ∆z = −1 es la ecuación de Euler-Ostrogradski del

funcional
 ∂z 2  ∂z 2 
v[ y ] = ∫∫   +   − 2 z dxdy.
   
D 
∂x ∂y 
La solución se busca en la forma
z1 ( x, y ) = (b 2 − y 2 )α ( x);
la función z1 ( x, y ) satisface las condiciones de contorno del problema en las

rectas y = ±b.
Sustituyendo este valor de z1 en el funcional v[ y ] anterior, obtenemos
a
 16 8 8 
v[ z1 ] = ∫  b5α '2 + b3α 2 − b3α  dx.
−a  
15 3 3
La ecuación de Euler de este funcional es
5
α ''− α = − 52
2b 2 4b
Se trata de una ecuación lineal homogénea con coeficientes constantes y su

solución general es
5 x 5 x 1
α ( x) = C1ch ⋅ + C2 sh ⋅ + .
2 b 2 b 2
Las constantes C1 , C2 se hallan a partir de las condiciones de contorno
α (− a ) = α (a ) = 0,
1
lo cual proporciona C2 = 0, C1 = − , de manera que
5 a
2ch ⋅
2 b

 5 x
⋅ 
1
ch
α ( x) =  1 − 2 b .
2 5 a
 ch ⋅ 
 2 b
Así pues,
 5 x
⋅ 
b −y 
2 2ch
z1 ( x, y ) = 1 − 2 b .
2  5 a
 ch ⋅ 
 2 b
Para obtener una aproximación más exacta se puede buscar la solución del
problema en la forma
z2 ( x, y ) = (b 2 − y 2 )α1 ( x) + (b 2 − y 2 ) 2 α 2 ( x).
■
7.4. Métodos variacionales de búsqueda de

valores propios y funciones propias
7.4.1 El problema de Sturm-Liouville.
La ecuación de Sturm-Liouville

− (() S ) + m() = ¡ (33)

con condiciones
() = 0, () = 0 (34)
donde () > 0 tiene derivada continua y m() es continua, siempre tiene la
solución trivial (nula) ≡ 0 para todo ¡ real o complejo.
Definición 14. La ecuación (33) con las condiciones de contorno (34) se
parámetro ¡ para los cuales el problema de contorno (33)-(34) tiene soluciones

denomina problema de contorno de Sturm-Liouville (33)-(34). Los valores del
no triviales () ≠ 0, se denominan valores propios y las propias soluciones,

funciones propias del problema de contorno dado.
La ecuación (33) es la ecuación de Euler del siguiente problema de

extremo condicionado: hallar el mínimo del funcional

= ( S + m 4 ) (35)
4

con las condiciones (34) y

()4 = 1 (36)

Si cierta función () es solución de este problema variacional, entonces,

en virtud de la condición (36), también es una solución no trivial del problema
(33)-(34). Por eso, los valores propios y las funciones propias del problema de
contorno de Sturm-Liouville se llaman también valores propios y funciones
propias del funcional (35) con las condiciones (34) y (36).
La función propia () se denomina función normalizada si

()4 = 1

Ejemplo 22. Hallar los valores propios y las funciones propias del funcional
K

=
(2 + 3)4 S − 4
4
con las condiciones

K
(0) = 0, (3) = 0, ()4 = 1.

Solución. La ecuación de Sturm-Liouville es

− − ((2 + 3)4 ′) = ¡,

o bien,
(2 + 3)4 SS + 4(2 + 3) S + (¡ + 1) = 0. (37)
Mediante el cambio de variables 2 + 3 = J { la ecuación anterior se reduce a la

ecuación lineal con coeficientes constantes
4
4 +4 + (¡ + 1) = 0. (38)
} 4 }
Su ecuación característica
4: 4 + 4: + ¡ + 1 = 0
tiene raíces

1 1
:,4 = − ± √−¡.
2 2
Consideremos tres casos:
1. ¡ < 0. Entonces la solución general de la ecuación (38) es
(}) = 1 J { + 14 J f { ,
donde : y :4 son números reales. La solucioón general de la ecuación

(37) es
() = 1 (2 + 3) + 14 (2 + 3)f .
1 = 0 14 = 0. Luego () ≡ 0.
Utilizando las condiciones de contorno del problema obtenemos que
2. ¡ = 0. Entonces,
{
(}) = (1 + 14 })J ^4 ,
es decir,
1
() = (1 + 14 ln(2 + 3)) .
√2 + 3
De las condiciones de contorno obtenemos que 1 = 0 14 = 0. Luego

() ≡ 0.
¡ > 0. Entonces :,4 = − 4 ± Y

√©
4
3. y la solución general de la
ecuación (38) es
{ √¡ √¡
(}) = J ^4 d1 cos } + 14 sin } e.
2 2
Regeresando a la variable , resulta
√¡ √¡
1 cos d 2 ln(2 + 3)e + 14 sin d 2 ln(2 + 3)e
() = . (39)
√2 + 3
Las condiciones de contorno proporcionan

1 cos d√¡ ln(3)e + 14 sin d√¡ ln(3)e = 0
2 2 , (40)
√¡ √¡
1 cos d ln(9)e + 14 sin d ln(9)e = 0
2 2
El sistema tiene soluciones no triviales si su determinante es igual a cero:
√¡ √¡
cos d ln(3)e sin d ln(3)e
¬ 2 2 ¬
= 0,
¬ √¡ √¡ ¬
cos d ln(9)e sin d ln(9)e
2 2
o bien
√¡
sin d√¡ ln 3 − ln 3e = 0,
2
es decir, sin 4 ln 3 = 0, de donde ln 3 = =B. Los valores propios son

√© √©
4
4=4 B 4
¡; = , = = 1,2, ….
ln4 3
¡ por ¡; obtenemos
Tomando cualquier ecuación del sistema (40) y sustituyendo el valor de
1 cos =B + 14 sin =B = 0,
de donde se tiene que 1 = 0. Haciendo en (39) 1 = 0 y ¡; = ,

_;f C f
®f K
obtenemos las funciones propias del problema dado:
=B ln(2 + 3)
sin ) *
ln 3
; () = 1; , = = 1,2, …
√2 + 3
Los coeficientes 1; se hallan de la condición de normalización

K
; ()4 = 1,

lo cual proporciona
2
1; = ± .
√ln 3

Entonces,
=B ln(2 + 3)
2 sin ) ln 3
*
; () = ± , = = 1,2, …
√ln 3 √2 + 3
Los valores propios y las funciones propias del problema variacional (35),
(34), (36) poseen una serie de propiedades importantes.
Si ¡h y ¡; son dos valores propios diferentes del funcional (35) con las
condiciones (34) y (36), e h () e ; () son las correspondientes funciones
1.
propias, entonces estas funciones son ortogonales, es decir,

h () ; () = 0, j ≠ =.

Todos los valores propios ¡; del funcional (35) son reales.

Si ¡; es un valor propio del funcional (35), e ; () es su función propia
2.
3.
normalizada correspondiente, entonces

; () = ¡; .
4. El menor de los valores propios coincide con el mínimo del funcional (35)
con las condiciones (34) y (36).
7.4.2 Principio de Rayleigh.
Consideremos el problema de valores propios

= − )() * + m() = ¡|(), (41)

() + ° S () = 0, 4 + °4 > 0,

¯ (42),
4 () + °4 S () = 0, 44 + °44 > 0,
donde (), S (), m(), |() son continuas en

, y () > 0 en
, .
Definición 15. La función () se denomina función admisible ( ∈ v) si es

dos veces diferenciable y satisface las condiciones de contorno (42).
Supongamos que para toda función admisible () se cumple la condición

≥ 0.

¡.
En este caso el problema de contorno (41)-(42) tiene sólo valores propios reales
siguiente problema variacional: entre todas las funciones admisibles () tales
Al problema de valores propios se le puede poner en correspondencia el
que

|() 4 > 0, (43)

hallar la función para la cual

= jY=.
|() 4

Sea = ² () la solución de este problema.
Si ¡ es el valor propio mínimo, es decir, si
´
² ´
²

¡ = min = ,
|() 4 | ²
!∈l 4
entonces ¡ es el menor valor propio positivo, y ² () es su función propia

correspondiente.
Si, además de las condiciones (43), sobre las funciones admisibles se

impone la condición

|² = 0

(condición de ortogonalidad), entonces el problema

= jY=
| 4

tiene nuevamente cierta solución ²4 ().
Si ¡4 es el valor propio mínimo correspondiente, entonces ¡4 es el

siguiente valor propio en magnitud (¡4 ≥ ¡ ), y ²4 () es su función propia
correspondiente, ortogonal a ² (). En general, si ya se conocen los : primeros
valores propios positivos

¡ ≤ ¡4 ≤ ⋯ ≤ ¡
y su sistema ortogonal de funciones correspondientes
² (), ²4 (), … , ² (),
entonces el siguiente valor propio es
´

¡ = min ,
| 4
!∈l
con la particularidad de que ahora se consideran las funciones admisibles ()

que, a excepción de (43), cumplen las condiciones adicionales

|()²] ()() = 0, = 1,2, … , :.

Si en la ecuación (41) la función |() > 0 en

, , entonces para estimar
superiormente el menor valor propio positivo ¡ , a menudo se utiliza la
desigualdad (principio de Rayleigh)

¡ ≤ .
| 4

Ejemplo 23. Utilizando el principio de Rayleigh, estimar ¡ en el problema de

contorno
− SS = ¡; S (0) = 0, (1) = 0.
Solución. En este caso

= − SS , es decir, () ≡ 1 > 0, m() ≡ 0 y |() ≡
1 > 0 en
0,1. Evidentemente, = 0, ° = 1, 4 = 1, °4 = 0, de modo que
4 + °4 = 1 > 0, 44 + °44 = 1 > 0. Tomemos como función admisible () =
1 − 4 . De acuerdo con el principio de Rayleigh tenemos

4µ
2(1 − 4 )

¡ ≤ = = 3 = 2,5
| (1 − ) 8µ
15
4 4
El valor exacto es ¡ = ≈ 2,4674.

Cf
_
■

8. Biografías
8.1. Leonhard Euler
(Basilea, Suiza, 1707-San

1707 San Petersburgo, 1783) Matemático suizo. Las
facultades
tades que desde temprana edad demostró para las matemáticas pronto le
ganaron la estima del patriarca de los Bernoulli, Johann, uno de los más
eminentes matemáticos de su tiempo y profesor de Euler en la Universidad de
Basilea. Tras graduarse en dicha institución
institución en 1723, cuatro años más tarde fue
invitado personalmente por Catalina I para convertirse en asociado de la
Academia de Ciencias de San Petersburgo, donde coincidió con otro miembro de
la familia Bernoulli, Daniel, a quien en 1733 relevó en la cátedra cáte de
matemáticas.
A causa de su extrema dedicación al trabajo, dos años más tarde perdió la

visión del ojo derecho, hecho que no afectó ni a la calidad ni al número de sus
hallazgos. Hasta 1741, año en que por invitación de Federico el Grande se
trasladó a la Academia de Berlín, refinó los métodos y las formas del cálculo
integral (no sólo gracias a resultados novedosos, sino también a un cambio en
los habituales métodos de demostración geométricos, que sustituyó por
métodos algebraicos), que convirtió en una herramienta de fácil aplicación a
problemas de física. Con ello configuró en buena parte las matemáticas
aplicadas de la centuria siguiente (a las que contribuiría luego con otros
resultados destacados en el campo de la teoría de las ecuaciones diferenciales
difere
lineales), además de desarrollar la teoría de las funciones trigonométricas y
logarítmicas (introduciendo de paso la notación e para definir la base de los
logaritmos naturales).

En 1748 publicó la obra Introductio in analysim infinitorum, en la que
expuso el concepto de función en el marco del análisis matemático, campo en el
que así mismo contribuyó de forma decisiva con resultados como el teorema
sobre las funciones homogéneas y la teoría de la convergencia. En el ámbito de
la geometría desarrolló conceptos básicos como los del ortocentro, el
circuncentro y el baricentro de un triángulo, y revolucionó el tratamiento de las
funciones trigonométricas al adoptar ratios numéricos y relacionarlos con los
números complejos mediante la denominada identidad de Euler; a él se debe la
moderna tendencia a representar cuestiones matemáticas y físicas en términos
aritméticos.
En el terreno del álgebra obtuvo así mismo resultados destacados, como el de

la reducción de una ecuación cúbica a una bicuadrada y el de la determinación
de la constante que lleva su nombre. A lo largo de sus innumerables obras,
tratados y publicaciones introdujo gran número de nuevas técnicas y contribuyó
sustancialmente a la moderna notación matemática de conceptos como función,
suma de los divisores de un número y expresión del número imaginario raíz de
menos uno. También se ocupó de la teoría de números, campo en el cual su
mayor aportación fue la ley de la reciprocidad cuadrática, enunciada en 1783.
A raíz de ciertas tensiones con su patrón Federico el Grande, regresó

nuevamente a Rusia en 1766, donde al poco de llegar perdió la visión del otro
ojo. A pesar de ello, su memoria privilegiada y su prodigiosa capacidad para el
tratamiento computacional de los problemas le permitieron continuar su
actividad científica; así, entre 1768 y 1772 escribió sus Lettres à une princesse
d’Allemagne, en las que expuso concisa y claramente los principios básicos de la
mecánica, la óptica, la acústica y la astrofísica de su tiempo.
De sus trabajos sobre mecánica destacan, entre los dedicados a la mecánica de

fluidos, la formulación de las ecuaciones que rigen su movimiento y su estudio
sobre la presión de una corriente líquida, y, en relación a la mecánica celeste, el
desarrollo de una solución parcial al problema de los tres cuerpos –resultado de
su interés por perfeccionar la teoría del movimiento lunar–, así como la
determinación precisa del centro de las órbitas elípticas planetarias, que
identificó con el centro de la masa solar. Tras su muerte, se inició un ambicioso
proyecto para publicar la totalidad de su obra científica, compuesta por más de
ochocientos tratados, lo cual lo convierte en el matemático más prolífico de la
historia.

8.2. Karl Gustav Jacobi
(Potsdam, actual Alemania, 1804-Berlín, 1851) Matemático alemán. Hijo de

una familia de banqueros de origen judío, estudió en la Universidad de Berlín,
donde se doctoró en 1825. Convertido al cristianismo, tuvo oportunidad de
acceder a un puesto de profesor en la Universidad de Königsberg.
Destacadísimo pedagogo, influyó en numerosas generaciones posteriores de
matemáticos alemanes. Sus trabajos más relevantes se produjeron en el campo
del álgebra, en el que introdujo y desarrolló el concepto de determinante,
aplicándolo así mismo al estudio de las funciones de variables múltiples. Entre
1826 y 1827 estableció, independientemente del noruego Niels Henrik Abel, los
principios fundamentales de la teoría de las funciones elípticas. En el ámbito de
la teoría de números, demostró el teorema de Bachet sobre el total de las
descomposiciones posibles de un entero, y en el de la mecánica física, trató con
profundidad y rigor el problema de los tres cuerpos.
Su obra más notable es sobre la formación y propiedades de los
determinantes (1841).

8.3. Adrien Marie Legendre
Auteuil, Francia, 1833) Matemático francés. Tras completar sus

(París, 1752-Auteuil,
estudios en el Collège Mazarin, entró a trabajar en la Escuela Militar, para la
que completó un estudio sobre la trayectoria de los proyectiles que le supuso
sup el
Premio de la Academia de Berlín en 1782. A partir de 1795 enseñó matemáticas
en la École Normale.
En sus primeros trabajos, centrados en la mecánica, introdujo conceptos
como la función que lleva su nombre o la primera demostración del método de
los mínimos cuadrados.
Tras los pasos de Euler y Lagrange, estudió las funciones elípticas y las redujo
a tres formas básicas. Fue el primero en dedicar una obra estrictamente a la
teoría de números, ámbito en el que obtuvo resultados fundamentales como la
demostración en 1830 de la ley de la reciprocidad cuadrática.
En 1794 publicó los Elementos de geometría, una versión reordenada y
simplificada de la obra original de Euclides, que fue traducida a más de treinta
idiomas.
8.4. Karl Weiertrass

Karl Weierstrass.. (Ostenfelde, actual Alemania, 1815- 1815-Berlín, 1897)
Matemático alemán, e hijo de un oficial a las órdenes de Napoleón, Karl era el
mayor de cuatro hermanos. Más tarde, su padre ingresó en en el servicio de
recaudación de impuestos en Prusia, lo que obligó a la familia a trasladarse
constantemente. Con catorce años, Karl fue aceptado en la escuela católica de
enseñanza secundaria de Paderborn. Ganó algunos premios antes de graduarse,
y en 1834,
4, siguiendo los deseos de su padre, ingresó en la Universidad de Bonn
para estudiar comercio y finanzas. Sin embargo, estas materias no le
interesaban y pasó la mayor parte del tiempo bebiendo, practicando esgrima y
leyendo libros de matemáticas.
En 1839 fue aceptado en la Academia de Teología y Filosofía de Münster,

donde encontró la inspiración matemática de manos de Christof Guderman.
Éste le introdujo en la teoría de las series de potencias, que más tarde serían la
base de todo su trabajo. Su primer escrito
escrito importante, publicado en 1841, fue un
ensayo sobre funciones elípticas.
Durante los quince años siguientes se dedicó a dar clase en una escuela de
enseñanza secundaria. En 1854 envió un trabajo sobre funciones abelianas a
una publicación matemática de prestigio, y sorprendió a la comunidad
matemática con su genio. Por este trabajo recibió el doctorado honorífico de la
Universidad de Königsberg y en 1856 fue aceptado como profesor asociado en la
Universidad de Berlín.
Abrumado por las enormes responsabilidades

responsabilidades de su nuevo cargo, sufrió una
crisis nerviosa en 1861, que le apartó de las aulas dos años. A pesar de ello, en
1864 fue ascendido a profesor, cargo que ostentó el resto de su vida.
Desafortunadamente, tras los ataques públicos de Kronecker por su apoyo a las
ideas de Cantor, y la muerte de su amiga Sonja Kovalevsky, se hundió
mentalmente y pasó el resto de su vida en una silla de ruedas hasta que murió
víctima de una neumonía.

9. Bibliografía
Las consultas realizadas para la elaboración de este trabajo han sido las
siguiente:
[1] Biografía y vidas http://www.biografiasyvidas.com
[2] Elsgoltz, L. Ecuaciones Diferenciales y Cálculo Variacional. Editorial MIR

(Moscú) Rusia. 1969.
[3] Galiano, G. Introducción al Cálculo Variacional. 2003.

http://orion.ciencias.uniovi.es/~galiano/docencia/main03.pdf
[4] Krasnov, M.L., Kiseliov, V, Makárenko, G.I. Cálculo Variacional. Editorial

URSS (Moscú). 2002.

Calculo Variacional PDF

Cargado por

Copyright:

Formatos disponibles

Calculo Variacional PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Calculo Variacional PDF

Cargado por

Copyright:

Formatos disponibles

Cálculo Variacional.

Condiciones necesarias y suficientes de extremo. La

6.1. Función (, , , ′). Condición de Weiertrass………….……….………..59

7. Anexo: Métodos directos en el cálculo variacional. …………………….………..72

Universidad de Almería. | Cálculo Variacional. 2

Universidad de Almería. | Cálculo Variacional. 3

1.1. Introducción histórica

Las raíces del cálculo variacional se extienden a tiempos anteriores a la

Hasta aquí hemos hablado de problemas de cálculo variacional, pero no

Se le atribuye a Pierre de Fermat, matemático francés del siglo XVII, el

Universidad de Almería. | Cálculo Variacional. 4

En 1696 Johann Bernoulli publica un desafío para los matemáticos de su

Johann Bernoulli también estudio geodésicas en varias superficies. Este es

Durante el siglo XIX los trabajos de Euler y Lagrange son formalizados y

El desarrollo del cálculo variacional está relacionado con el desarrollo de la

1.2. Definición de Cálculo Variacional

Las funciones definidas sobre un conjunto cuyos elementos son funciones, se

Universidad de Almería. | Cálculo Variacional. 5

El área S de cierta superficie es también un funcional, puesto que se

donde D es la proyección de la superficie en el plano OXY .

Los momentos de inercia, los momentos estáticos, las coordenadas del

En todos estos ejemplos se tiene una dependencia que es característica para

En el Cálculo Variacional se consideran los métodos para hallar el valor

Durante más de dos siglos, el Cálculo de Variaciones ha sido una de las

Es fácil captar el interés del tema si se toman en cuenta algunos de sus

1. Un espacio de funciones V , tal que u : Ω → ℜ q , donde Ω es un

2. Restricciones sobre el conjunto de soluciones, que pueden imponerse

3. Un funcional J : V → ℜ de la forma siguiente:

Universidad de Almería. | Cálculo Variacional. 6

La hipótesis sobre V y L deben asegurar la existencia de J sobre V , o al

El problema de optimización consiste en hallar el mínimo, u ∈ U , del

1.3. Problemas importantes

1.3.1. Problema isoperimétrico

1.3.2. Principio de acción mínima

En Mecánica Clásica, cuando una partícula se mueve bajo la acción de un

Una pregunta interesante es: ¿Podemos singularizar el movimiento real

Universidad de Almería. | Cálculo Variacional. 7

La respuesta a esta pregunta es un principio básico en Física, que en

La caracterización dada por el principio de Hamilton asocia una cantidad,

El principio de la mínima acción dice: entre todos los movimientos

¿Cuál es la relación entre este principio y la forma newtoniana de planear las

1.3.3. Problema de la Braquistocrona

El problema de la braquistocrona, o curva de descenso más rápido, es uno

donde y es una función regular. Como la curva pasa por A y B , la función y

0 = y (0) , y (x1 ) = y1 (3)

El movimiento de la masa puntual puede describirse por medio de la ley de

donde E > 0 es una constante y h A es la altura a la que se encuentra el punto A .

y tomando la coordenada vertical como x = h A − h , deducimos que la velocidad

siendo s una parametrización de la trayectoria del punto material. Deducimos

y como la longitud de arco de la curva viene dada por

Universidad de Almería. | Cálculo Variacional. 9

Hallar la braquistocrona es equivalente a resolver el siguiente problema de

1.3.4. El problema de las Geodésicas

1. Consideremos una superficie regular S ⊂ ℜ3 definida por la

x = x(u, v ) , y = y (u, v ) , z = z (u, v ) ,

con (u , v ) ∈ [u 0 , u1 ]× [v 0 , v1 ] . Cualquier curva contenida en S puede

El elemento de arco de las curvas contenidas en S está determinado

ds 2 := Eu ' 2 +2 Fu ' v'+Gv ' 2

De modo que la longitud del arco entre los puntos correspondientes a

6.1. Función (, , , ′). Condición de Weiertrass………….……….………..59

Además (, , ′ ) se considera derivable tres veces. Sabemos que la

el funcional. Supongamos que en la curva = (), derivando dos veces, se

también segunda derivada). Tomemos cierta curva admisible = () cercana a

la curva = () e incluyamos ambas curvas en la familia mono-paramétricas

(, ) = () + (() − ())

cuando = 0, se obtiene la curva = (); para = 1, se tiene = (). Lo

de la variación. De este modo, consideremos la familia = (, ), donde

(, ) = () + , que contiene para = 0 la curva en la cual se alcanza el

sólo en las curvas de la familia = (, ), el funcional se transforma en una

ya que el valor del parámetro determina una curva de la familia = (, )

= (, ) . La condición necesaria para que la función () tenga un extremo

en = 0, es la anulación de su derivada para = 0, es decir, ′ (0) = 0. Como

sustituyendo por = 0, obtenemos:

anulación de su variación: = 0. Para el funcional:

′ = ()′ obtenemos que:

donde el primer factor ! − . !′ es una función continua dada en la curva

= () que realiza el extremo, y el segundo factor es una función arbitraria

Lema 2. (Fundamental del Cálculo de Variaciones). Sea / ∈ 1 (

para todo ∈ 1 (

ya que el producto /() () conserva su signo en el segmento ̅ ≤ ≤ y se

extremo, el factor ! − !′ es función continua, y la variación es una

ya previstas en el lema fundamental. Por lo tanto, ! − !′ ≡ 0 en la curva

! − !′ − !!′ ′ − !′!′ ′′ = 0

variaciones. Las curvas integrales de la ecuación de Euler = (, 1 , 14 ) se

la solución general de esta ecuación, de las condiciones de frontera ( ) =

, ( ) = . Sólo en las extremales que satisfacen estas condiciones se puede