TEMA 4

ECUACIONES DIFERENCIALES PARCIALES

4.1. INTRODUCCIÓN

En este tema se verán procedimientos para resolver ecuaciones en derivadas

parciales que surgen con frecuencia en problemas donde aparecen vibraciones,
potenciales y distribuciones de temperatura. Se llaman problemas de valores en la
frontera y se describen mediante ecuaciones en derivadas parciales de primer o
segundo orden, que son relativamente simples. En el caso de las ecuaciones de
segundo orden, lo que se hace es hallar las soluciones particulares de una ecuación en
derivadas parciales reduciéndolas a dos o más ecuaciones diferenciales ordinarias. Se
comenzará con el método de separación de variables para ecuaciones en derivadas
parciales lineales. La aplicación de este método recordará conceptos de temas
anteriores como valores propios, funciones propias y el desarrollo de una función en
una serie infinita de funciones ortogonales.

4.2. ESPACIO VECTORIAL CON PRODUCTO INTERIOR

Definición 1. Un espacio vectorial real es un conjunto V (cuyos elementos se

llamarán vectores), con una operación binaria de adición y una multiplicación escalar
que satisfacen las siguientes propiedades: (α, β ∈ R)
a. Dado cualquier par de vectores x, y ∈ V , existe un único vector x + y ∈ V , llamado
la suma de x y y
b. ∀x , y ∈ V , x + (y + z) = (x + y) + z
c. ∀x , y ∈ V , x + y = y + x
d. ∃0 ∈ V ∀x ∈ V , x + 0 = 0 + x = x
e. ∃(− x) ∈ V ∀x ∈ V , x + (− x) = 0
f. Dado cualquier vector x ∈ V y cualquier número real α , existe un único vector
αx ∈ V llamado el producto de α y x
g. α(x + y) = αx + αy
h. (α + β)x = αx + βx
i. (αβ)x = α(βx)
j. 1.x = x

Definición 2. Se dice que un subconjunto S ≠ φ de un espacio vectorial V, es un

subespacio de V, si ∀x, y ∈ S y ∀α, β ∈ R el vector αx + β y ∈ S .

Observación 1. Es claro que {0} y V son subespacios de V.

Prof. José Luis Quintero 1

Definición 3. Se dice que un producto interior está definido en un espacio vectorial
real V, si con cada par de vectores x e y en V se asocia un número real x • y tal que
a. x•y = y•x
b. (αx) • y = α(x • y) para cada número real α
c. (x1 + x2 ) • y = x1 • y + x2 • y
d. x • x ≥ 0 y x • x = 0 si, y solamente si, x = 0 .

Definición 4. Un espacio vectorial con un producto interior se conoce como un

espacio euclidiano o espacio vectorial con producto interior.

Si se aplica (a) a (b) y a (c), se ve que un producto interior también satisface

las propiedades x • (αy) = α(x • y) y x • (y1 + y2 ) = x • y1 + x • y2 . En general
(α1x1 + α2 x2 ) • (β1y1 + β2 y2 ) = α1β1(x1 • y1 ) + α1β2 (x1 • y2 ) + α2β1(x2 • y1 ) + α2β2 (x2 • y2 )
donde α1 , α2 , β1 y β2 son escalares arbitrarios.

Ejemplo 1. Sean x = (x1 ,..., xn ) y y = (y1 ,..., yn ) vectores en Rn , y defínase a x • y

por x • y = x1y1 + ... + xnyn . Entonces Rn se convierte en un espacio euclidiano y, como
tal, se llama espacio euclidiano de dimensión n. En R 2 y R 3 este producto interior
no es otro que el “producto punto” usado en física, donde la definición generalmente
se expresa geométricamente como el producto de la longitud de x por la longitud de y
por el coseno del ángulo que forman ambos.

Ejemplo 2. Considere el espacio C a,b de todas las funciones continuas en a,b , con
f • g definido por
b
f •g =
∫ a
f(x)g(x)dx .

Pruebe que f • g define un producto interior.

Solución.
No es difícil demostrar que satisface todas las exigencias de un producto interior: Se
tiene que
b b

∫ a
αf(x)g(x)dx = α
∫ a
f(x)g(x)dx ,

siendo α un número real, y

b b b

∫ a
f1(x) + f2 (x) g(x)dx = ∫ a
f1(x)g(x)dx +
∫ a
f2 (x)g(x)dx .

Finalmente, al recordar que la integral de una función no negativa es no negativa, y

que la integral de una función continua no negativa es cero si, y solamente si, la
función es idénticamente cero se ve que
b
f•f =
∫ a
(f(x))2 dx ≥ 0 ,

y que f • f = 0 si, y solamente si, f = 0 .

Prof. José Luis Quintero 2

 En adelante, siempre que se haga referencia a Rn o a C a,b como espacios
euclidianos, se considerará que se están utilizando los productos interiores definidos en
los ejemplos anteriores, a menos que se mencione expresamente lo contrario.

Ejemplo 3. Sea r una función no negativa en C a,b que se anula, cuando más, en un
número finito de puntos en el intervalo a,b . Entonces
b
f•g=
∫ a
f(x)g(x)r(x)dx

también define un producto interior en C a,b .

La función r se llama la función peso para este producto interior, y se

observará que cuando r es idénticamente igual a 1 se tendrá la fórmula definida en el
ejemplo anterior. Se encontrará este producto interior nuevamente cuando se estudien
problemas de valor en la frontera, para ecuaciones diferenciales.

4.3. LONGITUD. MEDIDA ANGULAR. DISTANCIA

Se definirá longitud, medida angular y distancia, en términos del producto

interior de un espacio euclidiano. Cada uno de estos conceptos tiene un significado
bien definido en el espacio euclidiano, de dimensión 2, y es razonable exigir que
cualquier definición que se adopte, se reduzca a la ya usual en R2 .

Sea x = (x1 , x2 ) cualquier vector en R 2 . Entonces, la longitud de x indicada por

x , es el número real no negativo x = x12 + x22 . Esta expresión puede replantearse

en términos del producto interior en R 2 como x = x • x y se tiene la siguiente
definición:

Definición 5. La longitud (o norma) de un vector x en un espacio euclidiano, se

define como el número real no negativo x = x•x.

Así, la longitud de un vector x = (x1 ,..., xn ) en Rn , es x = x12 + ... + xn2 ,

mientras que la longitud de un vector f en C a,b es
1/2
 b 
f =
 ∫ a
(f(x))2 dx 

.

A continuación, se observa que si x y y son dos vectores no nulos cualesquiera,

2
en R , la fórmula
x•y
cos θ = , 0 ≤ θ ≤ π,
x y
es una consecuencia inmediata de la ley de los cosenos. La expresión

Prof. José Luis Quintero 3

 x•y
x y
también es significativa en un espacio euclidiano arbitrario, hecho que sugiere
considerarla como un candidato razonable para la definición general del cos θ . Sin
embargo, antes de proceder, debe establecerse la desigualdad
x•y
−1 ≤ ≤1
x y
para cada par de vectores no nulos en un espacio euclidiano, ya que, por supuesto,
cualquier definición de cos θ debe satisfacer la desigualdad −1 ≤ cos θ ≤ 1 . Este hecho
surgirá como una consecuencia del siguiente resultado importante, conocida como la
desigualdad de Schwarz o de Cauchy-Schwarz.

TEOREMA 1. (Desigualdad de Schwarz) Si x y y son dos vectores cualesquiera en

un espacio euclidiano, entonces (x • y)2 ≤ (x • x)(y • y) .
Demostración.
Primero se observa que esta desigualdad es inmediata si x o y, es el vector cero, ya
que entonces ambos miembros de la desigualdad son cero. Así, es suficiente
considerar el caso en que x y y son no nulos. Aquí se usan propiedades del producto
interior para desarrollar (αx − βy) • (αx − βy) donde α y β son números reales
arbitrarios. Por la propiedad d de la definición 3 se tiene
0 ≤ (αx − βy) • (αx − βy) = α2 (x • x) − 2αβ(x • y) + β2 (y • y) ,
de donde 2αβ(x • y) ≤ α2 (x • x) + β2 (y • y) . Ahora si α = y•y y β = x • x se obtiene
2 x • x y • y (x • y) ≤ 2(x • x)(y • y) que es equivalente a x • y ≤ x • x y • y . Ahora
bien, si se eleva al cuadrado, se obtiene lo que se quería demostrar. ■

En Rn la desigualdad de Schwarz toma la forma

2
 n   n  n 



∑
i =1
xi yi  ≤ 



 i =1
∑
xi2  
 ∑
  i =1
yi2 


y en C a,b se convierte en
2
 b   b  b 

 ∫ a
f(x)g(x)dx  ≤ 
  ∫ a
(f(x))2 dx  
 ∫ a
(g(x))2 dx  .


La primera de estas desigualdades es válida para cualquier colección x1 ,..., xn ,

y1 ,..., yn de números reales, y generalmente se le llama la desigualdad de Cauchy.
Es importante recordarla ya que es frecuentemente útil en la deducción de otras
desigualdades aritméticas.

En la notación de la definición 5, la desigualdad de Schwarz se convierte en

x • y ≤ x y , y afirma que el valor absoluto del producto interior de dos vectores no
excede el producto de las longitudes de los vectores. Así
x•y
≤ 1,
x y

Prof. José Luis Quintero 4

siempre que x y y sean no nulos.

Ejemplo 4. Sean a1 ,..., an números reales positivos. Use la desigualdad de Cauchy

para demostrar que
1 1  2
(a1 + ... + an )  + ... + ≥n .
a
 1 an 

Solución.
Sean
 1 1 
x = ( a1 ,..., an ) y =  ,..., .
 a an 
 1 
De modo que
1 1
x • x = a1 + ... + an , y • y = + ... + , x•y =n
a1 an
y aplicando la desigualdad de Cauchy se obtiene lo que se quería demostrar.

Ejemplo 5. Sean a, b y c, números reales positivos, tales que a + b + c = 1 . Utilice la

desigualdad de Cauchy para demostrar que
1 1  1 
 a − 1  b − 1  c − 1 ≥ 8 .
   
Solución.
Sean
 1 1 1 
x = ( a, b, c) y =  , , .
 a b c
De modo que
1 1 1
x • x = a + b + c = 1, y • y = + + , x•y =3
a b c
y aplicando la desigualdad de Cauchy
1 1 1
+ + ≥ 9.
a b c
Por lo tanto manipulando se tiene
bc + ac + ab bc + ac + ab bc + ac + ab − abc (1 − a)(1 − b)(1 − c)
≥9⇒ −1 ≥ 8 ⇒ ≥8⇒ ≥8
abc abc abc abc
1 − a 1 − b  1 − c  1 1  1 
⇒      ≥ 8 ⇒  − 1  − 1  − 1 ≥ 8
 a  b  c  a  b c 

Ejemplo 6. Demuestre que la desigualdad

2
 a1 + ... + an  a12 + ... + an2
  ≤
 n  n
persiste para cualquier colección de números reales a1 ,..., an .
Solución.
Sean
a a 
x = (1,...,1) y =  1 ,..., n  .
 n n
Entonces

Prof. José Luis Quintero 5

 a12 + ... + an2 a1 + ... + an
x • x = n, y•y = 2
, x•y =
n n
y aplicando la desigualdad de Cauchy se obtiene lo que se quería demostrar.

Ejemplo 7. Sea f ∈ C1 0,1 . Pruebe que

1/2
1  1 
∫ 0
( 12 − t)f(t)dt ≤ 1
4
2
15 
 ∫ 0
f '(t) 2dt 

.

Solución.
t2
Integrando por partes: u = f(t) ⇒ du = f '(t)dt , dv = 1
2
−t⇒v = t
2
− 2
.
1 1 1 1

∫ t2
∫ t2
∫
2
( 12 − t)f(t)dt = ( 2t − 2
)f(t) − ( 2t − 2
)f '(t)dt = ( t2 − 2t )f '(t)dt .
0 0 0 0

Aplicando la desigualdad de Schwarz:

1/2 1/2
1  1   1 
∫ ∫ ∫
2 2
( t2 − 2t )f '(t)dt ≤  ( t2 − 2t )2 dt   (f '(t))2 dt  .
0  0   0 
1 1 1

∫ ∫
2 4 t3 t2 5 t4 t3 6 −15 +10
( t2 − 2t )2 dt = ( t4 − 2
+ 4
)dt = ( 20
t
− 8
+ 12
) = 1
20
− 1
8
+ 1
12
= 120
= 1
120
.
0 0 0

Por lo tanto:
1/2
 1 
∫
2 t 2 1 1/2
 ( t2 − 2
) dt  = (120 ) = 1
= 2
= 1
4
2
15
.
2 2 15 4 15
 0 
En consecuencia se demuestra que:
1/2
1  1 
∫ ∫
2
( 12 − t)f(t)dt ≤ 1
4
2
15  f '(t) dt  .
0  0 

Definición 6. Si x y y son vectores no nulos en un espacio euclidiano, se define el

coseno del ángulo entre ellos como
x•y
cos θ = .
x y

Queda finalmente pendiente por definir la distancia entre dos puntos

cualesquiera (esto es, vectores) en un espacio euclidiano. Nuevamente esto se hace
sencillamente copiando la definición de R 2 donde la distancia entre x y y es la longitud
del vector x − y . Así se establece la definición de distancia:

Definición 7. La distancia entre dos vectores x y y en un espacio euclidiano es, por

definición, d(x, y) = x − y .

Pero ¿es ésta una definición razonable del término “distancia”? Con objeto de
contestar esta pregunta, primero debe decidirse qué propiedades se requieren de la
distancia, en general. En ese aspecto, los matemáticos están de acuerdo, y lo han
decidido como sigue: la distancia entre dos puntos debe ser un número real no
negativo que es cero si, y solamente si, los puntos coinciden. Debe ser independiente

Prof. José Luis Quintero 6

del orden en el que se consideren los puntos y, finalmente, la desigualdad del
triángulo, famosa en la geometría plana, debe satisfacerse.

Así, con el objeto de justificar el uso del término “distancia” en la definición 7

se debe demostrar que d(x,y) es un número real que satisface las siguientes
propiedades:
a. d(x, y) ≥ 0
b. d(x, y) = 0 si, y solamente si, x = y
c. d(x, y) = d(y, x)
d. d(x, y) + d(y, z) ≥ d(x, z) , x, y y z tres vectores cualesquiera.

Las tres primeras de estas propiedades se deducen inmediatamente de la

definición de longitud y de los axiomas que gobiernan un producto interior. La última,
sin embargo, no es tan obvia. Para demostrarla, primero se establecerá una
desigualdad que es de cierta importancia por sí misma.

LEMA 1. Si x y y son vectores arbitrarios en un espacio euclidiano, entonces

x+y ≤ x + y .
Demostración.
1 /2 1/2
x + y = (x + y) • (x + y) = (x • x) + 2(x • y) + (y • y)
1/2

( )
1/2 2
≤ (x • x) + 2 x • x y • y + (y • y) = x•x + y•y  = x•x + y•y
  
≤ x + y
■
La desigualdad del triángulo sigue inmediatamente de este resultado. En
efecto, x − z = (x − y) + (y − z) ≤ x − y + y − z , que precisamente es lo que se tenía
que demostrar. Finalmente, se observa que la función distancia definida antes también
goza de las propiedades notables siguientes:
a. d(αx, αy) = α d(x, y) para todo número real α
b. d(x + z, y + z) = d(x, y)

4.4. ORTOGONALIDAD. ORTOGONALIZACIÓN

Definición 8. Se dice que dos vectores en un espacio euclidiano son ortogonales o

perpendiculares, si el coseno del ángulo que forman, vale cero.

En seguida se generalizará algo más la noción de ortogonalidad, pero primero

se probará un teorema particularmente notable.

TEOREMA 2. (Pitágoras) Dos vectores x y y en un espacio euclidiano se dice que son

2 2 2
ortogonales si, y solamente si, x + y = x + y .
Demostración.

Prof. José Luis Quintero 7

 2 2 2
x+y = (x + y) • (x + y) = x • x + 2(x • y) + y • y = x + 2(x • y) + y
2 2 2
Así x + y = x + y si, y solamente si, x • y = 0 , como se afirmó. ■

Ejemplo 8. Sean x1 , ..., xn vectores mutuamente perpendiculares en un espacio

2 2 2
euclidiano. Demuestre que x1 + ... + xn = x1 + ... + xn .
Solución. Se sabe que
2
x1 + ... + xn = (x1 + ... + xn ) • (x1 + ... + xn )
n n n n n n
= ∑
i =1
(xi • xi ) + ∑∑i =1 j =1
(xi • x j ) = ∑
i =1
xi
2
+ ∑∑ (x • x )
i =1 j =1
i j

i≠ j i≠ j

Así
n
x1 + ... + xn
2
= ∑x
i =1
i
2
= x1
2
+ ... + xn
2

ya que xi • x j = 0, i ≠ j .

Dicho esto, se continúa la discusión, dando la siguiente definición:

Definición 9. Se dice que un conjunto de vectores x1 , x2 ,..., xi ,... en un espacio

euclidiano, es un conjunto ortogonal si xi ≠ 0 para todo i, y que xi • x j = 0 siempre
que i ≠ j . Si, además, xi • xi = 1 para cada i, se dice que el conjunto es ortonormal.

Así, un conjunto ortogonal es un conjunto de vectores no nulos perpendiculares

entre sí, mientras que un conjunto ortonormal es un conjunto ortogonal en el cual cada
uno de los vectores es de longitud unitaria. Por economía de notación, cuando se
discutan conjuntos ortonormales, las propiedades anteriores se combinarán
frecuentemente escribiendo
0 si i ≠ j
xi • x j = δij =  .
1 si i = j
El símbolo δij que se introduce aquí, se le llama la delta de Kronecker.

TEOREMA 3. Un conjunto ortogonal de vectores {x1 , x2 ,..., xn } es linealmente

independiente.
Demostración.
Se sabe de cursos anteriores de Álgebra Lineal que un conjunto de vectores
{x1 , x2 ,..., xn } es linealmente independiente si
α1x1 + α2 x2 + ... + αn xn = 0 ⇔ α1 = α2 = ... = αn = 0 .
Se tiene una combinación lineal de los vectores ortogonales:
α1x1 + α2 x2 + ... + αn xn = 0
Haciendo producto escalar en la expresión anterior por un vector del conjunto xi se
tiene
2
α1x1 • xi + α2 x2 • xi + ... + αixi • xi + ... + αn xn • xi = 0 ⇒ αi xi = 0.

Prof. José Luis Quintero 8

De modo que αi = 0 ∀i ∈ {1, 2,...,n} . Esto implica que el conjunto ortogonal de vectores
{x1 , x2 ,..., xn } es linealmente independiente (forman una base) por lo tanto forman un
sistema completo.

Ejemplo 9. En R3 los vectores (1, 0, 0), (0, 2, 0), (0, 0, −1 / 2) forman un conjunto
ortogonal, mientras que los vectores de base canónica (1,0,0), (0,1,0) y (0,0,1)
forman un conjunto ortonormal. En forma más general, el conjunto que consiste de los
vectores de la base en Rn es ortonormal.

Ejemplo 10. Se define un polinomio trigonométrico de grado 2n + 1 como una

expresión de la forma
a
f(x) = 0 + a1 cos(x) + a2 cos(2x) + ... + an cos(nx) + b1sen(x) + b2sen(2x) + ... + bnsen(nx)
2
donde a0 ,...,bn son números reales, y an ≠ 0 o bn ≠ 0 o ambos. Sea Tn el conjunto de
todos los polinomios trigonométricos de grado ≤ 2n + 1 , junto con el polinomio cero.
Haciendo que Tn sea un espacio euclidiano al definir la adición y la multiplicación
escalar de los polinomios trigonométricos, término a término como con los polinomios
ordinarios y un producto interior dado por
π
f •g =
∫ −π
f(x)g(x)dx ,

el conjunto de funciones 1, cos(x), sen(x), ..., cos(nx), sen(nx) es un conjunto ortogonal

en Tn ya que, para enteros m y n no negativos
π π

∫ −π
sen(mx)sen(nx)dx = 0,
∫ −π
sen(mx) cos(nx)dx = 0,
π

∫ −π
cos(mx) cos(nx)dx = 0, si m ≠ n .

Para normalizar este conjunto se observa que

π π π

∫ −π
1dx = 2π y
∫ −π
sen2 (mx)dx =
∫ −π
cos2 (mx)dx = π, si m > 0 .

Por tanto, las funciones

1 cos(x)
cos(nx) sen(x) sen(nx)
, , ,...,
,...,
2π π π π π
forman un conjunto ortonormal en Tn .

Ejemplo 11. Se deduce del ejemplo precedente que el conjunto (infinito) 1, cos(x),
sen(x), …, cos(nx), sen(nx), … es ortogonal en C −π, π .

En lugar de empezar con la situación más general, se introducirá el proceso de

ortogonalización con dos ejemplos. El primero es extraído de R2 donde se considera
un par de vectores linealmente independientes x1 y x2 . Entonces x1 y x2 forman
una base para R2 y, en este caso, el problema se convierte en la sustitución de x1 y
x2 por una base ortogonal e1 y e2 construida en términos de x1 y x2 de alguna
manera razonable. Una forma sencilla de resolver el problema es simplemente tomar
e1 = x1 y luego se expresa e2 en la forma e2 = x2 − αe1 , y luego se determina α de

Prof. José Luis Quintero 9

modo que la condición de ortogonalidad e2 • e1 = 0 quede satisfecha. Esto da la
ecuación x2 • e1 − α(e1 • e1 ) = 0 , y en consecuencia, el valor de α es
x2 • e1
α= .
e1 • e1
Con esto e2 se ha determinado en términos de x1 (= e1 ) y x2 , y la base x1 y x2 en
R 2 se ha ortogonalizado.

Ejemplo 12. Si x1 = (1,1) y x2 = (0,1) , entonces

(0,1) • (1,1) 1
α= = ,
(1,1) • (1,1) 2
y así e1 = (1,1), e2 = (− 12 , 12 ) .

Como otro ejemplo, se ortogonaliza una base arbitraria x1 , x2 y x3 en R 3 . El

procedimiento esencialmente es el mismo usado antes en R 2 , y se inicia eligiendo
e1 = x1 . El segundo paso consiste en determinar e2 de acuerdo con el par de
ecuaciones e2 • e1 = 0, e2 = x2 − αe1 , lo que da nuevamente
x2 • e1
α= .
e1 • e1

Está claro que e2 no es el vector cero, y también lo es que e1 y e2

pertenecen al subespacio de R 3 generado por x1 y x2 . En consecuencia δ(e1 , e2 ) es
un subespacio de δ(x1 , x2 ) . Más aún, ya que los vectores ortogonales e1 y e2 son
linealmente independientes, δ(e1 , e2 ) tiene la misma dimensión que δ(x1 , x2 ) . Así,
δ(e1 , e2 ) = δ(x1 , x2 ) .

Combinado con el hecho de que x1 , x2 y x3 forman una base para R 3 , esta

igualdad implica que x3 no pertenece al subespacio de R 3 generado por e1 y e2 . El
proceso de ortogonalización debe completarse aceptando que e3 sea la componente
de x3 , perpendicular al subespacio δ(e1 , e2 ) . Así, se hace e3 = x3 − α1e1 − α2e2 , y se
encuentran α1 y α2 por medio de las condiciones de la ortogonalidad e1 • e2 = e1 • e3
= e2 • e3 = 0 . Ellas dan el par de ecuaciones dadas por 0 = x3 • e1 − α1(e1 • e1 ),
0 = x3 • e2 − α2 (e2 • e2 ) de donde
x3 • e1 x • e2
α1 = , α2 = 3 .
e1 • e1 e2 • e2

Con esto se completa el proceso de ortogonalización de la base x1 , x2 y x3 .

Ahora se puede disponer de la situación general requerida para ortogonalizar un
conjunto arbitrario de vectores linealmente independientes x1 , x2 , ... en un espacio
euclidiano. Se va a presentar un teorema que puede resumir este proceso de
generalización.

TEOREMA 4. Sea x1 , x2 , ... un conjunto (finito o infinito) de vectores linealmente

independientes en un espacio euclidiano T. Entonces existe un conjunto ortogonal
e1 , e2 , ... en T tal que para cada entero n, δ(e1 ,..., en ) = δ(x1 ,..., xn ) . Más aún, el en

Prof. José Luis Quintero 10

puede escogerse de acuerdo con la regla e1 = x1 , y en +1 = xn +1 − α1e1 − ... − αnen ,
donde
xn +1 • e1 x • e2 x • en
α1 = , α2 = n +1 , ... , αn = n +1 .
e1 • e1 e2 • e2 en • en

El método de ortogonalización descrito en este teorema, es conocido como el

proceso de ortogonalización de Gram-Schmidt.

Ejemplo 13. En este ejemplo se aplica el proceso de ortogonalización a los conjuntos

infinitos de vectores, linealmente independientes, 1, x, x2 , ... en el espacio de
polinomios P con producto interior definido por
1
p•q =
∫ −1
p(x)q(x)dx .

La ortogonalización va como sigue:

1 1
e1 = 1 , e1 • e1 =
∫ −1
dx = 2 , e2 = x − α, donde α = 1
2 ∫ −1
xdx = 0 .

Así,
1
e2 = x , e2 • e2 =
∫ −1
x2dx = 2
3
, e3 = x2 − α1 − α2 x, donde
1 1
α1 = 1
2 ∫ −1
x2dx = 1
3
, α2 = 3
2 ∫ −1
x3dx = 0 .

Así,
1
e3 = x2 − 1
3
, e3 • e3 =
∫ −1
(x2 − 13 )2 dx = 8
45
.

Continuando de este modo, se obtiene la sucesión ortogonal

1, x, x2 − 13 , x3 − 3
5
x, x4 − 6
7
x2 + 3
35
, ...
en P. Cuando se presentan con constantes de multiplicación adecuadas, estos
polinomios se convierten en los famosos polinomios de Legendre del Análisis.

4.5. EL ESPACIO DE LAS FUNCIONES CONTINUAS POR

TRAMOS

Dada la importancia de las diferentes series que están por considerarse, reviste
cierto interés expresar la discusión siguiente de tal modo que, de una categoría de
funciones, se abarque tanto como sea posible. Desafortunadamente, la mayoría de las
generalizaciones que van más allá del espacio de funciones continuas, son sumamente
técnicas e inaccesibles en un curso de este nivel. Sin embargo, es posible extender las
conclusiones para incluir funciones continuas por tramos, y ya que bien vale la pena
hacer esta modesta generalización, esta sección se destina a demostrar cómo el
conjunto de todas estas funciones en un intervalo determinado, puede estar dentro de
un espacio euclidiano. Por conveniencia, se empieza por recordar la definición de
continuidad por tramos.

Prof. José Luis Quintero 11

Definición 10. Se dice que una función de valor real f es continua por tramos, en un
intervalo a,b si
(i) f es definida y continua en todos, excepto un número finito, de puntos de a,b y
(ii) los límites f(x0+ ) = lím f(x0 + h) , f(x0− ) = lím f(x0 − h) existen en cada punto x0 en
h → 0+ h → 0−

a,b . (Observar que sólo uno de estos límites es apropiado, si x0 es un punto

extremo de a,b .)

Cuando x0 es un punto de continuidad de f cada uno de estos límites es igual

al valor de f en x0 y, entonces, se tiene que f(x0+ ) = f(x0− ) = f(x0 ) . Con mayor
generalidad, la exigencia de que ambos límites sean finitos en todas partes de a,b ,
implica que las únicas discontinuidades de f son “discontinuidades por salto”. Más aún,
la diferencia f(x0+ ) − f(x0− ) mide la magnitud del salto de la función f en x0 .
1. Si f es continua por tramos en a,b , entonces
b

∫ a
f(x)dx

existe, y es independiente de los valores que f (si los toma) en sus puntos de
discontinuidad. En particular, si f y g son idénticas por todas partes en a,b , salvo
en sus puntos de discontinuidad, entonces
b b

∫ a
f(x)dx =
∫ a
g(x)dx .

2. Si f y g son continuas por tramos en a,b , también lo es su producto f.g.

3. Cada función continua en a,b es continua por tramos.

Dicho esto, ahora se fija la atención en el problema de convertir el conjunto de

funciones continuas por tramos en a,b en un espacio euclidiano. En vista del hecho
de que este conjunto incluye las funciones continuas en a,b solo es razonable exigir
que la construcción sea concebida de tal modo que el espacio euclidiano resultante
tenga a C a,b como un subespacio. Esto, a su vez, sugiere que se defina a f • g por la
fórmula
b
f •g =
∫ a
f(x)g(x)dx .

¿Pero la fórmula anterior produce realmente un producto interior en el conjunto de

funciones continuas por tramos en a,b ? La respuesta desagradable es, NO. Para ver
lo que está equivocado, sea n(x) una función que es cero dondequiera en a,b
excepto en un número finito de puntos.

Se dice que tal función es una función nula, y tiene la molesta propiedad de
que
b

∫ a
n(x)2 dx = 0

Prof. José Luis Quintero 12

a pesar del hecho de que n no es la función cero. Esto, por supuesto, invalida el uso de
la fórmula cuestionada como un producto interior, ya que, por definición, el producto
interior de un vector no cero por sí mismo, no puede ser cero.

Está perfectamente claro, sin embargo, que la dificultad anterior desaparecerá

si se pasa por alto el hecho de que una función nula no es idénticamente cero, y se la
trata como si lo fuera. Pero entonces, para que sea consistente, deben mirarse
también dos funciones continuas por tramos cualesquiera como idénticas, siempre
que difieran solamente en un número finito de puntos, y esto es exactamente lo que se
requiere para hacer que la fórmula cuestionada dé un producto interior.

4.6. CONVERGENCIA DE SUCESIONES

Definición 11. Se dice que una sucesión {xk } = {x1 , x2 ,...} de vectores, en un espacio
euclidiano V converge al vector x en V si, y solamente si, lim xk − x = 0 . En este
k →∞
caso se dice que x es el límite de {xk } , denotándose el hecho de que la sucesión {xk }
converge a x, escribiendo lim {xk } = x .
k →∞

El estudio de convergencia en espacios euclidianos de dimensión finita es

esencialmente el mismo que el estudio de convergencia de sucesiones de números
reales. Sin embargo, en espacios de dimensión infinita, tales como C a,b , esta
situación se vuelve mucho más compleja y, en consecuencia, más interesante, porque
entonces el tipo de convergencia definido antes es totalmente distinto del estudiado en
cálculo. En efecto se verá momentáneamente, que en un espacio de funciones con
producto interior integral, la afirmación
1/2
 b 
∫
2
lí m fk − f = lí m  fk (x) − f(x) dx  =0
k →∞ k →∞  a 
de ninguna manera es la misma que decir que la sucesión {fk } converge a la función f
en cada punto de a,b . En análisis, tal convergencia se conoce como convergencia
media para enfatizar que se calcula por integración, lo que, en un sentido, es un
proceso de promedio generalizado.

{ }
Ejemplo 14. La sucesión de funciones x, x2 , x3 ,... converge en la media en C −1,1

a la función cero, ya que

1/2 1/2
 1   2 
lí m xk − 0 = lí m 
k →∞ k →∞  ∫ −1
x2k dx 

= lí m  
k →∞  2k + 1 
= 0.

Aunque {x, x , x ,...}

2 3
no converge a cero en cada punto en el intervalo −1,1 . En

efecto, en x = 1 la sucesión converge a 1, mientras que en x = −1 no converge (ver

figura 1).

Prof. José Luis Quintero 13

 Figura 1. Sucesión de funciones del ejemplo 14

4.7. BASES EN ESPACIOS EUCLIDIANOS DE DIMENSIÓN

INFINITA

En secciones anteriores se utilizó el proceso de ortogonalización de Gram-

Schmidt para demostrar que cada espacio euclidiano de dimensión finita tiene una
base ortonormal e1 ,..., en y que cada vector en tal espacio puede expresarse
unívocamente en la forma
x = (x • e1 )e1 + ... + (x • en )en . (1)

Pero el proceso de Gram-Schmidt puede aplicarse igualmente bien en un

espacio euclidiano de dimensión infinita donde puede usarse para dar un conjunto
ortonormal infinito e1 , e2 ,... . Este hecho sugiere que una versión ampliada de los
resultados iniciales, podría aceptarse, y las secciones restantes del presente tema se
dedican a realizar esta sugerencia. La forma más obvia para intentar generalizar (1) en
presencia de un conjunto ortonormal e1 , e2 ,... en un espacio de dimensión infinita V,
es reemplazar su segundo miembro por la serie infinita
∞

∑ (x • e )e
k =1
k k . (2)

Sin embargo, en la ausencia de cualquier información adicional, está claro que

no hay razón a priori para suponer que esta serie converge y, mucho menos, que
converge a x* . No obstante, es conveniente tener una notación que exprese el hecho
de que (2) se deduce de x, aunque de modo puramente formal. La usada con más
frecuencia es
∞
x∼ ∑ (x • e )e
k =1
k k ,

Prof. José Luis Quintero 14

usándose el símbolo ∼ (el que, de paso, nada tiene que ver con el usado antes, para
relaciones equivalentes) para enfatizar que las series en cuestión podrían no converger
a x. Por supuesto, en caso de que suceda, se escribe
∞
x= ∑ (x • e )e
k =1
k k ,

y se dice que la serie converge en la media a x. En uno u otro caso, los productos
interiores x • ek son llamados las coordenadas o (generalizando) los coeficientes de
Fourier de x con respecto al conjunto ortonormal e1 , e2 ,...

Está claro que los coeficientes de Fourier de x dependen del conjunto

ortonormal con respecto al cual se calculan. No tan claro, pero igualmente ciero, es
que (1) puede converger en la media a x para un conjunto ortonormal, pero no para
otro. Los ejemplos siguientes ilustran ambos puntos.

Ejemplo 15. Calcule los coeficientes generalizados de Fourier de la función f(x) = x ,

−π ≤ x ≤ π en C −π, π con respecto al conjunto ortonormal
senx sen2x sen3x
, , ,...
π π π
Solución.
En este caso los coeficientes son
1 π
x • ek =
π ∫ −π
xsenkxdx, k = 1,2,...,

y usando integración por partes, se obtiene

2 π
π  , k = 1, 3,5,...,
1 2 π  k
π ∫ −π
xsenkxdx = −
k
cos kπ = 
 2 π
 − k k = 2, 4, 6,...

En consecuencia, los coeficientes de Fourier de x con respecto al conjunto ortonormal

dado son (−1)k −1(2 π / k) , y (2) adopta la forma dada por la expresión
 sen2x sen3x 
x ∼ 2  senx − + − ... .
 2 3 
En próximos temas se verá que esta serie realmente converge en la media a x, de
suerte que se justifica escribir
∞

∑ (−1) k −1 senkx
x=2 .
k =1
k

Ejemplo 16. Si en el ejemplo anterior se reemplaza al conjunto ortonormal por el

conjunto dado por
1 cos x cos 2x
, , ,...,
2π π π
entonces los coeficientes de Fourier son
1 π 1 π

2π − π ∫
xdx = 0 y
π ∫ −π
x cos kxdx = 0, k = 1,2,... .

Prof. José Luis Quintero 15

Pero todas estas integrales son cero, y en consecuencia (2) se convierte en x ∼ 0 .
Ahora queda claro que no se tiene una igualdad, ya que
1/2 1/2
 π   2π3 
x =
 ∫ −π
x2dx 

=
 3 
 
≠ 0.

Observación 2. Para desarrollar f en una serie de funciones ortogonales, es necesario

que no sea ortogonal a cada elemento del conjunto ortogonal. Para evitar este
problema se supondrá, en temas siguientes, que un conjunto ortogonal es completo.
Esto quiere decir que la única función ortogonal a cada miembro del conjunto es la
función cero.

4.8. DESIGUALDAD DE BESSEL. IGUALDAD DE PARSEVAL.

TEOREMA 5. Sea e1 , e2 ,... un conjunto ortonormal de vectores en un espacio

euclidiano V, de dimensión infinita, y sea x un vector arbitrario en V. Entonces
∞

∑ (x • e )
k =1
k
2
≤ x
2

(desigualdad de Bessel).
Más aún, e1 , e2 ,... es una base para V si, y solo si,
∞

∑ (x • e )
k =1
k
2
= x
2

(igualdad de Parseval).
Demostración. La demostración descansa en el cálculo del valor de la expresión
2
n
x− ∑ (x • e )e
k =1
k k

para cualquier x en V y cualquier entero n. Ahora

2
n  n   n 
x− ∑
k =1
(x • ek )ek = x −


∑
k =1
(x • ek )ek  •  x −
 
 
∑ (x • e )e
k =1
k k



n  n   n 
= x•x −2 ∑ (x • e )(x • e ) +  ∑ (x • e )e  •  ∑ (x • e )e
k k j j k k



k =1  j =1  k =1

(Se ha cambiado el índice de la sumatoria en el primer factor del último término, por
conveniencia de cálculo). Pero ya que ek es ortonormal, e j • ek = δ jk , se sigue que
 n   n  n n n

 ∑ (x • e j )e j  • 
  ∑ (x • ek )ek  =
 ∑∑ (x • e j )(x • ek )(e j • ek ) = ∑ (x • e ) .
k
2

 j =1   k =1  j =1 k =1 k =1

Entonces
2
n n
x− ∑
k =1
(x • ek )ek = x
2
− ∑ (x • e )
k =1
k
2
. (3)

Prof. José Luis Quintero 16

Para demostrar la desigualdad de Bessel, es suficiente con observar que
2
n n
0≤ x− ∑
k =1
(x • ek )ek = x
2
− ∑ (x • e )
k =1
k
2

para todo x en V y toda n. En consecuencia

n

∑ (x • e )
k =1
k
2
≤ x
2

para toda n, y las sumas parciales de la serie

∞

∑ (x • e )
k =1
k
2

forman una sucesión no decreciente (acotada) de números reales no negativos. Por un

teorema bien conocido del cálculo se concluye que esta serie converge y que
n

∑ (x • e )
k =1
k
2
≤ x
2
.

Terminando la demostración ahora se supone que e1 , e2 ,... es una base ortonormal

para V. Entonces, se sabe que
∞
x= ∑ (x • e )e
k =1
k k .

Por lo tanto
2
n
lím x −
n →∞ ∑ (x • e )e
k =1
k k = 0,

y se sigue de (3) que

 n 
lím  x
n →∞ 

2
− ∑
k =1
(x • ek )2  = 0 .

Por tanto
∞

∑ (x • e )
k =1
k
2
= x
2
,

que es la igualdad de Parseval. Finalmente, ya que los eslabones de esta cadena de

razonamientos son reversibles, se concluye que e1 , e2 ,... es una base siempre que la
igualdad de Parseval se satisfaga, y la demostración está completa.

4.9. PROBLEMA DE STURM-LIOUVILLE

Ejemplo 17. y ''+ y = 0, y(0) = 0, y(π) = 0 .

Solución.
Para resolver este ejemplo simplemente se aplican las condiciones de frontera a la
solución general c1senx + c2 cos x para deducir que c2 = 0 y que c1 es arbitraria. Así,
y = csenx donde c es una constante arbitraria, es una solución general de este

problema particular.

Prof. José Luis Quintero 17

 Se dirá que las condiciones de frontera dadas, son homogéneas siempre que
γ1 = γ2 = 0 . Las soluciones de un problema con valor en la frontera, que considera un
operador diferencial lineal L : δ → C a,b , están íntimamente ligadas a las soluciones

de la ecuación
Ly = λy , (4)

donde λ es un parámetro desconocido. En este marco, se requiere encontrar todos los

valores de λ para los cuales (4) admite soluciones no triviales en δ , y luego encontrar
las soluciones correspondientes a estas λ . Se puede observar que (4) puede volver a
escribirse (L − λI)y = 0 , donde I indica la transformación identidad que envía cada
función en C a,b  sobre sí misma, o como a2 (x)y ''+ a1(x)y '+ a0 (x) − λ  y = 0 si
2
L = a2 (x)D + a1(x)D + a0 (x) .

Así, para cada valor de λ , (4) es una ecuación diferencial lineal homogénea de
segundo orden, y el conjunto solución de cualquier problema (homogéneo) con valor
en la frontera que considere esta ecuación, es el espacio nulo, en δ , del operador
L − λI .

Ejemplo 18. Resuelve el problema con valor en la frontera

y ''+ λy = 0, y(0) = 0, y(π) = 0 . (5)

Solución.
Aquí, δ es el espacio de todas las funciones dos veces continuamente diferenciables en
0, π  que se anulan en los puntos extremos del intervalo, y L es el operador lineal de

segundo orden −D2 .(El signo menos se ha introducido únicamente para simplificar los
resultados finales. Sin él, los valores adecuados de λ serían negativos.) Se distinguen
tres casos, según que λ = 0, λ < 0 y λ > 0 .
Caso I. λ = 0 . Aquí la ecuación diferencial tiene a c1 + c2 x como su solución general, y
las condiciones de frontera implican que c1 = c2 = 0 . Así, (5) no tiene soluciones no

triviales cuando λ = 0 .
−λx
Caso 2. λ < 0 . Ahora, y = c1e + c2e− −λ x
, y las condiciones de frontera dan otra vez
y ≡ 0.
Caso 3. λ > 0 . Aquí, la solución general de y ''+ λy = 0 es
y = c1sen λ x + c2 cos λ x , (6)
y las condiciones de frontera dan el par de ecuaciones c2 = 0 , c1sen λπ = 0 .

Así, (5) admite soluciones no triviales si, y solamente si, sen λπ = 0 ; esto es,
si, y solamente si, λ asume uno de los valores λn = n2 , donde n = 1, 2,...

Prof. José Luis Quintero 18

 Para cada uno de estos valores de λ , la constante c1 en (6) permanece
arbitraria, siguiéndose que el espacio solución correspondiente a λn es el subespacio
unidimensional de C 0, π generado por la función sen(nx). Los números λn = n2 que

determinan los casos en que (5) tiene soluciones no triviales son llamados los valores
propios para este problema, y cada solución no trivial correspondiente al valor propio
λn es llamada vector propio o función propia perteneciente a λn . Esta terminología se

generalizará luego, para incluir una clase mucho más amplia de problemas y, por el
momento, sólo se pide que observe que cualquier conjunto de funciones propias, tal
como sen(x), sen(2x), sen(3x), …, uno para cada valor propio, es ortogonal en C 0, π .

Definición 12. Un problema de Sturm-Liouville es una ecuación diferencial ordinaria

de segundo orden homogénea de la forma
d
(r(x)y ') + (q(x) + λp(x))y = 0 , a < x < b (7)
dx
sujeta a las condiciones de frontera
α1y(a) + β1y '(a) = 0, α2 y(b) + β2 y '(b) = 0 , (8)
donde p(x), q(x), r(x) y r '(x) son funciones continuas en a,b  . Los coeficientes en las
ecuaciones (7) y (8) se suponen reales e independientes de λ . Además, α1 y β1 no
son cero ambas ni α2 y β2 no son cero ambas.

Resolver un problema de Sturm-Liouville es hallar los valores propios

(autovalores) y funciones propias (autofunciones) no triviales que satisfagan (7) y (8).
Al problema anterior se le suele llamar problema regular o problema normal de
Sturm-Liouville.

TEOREMA 6. (Propiedades del problema regular de Sturm-Liouville).

a. Existe un número infinito de valores propios reales que se pueden ordenar en
forma creciente, λ1 < λ2 < ... < λn < ... de modo que λn → ∞ cuando n → ∞ .

b. A cada valor propio corresponde solo una función propia .

c. Las funciones propias que corresponden a los diversos valores propios son
linealmente independientes.
d. El conjunto de funciones propias que corresponde al conjunto de los valores
propios es ortogonal con respecto a la función peso p(x) en el intervalo a,b  .

Demostración de la propiedad (d).

Sean ym y yn las funciones propias que corresponden a los valores propios λm
y λn , respectivamente. Entonces
d
(r(x)y'm ) + (q(x) + λmp(x))ym = 0 (9)
dx

Prof. José Luis Quintero 19

 d
(r(x)yn' ) + (q(x) + λnp(x))yn = 0 (10)
dx
Al multiplicar la ecuación (9) por yn y la ecuación (10) por ym y restar los
d d
resultados se obtiene (λm − λn )p(x)ymyn = ym (r(x)yn' ) − yn (r(x)y'm ) .
dx dx
Se integra por partes este resultado, desde x = a hasta x = b , y se llega a
b
(λm − λn )
∫ a
p(x)ymyndx = r(b) ym (b)yn' (b) − yn (b)y'm (b) − r(a)  ym (a)yn' (a) − yn(a)y'm (a)
   
Ahora bien, las funciones propias ym y yn deben satisfacer, ambas, las

condiciones en la frontera (8). Por ejemplo, de acuerdo a una de las condiciones se

tiene α1ym (a) + β1y'm (a) = 0, α1yn (a) + β1yn' (a) = 0 . Para que α1 y β1 satisfagan este

sistema, no nulos simultáneamente, el determinante de los coeficientes debe ser igual

a cero: ym (a)yn' (a) − yn (a)y'm (a) = 0 . Se aplica un argumento semejante a la otra
condición de (8) y se obtiene ym (b)yn' (b) − yn (b)y'm (b) = 0 . Como los dos miembros del

lado derecho de la ecuación son cero, se ha establecido la relación de ortogonalidad

b

∫ a
p(x)ym (x)yn(x)dx = 0, λm ≠ λn .

Ejemplo 19. Resuelva el problema de valores en la frontera

y ''+ λy = 0, y(0) = 0, y(1) + y '(1) = 0 . (11)

Solución.
Se puede comprobar que para λ < 0 y para λ = 0 , el problema planteado por la
ecuación (11) sólo posee la solución trivial y = 0 . Cuando λ > 0 , la solución general de
la ecuación diferencial es y = c1 cos λx + c2sen λ x . Ahora bien, la condición y(0) = 0
implica que en esta solución, c1 = 0 . Cuando y = c2sen λ x , se satisface la segunda
condición de frontera y(1) + y '(1) = 0 , siempre que c2sen λ + c2 λ cos λ = 0 . Al
escoger c2 ≠ 0 , se ve que esta última ecuación es equivalente a tg λ = − λ . Si se
hace x = λ , los valores propios del problema (11) son, entonces, λn = xn2 , donde
xn , n = 1, 2, 3,... son las raíces positivas consecutivas. En general, las funciones propias

{ }
del problema son sen( λn x , n = 1, 2,3,...

Al identificar r(x) = 1, q(x) = 0, p(x) = 1, α1 = 1, β1 = 0, α2 = 1, β2 = 1 se ve que

la ecuación (11) es un problema regular de Sturm-Liouville. Así, {sen( λn x , }
n = 1,2, 3,... es un conjunto ortogonal con respecto a la función peso p(x) = 1 en el
intervalo 0,1 .

A continuación se consideran cierto número de problemas de Sturm-Liouville

que se encontrarán repetidamente en trabajos posteriores y que, por conveniencia de

Prof. José Luis Quintero 20

referencia, se resuelven aquí, de una vez por todas. En cada caso se limitará a
encontrar valores propios y funciones propias, y se pospone cualquier discusión de
cómo usar esta información para resolver problemas específicos con valor en la
frontera que involucren el operador considerado.

Ejemplo 20. Resuelva el sistema de Sturm-Liouville

y ''+ λy = 0, y(0) = 0, y(L) = 0 , 0 < x < L

Solución.
Esta es una variante de un problema que se ha considerado muchas veces y, en esta
ocasión, se encuentra que las constantes dadas por
n2 π2
λn = 2 , n = 1, 2,...
L
y las funciones
 nπx 
yn (x) = sen   , n = 1,2,...,
 L 
son considerados conjuntos completos de valores propios y funciones propias. Se verá
en detalle su resolución: Se tienen tres opciones:
λ < 0 → λ = −k2 : y ''− k2 y = 0 ⇒ m2 − k2 = 0 ⇒ m = ±k
y(x) = c1ekx + c2 e−kx , y(0) = c1 + c2 = 0 ⇒ c2 = −c1 .

Sustituyendo:
y(x) = c1(ekx − e−kx ) , y(L) = c1(ekL − e−kL ) = 0 ⇒ c1 = 0, c2 = 0 ⇒ y(x) = 0

Conduce a solución trivial.

λ = 0 : y '' = 0 ⇒ y(x) = c1x + c2 ⇒ y(0) = c2 = 0 ⇒ y(x) = c1x
y(L) = c1L = 0 , pero L ≠ 0 , entonces conduce a solución trivial.
λ > 0 : λ = k2 : y ''+ k2 y = 0 ⇒ y(x) = c1 cos(kx) + c2sen(kx)
Al evaluar y(0) = c1 = 0 ⇒ y(x) = c2sen(kx)
nπ
y(L) = c2sen(kL) = 0 ⇒ sen(kL) = 0 ⇒ kL = nπ ⇒ k = .
L
Al sustituir se tiene:
 nπ 
yn = c2sen  x
 L 
(Solución no trivial).
Como
2
nπ  nπ 
k= ⇒ λn =   Autovalores.
L  L 

Ejemplo 21. Resuelva el sistema de Sturm-Liouville y ''+ λy = 0, y '(0) = 0, y '(L) = 0 .

Solución.

Prof. José Luis Quintero 21

Un cálculo similar en todos los aspectos al que se usó para resolver el ejemplo
anterior, revela que los valores propios para este problema son las constantes no
negativas
n2 π2
λn = , n = 0,1,2,...
L2
y que
 nπx 
yn (x) = cos   , n = 0,1, 2,...,
 L 
constituye un conjunto completo de funciones propias.

Ejemplo 22. Resuelva el sistema de Sturm-Liouville

y ''+ λy = 0, y(0) = 0, hy(L) + y '(L) = 0,

considerando que h y L son constantes positivas.

Solución.
Como es usual, se estudia por casos, dependiendo del signo algebraico de λ , y
nuevamente se encuentra que no hay valores propios ≤ 0 . Por otra parte, cuando
λ > 0 , y = c1sen λ x + c2 cos λ x, y la primera condición de frontera implica que c2 = 0 .
Así, queda pendiente (si es que queda) elegir λ en forma tal que la función
y = c1sen( λ x) , con c1 ≠ 0 , satisfaga la ecuación hy(L) + y '(L) = 0 . Esto, a su vez,
implica que λ debe escogerse en forma tal que
λL
sen( λL) = − cos( λL) ; (12)
hL
una ecuación que vuelve a escribirse como
µ
t gµ = − (13)
hL
haciendo µ = λL . Aunque es imposible resolver (12) explícitamente para λ , sus

soluciones pueden visualizarse conforme surgen, a través de (13), en los puntos de

intersección de las gráficas de las funciones tgµ y −µ / hL.

Hay una infinidad de tales puntos localizados simétricamente con respecto al

origen. Así, el problema considerado tiene un número infinito de valores propios
µ2
positivos λn = 2n , n = 1,2,... . A continuación se muestra la gráfica (figura 2) para el
L
caso h = 1, L = 1, µ = x ):

Prof. José Luis Quintero 22

 Figura 2. Problema del Sturm-Liouville del ejemplo 22

con λn < λn +1 para toda n, y lím λn = ∞ . (También observe, de la geometría del caso,
n →∞
se tiene que lím(λn +1 − λn ) = π ). Las funciones
n →∞

µ x
yn (x) = sen  n  = sen( λx), n = 1,2,...
 L 
constituyen un conjunto completo de funciones propias para este problema.

Ejemplo 23. Considere el problema con valores en la frontera

y ''+ 2y '+ (1 − λ)y = 0
 .
 y(0) = 0, y(1) = 0
Halle los valores propios y las funciones propias y determine una relación de
ortogonalidad entre las soluciones.
Solución.
Si λ < 0 se tiene que la solución es y(x) = e− x c1 cos(kx) + c2sen(kx) , donde k2 = −λ .
Con las condiciones iniciales se tiene 0 = c1 , 0 = e−1c2sen(k) . Si c2 ≠ 0, k = nπ . Por lo
tanto λ = −n2 π2 y yn (x) = e− x sen(nπx) . Sea la función de peso p(x) = e2x , se tiene la
siguiente relación de ortogonalidad entre las soluciones:
1

∫ 0
e2x e−x sen(mπx)e−x sen(nπx)dx = 0 .

4.10. FORMA AUTOADJUNTA DE PROBLEMAS DE

STURM-LIOUVILLE

Si los coeficientes son continuos y a(x) ≠ 0 para todo x en algún intervalo,

entonces cualquier ecuación diferencial de segundo orden

Prof. José Luis Quintero 23

 a(x)y ''+ b(x)y '+ (c(x) + λd(x))y = 0 (14)

se puede llevar a la forma autoadjunta, ecuación (7), si se multiplica por el factor

integrante
1
e∫
(b(x) a(x))dx
.
a(x)
Para comprobarlo, se observa que la ecuación diferencial
b(x) ∫ (b / a)dx  c(x) ∫ (b / a)dx d(x) ∫ (b / a)dx 
e∫
(b / a)dx
y ''+ e y '+  e +λ e y = 0
a(x)  a(x) a(x) 
es la misma que
d e ∫ (b / a)dx y ' +  c(x) e∫ (b / a)dx + λ d(x) e ∫ (b / a)dx  y = 0 . (15)
dx    a(x) a(x)



Entonces
c(x) ∫ (b / a)dx d(x) ∫ (b / a)dx
r(x) = e∫
(b / a)dx
, q(x) = e , p(x) = e .
a(x) a(x)

Desde luego que no es necesario llevar a una ecuación diferencial de segundo

orden como la (14) a la forma autoadjunta (7), para resolverla. Se usará la forma de
la ecuación (15) para determinar la función peso p(x) que se necesita en la relación de
ortogonalidad.

Ejemplo 24. Encuentre los valores propios y las funciones propias para el problema de
Sturm-Liouville y ''+ 4y '+ (4 − 9λ)y = 0, y(0) = 0, y(a) = 0 .

Solución.
En primer lugar, se escribirá en forma autoadjunta como
d  4x dy 
e + (4e4x − λ(9e4x ))y = 0, y(0) = 0, y(a) = 0 .
dx  dx 

Así, las funciones propias correspondientes a diferentes valores propios para

este problema serán mutuamente ortogonales en el espacio euclidiano C 0, a con
producto interior calculado con relación a la función peso p(x) = 9e4x .
Caso I. λ > 0 . Aquí la solución general de y ''+ 4y '+ (4 − 9λ)y = 0 viene dada por
y = c1e( −2 + 3 λ )x
+ c2e(−2 − 3 λ )x
,

y las condiciones de frontera implican que

c1 + c2 = 0, c1e(−2 + 3 λ )a
+ c2e(−2 − 3 λ )a
= 0.
Así c1 = c2 = 0 y el problema no tiene valores propios positivos.
Caso II. λ = 0 . Ahora la solución general de la ecuación es (c1 + c2 x)e−2a , y
nuevamente se encuentra que c1 = c2 = 0 .

Caso III. λ < 0 . Aquí

y = e−2x (c1sen(3 −λ x) + c2 cos(3 −λ x))

Prof. José Luis Quintero 24

y la exigencia de que y(0) = y(a) = 0 da c2 = 0, c1sen −λa = 0 . En consecuencia el

problema tiene soluciones no triviales si, y solamente si, λ satisface la ecuación

sen(3 −λa) = 0 , siguiéndose que los valores propios para este problema son
π2n2
λn = − , n = 1, 2,...
9a2
Para encontrar un conjunto de funciones propias correspondientes, ahora se hace
c1 = 1, c2 = 0 y λ = λn y se obtienen
 nπx 
ϕn (x) = e−2x sen   , n = 1, 2,... .
 a 

Ejemplo 25. Encuentre los valores propios y las funciones propias del problema de
Sturm-Liouville dado por y ''+ ky = 0, y(0) = y(2π), y '(0) = y '(2π) .
Solución.
Caso 1. k = 0 . y ''(x) = 0 ⇒ y(x) = c1 + c2 x .
Autovalor k = 0 . Autofunción y = c1 c1 ≠ 0
y(0) = y(2π) ⇒ c1 = c1 + 2πc2 ⇒ c2 = 0 ⇒ y(x) = c1 , y '(0) = y '(2π) ⇒ 0 = 0 ⇒ y(x) = c1 .
Caso 2. k = −λ2 , λ ≠ 0 . y ''− λ2 y = 0 ⇒ y(x) = c3 cosh(λx) + c4senh(λx) .
y(0) = y(2π) ⇒ c3 = c3 cosh(2πλ) + c4senh(2πλ) ,
y '(0) = y '(2π) ⇒ c4 λ = c3λsenh(2πλ) + c4 λ cos h(2πλ)
c3 (cosh(2πλ) − 1) + c4senh(2πλ) = 0 , c3senh(2πλ) + c4 (cos h(2πλ) − 1) = 0 . Por lo tanto se
tiene un sistema de ecuaciones lineales homogéneo. El valor del determinante del
sistema corresponde a
(cosh(2πλ) − 1)2 − sen2h(2πλ) = cos2 h(2πλ) − 2 cosh(2πλ) + 1 − sen2h(2πλ)
= 2 − 2 cosh(2πλ) ≠ 0
En consecuencia c3 = c4 = 0 es la única solución y se tiene y(x) = 0 .
Caso 2 (Alternativo). k = −λ2 , λ ≠ 0 . y ''− λ2 y = 0 ⇒ y(x) = c3eλx + c4e−λx .
y(0) = y(2π) ⇒ c3 + c4 = c3e2 πλ + c4e−2πλ , y '(0) = y '(2π) ⇒ λ(c3 − c4 ) = λ(c3e2πλ − c4e−2πλ )
c3 (e2 πλ − 1) + c4 (e−2πλ − 1) = 0 , c3 (e2πλ − 1) − c4 (e−2πλ − 1) = 0 . Por lo tanto se tiene un
sistema de ecuaciones lineales homogéneo. El valor del determinante del sistema
corresponde a −2(e2πλ − 1)(e−2πλ − 1) = −2(1 − e2 πλ − e−2πλ + 1) = −2(2 − (e2πλ + e−2πλ )) ≠ 0 .
En consecuencia c3 = c4 = 0 es la única solución y se tiene y(x) = 0 .
Caso 3. k = λ2 , λ ≠ 0 . y ''+ λ2 y = 0 ⇒ y(x) = c5 cos(λx) + c6 sen(λx) .
y(0) = y(2π) ⇒ c5 = c5 cos(2πλ) + c6 sen(2πλ) , y '(0) = y '(2π)
⇒ λc6 = λ(c6 cos(2πλ) − c5sen(2πλ))
c5 (cos(2πλ) − 1) + c6 sen(2πλ) = 0 , c5sen(2πλ) + c6 (1 − cos(2πλ)) = 0 . Por lo tanto se tiene
un sistema de ecuaciones lineales homogéneo. El valor del determinante del sistema
corresponde a
−(cos(2πλ) − 1)2 − sen2 (2πλ) = − cos2 (2πλ) + 2 cos(2πλ) − 1 − sen2 (2πλ) = 2(cos(2πλ) − 1) .
cos(2πλ) − 1 = 0 ⇒ cos(2πλ) = 1 ⇒ 2πλ = 2nπ ⇒ λ = n . Por tanto los autovalores vienen
dados por kn = λ2 = n2 n = 1,... y las autofunciones yn (x) = c5 cos(nx) + c6 sen(nx) .
SOLUCIÓN FINAL: k = n2 n = 0,... y las autofunciones yn (x) = c5 cos(nx) + c6 sen(nx) .

Prof. José Luis Quintero 25

4.11. FUNCIONES PARES Y FUNCIONES IMPARES

La tarea de calcular las integrales que surgen en el estudio de series

ortogonales, puede simplificarse a menudo, al sacar provecho de la simetría de las
funciones involucradas. Esta técnica, generalmente se formaliza por la introducción de
las nociones de funciones pares e impares, como sigue:

Definición 13. Se dice que una función f, definida en un intervalo centrado en el

origen, es par si f(−x) = f(x) para toda x en el dominio de f, e impar si f(−x) = −f(x) .

Esto, por supuesto, es justamente otra forma de decir que una función es par si
su gráfica es simétrica con respecto al eje vertical, e impar si su gráfica es simétrica
con respecto al origen. Así, para valores enteros de n, xn es par, si n es par; e impar,
si n es impar. La importancia de las funciones pares e impares para este trabajo, nace
de las igualdades
a a

∫ −a
f(x)dx = 2
∫ 0
f(x)dx

siempre que f sea par e integrable, y

a

∫ −a
f(x)dx = 0

siempre que f sea impar e integrable. Una propiedad un tanto menos obvia de las
funciones pares e impares, se establece en el siguiente

TEOREMA 7. Cada función en el intervalo −a, a puede expresarse justamente en

forma de una suma de una función par y una función impar.
Demostración. Sea f una función arbitraria en −a, a , y sea
f(x) + f(−x) f(x) − f(−x)
fE (x) = , f0 (x) = .
2 2

Es trivial verificar que fE es par, f0 impar, y que f = fE + f0 . Así, f tiene cuando

menos una descomposición de la forma deseada, y queda por demostrar que es la
única. Con este fin, se supone que también se tiene f = gE + g0 , con gE par y g0 impar.
Entonces, fE + f0 = gE + g0 , y fE − gE = g0 − f0 . Pero la diferencia de dos funciones pares
es par, y la diferencia de dos funciones impares es impar. Así, la función definida por la
igualdad anterior es simultáneamente par e impar, y así, debe ser la función cero. En
otras palabras, fE − gE = g0 − f0 = 0 , de donde se sigue que gE = fE y g0 = f0 , como se
deseaba. Las funciones fE y f0 son conocidas, respectivamente, como las partes par e
impar de f.

Ejemplo 26. Si f(x) = ex , entonces

ex + e− x ex − e − x
fE (x) = y f0 (x) = .
2 2
Así, las partes par e impar de la función exponencial son el coseno hiperbólico y el
seno hiperbólico, respectivamente.

Prof. José Luis Quintero 26

4.12. SERIES DE FOURIER

En esta sección se inicia el estudio de series ortogonales considerando el

desarrollo en serie relativo a las funciones 1, cos x, senx, cos(2x), sen(2x),... . Ya se
ha visto que estas funciones son mutuamente ortogonales en PC −π, π , y se probará
en forma breve que también son una base. Garantizando lo cierto de este hecho,
puede utilizarse la fórmula
∞
x • ek
x= ∑
k =1 ek
2
ek

para expresar cualquier función continua por tramos f en el intervalo −π, π en la
forma
f •1
∞  f • cos(kx) f • sen(kx) 
f(x) =
1
2
+ ∑ 

k =1  cos(kx)
2
cos(kx) +
sen(kx)
2
sen(kx) (media)


donde la notación “(media)” indica que la serie en cuestión converge en la media a f.
Pero ya que
π π π

∫ ∫ ∫
2 2 2
1 = dx = 2π , cos(kx) = cos2 (kx)dx = π , sen(kx) = sen2 (kx)dx = π ,
−π −π −π

entonces se tiene que

∞
a0
f(x) =
2
+ ∑ (a
k =1
k cos(kx) + bk sen(kx) (media) , (16)

donde
1 π 1 π
ak =
π ∫ −π
f(x) cos(kx)dx , k = 0,1,2,... , bk =
π ∫ −π
f(x)sen(kx)dx , k = 1,2,...

Esta representación particular es conocida como su desarrollo en serie de

Fourier en el intervalo −π, π , y las ak y bk son llamadas los coeficientes de
Fourier. Una vez más se enfatiza que (16) debe leerse como una afirmación de que la
serie en cuestión converge en la media a f, no que converja punto por punto en el
sentido de que
∞
a0
f(x0 ) =
2
+ ∑ (a
k =1
k cos(kx0 ) + bk sen(kx0 ))

para toda x0 en −π, π . En efecto ya que el valor de f en x0 puede cambiarse

arbitrariamente, sin cambiar los valores de sus coeficientes de Fourier, esto sería
demasiado esperar. Pero, sorprendentemente, siempre que f se comporte
razonablemente bien, converge a f(x) para toda x.

Ejemplo 27. Encuentre el desarrollo en serie de Fourier, de la función

−1 −π < x < 0
f(x) =  .
1 0 < x < π
Solución.

Prof. José Luis Quintero 27

En este caso f es una función impar en −π, π . En consecuencia, así es f(x)cos(kx), de
modo que ak = 0 , para toda k. Por otra parte, f(x)sen(kx) es par y, por lo tanto, se
tiene
4
2 π 2  k = 1, 3,5,...
bk =
π ∫ 0
sen(kx)dx =
kπ
(1 − cos(kπ)) = kπ
 0 k = 2, 4, 6,...
.

Por tanto, el desarrollo en serie de Fourier de f es

∞
4  4 sen(2k − 1)x
∑
sen(3x) sen(5x)
f(x) =  sen(x) + + + ... =
π 3 5  π k =1 2k − 1
A continuación se presentan algunas gráficas (ver figuras 3, 4 y 5):
4sen(x)
y=
π

Figura 3. Primer desarrollo de Fourier para el ejemplo 27

4 sen(3x) sen(5x) 
y=  sen(x) + + 
π 3 5 

Figura 4. Segundo desarrollo de Fourier para el ejemplo 27

Prof. José Luis Quintero 28

 4 sen(3x) sen(5x) sen(7x) sen(9x) 
y=  sen(x) + + + + 
π 3 5 7 9 

Figura 5. Tercer desarrollo de Fourier para el ejemplo 27

TEOREMA 8. Sea f una función suave por tramos en PC −π, π , por lo que se entiende
que f tiene una primera derivada continua por tramos en −π, π . Entonces el desarrollo
en serie de Fourier para f converge punto por punto dondequiera en −π, π , y tiene el
valor
f(x0+ ) + f(x0− )
2
en cada punto x0 en el interior del intervalo, y
f(−π+ ) + f(π− )
2
en ±π .

Cuando se aplican estos resultados al ejemplo anterior, permiten afirmar que la

serie
4 sen(3x) sen(5x) 
sen(x) + + + ...
π  3 5 
converge punto por punto en el intervalo −π, π a
−1 si − π < x < 0,
0 si x = −π, 0, π .
1 si 0<x<π

Así, por ejemplo, cuando x = π 2 , el valor de la serie es 1, y entonces

4 1 1 1 
1= 1 − + − + ... ,
π  3 5 7 
o

Prof. José Luis Quintero 29

 π 1 1 1
= 1 − + − + ...
4 3 5 7

TEOREMA 9. El desarrollo en serie de Fourier de una función f suave por tramos en

PC −π, π converge punto por punto en toda la recta real. Más aún, si F indica la
extensión periódica de f, entonces el valor de la serie es F(x0 ) cuando x0 es un punto
de continuidad de F, y
F(x0+ ) + F(x0− )
2
cuando x0 es una discontinuidad por salto de F.

Ejemplo 28. Halle el desarrollo en serie de Fourier de la función f(x) = x , − π < x < π .
Solución.
En este caso, f es una función par en −π, π . En consecuencia, bk = 0 para toda k,
mientras que, para k ≠ 0 ,
π
2 π 2  xsen(kx) 1 π 
ak =
π ∫ 0
x cos(kx)dx =
π


k 0
−
k ∫0
sen(kx)dx 

π  4
2 2 − k = 1, 3,5,...,
= cos(kx) = (cos(kπ) − 1) =  πk2
πk 2 πk 2

0
 0 k = 2, 4, 6,...
Finalmente, cuando k = 0 , se tiene
2 π
a0 =
π ∫ 0
xdx = π ,

siguiendo que
∞
π 4  π 4 cos(2k − 1)x
∑
cos(3x) cos(5x)
x = − cos(x) + + + ... = −
2 π  32
52
 2 π k =1
(2k − 1)2
Algunas gráficas (ver figuras 6 y 7):
π 4 cos(3x) 
y = −  cos(x) + 
2 π 9 

Figura 6. Primer desarrollo de Fourier para el ejemplo 28

Prof. José Luis Quintero 30

 π 4 cos(3x) cos(5x) cos(7x) 
y= −  cos(x) + + +
2 π 9 25 49 

Figura 7. Segundo desarrollo de Fourier para el ejemplo 28

Ejemplo 29. Sea g la función en PC −π, π definida por

0 −π < x < 0
g(x) =  .
1 0 < x < π
Entonces, con f como en el ejemplo 27,
1 1 1
g = (1 + f) = + f,
2 2 2
y se concluye que el desarrollo en serie de Fourier de g es
∞
sen(2k − 1)x
∑
1 2
f(x) = + .
2 π k =1 2k − 1

La moraleja de este ejemplo es que una serie de Fourier puede, algunas veces,
encontrarse sin recurrir a la integración. Otra forma de expresar la serie sería
∞
1 − (−1)k
∑
1 1
f(x) = + sen(kx)
2 π k =1 k
1 2
(ver figuras 8, 9, 10 y 11) . Gráficas: k = 1 ⇒ s1(x) = + sen(x)
2 π

Figura 8. Primer desarrollo de Fourier para el ejemplo 29

Prof. José Luis Quintero 31

 1 2 sen(3x) 
k = 3 ⇒ s3 (x) = + sen(x) + 
2 π 3 

Figura 9. Segundo desarrollo de Fourier para el ejemplo 29

1 2 sen(3x) sen(5x) sen(7x) 

k = 7 ⇒ s7 (x) = + sen(x) + + + 
2 π 3 5 7 

Figura 10. Tercer desarrollo de Fourier para el ejemplo 29

1 2 sen(3x) sen(5x) sen(7x) sen(9x) sen(11x) 

s11(x) = + sen(x) + + + + + 
2 π 3 5 7 9 11 

Prof. José Luis Quintero 32

 Figura 11. Cuarto desarrollo de Fourier para el ejemplo 29

4.13. SERIES DE SENOS Y COSENOS

En los ejemplos anteriores se abusó del hecho de que las funciones

consideradas eran pares o impares, para simplificar la tarea de encontrar su desarrollo
en serie de Fourier. Esta técnica puede ser explotada más frecuentemente de lo que se
pudiera esperar, y es de suma importancia para sacarla al descubierto.

Específicamente, si f es una función par en PC −π, π , entonces, para todos los
valores de k, f(x)cos(kx) es par, y f(x)sen(kx) es impar. Así, se tiene que

π π π

∫ −π
f(x) cos(kx)dx = 2
∫ 0
f(x) cos(kx)dx,
∫ −π
f(x)sen(kx)dx = 0,

siguiéndose que el desarrollo en serie de Fourier de una función par en PC −π, π
incluye solamente términos cosenos y pueden calcularse de acuerdo con la fórmula
∞
a0
f(x) =
2
+ ∑ a cos(kx) ,
k =1
k

donde
2 π
ak =
π ∫ 0
f(x) cos(kx)dx .

Un argumento similar muestra que el desarrollo en serie de Fourier de una función par
en PC −π, π incluye solamente términos senos, y se calcula con
∞
f(x) = ∑ b sen(kx) ,
k =1
k

donde
2 π
bk =
π ∫ 0
f(x)sen(kx)dx .

Prof. José Luis Quintero 33

 En las aplicaciones de la teoría de las series de Fourier, frecuentemente se
necesita obtener el desarrollo en serie para una función f continua por tramos, que
está definida solamente en el intervalo 0, π . Una forma de hacerlo es, extendiendo f
en todo el intervalo 0, π (donde con esto quiere decirse que una F est´s definida en
−π, π de tal modo, que F coincide con f en 0, π ), y entonces F se desarrolla como
una serie de Fourier.

Lo esencial de este método tiene relación con la forma en que f es extendida a

−π, π . Esto por supuesto, puede hacerse en la forma que sea (en tanto que la función
resultante pertenezca a PC −π, π ), pero las dos extensiones son las más convenientes
e importantes. La primera es la llamada extensión par de f, indicada por Ef , y definida
por
 f(x) 0 ≤ x ≤ π
Ef (x) = 
f(−x) −π ≤ x < 0
mientras la segunda es la extensión impar de f, indicada por Of , y definida por
 f(x) 0≤x≤π
Of (x) =  .
−f(−x) −π ≤ x < 0
Se tienen entonces los desarrollos
∞ π
a0
∑a
2
Ef (x) =
2
+
k =1
k cos(kx), ak =
π ∫ 0
f(x) cos(kx)dx

y
∞ π
∑ b sen(kx),
2
Of (x) =
k =1
k bk =
π ∫ 0
f(x)sen(kx)dx .

Estas series son llamadas respectivamente, los desarrollos en serie de

Fourier de coseno y de seno de f. Un término, un tanto extraviado “desarrollo semi-
intervalo” también es usado en este contexto.

Ejemplo 30. Halle el desarrollo de Fourier de

 0 −π < x ≤ 0
f(x) =  2
x 0≤x<π
y pruebe, usando el desarrollo anterior, que
∞
(−1)n +1 π2
∑
n =1
n2
=
12
.

Solución.
Sea
∞
a0
f(x) =
2
+ ∑ (a
n =1
n cos(nx) + bnsen(nx)) .

π
1 π 1 x3 π2
a0 =
π ∫ 0
x2dx = .
π 3
0
=
3
.

Prof. José Luis Quintero 34

  π 
1 π 1  x2sen(nx) 2 π

∫ [Link](nx)dx 
∫
2
an = x cos(nx)dx = −
π 0 π n 0 n 
 0 
 2 π π
1 x sen(nx) 2  [Link](nx) sen(nx)  
=  − − +  
π n n n n2 0
 0 
2 π 2 2  π 2 2
n impar : =   = − 2 . n par : = −  −  = 2 . Por tan to : an = (−1)n 2 .
πn  n  n πn  n  n n
 π π
π 1  x2 cos(nx) 2 x
1 1 
∫ +  sen(nx) + 2 cos(nx) 
2
bn = x sen(nx)dx = −
π 0 π n n n n 0 
 0 
1  π2 2  2   1  2 4 1  π2  π2
n impar : =  + − 2   =  π − 2 n par : = −  = −
π  n n  n   πn  n  π  n  πn

(−1)n +1  2 2(1 − (−1)n )  (−1)n +1 π 2((−1)n − 1)

Por tan to : bn = π − = +
πn  n2  n πn3
De modo que:
∞
π2   (−1)n +1 π 2((−1)n − 1)  
∑ (−1) n 2
f(x) = + cos(nx) +  +  sen(nx).
6 n2  n πn3 
n =1   
Si se hace x = 0 se tiene:
∞ ∞ ∞
π2 π2 (−1)n +1 π2
∑ ∑ ∑
2 2
0= + (−1)n ⇒− (−1)n = ⇒ = .
6 n =1
n2 n =1
n2 6 n =1
n2 12

Ejemplo 31. Encuentre la serie de Fourier de

0 −π < x < 0
f(x) =  .
1 0 ≤ x < π

∞
1 − (−1)n
∑
1 1
Rta. f(x) = + sen(nx) (ver figuras 12, 13, 14 y 15)
2 π n =1 n
1 2
n = 1 ⇒ s1(x) = + sen(x)
2 π

Figura 12. Primer desarrollo de Fourier para el ejemplo 31

Prof. José Luis Quintero 35

 1 2 sen(3x) 
n = 3 ⇒ s3 (x) = + sen(x) + 
2 π 3 

Figura 13. Segundo desarrollo de Fourier para el ejemplo 31

1 2 sen(3x) sen(5x) sen(7x) 

n = 7 ⇒ s7 (x) = + sen(x) + + + 
2 π 3 5 7 

Figura 14. Tercer desarrollo de Fourier para el ejemplo 31

1 2 sen(3x) sen(5x) sen(7x) sen(9x) sen(11x) 

s11(x) = + sen(x) + + + + + 
2 π 3 5 7 9 11 

Prof. José Luis Quintero 36

 Figura 15. Cuarto desarrollo de Fourier para el ejemplo 31

Ejemplo 32. Encuentre la serie de Fourier de

0 −π < x < 0
f(x) = 
x 0 ≤ x < π
y usela para calcular el valor de la serie
∞
 (−1)k 
∑
1
 + .
k =1 
(2k − 1)2 (2k − 1) 
Solución.
π
1 π x2 π2 π
a0 =
π ∫ 0
xdx =
2π
0
= = .
2π 2

1 π (−1) − 1k
1 π (−1)k +1
ak =
π ∫ 0
x cos(kx)dx =
πk 2
. bk =
π ∫ 0
xsen(kx)dx =
k
.

Se tiene entonces que

∞
π  (−1)k − 1 (−1)k +1 
f(x) =
4
+ ∑ 
k =1 
πk 2
cos(kx) +
k
sen(kx) 


Si x = 0 :
∞ ∞
π π2
∑ ∑ (2k − 1)
2 1 1
− =− ⇒ =
4 π k =1 (2k − 1)2 8 k =1
2

π
Si x = :
2
∞ ∞
π 2 1 π  (−1)k +1 π 
4
=− ∑
π k =1 (2k − 1)2
cos  (2k − 1)  +
 2 
∑
k =1
k
sen  k 
2 
∞ ∞ ∞
π (−1)k +1 π  π (−1)k +1 π (−1)k
4
= ∑
k =1
k
sen  k  ⇒ =
 2  4 ∑
k =1
2k − 1
⇒− =
4 ∑
k =1
2k − 1
Por lo tanto:
∞
 (−1)k  π2 π
∑  (2k − 1)
1
+ = − .
k =1
2 (2k − 1)  8 4

Prof. José Luis Quintero 37

Ejemplo 33. Desarrolle
 0, −π < x < 0
f(x) = 
 π − x, 0<x<π
en serie de Fourier.
Solución.
Con p = π , se tiene que
π
1 π 1 0 π  1 x2  π
a0 =
π ∫ −π
f(x)dx = 
π ∫ −π
0dx +
∫ 0
(π − x)dx  =  πx −
 π 
 = ,
2 
0
2

1 π 1 0 π 
an =
π ∫ −π
f(x) cos(nx)dx = 
π ∫ −π
odx +
∫ 0
(π − x) cos(nx)dx 

π π
1 sen(nx) 1 π  1 cos(nx) − cos(nπ) + 1 1 − (−1)n
= (π − x)
π

n 0
+
n ∫ 0
sen(nx)dx  = −
 nπ n 0
=
n2 π
=
n2 π
En forma análoga, se obtiene que
1 π 1
bn =
π ∫ 0
(π − x)sen(nx)dx =
n
.

Por lo tanto,
∞
π  1 − (−1)n 
∑ 
1
f(x) = + cos(nx) + sen(nx) 
4 n π 2 n 
n =1 

4.14. CAMBIO DE INTERVALO

Hasta ahora se ha tratado exclusivamente con funciones en los intervalos −π, π
y 0, π . Para muchas finalidades, sin embargo, este marco es demasiado restrictivo, y
ahora se propone generalizar los resultados a un intervalo arbitrario a,b . Pero más
bien empezar desde luego con el caso más general, resultará más sencillo si primero
se consideran intervalos de la forma −p,p  y sus espacios euclidianos asociados
PC −p,p . Porque así, la situación puede resolverse mejor. En realidad, casi es obvio
que las funciones
 πx   πx   2πx   2πx 
1, cos   , sen   , cos   , sen   , ...
 p   p   p   p 
son mutuamente ortogonales en PC −p,p . Más aún, justamente como en el caso
donde p = π , puede demostrarse que estas funciones son una base para este espacio
y, en consecuencia, que sus series ortogonales asociadas (las que, de paso, son aún
llamadas series de Fourier) convergen en la media. Y finalmente, teniendo en cuenta la
longitud del intervalo, todas las observaciones anteriores referentes a la convergencia
por puntos, son válidas en este caso. Entonces
∞
a0   kπx   kπx  
f(x) =
2
+ ∑  a
k =1
k cos 
 p
 + bk sen 
  p 
 ,

donde
1 p  kπx  1 p  kπx 
ak =
p ∫ −p
f(x) cos 
 p
 dx , bk =
 p ∫ −p
f(x)sen 
 p 
 dx

Prof. José Luis Quintero 38

para toda k. Y con esto se ha logrado.

La discusión anterior puede adaptarse fácilmente para manipular el espacio

euclidiano PC a,b  . En efecto, si se establece 2p = b − a , de manera que
a,b = a, a + 2p , las funciones anteriores también son una base para PC a, a + 2p .
Esto inmediatamente lleva a las fórmulas para calcular el desarrollo en serie de Fourier
de una función f en PC a,b  de la forma:
∞
a0   2kπx   2kπx  
f(x) =
2
+ ∑  a
k =1
k cos 
 b − a  + bk sen  b − a   ,
  
donde
2 b
 2kπx  2 b
 2kπx 
ak =
b−a ∫ a
f(x) cos 
b − a
 dx, bk =
b−a ∫ a
f(x)sen 
b − a
 dx,

para toda k.

Ejemplo 35. Halle el desarrollo en serie de Fourier, en PC 0,1 de la función f(x) = x .

Solución.
En este caso, b − a = 1 , de modo que
1 1
ak = 2
∫ 0
x cos(2kπx)dx, bk = 2
∫ 0
xsen(2kπx)dx .

1
Integrando por partes se tiene ahora: a0 = 1 , ak = 0, k ≠ 0 , bk = − .
kπ
Por lo tanto
1 1 sen(4πx) sen(6πx) 
f(x) = −  sen(2πx) + + + ... .
2 π 2 3 

Ejemplo 36. Encuentre el desarrollo en serie de Fourier de la función

x − 2 2 ≤ x ≤ 3
f(x) =  .
4 − x 3 ≤ x ≤ 4
Solución. Las fórmulas serían
4 3 4
ak =
∫2
f(x) cos(kπx)dx =
∫ 2
(x − 2) cos(kπx)dx +
∫ 3
(4 − x) cos(kπx)dx
4 3 4
bk =
∫ 2
f(x)sen(kπx)dx =
∫ 2
(x − 2)sen(kπx)dx +
∫ 3
(4 − x)sen(kπx)dx

Aunque estas integrales pueden calcularse directamente, los cálculos pueden

simplificarse considerablemente al tomar en cuenta el siguiente argumento:
Designando por F la extensión periódica de f en todo el eje x. Entonces, las funciones
F(x) cos(kπx) y F(x)sen(kπx) son periódicas con período 2, y se tiene
a+ 2 4 a+ 2 4

∫ a
F(x) cos(kπx)dx =
∫ 2
f(x) cos(kπx)dx,
∫ a
F(x)sen(kπx)dx =
∫ 2
f(x)sen(kπx)dx

para cualquier número real a. (En este punto se utiliza el hecho obvio de que si g es
continua por tramos en (−∞, ∞) con período 2p, entonces
a + 2p b + 2p

∫ a
g(x)dx =
∫ b
g(x)dx

Prof. José Luis Quintero 39

para cualquier par de números reales a y b.) Ahora se hace a = −1 y se obtiene
1 1
ak =
∫ −1
F(x) cos(kπx)dx , bk =
∫ −1
F(x)sen(kπx)dx

Pero en el intervalo −1,1 , F coincide con la función par x . Entonces bk = 0

para toda k, y
1
ak = 2
∫ 0
x cos(kπx)dx .

Así
 4
− k impar
a0 = 1, ak =  k2 π2 ,
 0 k par, k ≠ 0

y el desarrollo en serie de Fourier de f es
1 4  cos(3πx) cos(5πx) 
f(x) = − 2  cos(πx) + + + ... .
2 π  32 52 

Ejemplo 37. Encuentre la serie de Fourier de cosenos de

 1 0≤ x ≤1
f(x) =  .
2 − x 1 ≤ x ≤ 2
Solución.
x + 2 −2 ≤ x ≤ −1

f(x) =  1 −1 ≤ x ≤ 1 .
2 − x 1 ≤ x ≤ 2

2
1 2  x2  3 3
∫ ∫
1
a0 = dx + (2 − x)dx = x 0 +  2x −  = 1 + 2 − = .
 2  2 2
0 1  1
1
 nπ  2
 nπ  2
 nπ 
an =
∫ 0
cos  x  dx + 2
 2  ∫ 1
cos  x  dx −
 2  ∫ 1
x cos  x  dx
 2 
1
1
 nπ  2  nπ  2
n = 1,5, 9,13,...
∫ 0
cos 
 2
x  dx =
 n π
sen 
 2
x =
0 n π
.

1
1
 nπ  2  nπ  2
n = 3, 7,11,15,...
∫ 0
cos  x  dx =
 2  nπ
sen  x = −
 2 0 nπ
.

1
1
 nπ  2  nπ 
n = 2, 6,10,14,...
∫ 0
cos 
 2
x  dx =
 n π
sen  x  = 0.
 2 0
1
1
 nπ  2  nπ 
n = 4, 8,12,16,...
∫ 0
cos  x  dx =
 2  nπ
sen  x  = 0.
 2 0
2
2
 nπ  4  nπ  4
n = 1,5, 9,13,... an = 2
∫ 1
cos 
 2
x  dx =
 n π
sen 
 2
x = −
1 n π
.

2
2
 nπ  4  nπ  4
n = 3, 7,11,15,... an = 2
∫ 1
cos  x  dx =
 2  nπ
sen  x =
 2  1 nπ
.

2
2
 nπ  4  nπ 
n = 2, 6,10,14,... an = 2
∫ 1
cos 
 2
x  dx =
 n π
sen  x = 0.
 2 1

Prof. José Luis Quintero 40

 2
2
 nπ  4  nπ 
n = 4, 8,12,16,... an = 2
∫ 1
cos 
 2
x  dx =
 n π
sen  x = 0.
 2 1
2
2
 nπ   2x  nπ  4  nπ   4 2
n = 1,5, 9,13,...
∫ 1
x cos  x  dx = 
 2   n π
sen 
 2
x  + 2 2 cos 
 n π  2
x  = − 2 2 −
 1 n π n π
2
2
 nπ   2x  nπ  4  nπ   4 2
n = 3, 7,11,15,...
∫ 1
x cos  x  dx = 
 2   n π
sen 
 2
x  + 2 2 cos 
 n π  2
x  = − 2 2 +
 1 n π n π
2
2
 nπ   2x  nπ  4  nπ   8
n = 2, 6,10,14,...
∫ 1
x cos  x  dx = 
 2   n π
sen 
 2
x  + 2 2 cos 
 n π  2
x  = 2 2
 1 n π
2
2
 nπ   2x  nπ  4  nπ  
n = 4, 8,12,16,...
∫ 1
x cos  x  dx = 
 2   n π
sen 
 2
x  + 2 2 cos 
 n π
x  = 0
 2  1
2 4 4 2 4
n = 1,5, 9,13,... an = − + 2 2 + = 2 2.
nπ nπ n π nπ n π
2 4 4 2 4
n = 3, 7,11,15,... an = − + + 2 2 − = 2 2.
nπ nπ n π nπ n π
8
n = 2, 6,10,14,... an = − 2 2 . n = 4, 8,12,16,... an = 0.
n π
∞
  (2n − 1)π  (−1)n − 1  (2n)π  
∑
3 4 1
f(x) = +  cos  x + cos  x
4 π2 n =1 
(2n − 1)2
 2  (2n)2
 2  
(ver figuras 16, 17 y 18)

3 4  π  1 
s2 (x) = +  cos  x  − cos(πx) 
4 π2   2  2 

Figura 16. Primer desarrollo de Fourier para el ejemplo 37

Prof. José Luis Quintero 41

 3 4  π  1 1  3π  
s4 (x) = + 2  cos  x  − cos(πx) + cos  x
4 π  2  2 9  2 

Figura 17. Segundo desarrollo de Fourier para el ejemplo 37

3 4  π  1 1  3π  1  5π  1 
s6 (x) = + 2  cos  x  − cos(πx) + cos  x + cos  x − cos(3πx) 
4 π  2  2 9  2  25  2  18 

Figura 18. Tercer desarrollo de Fourier para el ejemplo 37

Ejemplo 38. Sin recurrir a la integración, encuentre el desarrollo en serie de Fourier

de la función f(x) = cos3 (x).sen(2x) en el intervalo (−π, π) .
Solución.
Usando las identidades trigonométricas, se puede demostrar que

Prof. José Luis Quintero 42

 senα cos β =
sen(α + β) + sen(α − β) . Por tanto:
1
2
1 + cos(2x)
cos3 (x)sen(2x) = cos(x). .sen(2x) = 12 cos(x) sen(2x) + cos(2x)sen(2x)
2
= 12 cos(x)sen(2x) + 14 cos(x)sen(4x)
= 1
4
sen(3x) + sen(x) + 1
8
sen(5x) + sen(3x)
= 1
4
sen(x) + 3
8
sen(3x) + 1
8
sen(5x)

Ejemplo 39. Sea

 1
x 0 ≤ x ≤ 10
f(x) =  1 10 .
10 (20 − x) 10 ≤ x ≤ 20
Encuentre el desarrollo en serie de Fourier de:
a. cosenos de f(x) y úsela para demostrar que
∞
π2
∑ (2k + 1)
1
2
= .
k =0
8
Solución.
 10 20   2 10 20 
a0 = 2
200


∫ 0
xdx +
∫ 10
(20 − x)dx  =

1
100
 x2
 0
+ (20x − x2 )
2
=
10 
1 (50
100
+ 200 − 200 + 50) = 1

 10 20 
ak = 2
200


∫ 0
πx )dx +
x cos(k20
∫ 10
πx )dx . k = 1,...
(20 − x) cos(k20 

 20 20 20 
= 100
1


∫ 0
πx )dx + 20
x cos(k20
∫ 10
πx )dx −
cos(k20
∫ 10
πx )dx
x cos(k20 

 10 20 20 
πx πx πx πx πx
= ( 20x
1
100 kπ
sen(k20 ) + 400
2 2
cos(k20 )) + 400
k π
sen(k20 ) − (20xk π
sen(k20 ) + 400
2 2
cos(k20 )) 
 k π 0 10 k π 10 

1  200 sen( kπ ) + 400 cos( kπ ) − 400 − 400 sen( kπ ) − 400 cos(kπ) + 200 sen( kπ ) + 400 cos(kπ )
= 100  kπ 2 k2 π2 2 k2 π2 kπ 2 k2 π2 kπ 2 k2 π2 2 

= 24 2 2 cos(k2π ) − 1 − (−1)k  . Si k = 2n,

k π  
= 4 2 cos( 2nπ ) − 1 − (−1)2n  = 1
2 cos(nπ) − 2 =
2 (−1)n − 1
4n2 π2  2  n2 π2 n2 π2  
Por lo tanto
∞
f(x) = 1
2
+ 2
π2 ∑ 
k =1
( −1)k −1 
k2
kπx
 cos( 10 ) .

En particular, para x = 0 se obtiene
∞ ∞ ∞ ∞
0= 1
2
+ 2
π2 ∑
k =1
 (−1)k −1 
 k2  ⇒ 0 =
1
2
− 4
π2 ∑k =1
(2k +1)2
1
⇒ 4
π2 ∑
k =0
1
(2k + 1)2
= 1
2
⇒ ∑
k =0
1
(2k +1)2
= π2
8
.

b. senos de f(x) y use la identidad de Parseval para hallar el valor de

∞

∑ (2k + 1)
1
4
.
k =0

Solución.

Prof. José Luis Quintero 43

  10 20 
bk = 2 
200 
 0
∫xsen(kπx )dx +
20
10 ∫
(20 − x)sen(kπx )dx  . k = 1,...
20 

 20 20 20 
= 1 
100 
 0
∫ xsen(kπx )dx + 20
20
10 ∫
sen(kπx )dx −
20
10
20 ∫
xsen(kπx )dx 


 10 20 20 
= 1 (− 20x cos(kπx ) + 400 sen(kπx )) − 400 cos(kπx ) − (− 20x cos(kπx ) + 400 sen(kπx )) 
100  kπ 20 2 2 20 k π 20 k π 20 2 2 20
 k π 0 10 k π 10 

 
= 1  − 200 cos(kπ ) + 400 sen(kπ ) − 400 cos(kπ) + 400 cos(kπ ) + 400 cos(kπ) − 200 cos(kπ ) + 400 sen(kπ )
100  kπ 2 k2 π2 2 kπ kπ 2 kπ kπ 2 k2 π2 2 

 
= 1  400 sen(kπ ) + 400 sen(kπ ) = 8 sen(kπ ).
100  k2 π2 2 k2 π2 2  k2 π2 2

Por lo tanto
∞
 sen(k2π ) 
f(x) = 8
π2 ∑ 
k =1
k2

πx
 sen(k20 ).

Al calcular
20

∫
2
f(x) dx
0

se obtiene:
20 10 20  3 10 (20 − x)3
20 

∫ ∫ ∫
2
f(x) dx =
1
100
x2dx + 1
100
(20 − x)2 dx = 1
100
 x3 − 3
= 20
3
.
0 0 10
 0 
10 

Aplicando la igualdad de Parseval:

∞ ∞ ∞
sen2 (kπ ) sen2 ( kπ )
64
π4 ∑
k =1
k4
2
= 20
3
⇒ ∑
k =1
k4
2
= π4 20
64.3
⇒ ∑
k =0
1
(2k +1)4
= 5π4
48
.

4.15. ALGUNAS DEFINICIONES

Definición 14. Una ecuación diferencial parcial (EDP) es una relación funcional
entre una función incógnita u de varias variables independientes, dichas variables
independientes y las derivadas parciales de u.

Ejemplo 40.
∂z ∂2 t ∂2 t ∂2 t
= −2 , + 2 + 2 =0
∂y ∂x2
∂y ∂z
son ecuaciones diferenciales parciales o en derivadas parciales.

Ejemplo 41. Se tienen las siguientes ecuaciones diferenciales parciales:

a. utt − a2uxx = 0 (Ecuación de onda unidimensional que caracteriza la propagación de
ondas en algunos medios y las vibraciones mecánicas de una cuerda vibrante)
b. uxx + uyy + uzz = 0 (Ecuación de Laplace que se presenta en el estudio de
potenciales magnético, eléctrico, gravitatorio y en el flujo de calor)

Definición 15. El orden de una EDP es el mayor orden de todas las derivadas
parciales que figuran en la ecuación.

Prof. José Luis Quintero 44

Ejemplo 42. La ecuación
∂ 4t ∂2 t
+ =0
∂x 4 ∂x2
es una ecuación diferencial parcial de orden cuatro.

Definición 16. El grado de una EDP lo da la potencia de la derivada de mayor orden

en dicha ecuación.

Ejemplo 43.
∂2 z ∂z
− 3x =6
∂x 2 ∂x
es una ecuación diferencial parcial de orden 2 y grado 1.

Observación 3. Si se utiliza la notación

ux = ∂u / ∂x , uy = ∂u / ∂y , uxx = ∂2u / ∂x2 , uxy = ∂2u / ∂x∂y ,
se puede escribir una ecuación diferencial parcial general de segundo orden para una
función u de dos variables independientes x e y definida en un dominio D, en la forma

F(x, y,u,ux ,uy ,uxx ,uxy ,uyy ) = 0 , (17)

relación que debe verificarse para todo punto (x,y) del interior R del dominio D de la
función u.

Definición 17. Una solución de (17) en R es cualquier función u definida y dos veces
continuamente diferenciable en R que verifica (17) idénticamente para todo (x,y) en
R2 .

4.16. DIFERENCIAS ENTRE LAS ECUACIONES DIFERENCIALES

ORDINARIAS (EDO) Y LAS ECUACIONES DIFERENCIALES
PARCIALES (EDP)

a. En cuanto al tipo de problemas donde se presentan, se puede decir que las

EDO describen problemas unidimensionales, mientras que las EDP describen
problemas de dos o más dimensiones.
b. En cuanto al espacio de soluciones, se puede decir que en el caso de las EDO el
espacio de soluciones es de dimensión finita, y además la dimensión de dicho
espacio coincide con el orden de la ecuación diferencial, mientras que en el caso de
las EDP no necesariamente el espacio de soluciones es de dimensión finita y
además no depende del orden de la ecuación. Es bueno hacer notar que en la
mayoría de los casos los espacios de soluciones para EDP son de dimensión
infinita.
c. En cuanto a la cantidad de condiciones suplementarias que hay que
agregar a la ecuación para determinar en forma única una solución, en el
caso de las EDO, se necesitarán tantas como el número de constantes arbitrarias

Prof. José Luis Quintero 45

 que aparezcan en la solución, o sea, tantas como el orden de la ecuación. En el
caso de las EDP, aún en el caso comparativamente más simple de las ecuaciones
lineales la situación es mucho más complicada. Las soluciones de una EDP
dependen no ya de constantes arbitrarias, sino de funciones arbitrarias.

4.17. SOLUCIONES DE ALGUNAS ECUACIONES DIFERENCIALES

PARCIALES SENCILLAS

Con el objeto de obtener algunas ideas relacionadas con la naturaleza de las

soluciones de ecuaciones diferenciales parciales, considere los siguientes ejemplos:

Ejemplo 44. La ecuación de segundo orden

∂2u
=0
∂y2
puede ser resuelta integrando dos veces respecto de
∂u
y: = f(x) , u = yf(x) + g(x)
∂y
donde f y g son funciones arbitrarias.

Ejemplo 45. Obtenga soluciones de la ecuación diferencial parcial

∂2U
= 6x + 12y2 . (18)
∂x∂y
Solución. Aquí la variable dependiente U depende de dos variables independientes x y
y. Para hallar soluciones, se busca determinar U en términos de x y y, esto es, U(x,y)
Si se escribe (18) como
∂  ∂U  2
  = 6x + 12y ,
∂x  ∂y 
se puede integrar con respecto a x, manteniendo y constante, para hallar
∂U
= 3x2 + 12xy2 + F(y) (19)
∂y
donde se ha añadido la “constante” de integración, la cual depende de y y así
realmente es una función de y denotada por F(y). Ahora se integra (19) con respecto a
y manteniendo x constante para hallar
U = 3x2 y + 4xy3 +
∫ F(y)dy + G(x) (20)

esta vez añadiendo una función arbitraria de x dada por G(x). Puesto que la integral de
una función arbitraria de y es otra función arbitraria de y, se puede escribir (20) como
U = 3x2 y + 4xy3 + H(y) + G(x) . (21)

Esto se puede chequear al sustituirla en (18) y obtener una identidad. Puesto

que (18) es una ecuación diferencial parcial de segundo orden, y (21) tiene dos
funciones arbitrarias, por analogía con las ecuaciones diferenciales ordinarias se llama
a (21) la solución general de (18). Usando la misma analogía, se llamaría por
supuesto a cualquier solución obtenida de la solución general (21) por selecciones

Prof. José Luis Quintero 46

particulares de las funciones arbitrarias, tal como por ejemplo H(y) = y3 , G(x) = sen2x ,
una solución particular. Se construye así la siguiente definición:

Definición 18. Dada una ecuación diferencial parcial de orden n, una solución que
contenga n funciones arbitrarias se llama la solución general y cualquier solución
obtenida de esta solución general por selecciones particulares de las funciones
arbitrarias se llama una solución particular.

Observación 4. Como en el caso de ecuaciones diferenciales ordinarias, puede

suceder que existan soluciones singulares las cuales no se puedan obtener de la
solución general por cualquier selección de las funciones arbitrarias.

Como en el caso de ecuaciones diferenciales ordinarias, con frecuencia se

necesitan determinar soluciones de ecuaciones diferenciales parciales que satisfagan
condiciones dadas. Por ejemplo, suponga que se desea resolver la ecuación diferencial
(18) sujeta a las dos condiciones
U(1, y) = y2 − 2y , U(x, 2) = 5x − 5 . (22)

Entonces de la solución general (21) y la primera condición (22) se tiene

U(1, y) = 3(1)2 y + 4(1)y3 + H(y) + G(1) = y2 − 2y o H(y) = y2 − 5y − 4y3 − G(1)
de modo que
U = 3x2 y + 4xy3 + y2 − 5y − 4y3 − G(1) + G(x) . (23)

Si ahora se usa la segunda condición en (22), se tiene

U(x, 2) = 3x2 (2) + 4x(2)3 + (2)2 − 5(2) − 4(2)3 − G(1) + G(x) = 5x − 5
de donde G(x) = 33 − 27x − 6x2 + G(1) . Usando esto en (23), se obtiene la solución
U = 3x2 y + 4xy3 + y2 − 5y − 4y3 − 27x − 6x2 + 33 .

Se podría usar la misma terminología de problemas de valor inicial y de

frontera para las ecuaciones diferenciales parciales como se hizo para las ecuaciones
diferenciales ordinarias. Sin embargo, debido a que generalmente hay una
combinación de condiciones de frontera e iniciales, con frecuencia se refiere a tales
problemas como problemas de valor de frontera.

Ejemplo 46. Hallar una solución al problema de valor de frontera

∂2U ∂U
= +2 , U(0, y) = 0 , Ux (x, 0) = x2
∂x∂y ∂x
Solución.
Escribiendo la ecuación como
∂  ∂U 
 − U = 2
∂x  ∂y 
e integrando con respecto a x se tiene

Prof. José Luis Quintero 47

 ∂U
− U = 2x + F(y) ,
∂y
la cual es una ecuación lineal con factor integrante e− y . Por tanto,
∂
(e− yU) = 2xe− y + e− yF(y)
∂y
o

∫
U(x, y) = −2x + ey e− yF(y)dy + ey G(x) ,

donde G(x) es arbitraria. Escribiendo

∫
H(y) = ey e− yF(y)dy ,

se tiene
U(x, y) = −2x + H(y) + ey G(x) . (24)

De U(0, y) = 0 se encuentra H(y) = −G(0)ey , así (24) llega a ser

U(x, y) = −2x − G(0)ey + ey G(x) .

Diferenciando con respecto a x y colocando y = 0 , se encuentra

x3
Ux (x, 0) = −2 + G '(x) = x2 o G(x) = + 2x + c ,
3
de donde,
 x3 
U(x, y) = −2x − G(0)ey + ey  + 2x + c  .
 3 
 
Puesto que c = G(0), entonces se tiene que
x3ey
U(x, y) = + 2xey − 2x .
3

4.18. SIGNIFICADO GEOMÉTRICO DE LAS SOLUCIONES

GENERAL Y PARTICULAR

En el ejemplo 45, se obtuvo la solución general

U = 3x2 y + 4xy3 + H(y) + G(x) . (25)

Suponga ahora que se escogen funciones particulares para H(y) y G(x), y se

reemplaza U por z. Entonces (25) toma la forma z = f(x, y) , la cual se interpreta como
una superficie S en un sistema de coordenadas rectangular o xyz. La superficie está
formada por los puntos con coordenadas (x,y,z) que satisfacen z = f(x, y) .

Para funciones arbitrarias H(y) y G(x), se obtiene una familia de superficies

cada miembro de la cual corresponde a una selección particular de H(y) y G(x), esto
es, una solución particular. La ecuación diferencial que tenga esto como una solución
se llama entonces la ecuación diferencial de la familia de superficies.

Prof. José Luis Quintero 48

4.19. ECUACIONES DIFERENCIALES PARCIALES QUE SURGEN
DE LA ELIMINACIÓN DE FUNCIONES ARBITRARIAS

Puesto que las soluciones de ecuaciones diferenciales parciales involucran

funciones arbitrarias, parece lógico que se debiera obtener ecuaciones diferenciales
parciales por el proceso inverso de eliminar tales funciones. Esta idea es útil porque
ayuda a enriquecer el conocimiento de cómo se pueden resolver ecuaciones
diferenciales parciales. Se van a considerar algunos ejemplos.

Ejemplo 47. Encuentre una ecuación diferencial parcial de primer orden que tenga
como su solución general
U = y2F(x) − 3x + 4y , (26)
donde F(x) es una función arbitraria de x.
Solución.
Si se diferencia a (26) con respecto a y, se obtiene
∂U
= 2yF(x) + 4 . (27)
∂y
Entonces eliminando F(x) entre (26) y (27), se encuentra la ecuación deseada
∂U
y − 2U = 6x − 4y .
∂y
∂U
Chequeo. y − 2U = y[2yF(x) + 4] − 2[y2F(x) − 3x + 4y] = 6x − 4y .
∂y

Ejemplo 48. Encuentre una ecuación diferencial parcial de primer orden que tenga
como solución general
z = F(3x − 4y) , (28)
donde F es una función arbitraria.
Solución.
Sea u = 3x − 4y . Entonces (28) llega a ser
z = f(u) (29)
Diferenciando (29) con respecto a x, se tiene
∂z ∂z ∂u
= . = F '(u).3 = 3F '(u) . (30)
∂x ∂u ∂x
Diferenciando (29) con respecto a y, se tiene
∂z ∂z ∂u
= . = F '(u).(−4) = −4F '(u) . (31)
∂y ∂u ∂y
Eliminando F '(u) entre (30) y (31) produce la ecuación
∂z ∂z
4 +3 = 0.
∂x ∂y

Ejemplo 49. Encuentre una ecuación diferencial parcial de segundo orden que tenga
como su solución general
U = xF(y) + yG(x) , (32)
donde F y G son funciones arbitrarias.
Solución.

Prof. José Luis Quintero 49

Se puede eliminar F(y) en (32) dividiendo ambos lados de (32) por x y diferenciando el
resultado con respecto a x. Entonces se encuentra
∂ U  ∂  y  ∂U
  = F(y) + G(x) , esto es, x − U = xyG '(x) − yG(x) ,
∂x  x  ∂x  x  ∂x
la cual se puede escribir
∂U
x − U = y[xG '(x) − G(x)] . (33)
∂x
Si ahora se dividen ambos lados de (33) por y y se diferencia con respecto a y, se
encuentra
∂  1  ∂U  ∂2U ∂U ∂U
 x − U   = 0 ó xy −x −y +U = 0,
∂y y
  ∂x  ∂y ∂ x ∂ x ∂y
la cual construye la ecuación de segundo orden deseada. Note que la segunda
ecuación también se puede escribir como
∂2U ∂U ∂U ∂2U ∂2U
xy −x −y + U = 0 puesto que = .
∂x∂y ∂x ∂y ∂x∂y ∂y∂x

Se deberían hacer algunas observaciones acerca de los resultados anteriores.

Observación 5. Al diferenciar las funciones arbitrarias se asume por supuesto que

ellas son diferenciables. De otra manera, no se tiene derecho a diferenciar.

Observación 6. La ecuación diferencial obtenida en cada ejemplo representa la

ecuación diferencial de la familia representada por la solución general.

Observación 7. Si una solución tiene un número dado de n funciones arbitrarias, con

frecuencia es fácil escribir una ecuación diferencial de orden mayor que n teniendo
esta solución. Por ejemplo, es fácil ver que (32) es una solución de
∂ 4U
=0
∂x2 ∂y2
pero ésta es de orden 4 y no 2. Cuando se busca la ecuación diferencial, se busca
aquella del menor orden.

4.20. EL MÉTODO DE SEPARACIÓN DE VARIABLES

A veces, para una ecuación diferencial en derivadas parciales lineal

homogénea, es posible obtener soluciones particulares en forma de producto
u(x, y) = X(x)Y(y) . El uso de este producto, llamado método de separación de
variables, permite reducir la ecuación diferencial en derivadas parciales a varias
ecuaciones diferenciales ordinarias. Con este propósito, se hace notar que
∂u ∂u ∂2u ∂2u
= X'Y , = XY ' , = X '' Y , = XY '' ,
∂x ∂y ∂x2 ∂y2
en donde las primas indican diferenciales ordinarias.

Prof. José Luis Quintero 50

Ejemplo 50. Halle soluciones en forma de producto de la ecuación
∂2u ∂u
=4 . (34)
∂x 2 ∂y
Solución.
Si u = X(x)Y(y), entonces (34) se transforma en X '' Y = 4XY ' . Después de dividir
ambos miembros entre 4XY, se logra separar las variables:
X '' Y '
= .
4X Y
Puesto que el lado izquierdo de esta ecuación es independiente de y y es idéntico al
lado derecho, el cual es independiente de x, se concluye que ambos miembros deben
ser iguales a una constante. En la práctica es conveniente escribir esta constante real
como λ2 , o bien como −λ2 . Se distinguen los casos siguientes:
CASO I: Si ctte = λ2 , las igualdades
X '' Y '
= = λ2
4X Y
conducen a las ecuaciones X ''− 4λ2 X = 0 y Y '− λ2 Y = 0 . Estas últimas tienen las
2y
soluciones X = c1e−2λx + c2e2λx y Y = c3eλ , respectivamente. Así, una solución
particular de (34) es
2y 2y − 2 λx 2y + 2 λx
u = XY = (c1e−2λx + c2e2λx )c3eλ = A1eλ + B1eλ ,
en donde A1 = c1c3 y B1 = c2c3 .
X '' Y '
CASO II. Si ctte = −λ2 , las igualdades = = −λ2 conducen a X ''+ 4λ2 X = 0 y
4X Y
Y '+ λ2 Y = 0 . Puesto que las soluciones de estas ecuaciones son
2y
X = c4 cos(2λx) + c5sen(2λx) y Y = c6 e−λ ,
2y 2
respectivamente, otra solución de (34) es u = A2e−λ cos(2λx) + B2e−λ y sen(2λx) , en
donde A2 = c4c6 y B2 = c5c6 .
CASO III. Si ctte = 0 , se tiene que X '' = 0 y Y ' = 0 . En este caso X = c7 x + c8 y
Y = c9 , de modo que u = A3x + B3 en donde A3 = c7c9 y B2 = c8c9 .

Se deja como ejercicio verificar que las tres soluciones obtenidas satisfacen la
ecuación dada.

Ejemplo 51. Halle una solución en forma de producto de

∂2u ∂u
k 2 = , k>0
∂x ∂t
que satisfaga las condiciones u(0, t) = 0, u(L, t) = 0 .
Solución.
Si u = XT , se puede escribir la ecuación dada como
X '' T '
= = −λ2
X kT
lo que conduce a X ''+ λ2 X = 0 y T '+ kλ2 T = 0 , de modo que se tiene
2
X = c1 cos(λx) + c2sen(λx) y T = c3e−kλ t ,

Prof. José Luis Quintero 51

respectivamente. Ahora bien, puesto que u(0, t) = X(0)T(t) = 0 , u(L, t) = X(L)T(t) = 0 se
debe tener X(0) = 0 y X(L) = 0 . Estas son condiciones de frontera para la ecuación
diferencial ordinaria. Aplicando la primera de tales condiciones, resulta de inmediato
c1 = 0 . Por lo tanto, X = c2sen(λx) . La segunda condición de frontera implica ahora que
X(L) = c2sen(λL) = 0 . Si c2 = 0 , entonces X = 0 , de modo que u = 0 . Para obtener una
solución no trivial u, se debe tener c2 ≠ 0 y entonces la última ecuación se satisface
cuando sen(λL) = 0 . Entonces λL = nπ o
nπ
λ= ,
L
siendo n = 1,2, 3,...
Por tanto,
2  nπ  2 π2 / L2 )t
u = (c2sen(λx)(c3e−kλ t ) = Ane−k(n
sen  x
 L 
satisface la ecuación dada y ambas condiciones adicionales. El coeficiente c2c3 se
reescribe como An para recalcar que se obtiene una solución diferente para cada n. Se
deja como ejercicio verificar que si se usa λ2 ≥ 0 , no se llega a una solución que
satisfaga u(0, t) = 0 y u(L, t) = 0 .

4.21. PRINCIPIO DE SUPERPOSICIÓN

TEOREMA 10. Si u1 , u2 , ..., un son soluciones de una ecuación diferencial paracial

lineal homogénea, entonces la combinación lineal u = c1u1 + c2u2 + ... + ckuk , donde los
ci , i = 1, 2,...,k son constantes, también es una solución.

En la próxima sección se supondrá formalmente que cada vez que se tenga un

conjunto infinito u1 , u2 , ... de soluciones de una ecuación lineal homogénea, aún se
puede obtener otra solución u formando la serie infinita
∞
u= ∑u
k =1
k .

Ejemplo 52. En virtud del principio de superposición, la función definida mediante la

serie
∞
 nπ 
u= ∑A e
n =1
n
−k(n2 π2 / L2 )t
sen  x
 L 
debe satisfacer también, aunque sea formalmente, a la ecuación del ejemplo anterior,
junto con las condiciones dadas.

4.22. PROBLEMAS DE CONDICIÓN DE FRONTERA

Las siguientes ecuaciones diferenciales parciales lineales

∂2u ∂u ∂2u ∂2u ∂2u ∂2u
k 2 = , k > 0, a2 2 = 2 , + = 0,
∂x ∂t ∂x ∂t ∂x2 ∂v2

Prof. José Luis Quintero 52

desempeñan un papel importante en muchas áreas de física e ingeniería. Las primeras
dos ecuaciones son conocidas como ecuación del calor en una dimensión y
ecuación de onda en una dimensión respectivamente. “Una dimensión” se refiere al
hecho de que x denota una dimensión espacial, en tanto que t generalmente
representa tiempo. La última ecuación se llama ecuación de Laplace. Se usará el
método de separación de variables para resolver varios problemas aplicados, cada uno
de los cuales es descrito por una de las ecuaciones anteriores además de ciertas
condiciones adicionales. Estas condiciones adicionales consisten en:
a. Condiciones de frontera: u o ∂u ∂x especificada para x = cons tan te ; u o ∂u ∂y
especificada para y = cons tan te ,
b. Condiciones iniciales: u en t = 0 para la primera ecuación, o bien, u y ∂u ∂t en
t = 0 para la segunda ecuación.

La descripción matemática colectiva de un problema de esta naturaleza es

conocida como problema de condición en la frontera.

En una frontera se pueden especificar los valores de una de las siguientes

funciones:
(i) u , (ii) ux , o bien (iii) ux + hu , h cons tan te .

Una condición a la frontera del primer tipo (i), se llama condición de

Dirichlet; del segundo tipo (ii), condición de Neuman, y del tercer tipo (iii),
condición de Robin.

4.23. ECUACIONES DIFERENCIALES PARCIALES LINEALES DE

SEGUNDO ORDEN

Ya se han visto cómo se resuelven algunos problemas de frontera que

involucran ecuaciones diferenciales parciales. Se enfocará ahora el estudio sobre cierto
tipo de ecuaciones diferenciales parciales importantes que surgen en varios problemas
aplicados. Hay tres tipos importantes de problemas que involucran:
a. Conducción o difusión de calor
b. Vibraciones u oscilaciones
c. Temperatura estacionaria, potencial eléctrico o gravitacional

Definición 19. Si u representa la variable dependiente, x e y las variables

independientes, entonces la forma general de una ecuación diferencial parcial
lineal de segundo orden con dos variables independientes x e y, es
∂2u ∂2u ∂2u ∂u ∂u
A 2 +B +C 2 +D +E + Fu = G , (35)
∂x ∂x∂y ∂y ∂x ∂y
en donde A, B, C, D, E, F y G son funciones de x e y.

Observación 8. Cuando G(x, y) = 0 , la ecuación se llama homogénea; en cualquier

otro caso es no homogénea.

Prof. José Luis Quintero 53

 ∂2u ∂2u ∂2u ∂u
Ejemplo 53. + = 0 es homogénea. − = xy no es homogénea.
∂x2
∂y 2
∂x 2 ∂y

Definición 20. Una solución de una ecuación (35) en derivadas parciales con dos
variables independientes x e y es una función u(x,y) que posee todas las derivadas
parciales que indica la ecuación y que la satisface en alguna región del plano xy.

Una ecuación en derivadas parciales, lineal de segundo orden con dos variables
independientes y con coeficientes constantes, puede pertenecer a uno de tres tipos
generales. Esta clasificación sólo depende de los coeficientes de las derivadas de
segundo orden. Naturalmente, se supone que al menos uno de los coeficientes A, B y C
no es cero.

Definición 21. La ecuación diferencial parcial lineal y de segundo orden

∂2u ∂2u ∂2u ∂u ∂u
A 2 +B +C 2 +D +E + Fu = 0 ,
∂x ∂x∂y ∂y ∂x ∂y
en donde A, B, C, D, E, F y G son constantes reales es

a. parabólica, si B2 − 4AC = 0 .
∂2u ∂u
Ejemplo 54. 3 − = 0 (A = 3,B = 0, C = 0) . Esta ecuación modela la distribución
∂x 2 ∂y
de temperatura en una varilla (ecuación de calor unidimensional).
b. hiperbólica, si B2 − 4AC > 0 .
∂2u ∂2u
Ejemplo 55. − = 0 (A = 1,B = 0, C = −1) . Esta ecuación modela las vibraciones
∂x2 ∂y2
en una cuerda (ecuación de onda unidimensional).
c. elíptica, si B2 − 4AC < 0 .
∂2u ∂2u
Ejemplo 56. + = 0 (A = 1,B = 0, C = 1) . Esta ecuación modela el potencial o
∂x2 ∂y2
distribución de temperatura estacionaria no dependiente del tiempo (ecuación de
Laplace).

Estas ecuaciones tienen métodos para hallar su solución, entre los que
destacan:
a. Método de separación de variables
b. Método de la transformada de Laplace
c. Método de la integral de Fourier

4.24. ECUACIÓN DE CALOR UNIDIMENSIONAL

A continuación se empezará formulando las ecuaciones del flujo de calor que

describen la transferencia de energía térmica (o calorífica). La energía térmica está
causada por la agitación de las moléculas. Existen dos procesos básicos que

Prof. José Luis Quintero 54

intervienen en la transferencia de la energía térmica: conducción y convección. La
conducción resulta de las colisiones de moléculas colindantes, por las que la energía
cinética de vibración de una molécula se transfiere a su vecina más cercana. La
energía térmica se transporta de esta manera por conducción incluso si las moléculas
mismas no se mueven de forma apreciable de su posición. Por otro lado, una molécula
vibrante que se traslada de una región a otra, lleva su energía térmica con ella. Este
segundo fenómeno se conoce con el nombre de convección. Con objeto de comenzar el
estudio con problemas relativamente simples, se estudiará el flujo de calor sólo en los
casos en que la conducción de energía térmica predomina sobre la convección. Éste es
el caso del flujo de calor en sólidos, o incluso en fluidos (líquidos y gases) en los que la
velocidad es suficientemente pequeña.

4.25. DEDUCCIÓN DE LA CONDUCCIÓN DEL CALOR EN UNA

VARILLA UNIDIMENSIONAL

Suponga que una varilla circular delgada de longitud L tiene un área transversal
A y coincide con el eje x en el intervalo 0,L  (ver figura 19). Suponga lo siguiente:

x
Q
x=0 ∆x
x x+∆ x=L

Figura 19. Flujo de calor en una varilla delgada

a. El flujo de calor dentro de la varilla toma lugar solo en la dirección x.

b. La superficie lateral, o curva, de la varilla está aislada; es decir, el calor no escapa
de su superficie.
c. No se genera calor dentro de la varilla.
d. La varilla es homogénea; es decir, su masa por unidad de volumen ρ es una
constante.
e. El calor específico γ y la conductividad térmica K del material de la varilla, son
constantes.

Prof. José Luis Quintero 55

 A fin de obtener la ecuación diferencial parcial que se cumple para la
temperatura u(x,t), son necesarias dos leyes empíricas de conducción del calor:
a. La cantidad de calor Q en un elemento de masa m es
Q = γmu , (36)
donde u es la temperatura del elemento.
b. La rapidez de flujo de calor Qt por la sección transversal indicada en la figura, es
proporcional al área A de la sección transversal y la derivada parcial con respecto a
x de la temperatura:
Qt = −KAux . (37)

Puesto que el calor fluye en la dirección de temperatura decreciente, el signo

menos en (37) se usa para asegurar que Qt sea positiva para ux < 0 (flujo de calor a
la derecha) y negativa para ux > 0 (flujo de calor a la izquierda). Si la parte circular de
la varilla mostrada en la figura entre x y x + ∆x es muy delgada, entonces u(x,t) se
toma como la temperatura aproximada en cada punto del intervalo. Ahora la masa de
la porción es m = ρA∆x , y de (36) se deduce que la cantidad de calor en ésta es
Q = γρA∆xu . (38)
Además, cuando el calor fluye en la dirección positiva de x, se ve de (37) que el
calor se acumula en el corte con razón neta
−KAux (x, t) − −KAux (x + ∆x, t) = KA ux (x + ∆x, t) − ux (x, t) . (39)

Al diferenciar (38) con respecto a t, se ve que esta razón neta se expresa

también mediante
Qt = γρ∆xut . (40)

Al igualar (39) y (40), se obtiene

K ux (x + ∆x, t) − ux (x, t)
= ut . (41)
γρ ∆x

Por último, al tomar el límite de (41) cuando ∆x → 0 , se obtiene la ecuación

(K / γρ)uxx = ut (La definición de la segunda derivada parcial viene dada por la
ux (x + ∆x, t) − ux (x, t)
expresión uxx = lím . Normalmente se establece que k = K / γρ , y
∆x → 0 ∆x
esta constante se llama difusividad térmica.

4.26. PROBLEMAS ILUSTRATIVOS DE LA ECUACIÓN DE CALOR

UNIDIMENSIONAL

Ejemplo 57. Considere una barra delgada o varilla de largo L con una distribución
longitudinal de temperatura f(x) y cuyos extremos se mantienen a una temperatura
constante de cero grados en todo instante. Si la varilla satisface las suposiciones
enunciadas en la sección anterior, entonces la temperatura u(x,t) en la varilla se
determina a partir del problema de valores en la frontera

Prof. José Luis Quintero 56

 ∂2u ∂u
k = , k > 0, 0 < x < L, t > 0 , (42)
∂x 2 ∂t
u(0, t) = 0, u(L, t) = 0, t > 0 , (43)
u(x, 0) = f(x), 0 < x < L . (44)

La constante k es proporcional a la conductividad térmica y se llama

difusividad térmica. Este problema representa el caso de una ecuación de calor
homogénea con condición inicial no nula y condiciones de frontera nulas.
Solución.
Usando el producto u = XT y −λ2 como constante de separación, se llega a
X '' T '
= = −λ2
X kT
y
X ''+ λ2 X = 0, X(0) = 0, X(L) = 0, (45)
2
T '+ kλ T = 0 (46)
Ya se ha obtenido resolviendo (45) la solución de la ecuación (42) sujeta a las
condiciones de frontera (43). Se vió que (45) tenía una solución no trivial sólo si el
parámetro λ adoptaba los valores
nπ
λ= , n = 1,2, 3,... (47)
L
Las correspondientes soluciones de (45) eran entonces
 nπ 
X = c1sen  x, n = 1, 2, 3,... (48)
 L 

Los valores (47) para los cuales (45) tiene una solución no trivial son conocidos
como valores característicos o como valores propios. Las soluciones (48) se
denominan funciones características o funciones propias. Se hace notar que para
un valor de λ distinto de los dados en (47), la única solución de (45) es la función
X = 0 . A su vez, esto implicaría que una función que satisface (42) y (43) en u = 0 no
es una solución del problema original de condición en la frontera cuando f(x) ≠ 0 .
Naturalmente se supone que se cumple esta última condición. Puesto que la solución
2t 2 π2 / L2 )t
de (46) es T = c3e−kλ = c3e−k(n , los productos
2 π2 / L2 )t  nπ 
un = XT = Ane−k(n sen  x (49)
 L 
satisfacen la ecuación diferencial (42) y las condiciones de frontera (43) para cada
valor del entero positivo n. Por conveniencia se ha reemplazado la constante c1c3 por
An . Para que las funciones dadas en (49) satisfagan la condición inicial (44), se
tendría que elegir coeficientes constantes An de modo que
 nπ 
un (x, 0) = f(x) = Ansen  x . (50)
 L 

En general, la condición (50) no se satisface para cualquier f. Por lo tanto, hay

que admitir que (47) no es una solución del problema dado. Sin embargo, por el
principio de superposición, la función

Prof. José Luis Quintero 57

 ∞ ∞
 nπ 
u(x, t) = ∑n =1
un = ∑A e
n =1
n
−k(n2 π2 / L2 )t
sen  x
 L 
(51)

también satisface (42) y (43). Sustituyendo t = 0 en (51) se obtiene

∞
nπ 
u(x, 0) = f(x) = ∑ A sen  L
n =1
n x .


Se ve que esta última expresión es el desarrollo en serie de senos en medio

intervalo. Por consiguiente si se hace la identificación An = bn , n = 1,2, 3,... , se tiene
que
2 L
 nπ 
An =
L ∫ 0
f(x)sen  x  dx .
 L 
Se concluye que la solución del problema de condición de frontera descrito en (42),
(43) y (44) está dada por la serie infinita
∞
 L
 nπ    nπ 
∑∫
2 2 2 2
u(x, t) =  f(x)sen  x  dx  e−k(n π / L )t sen  x .
L n =1  0  L    L 

Ejemplo 58. Ecuación de calor homogénea con condición inicial no nula y condiciones
de frontera no nulas y constantes.
∂2U ∂U
k 2 = , 0 < x < L, k > 0, t > 0 ,
∂x ∂t
U(0, t) = K1 , U(L, t) = K2 , t > 0 ,
U(x, 0) = f(x), 0 < x < L .
Solución.
Sea U(x, t) = W(x, t) + ϕ(x) . Se tiene Ut (x, t) = Wt (x, t) , Uxx (x, t) = Wxx (x, t) + ϕ ''(x) .
El problema se convierte en
 ∂2 W  ∂W
k  2 + ϕ ''(x)  = , 0 < x < L, k > 0, t > 0 ,
  ∂t
 ∂x 
W(0, t) = K1 − ϕ(0), W(L, t) = K2 − ϕ(L), t > 0 ,
W(x, 0) = f(x) − ϕ(x), 0 < x < L .
Condiciones sobre ϕ(x) : ϕ ''(x) = 0, ϕ(0) = K1 , ϕ(L) = K2 .
Entonces:
K2 − K1
ϕ(x) = αx + β , ϕ(0) = K1 , ϕ(L) = K2 ⇒ β = K1 , α = .
L
El problema es:
 ∂2W  ∂W
k 2  = , 0 < x < L, k > 0, t > 0 ,
  ∂t
 ∂x 
W(0, t) = 0, W(L, t) = 0, t > 0 ,
K2 − K1
W(x, 0) = f(x) − L
x − K1 , 0 < x < L .
Solución:
∞
K2 − K1 −k( nπ )2 t
U(x, t) =
L
x + K1 + ∑
n =1
ansen(nLπx )e L ,

donde

Prof. José Luis Quintero 58

 2 L
an =
L ∫ 0
nπx
f(x) − β(x) sen( L )dx

Ejemplo 59. La temperatura u(x,t) en una varilla no aislada de longitud L, con

condiciones de frontera no nulas se determina a partir del problema de valores en la
frontera
kuxx + f(x, t) = ut , k > 0, 0 < x < L, t > 0 ,
u(0, t) = α1(t), u(L, t) = α2 (t), t > 0,
u(x, 0) = ϕ(x), 0 < x < L .
Solución.
En primer lugar, se puede escribir u(x, t) = F(x, t) + w(x, t) , donde F(x,t) es una función
que debe resolver el problema Fxx (x, t) = 0, F(0, t) = α1(t), F(L, t) = α2 (t) . De modo que
Fxx (x, t) = 0 ⇒ F(x, t) = g1(t) + xg2 (t) . Al imponer la primera condición de borde se tiene
F(0, t) = α1(t) = g1(t) . Al imponer la segunda condición de borde se tiene
α2 (t) −α1(t)
F(x, t) = α1(t) + xg2 (t) ⇒ F(L, t) = α2 (t) = α1(t) + Lg2 (t) ⇒ g2 (t) = L
.
Por lo tanto se obtiene la función definida por
F(x, t) = α1(t) + Lx α2 (t) − α1(t) .

Por otro lado w(x,t) es una función incógnita auxiliar cuya determinación se estudiará a
continuación. Como F(0, t) = α1(t) y F(L, t) = α2 (t) por construcción, resulta que
w(0, t) = 0 y w(L, t) = 0 . Además w(x, 0) = u(x, 0) − F(x, 0) = ϕ(x) − F(x, 0) en 0 < x < L , y
kwxx + f(x, t) = wt + Ft (x, t) ⇒ kwxx = wt + Ft (x, t) − f(x, t) .
Si se denotan H(x) = ϕ(x) − F(x, 0) y G(x, t) = f(x, t) − Ft (x, t) , el problema original queda
convertido como
kwxx + G(x, t) = wt , k > 0, 0 < x < L, t > 0 ,
w(0, t) = 0, w(L, t) = 0, t > 0 ,
w(x, 0) = H(x), 0 < x < L .
Este problema representa el caso de una ecuación de calor no homogénea con
condición inicial no nula y condiciones de frontera nulas.
Si se ensaya con una solución de la forma
∞
w(x, t) = ∑ w (t)sen(
n =1
n
nπx
L
),

y se expresa
∞
G(x, t) = ∑ G (t)sen(
n =1
n
nπx
L
),

con
2 L
Gn (t) =
L ∫ 0
G(x, t)sen(nLπx )dx, n = 1,... ,

entonces
∞ ∞
wxx (x, t) = − ∑
n =1
wn (t)(nLπ )2 sen(nLπx ) , wt (x, t) = ∑ w (t)sen(
n =1
'
n
nπx
L
).

Sustituyendo en la ecuación kwxx + G(x, t) = wt se tiene

Prof. José Luis Quintero 59

 ∞ ∞ ∞
−k ∑
n =1
wn(t)( nLπ )2 sen( nLπx ) + ∑
n =1
Gn(t)sen(nLπx ) = ∑ w (t)sen(
n =1
'
n
nπx
L
).

Reacomodando términos:
∞

∑ w (t) + k(
n =1
'
n
nπ 2
L
) wn (t) − Gn (t)sen(nLπx ) = 0 ⇒ wn' (t) + k(nLπ )2 wn (t) = Gn (t) .
 
Ecuación lineal de primer orden

Como
∞
w(x, t) = ∑ w (t)sen(
n =1
n
nπx
L
),

entonces de la condición inicial del problema se tiene

∞
H(x) = ∑ w (0)sen(
n =1
n
nπx
L
)

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente wn (0) como
2 L
Wn (0) =
L ∫ 0
H(x)sen(nLπx )dx .

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de
primer orden sujeta a una condición:
2 L
wn' (t) + k(nLπ )2 wn (t) = Gn (t) , Wn(0) =
L ∫ 0
H(x)sen(nLπx )dx

Método 1.
Solución general:
−k(nπ )2 t
 k( nπ )2 t  −k( nπ )2 t
wn(t) = e

L
∫
 Gn (t)e L dt + C  = e L

Mn(t) + C .

Valor C: C = wn (0) − Mn (0) . Solución particular:

−k( nπ )2 t
wn(t) = e L Mn (t) + wn (0) − Mn(0) .
Resumen de pasos a seguir:
1. Construcción de G(x,t). G(x, t) = f(x, t) − Ft (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)
5. Cálculo de Mn (t) y Mn (0) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

Método 2.
L{wn' (t)} + k(nLπ )2 L{wn (t)} = L{Gn (t)} ⇒ swn (s) − wn (0) + k(nLπ )2 wn(s) = Gn(s) .
Gn (s) + wn (0) Gn(s) wn (0)
s + k(nπ )2  wn(s) = Gn (s) + wn (0) ⇒ wn(s) = = + .
 L  s+ k( nLπ )2 s+ k(nLπ )2 s + k(nLπ )2

Prof. José Luis Quintero 60

  G (s)  −k( nπ )2 t −k( nπ )2 t
wn(t) = L-1{wn (s)} = L-1  n
 + wn (0)e L = Pn (t) + wn (0)e L .
nπ 2
 s + k( L ) 

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = f(x, t) − Ft (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)
5. Cálculo de L {Gn (t)} y Pn (t) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

Solución final:
∞
u(x, t) = F(x, t) + w(x, t) = α1(t) + x
L α2 (t) − α1(t) + ∑ w (t)sen(
n =1
n
nπx
L
).

Ejemplo 60. (Caso particular del ejemplo 59) Ecuación de calor homogénea con
condición inicial no nula y condiciones de frontera no nulas.
Solución.
En este caso G(x, t) = −Ft (x, t) y se aplican cualquiera de los métodos descritos en el
ejemplo 59.

Ejemplo 61. La temperatura u(x,t) en una varilla de longitud L con extremos aislados.
Paso 1. Planteamiento del problema.
kuxx = ut , k > 0, 0 < x < L, t > 0 ,
ux (0, t) = 0, ux (L, t) = 0, t > 0 ,
u(x, 0) = f(x), 0 < x < L .
Paso 2. Planteamiento y resolución de las ecuaciones diferenciales ordinarias.
Si se ensaya con soluciones del tipo u(x, t) = X(x)T(t) con X(x) ≠ 0 , T(t) ≠ 0 se tiene
kX '' T = XT ' . Al dividir cada miembro entre kXT e igualando a una constante negativa
se obtiene
X '' T '
= = −λ2 ,
X kT
lo que origina las dos ecuaciones diferenciales ordinarias:
2t
X ''+ λ2 X = 0 ⇒ X(x) = c1 cos(λx) + c2sen(λx) , T '+ kλ2 T = 0 ⇒ T(t) = c3e−kλ
Paso 3. Utilización de las condiciones de borde.
X ''+ λ2 X = 0 , X '(0) = 0 , X '(L) = 0 (Problema de Sturm-Liouville clásico)
Solución. Xn (x) = an cos( nLπx ) n = 0,1, 2,...
Paso 4. Aplicación del principio de superposición.
−k(nπ )2 t
un (x, t) = bn cos(nLπx )e L ,

entonces

Prof. José Luis Quintero 61

 ∞ ∞
−k(nπ )2 t
u(x, t) = ∑
n =0
un (x, t) = ∑b
n =0
n cos(nLπx )e L .

Paso 5. Utilización de la condición inicial.

∞ ∞ ∞
A0
u(x, 0) = f(x) = ∑
n=0
bn cos( nLπx ) = b0 + ∑
n =1
bn cos(nLπx ) =
2
+ ∑A
n =1
n cos(nLπx ) ,

donde
2 L 2 L
A0 =
L ∫ 0
f(x)dx , An =
L ∫ 0
f(x) cos(nLπx )dx .

Serie de cosenos de Fourier para f(x)

Paso 6. Solución del problema planteado.
∞
 L  −k(nπ )2 t
∑∫
2
u(x, t) =  f(x) cos(nLπx )dx  cos(nLπx )e L .
L n=0  0 

Ejemplo 62.
kuxx + f(x, t) = ut , k > 0, 0 < x < L, t > 0 ,
ux (0, t) = ux (L, t) = α(t), t > 0,
u(x, 0) = ϕ(x), 0 < x < L .
Solución.
En primer lugar, se puede escribir u(x, t) = F(x, t) + w(x, t) , donde F(x,t) es una función
que debe resolver el problema Fxx (x, t) = 0, Fx (0, t) = Fx (L, t) = α(t) . De modo que
Fxx (x, t) = 0 ⇒ Fx (x, t) = g2 (t) ⇒ F(x, t) = g1(t) + xg2 (t) . Al imponer la primera condición
de borde se tiene Fx (0, t) = α(t) = g2 (t) . Al imponer la segunda condición de borde se
tiene F(x, t) = g1(t) + xα(t) ⇒ Fx (L, t) = α(t) = α(t) . Por lo tanto se obtiene la función
definida por F(x, t) = g1(t) + xα(t) , donde g1(t) es una función arbitraria que por
comodidad se asumirá g1(t) = 0 , de modo que F(x, t) = α(t)x . Por otro lado w(x,t) es
una función incógnita auxiliar cuya determinación se estudiará a continuación.

Como Fx (0, t) = Fx (L, t) = α(t) , resulta que wx (0, t) = 0 y wx (L, t) = 0 . Además se tiene
que w(x, 0) = u(x, 0) − F(x, 0) = ϕ(x) − F(x, 0) en 0 < x < L , y
kwxx + f(x, t) = wt + Ft (x, t) ⇒ kwxx = wt + Ft (x, t) − f(x, t) .

Si se denotan H(x) = ϕ(x) − F(x, 0) y G(x, t) = f(x, t) − Ft (x, t) , el problema original queda
convertido como
kwxx + G(x, t) = wt , k > 0, 0 < x < L, t > 0 ,
wx (0, t) = 0, wx (L, t) = 0, t > 0 ,
w(x, 0) = H(x), 0 < x < L .
Este problema representa el caso de una ecuación de calor no homogénea con
condición inicial no nula y extremos aislados.

Si se ensaya con una solución de la forma

∞
w(x, t) = ∑ w (t) cos(
n=0
n
nπx
L
),

y se expresa

Prof. José Luis Quintero 62

 ∞
G(x, t) = ∑ G (t) cos(
n =0
n
nπx
L
),

con
2 L
Gn (t) =
L ∫ 0
G(x, t) cos(nLπx )dx, n = 0,... ,

se tiene
∞ ∞
wxx (x, t) = − ∑
n=0
wn (t)(nLπ )2 cos(nLπx ) , wt (x, t) = ∑ w (t) cos(
n=0
'
n
nπx
L
).

Sustituyendo en la ecuación kwxx + G(x, t) = wt se tiene

∞ ∞ ∞
−k ∑
n =0
wn(t)(nLπ )2 cos(nLπx ) + ∑
n=0
Gn (t) cos(nLπx ) = ∑ w (t) cos(
n =0
'
n
nπx
L
).

Reacomodando términos:
∞

∑ w (t) + k(
n =1
'
n
nπ 2
L
) wn (t) − Gn (t) cos(nLπx ) = 0 ⇒ wn' (t) + k( nLπ )2 wn (t) = Gn (t) .
 
Ecuación lineal de primer orden

Como
∞
w(x, t) = ∑ w (t) cos(
n=0
n
nπx
L
),

entonces de la condición inicial del problema se tiene

∞
H(x) = ∑ w (0) cos(
n=0
n
nπx
L
)

(serie de cosenos de Fourier)

y por tanto se obtiene el coeficiente wn (0) como
2 L
Wn (0) =
L ∫ 0
H(x) cos(nLπx )dx .

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de
primer orden sujeta a una condición:
2 L
wn' (t) + k(nLπ )2 wn(t) = Gn (t) , Wn (0) =
L 0
H(x) cos(nLπx )dx
∫
Método 1.
Solución general:
−k(nπ )2 t  k( nπ )2 t  −k( nπ )2 t
wn(t) = e L

∫
 Gn (t)e L dt + C  = e L

Mn(t) + C .

−k( nπ )2 t
Valor C: C = wn (0) − Mn (0) . Solución particular: wn(t) = e L
Mn (t) + wn (0) − Mn(0) .

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = f(x, t) − Ft (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)

Prof. José Luis Quintero 63

5. Cálculo de Mn (t) y Mn (0) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

Método 2.
L{wn' (t)} + k(nLπ )2 L{wn (t)} = L{Gn (t)} ⇒ swn (s) − wn (0) + k(nLπ )2 wn(s) = Gn(s) .
Gn (s) + wn (0) Gn(s) wn (0)
s + k(nπ )2  wn(s) = Gn (s) + wn (0) ⇒ wn(s) = = + .
 L  s + k( nLπ )2 s + k(nLπ )2 s + k(nLπ )2
 G (s)  −k(nπ )2 t −k( nπ )2 t
wn(t) = L-1{wn (s)} = L-1  n
 + wn (0)e L = Pn (t) + wn (0)e L .
nπ 2
 s + k( L ) 

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = f(x, t) − Ft (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)
5. Cálculo de L {Gn (t)} y Pn (t) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

∞
Solución final. u(x, t) = F(x, t) + w(x, t) = α1(t)x + w0 (t) + ∑ w (t) cos(
n =1
n
nπx
L
).

Ejemplo 63.
Paso 1. Planteamiento del problema.
kuxx = ut , k > 0, 0 < x < L, t > 0 ,
u(0, t) = 0, ux (L, t) = 0, t > 0 ,
u(x, 0) = f(x), 0 < x < L .
Paso 2. Planteamiento y resolución de las ecuaciones diferenciales ordinarias.
Si se ensaya con soluciones del tipo u(x, t) = X(x)T(t) con X(x) ≠ 0 , T(t) ≠ 0 se tiene
kX '' T = XT ' . Al dividir cada miembro entre kXT e igualando a una constante negativa
se obtiene
X '' T '
= = −λ2 ,
X kT
lo que origina las dos ecuaciones diferenciales ordinarias:
2t
X ''+ λ2 X = 0 ⇒ X(x) = c1 cos(λx) + c2sen(λx) , T '+ kλ2 T = 0 ⇒ T(t) = c3e−kλ
Paso 3. Utilización de las condiciones de borde.
X ''+ λ2 X = 0 , X(0) = 0 , X '(L) = 0 (Problema de Sturm-Liouville clásico)
(2n + 1)πx
Solución. Xn (x) = ansen( 2L
) n = 0,1, 2,...
Paso 4. Aplicación del principio de superposición.

Prof. José Luis Quintero 64

 (2n +1)π 2 ∞ ∞ (2n +1)π 2
−k(
∑ ∑ b sen(
(2n + 1)πx ) t (2n +1)πx −k( ) t
un (x, t) = bnsen( 2L
)e 2L , u(x, t) = un (x, t) = n 2L
)e 2L .
n =0 n =0

Paso 5. Utilización de la condición inicial.

∞ L
∑
2
∫
+1)πx (2n +1)πx
u(x, 0) = f(x) = bnsen((2n2L ) , donde bn = f(x)sen( 2L
)dx .
n=0
L 0

Serie de senos de Fourier para f(x).

Paso 6. Solución del problema planteado.
∞ (2n +1)π 2
 L 
∑∫
2 +1)πx + 1)πx −k( ) t
u(x, t) =  f(x)sen((2n2L )dx  sen((2n2L )e 2L .
L n=0  0 

Ejemplo 64.
kuxx + f(x, t) = ut , k > 0, 0 < x < L, t > 0 ,
u(0, t) = α1(t), ux (L, t) = α2 (t), t > 0,
u(x, 0) = ϕ(x), 0 < x < L .
Solución.
En primer lugar, se puede escribir u(x, t) = F(x, t) + w(x, t) , donde F(x,t) es una función
que debe resolver el problema Fxx (x, t) = 0, F(0, t) = α1(t), Fx (L, t) = α2 (t) . De modo que
Fxx (x, t) = 0 ⇒ F(x, t) = g1(t) + xg2 (t) . Al imponer la primera condición de borde se tiene
F(0, t) = α1(t) = g1(t) . Al imponer la segunda condición de borde se tiene
F(x, t) = α1(t) + xg2 (t) ⇒ Fx (L, t) = α2 (t) = g2 (t) . Por tanto se obtiene la función definida
por F(x, t) = α1(t) + α2 (t)x . Por otro lado w(x,t) es una función incógnita auxiliar cuya
determinación se estudiará a continuación.

Como F(0, t) = α1(t) y Fx (L, t) = α2 (t) por construcción, resulta que w(0, t) = 0 y
wx (L, t) = 0 . Además w(x, 0) = u(x, 0) − F(x, 0) = ϕ(x) − F(x, 0) en 0 < x < L , y
kwxx + f(x, t) = wt + Ft (x, t) ⇒ kwxx = wt + Ft (x, t) − f(x, t) .

Si se denotan H(x) = ϕ(x) − F(x, 0) y G(x, t) = f(x, t) − Ft (x, t) , el problema original queda
convertido como
kwxx + G(x, t) = wt , k > 0, 0 < x < L, t > 0 ,
w(0, t) = 0, wx (L, t) = 0, t > 0 ,
w(x, 0) = H(x), 0 < x < L .

Este problema representa el caso de una ecuación de calor no homogénea con

condición inicial no nula, condición de frontera nula y extremo aislado.

Si se ensaya con una solución de la forma

∞
w(x, t) = ∑ w (t)sen(
n=0
n
(2n +1)πx
2L
),

y se expresa
∞
G(x, t) = ∑ G (t)sen(
n =0
n
(2n +1)πx
2L
),

con

Prof. José Luis Quintero 65

 2 L

∫
(2n +1)πx
Gn (t) = G(x, t)sen( 2L
)dx, n = 0,... ,
L 0
∞ ∞
wxx (x, t) = − ∑n=0
+1)π 2
wn (t)((2n2L +1)πx
) sen((2n2L ) , wt (x, t) = ∑ w (t)sen(
n=0
'
n
(2n +1)πx
2L
).

Sustituyendo en la ecuación kwxx + G(x, t) = wt se tiene

∞ ∞ ∞
−k ∑
n =0
+ 1)π 2
wn(t)((2n2L +1)πx
) sen((2n2L )+ ∑
n=0
+ 1)πx
Gn(t)sen((2n2L )= ∑ w (t)sen(
n= 0
'
n
(2n +1)πx
2L
).

Reacomodando términos:
∞

∑ w (t) + k(
n=0
'
n
nπ 2
L
) wn(t) − Gn (t)sen(

(2n +1)πx
2L
) = 0 ⇒ wn' (t) + k(nLπ )2 wn(t) = Gn (t) .

Ecuación lineal de primer orden

Como
∞
w(x, t) = ∑ w (t)sen(
n=0
n
(2n +1)πx
2L
),

entonces de la condición inicial dada en el problema se tiene

∞
H(x) = ∑ w (0)sen(
n=0
n
(2n +1)πx
2L
)

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente wn (0) como
2 L

∫
(2n +1)πx
Wn (0) = H(x)sen( 2L
)dx .
L 0

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de
primer orden sujeta a una condición:
2 L
∫
(2n +1)πx
wn' (t) + k(nLπ )2 wn(t) = Gn (t) , Wn (0) = H(x)sen( 2L )dx .
L 0

Método 1.
Solución general:
(2n + 1)π 2  (2n +1)π 2  (2n +1)π 2
−k( ) t k( ) t −k( ) t
wn(t) = e 2L
∫
 Gn (t)e

2L dt + C  = e

2L Mn (t) + C .

(2n +1)π 2
−k( ) t
Valor C: C = wn (0) − Mn (0) . Sol. particular: wn(t) = e 2L
Mn (t) + wn (0) − Mn (0) .

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = f(x, t) − Ft (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)
5. Cálculo de Mn (t) y Mn (0) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

Prof. José Luis Quintero 66

Método 2.
+1)π 2 + 1)π 2
L{wn' (t)} + k((2n2L ) L{wn (t)} = L{Gn (t)} ⇒ swn (s) − wn (0) + k((2n2L ) wn (s) = Gn (s) .
s + k((2n +1)π )2  w (s) = G (s) + w (0)
 2L  n n n

G (s) + wn (0) Gn (s) wn (0)

⇒ wn(s) = n (2n +1) π 2
= (2n +1)π 2
+ + 1)π 2
s + k( 2L ) s + k( 2L ) s + k((2n2L )
 Gn (s)  −k(
(2n +1)π 2
) t −k(
(2n + 1)π 2
) t
wn(t) = L-1{wn (s)} = L-1  (2n + 1)π 2  + wn (0)e 2L = Pn (t) + wn (0)e 2L .
 s + k( 2L ) 

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = f(x, t) − Ft (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)
5. Cálculo de L {Gn (t)} y Pn (t) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).
∞
Solución final: u(x, t) = F(x, t) + w(x, t) = α1(t) + α2 (t)x + ∑ w (t)sen(
n =0
n
(2n +1)πx
2L
).

Ejemplo 65.
Paso 1. Planteamiento del problema.
kuxx = ut , k > 0, 0 < x < L, t > 0 ,
ux (0, t) = 0, u(L, t) = 0, t > 0,
u(x, 0) = f(x), 0 < x < L .
Paso 2. Planteamiento y resolución de las ecuaciones diferenciales ordinarias.
Si se ensaya con soluciones del tipo u(x, t) = X(x)T(t) con X(x) ≠ 0 , T(t) ≠ 0 se tiene
kX '' T = XT ' . Al dividir cada miembro entre kXT e igualando a una constante negativa
se obtiene
X '' T '
= = −λ2 ,
X kT
lo que origina las dos ecuaciones diferenciales ordinarias:
2t
X ''+ λ2 X = 0 ⇒ X(x) = c1 cos(λx) + c2sen(λx) , T '+ kλ2 T = 0 ⇒ T(t) = c3e−kλ
Paso 3. Utilización de las condiciones de borde.
X ''+ λ2 X = 0 , X '(0) = 0 , X(L) = 0 (Problema de Sturm-Liouville clásico)
(2n +1)πx
Solución. Xn (x) = bn cos( 2L
) n = 0,1, 2,...
Paso 4. Aplicación del principio de superposición.
(2n +1)π 2 ∞ ∞ (2n + 1)π 2
−k(
∑ u (x, t) = ∑ a
(2n +1)πx ) t +1)πx −k( ) t
un (x, t) = an cos( 2L )e 2L , u(x, t) = n n cos((2n2L )e 2L .
n=0 n=0

Paso 5. Utilización de la condición inicial.

Prof. José Luis Quintero 67

 ∞ L
∑a
2
∫
+1)πx (2n +1)πx
u(x, 0) = f(x) = n cos((2n2L ) , donde an = f(x) cos( 2L
)dx .
n =0
L 0

Serie de cosenos de Fourier para f(x).

Paso 6. Solución del problema planteado.
∞ (2n +1)π 2
 L 
∑∫
2 +1)πx +1)πx −k( ) t
u(x, t) =  f(x) cos((2n2L )dx  cos((2n2L )e 2L .
L n=0  0 

Ejemplo 66.
kuxx + f(x, t) = ut , k > 0, 0 < x < L, t > 0 ,
ux (0, t) = α1(t), u(L, t) = α2 (t), t > 0,
u(x, 0) = ϕ(x), 0 < x < L .
Solución.
En primer lugar, se puede escribir u(x, t) = F(x, t) + w(x, t) , donde F(x,t) es una función
que debe resolver el problema Fxx (x, t) = 0, Fx (0, t) = α1(t), F(L, t) = α2 (t) . De modo que
Fxx (x, t) = 0 ⇒ F(x, t) = g1(t) + xg2 (t) . Al imponer la primera condición de borde se tiene
Fx (0, t) = α1(t) = g2 (t) . Al imponer la segunda condición de borde se tiene
F(x, t) = g1(t) + xα1(t) ⇒ F(L, t) = α2 (t) = g1(t) + Lα1(t) ⇒ g1(t) = α2 (t) − Lα1(t) . Por lo tanto
se obtiene la función definida por F(x, t) = α2 (t) + (x − L)α1(t) . Por otro lado w(x,t) es
una función incógnita auxiliar cuya determinación se estudiará a continuación.

Como Fx (0, t) = α1(t) y F(L, t) = α2 (t) por construcción, resulta que wx (0, t) = 0 y
w(L, t) = 0 . Además w(x, 0) = u(x, 0) − F(x, 0) = ϕ(x) − F(x, 0) en 0 < x < L , y
kwxx + f(x, t) = wt + Ft (x, t) ⇒ kwxx = wt + Ft (x, t) − f(x, t) .

Si se denotan H(x) = ϕ(x) − F(x, 0) y G(x, t) = f(x, t) − Ft (x, t) , el problema original queda
convertido como
kwxx + G(x, t) = wt , k > 0, 0 < x < L, t > 0 ,
wx (0, t) = 0, w(L, t) = 0, t > 0 ,
w(x, 0) = H(x), 0 < x < L .

Este problema representa el caso de una ecuación de calor no homogénea con

condición inicial no nula, condición de frontera nula y extremo aislado.

Si se ensaya con una solución de la forma

∞
w(x, t) = ∑ w (t) cos(
n=0
n
(2n + 1)πx
2L
),

y se expresa
∞ L
∑ G (t) cos(
2
∫
(2n +1)πx (2n + 1)πx
G(x, t) = n 2L
), con Gn (t) = G(x, t) cos( 2L
)dx, n = 0,... ,
n =0
L 0

∞ ∞
wxx (x, t) = − ∑
n=0
+1)π 2
wn (t)((2n2L +1)πx
) cos((2n2L ) , wt (x, t) = ∑ w (t) cos(
n=0
'
n
(2n +1)πx
2L
).

Sustituyendo en la ecuación kwxx + G(x, t) = wt se tiene

Prof. José Luis Quintero 68

 ∞ ∞ ∞
−k ∑
n=0
+ 1)π 2
wn(t)((2n2L +1)πx
) cos((2n2L )+ ∑
n=0
+ 1)πx
Gn (t) cos((2n2L )= ∑ w (t) cos(
n =0
'
n
(2n +1)πx
2L
).

Reacomodando términos:
∞

∑ w (t) + k(
n=0
'
n
nπ 2
L
) wn (t)

(2n + 1)πx
− Gn (t) cos( 2L ) = 0 ⇒ wn' (t) + k(nLπ )2 wn (t) = Gn (t) .

Ecuación lineal de primer orden

Como
∞
w(x, t) = ∑ w (t) cos(
n=0
n
(2n + 1)πx
2L
),

entonces de la condición inicial dada en el problema se tiene

∞
H(x) = ∑ w (0) cos(
n=0
n
(2n +1)πx
2L
)

(serie de cosenos de Fourier)

y por tanto se obtiene el coeficiente wn (0) como
2 L

∫
(2n +1)πx
Wn (0) = H(x) cos( 2L
)dx .
L 0

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de
primer orden sujeta a una condición:
2 L

∫
(2n +1)πx
wn' (t) + k(nLπ )2 wn(t) = Gn (t) , Wn (0) = H(x) cos( 2L
)dx
L 0

Método 1.
Solución general:
(2n + 1)π 2  (2n +1)π 2  (2n +1)π 2
−k( ) t k( ) t −k( ) t
wn(t) = e 2L
∫
 Gn (t)e

2L dt + C  = e

2L
Mn (t) + C .

(2n + 1)π 2
−k( ) t
Valor C: C = wn (0) − Mn (0) . Sol. particular: wn(t) = e 2L Mn (t) + wn (0) − Mn (0) .

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = f(x, t) − Ft (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)
5. Cálculo de Mn (t) y Mn (0) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

Método 2.
(2n +1)π 2 (2n + 1)π 2
L{wn' (t)} + k( 2L
) L{wn (t)} = L{Gn (t)} ⇒ swn (s) − wn (0) + k( 2L
) wn (s) = Gn (s) .

Prof. José Luis Quintero 69

 s + k((2n +1)π )2  w (s) = G (s) + w (0)
 2L  n n n

G (s) + wn (0) Gn (s) wn (0)

⇒ wn(s) = n (2n +1) π 2
= (2n +1)π 2
+ + 1)π 2
s + k( 2L ) s + k( 2L ) s + k((2n2L )
 Gn (s)  −k(
(2n +1)π 2
) t −k(
(2n + 1)π 2
) t
wn(t) = L-1{wn (s)} = L-1  (2n + 1)π 2  + wn (0)e 2L = Pn (t) + wn (0)e 2L .
 s + k( 2L ) 

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = f(x, t) − Ft (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)
5. Cálculo de L {Gn (t)} y Pn (t) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).
∞
Solución final: u(x, t) = F(x, t) + w(x, t) = α2 (t) + (x − L)α1(t) + ∑ w (t) cos(
n=0
n
(2n +1)πx
2L
).

Ejemplo 67. Sea una varilla de aluminio de longitud L que inicialmente está a la
temperatura uniforme de 25 C . Suponga que, en el instante t = 0 , el extremo x = 0
se enfría hasta 0 C , mientras que el extremo x = L se calienta hasta 60 C y que, a
partir de entonces, los dos se mantienen a esas temperaturas. Demuestre que la
distribución de temperatura u(x,t) viene dada como
∞
 nπx 
∑  nπ + nπ (−1)  e
60x 50 70 n −0.86n2 π2t / L2
u(x, t) = + sen  .
L n =1  L 
Solución.
Planteamiento del problema:
 0.86uxx = ut , 0 < x < L , t > 0

u(0, t) = 0 , u(L, t) = 60 , t > 0 .
 u(x, 0) = 25 , 0 < x < L

Sea u(x, t) = v(x, t) + w(x, t) . Se tiene entonces el siguiente problema:
 0.86wxx + 0.86vxx = wt + vt , 0 < x < L , t > 0

w(0, t) = 0 − v(0, t) , w(L, t) = 60 − v(L, t) , t > 0 .
 w(x, 0) = 25 − v(x, 0) , 0 < x < L

Exigencias sobre v(x,t): v '' = 0, v(0) = 0, v(L) = 60 . v(x) = 60x L

El nuevo problema sería:

0.86wxx = wt , 0 < x < L , t > 0

 w(0, t) = 0, w(L, t) = 0 , t > 0
 w(x, 0) = 25 − 60x , 0 < x < L
 L

cuya solución viene dada por

Prof. José Luis Quintero 70

 ∞
w(x, t) = ∑c e
n =1
n
−0.86n2 π2t / L2
sen(nLπx ) ,

L
cn = 2
L ∫ 0
f(x)sen(nLπx )dx .

En este caso
L
cn = 2
L ∫ 0
(25 − 60x
L
)sen(nLπx )dx .

Al calcular cn se tiene
L L L
cn = 2
L ∫ 0
(25 − 60x
L
)sen(nLπx )dx = 50
L ∫ 0
sen(nLπx )dx − 120
L ∫ 0
xsen(nLπx )dx

cn = − 50
nπ
cos(nπ) + 50
nπ
+ 120L
nπ
cos(nπ) = 50
nπ
+ 70L
nπ
(−1)n .
Por tanto
∞
 nπx 
∑  nπ + nπ (−1)  e
60x 50 70 −0.86n2 π2t / L2
u(x, t) = + n
sen  .
L n =1  L 

Ejemplo 68. Dado el problema de valor de frontera

ut = kuxx − hu , 0 < x < π , t > 0 , k > 0 , h > 0

u(0, t) = u(π, t) = 0 , t ≥ 0,
 x
u(x, 0) = e , 0 < x < π
haga u(x, t) = e−ht w(x, t) y obtenga su solución.
Solución.
u = e−ht w , ut = −he−ht w + e−ht wt , uxx = e−ht wxx , u(x, 0) = w(x, 0) = ex ,
u(0, t) = e−ht w(0, t) = 0 ⇒ w(0, t) = 0 , u(L, t) = e−ht w(L, t) = 0 ⇒ w(L, t) = 0 .
En consecuencia: −he−ht w + e−ht wt = ke−ht wxx − he−ht w ⇒ wt = kwxx .
Problema de valor en la frontera equivalente:
wt = kwxx , 0 < x < π , t > 0 , k > 0

w(0, t) = w(π, t) = 0 , t ≥ 0,
 x
w(x, 0) = e , 0 < x < π
Sea w(x, t) = X(x)T(t) . Por tanto wxx = X '' T y wt = XT ' . De modo que kX '' T = XT ' . Al
dividir por kXT ≠ 0 se obtiene
X '' T '
= .
X kT
Igualando a una constante de separación:
X '' T '
= =K,
X kT
se tienen el par de ecuaciones ordinarias: X ''− KX = 0 , T '+ kKT = 0 . Se estudian 3
casos:
Caso I: K > 0 , K = λ2 .
Ecuaciones ordinarias:
2
X ''− λ2 X = 0 ⇒ X(x) = Aeλx + Be−λx , T '+ kλ2 T = 0 ⇒ T(t) = Ce−kλ t .
Condiciones de borde homogéneas:

Prof. José Luis Quintero 71

w(0, t) = 0 ⇒ X(0) = 0 ⇒ A + B = 0 ⇒ A = −B
w(π, t) = 0 ⇒ X(π) = 0 ⇒ Aeπλ + Be−πλ = A(eπλ − e−πλ ) = 0 ⇒ A = 0 ⇒ B = 0
Solo se obtiene solución trivial.
Caso II: K = 0 .
Ecuaciones ordinarias: X '' = 0 ⇒ X(x) = A + Bx , T ' = 0 ⇒ T(t) = C .
Condiciones de borde homogéneas:
w(0, t) = 0 ⇒ X(0) = 0 ⇒ A = 0 , w(π, t) = 0 ⇒ X(π) = 0 ⇒ Bπ = 0 ⇒ B = 0 .
Solo se obtiene solución trivial.
Caso III: K < 0 , K = −λ2 .
Ecuaciones ordinarias:
2
X ''+ λ2 X = 0 ⇒ X(x) = A cos(λx) + Bsen(λx) , T '− kλ2 T = 0 ⇒ T(t) = Ce−kλ t .
Condiciones de borde homogéneas:
w(0, y) = 0 ⇒ X(0) = 0 ⇒ A = 0
w(π, t) = 0 ⇒ X(π) = 0 ⇒ Bnsen(λπ) = 0 ⇒ λ = n , n = 1,2, 3,...
Se tiene:
∞ ∞

∑ ∑D e
2kt −n2kt
w(x, t) = Bnsen(nx)Cne−n , w(x, t) = n sen(nx)
n =1 n =1

Condición inicial:
∞
w(x, 0) = ex = ∑ D sen(nx) .
n =1
n

En consecuencia:
π π
2 2 n  sen(nx) 
Dn =
π ∫ 0
ex sen(nx)dx = .
π n2 + 1
ex 
 n
− cos(nx)
0
2 n 2 n
n impar = . 2 (1 + eπ ) n par = . 2 (1 − eπ ).
π n +1 π n +1
Por tanto
2 n
Dn = . 2 (1 + (−1)n +1 eπ )
π n +1
Se obtiene:
∞

∑
2 n 2
w(x, t) = (1 + (−1)n +1 eπ )e−n kt sen(nx)
π n =1 n2 + 1
∞

∑
2 −ht n 2
u(x, t) = e (1 + (−1)n +1 eπ )e−n kt sen(nx) .
π n =1
2
n +1

Ejemplo 69. La temperatura u(x,t) en una varilla no aislada de longitud L, con

condiciones de frontera no nulas, con flujo de calor en la superficie de la varilla en
estudio hacia el medio que la rodea cuya temperatura se mantiene constante, digamos
um , se determina a partir del problema de valores en la frontera
kuxx − h(u − um ) = ut , k > 0, h ctte, 0 < x < L, t > 0 ,
u(0, t) = α1(t), u(L, t) = α2 (t), t > 0,
u(x, 0) = ϕ(x), 0 < x < L .
Solución.

Prof. José Luis Quintero 72

En primer lugar, se puede escribir u(x, t) = F(x, t) + w(x, t) , donde F(x,t) es una función
que debe resolver el problema Fxx (x, t) = 0, F(0, t) = α1(t), F(L, t) = α2 (t) . De modo que
Fxx (x, t) = 0 ⇒ F(x, t) = g1(t) + xg2 (t) . Al imponer la primera condición de borde se tiene
F(0, t) = α1(t) = g1(t) . Al imponer la segunda condición de borde se tiene
α2 (t) −α1(t)
F(x, t) = α1(t) + xg2 (t) ⇒ F(L, t) = α2 (t) = α1(t) + Lg2 (t) ⇒ g2 (t) = L
.
Por lo tanto se obtiene la función definida por F(x, t) = α1(t) + x
L
α2 (t) − α1(t) . Por otro
lado w(x,t) es una función incógnita auxiliar cuya determinación se estudiará a
continuación.

Como F(0, t) = α1(t) y F(L, t) = α2 (t) por construcción, resulta que w(0, t) = 0 y
w(L, t) = 0 . Además w(x, 0) = u(x, 0) − F(x, 0) = ϕ(x) − F(x, 0) en 0 < x < L , y
kwxx + f(x, t) = wt + Ft (x, t) ⇒ kwxx = wt + Ft (x, t) − f(x, t) .

Si se denotan H(x) = ϕ(x) − F(x, 0) y G(x, t) = −(h(F(x, t) − um ) + Ft (x, t)) , se tiene que el
problema original queda convertido como
kwxx + G(x, t) − hw = wt , k > 0, 0 < x < L, t > 0 ,
w(0, t) = 0, w(L, t) = 0, t > 0 ,
w(x, 0) = H(x), 0 < x < L .

Este problema representa el caso de una ecuación de calor no homogénea con

condición inicial no nula y condiciones de frontera nulas.

Si se ensaya con una solución de la forma

∞
w(x, t) = ∑ w (t)sen(
n =1
n
nπx
L
),

y se expresa
∞
G(x, t) = ∑ G (t)sen(
n =1
n
nπx
L
),

con
2 L
Gn (t) =
L ∫ 0
G(x, t)sen(nLπx )dx, n = 1,... ,

entonces
∞ ∞
wxx (x, t) = − ∑
n =1
wn (t)(nLπ )2 sen(nLπx ) , wt (x, t) = ∑ w (t)sen(
n =1
'
n
nπx
L
).

Sustituyendo en la ecuación kwxx + G(x, t) − hw = wt , se tiene

∞ ∞ ∞ ∞
−k ∑ n =1
wn(t)(nLπ )2 sen(nLπx ) + ∑
n =1
Gn (t)sen(nLπx ) − h ∑
n =1
wn(t)sen(nLπx ) = ∑ w (t)sen(
n =1
'
n
nπx
L
).

Reacomodando términos:
∞

∑ w (t) + (k(

n =1
'
n
nπ 2
L
) + h)wn (t) − Gn(t)sen(nLπx ) = 0 ⇒ wn' (t) + (k(nLπ )2 + h)wn (t) = Gn (t) .
 
Ecuación lineal de primer orden

Como

Prof. José Luis Quintero 73

 ∞
w(x, t) = ∑ w (t)sen(
n =1
n
nπx
L
),

entonces de la condición inicial del problema se tiene

∞
H(x) = ∑ w (0)sen(
n =1
n
nπx
L
)

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente wn (0) como
2 L
Wn (0) =
L ∫ 0
H(x)sen(nLπx )dx .

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de
primer orden sujeta a una condición:
2 L
wn' (t) + (k(nLπ )2 + h)wn(t) = Gn (t) , Wn (0) =
L ∫ 0
H(x)sen(nLπx )dx

Método 1.
−(k( nπ )2 + h)t  (k(nπ )2 + h)t  −(k( nπ )2 + h)t
Solución general: wn(t) = e L
∫
 Gn (t)e L

dt + C  = e

L Mn (t) + C .

−(k( nπ )2 + h)t
Valor C: C = wn (0) − Mn (0) . Sol. particular: wn(t) = e L Mn (t) + wn (0) − Mn (0) .
Resumen de pasos a seguir:
1. Construcción de G(x,t). G(x, t) = −(h(F(x, t) − um ) + Ft (x, t))
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)
5. Cálculo de Mn (t) y Mn (0) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

Método 2.
L{wn' (t)} + (k(nLπ )2+h)L{wn (t)} = L{Gn (t)} ⇒ swn (s) − wn (0) + (k( nLπ )2 + h)wn (s) = Gn (s) .
s + k(nπ )2 + h wn (s) = Gn (s) + wn (0)
 L 
G (s) + wn (0) Gn (s) wn (0)
⇒ wn(s) = n = +
s + k(nLπ )2 + h s + k(nLπ )2 + h s + k(nLπ )2 + h
 Gn(s)  −(k(nπ )2 + h)t −(k(nπ )2 + h)t
wn(t) = L-1{wn (s)} = L-1   + wn (0)e L = Pn(t) + wn(0)e L .
nπ 2
 s + k( L
) + h 

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = −(h(F(x, t) − um ) + Ft (x, t))
2. Cálculo de los Gn (t) .
3. Construcción de H(x). H(x) = ϕ(x) − F(x, 0)
4. Cálculo de wn (0)
5. Cálculo de L {Gn (t)} y Pn (t) .

Prof. José Luis Quintero 74

6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

∞
Solución final: u(x, t) = F(x, t) + w(x, t) = α1(t) + x
L
α2 (t) − α1(t) + ∑ w (t)sen(
n =1
n
nπx
L
).

Ejemplo 70. Encuentre la solución u(x,t) de la ecuación de calor unidimensional que

satisface las condiciones de frontera dadas por
 uxx + 1 = ut , 0 < x < π , t > 0

u(0, t) = t , u(π, t) = 0 , t > 0 .
 u(x, 0) = x , 0 < x < π

Solución.
Paso 1. Construcción de G(x,t). F(x, t) = t − x
π
t = t(1 − xπ ) . G(x, t) = 1 − 1 + x
π
= x
π
.
Paso 2. Cálculo de los Gn (t) .
π
 x π
Gn (t) =
0
2
π2 ∫ 
xsen(nx)dx =
 − n cos(nx) +
2
π2
1
n2
sen(nx)  =
0
2
πn
(−1)n +1 .

Paso 3. Construcción de H(x). H(x) = x − 0 = x .

π
Paso 4. Cálculo de wn (0) . wn(0) = 2
π ∫ 0
xsen(nx)dx = 2
n
(−1)n +1 .

Paso 5. Cálculo de Mn (t) y Mn (0) . (Método 1)

2 2t
Mn (t) = 2
πn
(−1)n +1
∫ en t dt = 2
πn3
(−1)n +1 en , Mn (0) =
πn3
2
(−1)n +1 .
Paso 6. Cálculo de los wn(t) . (Método 1)

wn(t) = e−n  2 (−1)n +1 en2t + 2 (−1)n +1 − 2 (−1)n +1  =

2t
2
(−1)n +1 e−n
2t en2t + πn2 − 1 .
 πn3 n πn3  πn3  
Paso 5. Cálculo de L {Gn (t)} y Pn (t) . (Método 2)
Gn (s) = L {Gn (t)} =L { 2
πn
(−1)n +1 }= 2
πn
(−1)n +1 1s .

 2 (−1)n +1 1     
 1   1 
Pn (t) = L-1  πn 2
s
=
2
πn
(−1)n +1 L-1  2  = 2
πn
(−1)n +1 L-1  
s + n 2 n4 4 2 n4
   s + n s + 4
− n4  (s + n2 ) − 
4 

 n2  n2
 n +1 -1  n +1 − 2 t 2
Pn (t) = (−1) L 
4 2
 = 4
( −1) e senh(n2 t )
πn3 2 4 3
(s + n2 ) − n4  πn
Paso 6. Cálculo de los wn(t) . (Método 2)
2 2  n2t − n2t 
−n t 2 −n t
wn(t) = 4 2
(−1)n +1 e 2 senh(n2 t ) + 2
(−1)n +1 e−n t = 4
(−1)n +1 e 2  e 2 − e 2  + 2 (−1)n +1 e−n2t
πn3 n πn3  2  n
 
= (−1)n +1 − (−1)n +1 e−n
2t
+ 2n (−1)n +1 e−n
2t
= (−1)n +1 e−n
2t  n2t 2 
e + πn − 1
2 2 2
πn3 πn3 πn3

Paso 7. Construcción de w(x,t).

∞
( −1)n +1
∑ en2t + πn2 − 1 sen(nx) .
2t
w(x, t) = 2
π
e−n  
n3
n =1

Paso 8. Construcción de u(x,t).

Prof. José Luis Quintero 75

 ∞
( −1)n + 1
∑ en2t + πn2 − 1 sen(nx) .
2t
u(x, t) = t(1 − xπ ) + 2
e−n
π
n =1
n3  

Ejemplo 71. Si hay flujo de calor en la superficie de la varilla en estudio hacia el

medio que la rodea cuya temperatura se mantiene constante igual a cero, la ecuación
de calor se convierte en la expresión kuxx − hu = ut . Considere la temperatura u(x,t) en
una varilla delgada con condiciones u(0, t) = u(L, t) = 0 y u(x, 0) = sen( 3Lπx ) . Halle la

solución de la ecuación sujeta a las condiciones.

Solución.
Paso 1. Construcción de G(x,t). F(x, t) = 0 . G(x, t) = 0 .
Paso 2. Cálculo de los Gn (t) . Gn (t) = 0 .
Paso 3. Construcción de H(x). H(x) = sen( 3Lπx ) .
Paso 4. Cálculo de wn (0) .
L L

∫ ∫
(3 − n)πx (3 + n)πx
Wn (0) = 2
L
sen( 3Lπx )sen(nLπx )dx = 1
L
(cos( L
) − cos( L
))dx
0 0
L
= 1
( L
L 3 −n
sen((3 −Ln)πx ) − L
3 +n
sen((3 +Ln)πx ) =0 , n≠3
0
L L L
W3 (0) = 2
L ∫ 0
sen2 ( 3Lπx )dx = 1
L ∫ 0
(1 − cos( 6Lπx ))dx = 1
L
(x − L
6π
sen( 6Lπx ))
0
= 1.

Paso 5. Cálculo de Mn (t) y Mn (0) . (Método 1) Mn (t) = Mn (0) = 0 .

−(k( 3 π )2 + h)t
Paso 6. Cálculo de los wn(t) . (Método 1) w3 (t) = e L .
Paso 5. Cálculo de L {Gn (t)} y Pn (t) . (Método 2) L {Gn (t)} = Pn (t) = 0 .
−(k( 3 π )2 + h)t
Paso 6. Cálculo de los wn(t) . (Método 2) w3 (t) = e L .
−(k( 3π )2 + h)t
Paso 7. Construcción de w(x,t). w(x, t) = e L sen( 3Lπx ) .
−(k( 3 π )2 + h)t
Paso 8. Construcción de u(x,t). u(x, t) = e L sen( 3Lπx ) .

4.27. ECUACIÓN DE ONDA UNIDIMENSIONAL

Hasta el momento, en el estudio de ecuaciones en derivadas parciales, el único

problema físico que se ha introducido es la conducción del calor. Para ampliar la
discusión se estudiará ahora las vibraciones de cuerdas y membranas perfectamente
elásticas. Se comenzará formulando las ecuaciones que rigen las vibraciones de una
cuerda a partir de principios físicos y se demostrará que las condiciones de contorno
adecuadas son similares en sentido matemático a las condiciones de contorno para la
ecuación de calor. Finalmente, se resolverán algunos ejemplos por el método de
separación de variables.

Prof. José Luis Quintero 76

4.28. DEDUCCIÓN DE LA ECUACIÓN DE ONDA PARA LA
CUERDA VIBRANTE

Considere una cuerda de longitud L tal como una cuerda de guitarra tensa ante
dos puntos sobre el eje x, por ejemplo x = 0 y x = L . Cuando comienza a vibrar,
suponga que el movimiento se lleva a cabo en el plano xu, de tal modo que cada punto
de la cuerda se mueve en dirección perpendicular al eje x (vibraciones transversales).
De modo que u(x,t) representa el desplazamiento vertical de cualquier punto de la
cuerda, medido a partir del eje x, cuando t > 0 . Además se supone que:

T2 ( x + ∆ x ,t )
θ1 + ∆θ1 = θ 2
T1 ( x ,t )
θ1

x x x + ∆x

Figura 20. Vibraciones transversales en una cuerda delgada

• La cuerda es perfectamente flexible.

• La cuerda es homogénea, esto es, ρ , su masa por unidad de longitud es
constante.
• Los desplazamientos u son pequeños comparados con el largo de la cuerda.
• La pendiente de la curva es pequeña en todos sus puntos.
• La tensión T actúa tangente a la cuerda y su magnitud T es la misma en todos sus
puntos.
• La tensión es grande en comparación con la fuerza de gravedad.
• No hay otras fuerzas externas actuando sobre la cuerda.

En la figura 20, las tensiones T1 y T2 son tangentes a los extremos de la curva

en el intervalo x, x + ∆x  . Para θ1 y θ2 pequeños, la fuerza vertical neta que actúa
sobre el elemento ∆s correspondiente a la cuerda es, por consiguiente,
Tsenθ2 − Tsenθ1 = Ttgθ2 − Ttgθ1 = T ux (x + ∆x, t) − ux (x, t) ,

en donde T = T1 = T2 . Ahora, ρ∆s ≈ ρ∆x es la masa de la cuerda en x, x + ∆x  y al

aplicar la segunda ley de Newton se obtiene T ux (x + ∆x, t) − ux (x, t) = ρ∆xutt o sea
ux (x + ∆x, t) − ux (x, t) ρ
= utt .
∆x T

Prof. José Luis Quintero 77

 Si se toma el límite cuando ∆x → 0 , esta última ecuación se transforma en
uxx = (ρ / T)utt . Si se toma a2 = ρ / T , se obtiene uxx = a2utt .

4.29. PROBLEMAS ILUSTRATIVOS DE LA ECUACIÓN DE ONDA

UNIDIMENSIONAL

Ejemplo 72. Problema típico de condición en la frontera es

∂2u ∂2u
a2 2 = 2 , 0 < x < L, a > 0, t > 0 , (52)
∂x ∂t
u(0, t) = 0, u(L, t) = 0, t ≥ 0, (53)
∂u
u(x, 0) = f(x), = g(x), 0 < x <L. (54)
∂t t = 0
Solución.
Las condiciones de frontera (53) simplemente dicen que los extremos de la cuerda
permanecen fijos en todo instante. En t = 0 , las funciones f y g dadas en (54)
especifican la configuración inicial y la velocidad inicial de cada punto de la cuerda,
respectivamente. En el contexto está implícito que f es continua y que f(0) = 0 y
f(L) = 0 . Sea u(x, t) = X(x)T(t) . La ecuación se convierte en a2 X '' T = T '' X . De modo
que
a2 X '' T T '' X X '' T ''
2
= 2
, XT ≠ 0 , = 2 = −λ2 .
a XT a TX X a T
Se tienen
X ''+ λ2 X = 0 ⇒ X = A cos(λx) + Bsen(λx) , T ''+ a2 λ2 T = 0 ⇒ T = C cos(aλt) + Dsen(aλt) .
Condiciones de borde: u(0, t) = 0 ⇒ X(0)T(t) = 0 ⇒ X(0) = 0.
Se tiene: X(0) = A cos(0) + Bsen(0) = 0 ⇒ A = 0 , u(L, t) = 0 ⇒ X(L)T(t) = 0 ⇒ X(L) = 0.
nπ
Entonces: X(L) = Bsen(λL) = 0 ⇒ sen(λL) = 0 ⇒ λ = L
, n = 1, 2, 3,... . Por lo tanto:
  anπt   anπt    nπx 
un (x, t) =  Cn cos   + Dnsen    Bnsen  .
  L   L   L 
∞
  anπt   anπt    nπx 
u(x, t) = ∑  E
n =1
n cos 
 L 
 + Fnsen 
 L 
  sen  L  .
 
Usando condiciones iniciales:
∞
nπx  L
 nπx 
∑ E sen 
2
u(x, 0) = f(x) =
n =1
n
L 
⇒ En =
L ∫ 0
f(x)sen 
 L 
 dx

∞
 anπ  anπt  anπ  anπt    nπx 
ut (x, t) = ∑  −
n =1
L
Ensen 
 L  + L Fn cos  L   sen  L 
    
∞
anπ  nπx  anπ L
 nπx 
∑
2
ut (x, 0) = g(x) =
n =1
L
Fnsen 
 L 
 ⇒
L
Fn =
L ∫ 0
g(x)sen 
 L 
 dx

∞
anπ  nπx  2 L
 nπx 
ut (x, 0) = g(x) = ∑n =1
L
Fnsen 
 L  ⇒ Fn = anπ
 ∫ 0
g(x)sen 
 L 
 dx

Prof. José Luis Quintero 78

Observación 9. Si la cuerda se suelta en reposo (g(x) = 0) , entonces Fn = 0 y la
solución se expresa como
∞
 anπt   nπx 
u(x, t) = ∑E
n =1
n cos 
 L 
 sen 
 L 
.

Usando la identidad trigonométrica

sen(A + B) + sen(A − B)
sen(A) cos(B) = ,
2
la solución obtenida es
∞
1   nπ   nπ 
u(x, t) =
2 ∑E
n =1
n sen 
  L
(x + at)  + sen 
  L
(x − at)  
 
∞ ∞
1  nπ   nπ 
= 
2  n =1

∑
Ensen 
 L
(x + at)  +
 n =1
Ensen 
 L
∑
(x − at)  

1
= f(x + at) + f(x − at) (Ecuación de d' Alembert)
2

∞
 anπt   nπx 
Observación 10. u(x, t) = ∑  G sen 
n =1
n
L
+ φn   sen 
  L 
2 2
 donde Gn = En + Fn y

En Fn
sen(φn ) = , cos(φn ) = .
Gn Gn

Ejemplo 73.
4uxx = utt , 0 < x < 2 , t > 0

u(0, t) = 0 , u(2, t) = 0, t ≥ 0

ut (x, 0) = 0 , u(x, 0) = x , 0 < x < 2
Solución. Sea u(x, t) = X(x).T(t) . Por tanto: uxx = X '' T y utt = XT '' . De modo que:
4X '' T = XT '' .
Al dividir por 4XT ≠ 0 se obtiene
X '' T ''
= .
X 4T
Igualando a una constante de separación:
X '' T ''
= = K,
X 4T
se tienen el par de ecuaciones ordinarias: X ''− KX = 0 , T ''− 4KT = 0 . Se estudian 3
casos:
Caso I: K > 0 , K = λ2 . Ecuaciones ordinarias:
X ''− λ2 X = 0 ⇒ X(x) = Aeλx + Be−λx , T ''− 4λ2 T = 0 ⇒ T(t) = Ce2λt + De−2λt
Condiciones de borde homogéneas:
u(0, t) = 0 ⇒ X(0) = 0 ⇒ A + B = 0 ⇒ A = −B
u(2, t) = 0 ⇒ X(2) = 0 ⇒ Ae2λ + Be−2λ = A(e2λ − e−2λ ) = 0 ⇒ A = 0 ⇒ B = 0
Solo se obtiene solución trivial.
Caso II: K = 0 . Ecuaciones ordinarias:
X '' = 0 ⇒ X(x) = A + Bx , T '' = 0 ⇒ T(t) = C + Dt .
Condiciones de borde homogéneas:

Prof. José Luis Quintero 79

 u(0, t) = 0 ⇒ X(0) = 0 ⇒ A = 0 , u(2, t) = 0 ⇒ X(2) = 0 ⇒ 2B = 0 ⇒ B = 0 .
Solo se obtiene solución trivial.
Caso III: K < 0 , K = −λ2 . Ecuaciones ordinarias:
X ''+ λ2 X = 0 ⇒ X(x) = A cos(λx) + Bsen(λx) , T ''+ 4λ2 T = 0 ⇒ T(t) = C cos(2λt) + Dsen(2λt)
Condiciones de borde homogéneas:
u(0, t) = 0 ⇒ X(0) = 0 ⇒ A = 0
nπ
u(2, t) = 0 ⇒ X(2) = 0 ⇒ Bnsen(2λ) = 0 ⇒ λ = 2
, n = 1, 2,3,...
T '(t) = 2λ(−Csen(2λt) + D cos(2λt)) ut (x, 0) = 0 ⇒ T '(0) = 0 ⇒ D = 0.
Se tiene:
∞ ∞
u(x, t) = ∑
n =1
Bnsen(n2πx )Cn cos(nπt) ⇒ u(x, t) = ∑ E sen(
n =1
n
nπx
2
) cos(nπt)

Condición de borde no homogénea:

∞
u(x, 0) = f(x) = ∑ E sen(
n =1
n
nπx
2
).

En consecuencia:
2 2
En =
∫ 0
xsen ( n2πx ) dx = n2π n2π sen ( n2πx ) − x cos ( n2πx )
0
2 4 2 4 4
n impar = 2 = . n par = −2 = − . Por tan to : En = (−1)n +1 .
nπ   nπ nπ   nπ nπ
Se obtiene:
∞
(−1)n +1
∑
4
u(x, t) = cos(nπt)sen(n2πx ) = x .
π n =1
n

Ejemplo 74. Problema típico de condición en la frontera es

∂2u ∂2u
a2 2 = 2 , 0 < x < L, a > 0, t > 0 , (55)
∂x ∂t
ux (0, t) = 0, ux (L, t) = 0, t ≥ 0, (56)
∂u
u(x, 0) = f(x), = g(x), 0 < x <L. (57)
∂t t = 0
Solución.
En t = 0 , las funciones f y g dadas en (57) especifican la configuración inicial y la
velocidad inicial de cada punto de la cuerda, respectivamente. En el contexto está
implícito que f es continua y que f(0) = 0 y f(L) = 0 . Sea u(x, t) = X(x)T(t) . La ecuación
se convierte en a2 X '' T = T '' X .
a2 X '' T T '' X X '' T ''
De modo que 2
= 2
, XT ≠ 0 , = 2 = −λ2 . Se tienen
a XT a TX X a T
X ''+ λ X = 0 ⇒ X = A cos(λx) + Bsen(λx) , T ''+ a2 λ2 T = 0 ⇒ T = C cos(aλt) + Dsen(aλt) .
2

Condiciones de borde:
ux (0, t) = 0 ⇒ X '(0)T(t) = 0 ⇒ X '(0) = 0 , ux (L, t) = 0 ⇒ X '(L)T(t) = 0 ⇒ X '(L) = 0.
  anπt   anπt    nπx 
Por lo tanto: un (x, t) =  Cn cos   + Dnsen  L   Bn cos  L  , n = 0,1,...
  L     

Prof. José Luis Quintero 80

 ∞
  anπt   anπt    nπx 
u(x, t) = E0 + ∑  E
n =1
n cos 
 L  + Fnsen  L   cos  L  .
    
Usando condiciones iniciales:
∞ ∞
 nπx  2E0  nπx 
u(x, 0) = f(x) = E0 + ∑
n =1
En cos 
 L 
 ⇒
2
+ ∑E
n =1
n cos 
 L 


1 L 2 L
 nπx 
donde E0 =
L ∫ 0
f(x)dx , En =
L ∫ 0
f(x) cos 
 L 
 dx
∞
 anπ  anπt  anπ  anπt    nπx 
ut (x, t) = ∑  −
n =1
L
Ensen 
 L 

+
L
Fn cos 
 L   cos  L 
   
∞
anπ  nπx  anπ L
 nπx 
∑
2
ut (x, 0) = g(x) =
n =1
L
Fn cos 
 L 
 ⇒
L
Fn =
L ∫ 0
g(x) cos 
 L 
 dx

∞
anπ  nπx  2 L
 nπx 
ut (x, 0) = g(x) = ∑
n =1
L
Fnsen 
 L  ⇒ Fn = anπ
 ∫ 0
g(x) cos 
 L 
 dx

Ejemplo 75.
a2uxx + f(x, t) = utt , a > 0, 0 < x < L, t > 0 ,
u(0, t) = α1(t), u(L, t) = α2 (t), t > 0,
u(x, 0) = ϕ(x), ut (x, 0) = β(x), 0 < x < L.
Solución.
En primer lugar, se puede escribir u(x, t) = F(x, t) + w(x, t) , donde F(x,t) es una función
que debe resolver el problema Fxx (x, t) = 0, F(0, t) = α1(t), F(L, t) = α2 (t) . De modo que
Fxx (x, t) = 0 ⇒ F(x, t) = g1(t) + xg2 (t) . Al imponer la primera condición de borde se tiene
F(0, t) = α1(t) = g1(t) . Al imponer la segunda condición de borde se tiene
α2 (t) −α1(t)
F(x, t) = α1(t) + xg2 (t) ⇒ F(L, t) = α2 (t) = α1(t) + Lg2 (t) ⇒ g2 (t) = L
.
Por lo tanto se obtiene la función definida por F(x, t) = α1(t) + x
L
α2 (t) − α1(t) . Por otro
lado w(x,t) es una función incógnita auxiliar cuya determinación se estudiará a
continuación.

Como F(0, t) = α1(t) y F(L, t) = α2 (t) por construcción, resulta que w(0, t) = 0 y
w(L, t) = 0 . Por otro lado se tiene que
w(x, 0) = u(x, 0) − F(x, 0) = ϕ(x) − F(x, 0) , wt (x, 0) = ut (x, 0) − Ft (x, 0) = β(x) − Ft (x, 0)
en 0 < x < L , y a2wxx + f(x, t) = wtt + Ftt (x, t) ⇒ a2wxx = wtt + Ftt (x, t) − f(x, t) .
Si se denotan
H(x) = ϕ(x) − F(x, 0) , I(x) = β(x) − Ft (x, 0) , G(x, t) = f(x, t) − Ftt (x, t) ,
el problema original queda convertido como
a2wxx + G(x, t) = wtt , a > 0, 0 < x < L, t > 0 ,
w(0, t) = 0, w(L, t) = 0, t > 0 ,
w(x, 0) = H(x), wt (x, 0) = I(x), 0 < x < L .

Este problema representa el caso de una ecuación de onda no homogénea con

condiciones iniciales no nulas y condiciones de frontera nulas.

Prof. José Luis Quintero 81

Si se ensaya con una solución de la forma
∞
w(x, t) = ∑ w (t)sen(
n =1
n
nπx
L
),

y se expresa
∞
G(x, t) = ∑ G (t)sen(
n =1
n
nπx
L
),

con
2 L
Gn (t) =
L ∫ 0
G(x, t)sen(nLπx )dx, n = 1,... ,

entonces
∞ ∞
wxx (x, t) = − ∑
n =1
wn (t)(nLπ )2 sen(nLπx ) , wtt (x, t) = ∑ w (t)sen(
n =1
''
n
nπx
L
).

Sustituyendo en la ecuación a2wxx + G(x, t) = wtt , se tiene

∞ ∞ ∞
−a2 ∑
n =1
wn(t)(nLπ )2 sen(nLπx ) + ∑
n =1
Gn (t)sen( nLπx ) = ∑ w (t)sen(
n =1
''
n
nπx
L
).

Reacomodando términos:
∞

∑ w (t) + (
n =1
''
n
anπ 2
L
) wn (t) − Gn(t)sen(nLπx ) = 0 ⇒ wn'' (t) + ( anL π )2 wn (t) = Gn (t) .
 
Ecuación lineal de segundo orden

Como
∞
w(x, t) = ∑ w (t)sen(
n =1
n
nπx
L
),

entonces de la primera condición inicial del problema se tiene

∞
H(x) = ∑ w (0)sen(
n =1
n
nπx
L
)

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente wn (0) como
2 L
wn(0) =
L ∫ 0
H(x)sen(nLπx )dx .

De la segunda condición inicial del problema se tiene

∞
I(x) = ∑ w (0)sen(
n =1
'
n
nπx
L
)

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente wn' (0) como
2 L
wn' (0) =
L ∫ 0
I(x)sen(nLπx )dx .

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de
segundo orden sujeta a dos condiciones:
2 L 2 L
wn'' (t) + ( anL π )2 wn (t) = Gn(t) , wn (0) =
L 0 ∫
H(x)sen(nLπx )dx , wn' (0) =
L 0
I(x)sen(nLπx )dx
∫

Prof. José Luis Quintero 82

Método 1.
Solución general:
sen( anπt ) G (t) cos( anπt )dt − cos( anπt ) G (t)sen( anπt )dt 
wn(t) = C1 cos( anLπt ) + C2sen( anLπt ) + L
anπ  L ∫
n L L n L ∫
wn(t) = C1 cos( anLπt ) + C2sen( anLπt ) + L
anπ
Mn(t)sen( anLπt )
 − Nn (t) cos( anLπt ) .

Nn (0)L L  ' [Link]' (0) 

C1 = wn (0) + . C2 =  wn (0) + − Mn (0)  .
anπ 
anπ  anπ 

Resumen de pasos a seguir:
1. Construcción de G(x,t). G(x, t) = f(x, t) − Ftt (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x) y I(x). H(x) = ϕ(x) − F(x, 0) I(x) = β(x) − Ft (x, 0)
4. Cálculo de wn (0) y wn' (0) .
5. Cálculo de Mn (t), Mn (0), Nn (t), Nn (0) y Nn' (0) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

Método 2.
L{wn'' (t)} + ( anL π )2 L{wn (t)} = L{Gn (t)} ⇒ s2wn (s) − swn (0) − wn' (0) + ( anL π )2 wn(s) = Gn(s) .
Gn (s) + swn (0) + wn' (0)
s2 + ( anπ )2  wn (s) = Gn (s) + swn (0) + wn' (0) ⇒ wn (s) =
 L  s2 + ( anLπ )2
Gn(s) swn (0) wn' (0)
= + + .
s2 + ( anLπ )2 s2 + ( anL π )2 s2 + ( anL π )2
 G (s) 
anπt
wn(t) = L-1{wn (s)} = L-1  2 n  + wn(0) cos( L ) +
L
anπ
wn' (0)sen( anLπt )
anπ 2
 s + ( L
) 
= Pn (t) + wn (0) cos( anLπt ) + L
anπ
wn' (0)sen( anLπt ).

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = f(x, t) − Ftt (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x) y I(x). H(x) = ϕ(x) − F(x, 0) I(x) = β(x) − Ft (x, 0)
4. Cálculo de wn (0) y wn' (0) .
5. Cálculo de L {Gn (t)} y Pn (t) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

∞
Solución final: u(x, t) = F(x, t) + w(x, t) = α1(t) + x
L
α2 (t) − α1(t) + ∑ w (t)sen(
n =1
n
nπx
L
).

Prof. José Luis Quintero 83

Ejemplo 76. Resuelva la ecuación diferencial
a2uxx + 2ut − u = utt , a > 0 , 0 < x < L , t > 0

 u(0, t) = α1(t) , u(L, t) = α2 (t) , t > 0 .
 u(x, 0) = ϕ(x) , u (x, 0) = β(x) , 0 < x < L
 t

Solución.
En primer lugar, se puede escribir u(x, t) = F(x, t) + w(x, t) , donde F(x,t) es una función
que debe resolver la ecuación con condiciones
Fxx (x, t) = 0, F(0, t) = α1(t), F(L, t) = α2 (t) .
De modo que Fxx (x, t) = 0 ⇒ F(x, t) = g1(t) + xg2 (t) . Al imponer la primera condición de
borde se tiene F(0, t) = α1(t) = g1(t) . Al imponer la segunda condición de borde se tiene
α2 (t) − α1 (t)
F(x, t) = α1(t) + xg2 (t) ⇒ F(L, t) = α2 (t) = α1(t) + Lg2 (t) ⇒ g2 (t) = L
. Por lo tanto
F(x, t) = α1(t) +
α2 (t) − α1(t) . Por otro lado w(x,t) es una función auxiliar cuya
x
L
determinación se estudiará a continuación. Como F(0, t) = α1(t) y F(L, t) = α2 (t) por
construcción, resulta que w(0, t) = 0 y w(L, t) = 0 . Por otro lado se tiene que
w(x, 0) = u(x, 0) − F(x, 0) = ϕ(x) − F(x, 0) , wt (x, 0) = ut (x, 0) − Ft (x, 0) = β(x) − Ft (x, 0) en
0 < x < L , y a wxx + 2wt + 2Ft − w − F = wtt + Ftt ⇒ a2wxx + 2wt − w = wtt + Ftt − 2Ft + F  .
2

Si se denotan
H(x) = ϕ(x) − F(x, 0) , I(x) = β(x) − Ft (x, 0) , G(x, t) = − Ftt (x, t) − 2Ft (x, t) + F(x, t) ,
el problema original queda convertido en
a2wxx + 2wt − w + G(x, t) = wtt , a > 0 , 0 < x < L , t > 0

 w(0, t) = 0 , w(L, t) = 0 , t > 0 .
 w(x, 0) = H(x) , wt (x, 0) = I(x) , 0 < x < L


Este problema representa el caso de una ecuación de onda no homogénea con

condiciones iniciales no nulas y condiciones de frontera nulas. Si se ensaya con una
∞ ∞
solución de la forma w(x, t) = ∑
n =1
wn (t)sen(nLπx ) , y G(x, t) = ∑ G (t)sen(
n =1
n
nπx
L
), con

2 L
Gn (t) =
L ∫ 0
G(x, t)sen(nLπx )dx, n = 1,... ,

entonces
∞ ∞ ∞
wxx (x, t) = − ∑
n =1
wn (t)(nLπ )2 sen(nLπx ) , wt (x, t) = ∑
n =1
wn' (t)sen(nLπx ) , wtt (x, t) = ∑ w (t)sen(
n =1
''
n
nπx
L
).

Sustituyendo en la ecuación a2wxx + 2wt − w + G(x, t) = wtt , se tiene

∞ ∞ ∞ ∞
−a2 ∑
n =1
wn(t)(nLπ )2 sen(nLπx ) + 2 ∑
n =1
wn' (t)sen(nLπx ) − ∑
n =1
wn (t)sen(nLπx ) + ∑ G (t)sen(
n =1
n
nπx
L
)

∞
= ∑ w (t)sen(
n =1
''
n
nπx
L
)

Reacomodando términos:

Prof. José Luis Quintero 84

 ∞

∑ w (t) − 2w (t) + ((

n =1
''
n
'
n
anπ 2
L
) + 1)wn (t) − Gn (t)sen(nLπx ) = 0


⇒ wn'' (t) − 2wn' (t) + (( anL π )2 + 1)wn (t) = Gn(t)


Ecuación lineal de segundo orden

Como
∞
w(x, t) = ∑ w (t)sen(
n =1
n
nπx
L
),

entonces de la primera condición inicial del problema se tiene

∞
H(x) = ∑ w (0)sen(
n =1
n
nπx
L
)

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente wn (0) como
2 L
Wn (0) =
L ∫ 0
H(x)sen(nLπx )dx .

De la segunda condición inicial del problema se tiene

∞
I(x) = ∑ w (0)sen(
n =1
'
n
nπx
L
)

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente wn' (0) como
2 L
Wn' (0) =
L ∫ 0
I(x)sen(nLπx )dx .

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de
segundo orden sujeta a dos condiciones:
wn'' (t) − 2wn' (t) + (( anL π )2 + 1)wn (t) = Gn(t) ,
2 L 2 L
Wn (0) =
L ∫ 0
H(x)sen(nLπx )dx , Wn' (0) =
L ∫ 0
I(x)sen(nLπx )dx

Método 1. Solución general:

wn(t) = C1et cos( anLπt ) + C2et sen( anLπt ) +
sen( anπt ) G (t)e−2t cos( anπt )dt − cos( anπt ) G (t)e−2t sen( anπt )dt 
Let
anπ  L n ∫ L L n L  ∫
Let
wn(t) = C1e cos( anLπt ) + C2 e sen( anLπt ) +
t t
anπ
sen( anπt )Mn (t) − cos( anπt )Nn(t)
 L L 
Nn (0)L L  ' L(2Nn (0) + Nn' (0)) 
C1 = wn (0) + . C2 =  wn(0) − wn (0) − − Mn (0)  .
anπ anπ  anπ 

Resumen de pasos a seguir:
1. Construcción de G(x,t). G(x, t) = f(x, t) − Ftt (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x) y I(x). H(x) = ϕ(x) − F(x, 0) I(x) = β(x) − Ft (x, 0)
4. Cálculo de wn (0) y wn' (0) .
5. Cálculo de Mn (t), Mn (0), Nn (t), Nn (0) y Nn' (0) .
6. Cálculo de los wn(t) .

Prof. José Luis Quintero 85

7. Construcción de w(x,t).
8. Construcción de u(x,t).

Método 2.
L{wn'' (t)} − 2L{wn' (t)} + (( anL π )2 + 1)L{wn (t)} = L{Gn (t)}
⇒ s2wn(s) − swn(0) − wn' (0) − 2swn(s) + 2wn (0) + (( anL π )2 + 1)wn (s) = Gn (s)
s2 − 2s + 1 + ( anπ )2  wn (s) = Gn (s) + (s − 2)wn(0) + wn' (0)
 L 
Gn (s) + (s − 2)wn(0) + wn' (0) Gn (s) (s − 1)wn (0) wn' (0) − wn (0)
⇒ wn(s) = = + + .
2
(s − 1) + ( anL π )2 2
(s − 1) + ( anL π )2 2
(s − 1) + ( anL π )2 (s − 1)2 + ( anL π )2
 Gn(s) 
anπt
wn(t) = L-1{wn (s)} = L-1  t
 + wn (0)e cos( L ) +
L
anπ
(wn' (0) − wn (0))et sen( anLπt )
2 anπ 2
 (s − 1) + ( L ) 
= Pn (t) + wn (0)et cos( anLπt ) + L
anπ
(wn' (0) − wn (0))et sen( anLπt )

Resumen de pasos a seguir:

1. Construcción de G(x,t). G(x, t) = f(x, t) − Ftt (x, t)
2. Cálculo de los Gn (t) .
3. Construcción de H(x) y I(x). H(x) = ϕ(x) − F(x, 0) I(x) = β(x) − Ft (x, 0)
4. Cálculo de wn (0) y wn' (0) .
5. Cálculo de L {Gn (t)} y Pn (t) .
6. Cálculo de los wn(t) .
7. Construcción de w(x,t).
8. Construcción de u(x,t).

Solución final:
∞
u(x, t) = F(x, t) + w(x, t) = α1(t) + x
L
α2 (t) − α1(t) + ∑ w (t)sen(
n =1
n
nπx
L
).

4.30. ECUACIÓN DE LAPLACE

A continuación, se aplicarán las técnicas aprendidas anteriormente a un tipo

diferente de problemas. Se considera el problema estacionario del calor en un dominio
bidimensional, con condiciones no homogéneas todas. En ese caso, la ecuación en
derivadas parciales es lineal y homogénea, pero las condiciones de contorno, aunque
lineales, son todas no homogéneas. Si todas las condiciones de contorno no son
homogéneas, no se podrá aplicar el método de separación de variables a este
problema en su forma actual. Se podrá superar esta dificultad mediante la
descomposición en problemas donde el método de separación de variables, sea
apropiado.

Prof. José Luis Quintero 86

4.31. PROBLEMAS DONDE INTERVIENE LA ECUACIÓN DE
LAPLACE EN COORDENADAS CARTESIANAS

Aunque no se presentará su deducción, esta ecuación en dos o tres dimensiones

surge en problemas independientes del tiempo que conciernen a potenciales como el
electrostático, el gravitacional y la velocidad en mecánica de fluidos. Además, una
solución de la ecuación de Laplace se puede interpretar como la distribución de
temperatura en un estado estable. De manera que, una solución u(x,y) podría
representar la temperatura que varía de un punto a otro, pero no con el tiempo, en
una placa rectangular. La ecuación de Laplace, en dos y en tres dimensiones, se
abrevia en la forma ∇2u = 0 , donde
∂2u ∂2u ∂2u ∂2u ∂2u
∇2u = 2
+ 2
y ∇2u = 2
+ 2
+
∂x ∂y ∂x ∂y ∂z2
se llaman, respectivamente, laplaciano en dos dimensiones y laplaciano en tres
dimensiones, de una función u.

Ejemplo 77. Problema de Dirichlet para un rectángulo.

∂2u ∂2u
+ = 0, 0 < x < a, 0 < y < b
∂x2 ∂y2
u(x, 0) = 0, u(x,b) = 0, 0 < x < a
u(0, y) = 0, u(a, y) = f(y), 0 ≤ y ≤ b .
Solución.
Sea u(x, y) = X(x)Y(y) . La ecuación se convierte en X '' Y = − Y '' X . De modo que
X '' Y Y '' X X '' Y ''
=− , XY ≠ 0 , =− = λ2 .
XY XY X Y
Se tienen
X ''− λ2 X = 0 ⇒ X = Aeλx + Be−λx , Y ''+ λ2 Y = 0 ⇒ Y = C cos(λy) + Dsen(λy) .
Primeras condiciones de borde homogéneas:
u(x, 0) = 0 ⇒ X(x)Y(0) = 0 ⇒ Y(0) = 0. u(x,b) = 0 ⇒ X(x)Y(b) = 0 ⇒ Y(b) = 0.
Se tiene: Y(0) = C cos(0) + Dsen(0) = 0 ⇒ C = 0
nπ
Se tiene: Y(b) = Dsen(λb) = 0 ⇒ sen(λb) = 0 ⇒ λ = b
, n = 1, 2, 3,... .

Por lo tanto:
 nπx − nπx   nπy 
un (x, y) =  Ane b + Bne b  Dnsen  .
   b 
Otra condición de borde homogénea: u(0, y) = 0 ⇒ X(0)Y(y) = 0 ⇒ X(0) = 0.
Se tiene: X(0) = Ae0 + Be0 = 0 ⇒ A = −B . Por lo tanto:
 nπx − nπx   nπx − nπx 
 nπy   nπy 
un (x, y) =  Ane b − Ane b  Dnsen   ⇒ En  e b − e b  sen  
   b     b 
 nπx   nπy   nπx   nπy 
⇒ 2Ensenh   sen  b  ⇒ Fnsenh  b  sen  b 
 b       
∞
nπx   nπy 
u(x, y) = ∑ F senh 
n =1
n
b 
sen 
 b 


Condición de borde no homogénea:

Prof. José Luis Quintero 87

 ∞
anπ   nπy   anπ  2 b
 nπy 
U(a, y) = f(y) = ∑ F senh 
n =1
n
b 
sen 
 b  ⇒ Fnsenh  b  = b
   ∫ 0
f(y)sen 
 b 
 dy

b
 nπy 
⇒ Fn = 2
 anπ 
[Link] 
 b 

∫ 0
f(y)sen 
 b 
 dy

Ejemplo 78. Placa rectangular con caras laterales aisladas.

∂2u ∂2u
+ = 0, 0 < x < a, 0 < y < b
∂x2 ∂y2
u(x, 0) = 0, u(x,b) = f(x), 0 < x < a
ux (0, y) = 0, ux (a, y) = 0, 0 ≤ y ≤ b .
Solución.
Sea u(x, y) = X(x)Y(y) . La ecuación se convierte en X '' Y = − Y '' X . De modo que
X '' Y Y '' X X '' Y ''
=− , XY ≠ 0 , =− = −λ2 .
XY XY X Y
Se tienen
X ''+ λ2 X = 0 ⇒ X = A cos(λx) + Bsen(λx) , Y ''− λ2 Y = 0 ⇒ Y = Ceλy + De−λy .
Primeras condiciones de borde homogéneas:
ux (0, y) = 0 ⇒ X '(0)Y(y) = 0 ⇒ X '(0) = 0. ux (a, y) = 0 ⇒ X '(a)Y(y) = 0 ⇒ X '(a) = 0.
Se tiene: X '(0) = − Aλsen(0) + Bλ cos(0) = 0 ⇒ Bλ = 0 ⇒ B = 0
nπ
Se tiene: X '(a) = −Aλsen(λa) = 0 ⇒ λ = 0 ó λ = a
.
Si λ = 0 , u0 (x, y) = (A0 + B0 x)(C0 + D0 y)
uox (0, y) = 0 ⇒ X '(0)Y(y) = 0 ⇒ X '(0) = 0 ⇒ B = 0
uox (a, y) = 0 ⇒ X '(a)Y(y) = 0 ⇒ X '(a) = 0 ⇒ 0 = 0
Otra condición de borde homogénea: u0 (x, 0) = 0 ⇒ X(x)Y(0) = 0 ⇒ Y(0) = 0 ⇒ C = 0
nπ
Por lo tanto: u0 (x, y) = E0 y . Si λ = a
,
 nπy nπy 
−  nπx 
un (x, y) =  Cne a + Dne a  An cos  
   a 
 
Otra condición de borde homogénea: un (x, 0) = 0 ⇒ X(x)Y(0) = 0 ⇒ Y(0) = 0 ⇒ Cn = −Dn
Por lo tanto:
 nπy −
nπy 
 nπx   nπy   nπx 
un (x, y) =  e a − e a  En cos   ⇒ senh   Fn cos  
   a   a   a 
 
∞
nπy   nπx 
∑ F senh 
F0 y
u(x, y) = + n  cos  
2 n =1
a   a 
Condición de borde no homogénea:
∞
bnπ   nπx 
∑ F senh 
F0b
u(x,b) = f(x) = + cos   donde
a 
n
2 n =1  a 
2 a 2 a
F0b =
a 0 ∫ f(x)dx ⇒ F0 =
ab 0
f(x)dx
∫
 bnπ  2 a  nπx  a
 nπx 
Fnsenh 
 a = a
 0 ∫
f(x) cos 
 a 
 dx ⇒ Fn =
2
 bnπ 
[Link]
 a 

∫ 0
f(x) cos 
 a 
 dx

Prof. José Luis Quintero 88

Ejemplo 79.
∂2u ∂2u
+ = 0, 0 < x < a, 0 < y < b
∂x2 ∂y2
u(x, 0) = f(x), u(x,b) = g(x), 0 < x < a
u(0, y) = F(y), u(a, y) = G(y), 0 < y < b .
Solución.
PRINCIPIO DE SUPERPOSICIÓN. Un problema de Dirichlet para un rectángulo se
puede resolver con facilidad separando variables cuando se especifican condiciones
homogéneas para dos fronteras paralelas, sin embargo, el método de separación de
variables no se aplica a un problema de Dirichlet cuando las condiciones en la frontera
en los cuatro lados del rectángulo son no homogéneas. Para salvar esta dificultad se
descomponen el problema en dos problemas llamados PROBLEMA 1 y PROBLEMA 2.
Estos problemas vienen dados como

PROBLEMA 1
2 2
∂ u ∂ u
2
+ = 0, 0 < x < a, 0 < y < b
∂x ∂y2
u(x, 0) = f(x), u(x,b) = g(x), 0 < x < a
u(0, y) = 0, u(a, y) = 0, 0 < y < b .

PROBLEMA 2
∂2u ∂2u
+ = 0, 0 < x < a, 0 < y < b
∂x2 ∂y2
u(x, 0) = 0, u(x,b) = 0, 0 < x < a
u(0, y) = F(y), u(a, y) = G(y), 0 < y < b .

De modo que:
SOLUCIÓN PROBLEMA = SOLUCIÓN PROBLEMA 1 + SOLUCIÓN PROBLEMA 2.

Ejemplo 80. Ecuación de Poisson.

uxx + uyy = G(x, y), 0 < x < a, 0 < y < b ,
u(x, 0) = 0, u(x,b) = 0, 0 < x < a,
u(0, y) = f(y), u(a, y) = g(y), 0 < y < b.
Solución.
Si se ensaya con una solución de la forma
∞
u(x, y) = ∑ u (x)sen(
n =1
n
nπy
b
),

y se expresa
∞ L
∑ G (x)sen(
2
∫
nπy nπy
G(x, y) = n b
), con Gn (x) = G(x, y)sen( b
)dy, n = 1,... ,
n =1
L 0

se tiene que
∞ ∞
uxx (x, y) = ∑
n =1
un'' (x)sen(nbπy ) y uyy (x, y) = − ∑ u (x)(
n =1
n L
) sen(nbπy ) .
nπ 2

Prof. José Luis Quintero 89

Sustituyendo en la ecuación uxx + uyy = G(x, y) se tiene
∞ ∞ ∞

∑
n =1
un'' (x)sen(nbπy ) − ∑
n =1
un (x)(nbπ )2 sen(nbπy ) = ∑ G (x)sen(
n =1
n
nπy
b
).

Reacomodando términos:
∞

∑ u (x) − (
n =1
''
n
nπ 2
b
) un(x)

nπy
− Gn (x)sen( b ) = 0 ⇒ un'' (x) − (nbπ )2 un(x) = Gn (x) .

Ecuación lineal de segundo orden

Como
∞
u(x, y) = ∑ u (x)sen(
n =1
n
nπy
b
),

entonces de la condición de borde del problema se tiene

∞
f(y) = ∑ u (0)sen(
n =1
n
nπy
b
)

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente un(0) como
2 b

∫
nπy
un (0) = f(y)sen( b
)dy .
b 0

De la segunda condición de borde del problema se tiene

∞
g(y) = ∑ u (a)sen(
n =1
n
nπy
b
)

(serie de senos de Fourier)

y por tanto se obtiene el coeficiente un (a) como
2 b

∫
nπy
un (a) = g(y)sen( b
)dy .
b 0

De manera que hay que resolver una ecuación lineal ordinaria no homogénea de
segundo orden sujeta a dos condiciones:
2 b 2 b

∫ ∫
nπy nπy
un'' (x) + (nbπ )2 un(x) = Gn (x) , un (0) = f(y)sen( b )dy , un (a) = g(y)sen( b
)dy
b 0 b 0

Método 1.
Solución general:
un (x) = C1 cosh(nbπx ) + C2senh(nbπx ) +
senh(nπx ) G (x) cosh(nπx )dx − cosh(nπx ) G (x)senh( nπx )dx 
b
nπ  b n ∫ b b n b  ∫
un (x) = C1 cosh(nbπx ) + C2senh(nbπx ) + b
nπ
Mn (x)senh(nbπx ) − Nn(x) cosh(nbπx ) .
 

Nn (0)b
C1 = un (0) + .
nπ
C2 = 1 un (a) − C1 cosh( nπa) − b Mn (a)senh(nπa) − Nn (a) cosh(nπa)  .
senh( nπa)  b nπ  b b 
b

Prof. José Luis Quintero 90

Resumen de pasos a seguir:
1. Cálculo de los Gn (x) .
2. Cálculo de un (0) y un (a) .
3. Cálculo de Mn (x), Mn (a), Nn (x), Nn (0) y Nn (a) .
4. Cálculo de los un (x) .
5. Construcción de u(x,y).

Método 2.
L{un'' (x)} − (nbπ )2 L{un (x)} = L{Gn (x)} ⇒ s2un(s) − sun (0) − un' (0) − (nbπ )2 un (s) = Gn (s) .
Gn (s) + sun(0) + un (a)
s2 − (nπ )2  un(s) = Gn(s) + sun (0) + un' (0) ⇒ un (s) =
 b  s2 − (nbπ )2
Gn (s) sun (0) un' (0)
= + + .
s2 − (nbπ )2 s2 − (nbπ )2 s2 − (nbπ )2
 G (s) 
nπx
un (x) = L-1{un (s)} = L-1  2 n  + un(0) cosh( b ) +
b
nπ
un' (0)senh(nbπx )
s − ( nπ )2
 b 
= Pn(x) + un (0) cosh(nbπx ) + u (0)senh(nbπx ).
b '
nπ n

un(a) = Pn(a) + un(0) cosh(nbπa) + b u' (0)senh(nπa) ⇒ u' (0) = nπ un(a) − Pn(a) − un(0) cosh(nπa)
nπ n b n
[Link](nπa)  b 
b

Resumen de pasos a seguir:

1. Cálculo de los Gn (x) .
2. Cálculo de un (0) y un (a) .
3. Cálculo de L {Gn (x)} , Pn (x) y un' (0) .
4. Cálculo de los un (x) .
5. Construcción de u(x,y).

∞
Solución final: u(x, y) = ∑ u (x)sen(
n =1
n
nπy
b
).

Ejemplo 81. Placa semiinfinita.

∂2u ∂2u
+ = 0, 0 < x < L, y > 0
∂x2 ∂y2
u(0, y) = 0, u(L, y) = 0, y > 0
u(x, 0) = u0 , 0 < x < L , u(x, y) < M .
Solución.
X '' Y ''
X '' Y = −XY '' ⇒ =− = ctte ⇒ X ''− ctteX = 0 , Y ''+ ctteY = 0 .
X Y
X ''− ctteX = 0 , X(0) = X(L) = 0 ⇒ Xn (x) = C1sen(nLπx ) n = 1,...
nπy nπy
−
Yn (y) = C2e L + C3e L n = 1,...
∞ nπy nπy
 
∑ A
−
u(x, y) = n e L + Bn e L  sen(nLπx ) .
n =1 

Prof. José Luis Quintero 91

 nπy
Si y → +∞ entonces e L → +∞ por tanto para garantizar que la temperatura
esté acotada se tiene que hacer An = 0 . Usando la condición de borde se tiene
∞
2u0 L 2u0 1 + (−1)n +1 
u0 = ∑
n =1
Bnsen(nLπx ) ⇒ Bn =
L ∫ 0
sen( nLπx )dx = 
nπ
,

de modo que la solución se expresa como

∞ nπy ∞ 1 + ( −1)n +1  nπy
u(x, y) = ∑
n =1
Bne− L sen( nπx )
L
=
2u0
π ∑
n =1

n

e− L sen( nπx ) .
L

4.32. PROBLEMAS DONDE INTERVIENE LA ECUACIÓN DE

LAPLACE EN COORDENADAS POLARES

Pregunta.
Al usar coordenadas polares, ¿cómo se expresa la ecuación de Laplace dada en
coordenadas rectangulares como
∂2u ∂2u
+ = 0?
∂x2 ∂y2
Respuesta.
Sigamos los siguientes pasos:
a. Teniendo en cuenta que x = r cos θ , y = rsenθ , pruebe que:
∂x ∂x
r = x : Se tiene que r = [Link] θ = x .
∂r ∂r
∂y ∂y
r = y : Se tiene que r = [Link]θ = y .
∂r ∂r
∂x ∂x
= −y : Se tiene que = −[Link]θ = −y .
∂θ ∂θ
∂y ∂y
= x : Se tiene que = [Link] θ = x .
∂θ ∂θ
b. Emplee la regla de la cadena para funciones de dos variables y pruebe que:
∂u ∂u ∂u
r =x +y .
∂r ∂x ∂y
Prueba.
∂u ∂u ∂x ∂u ∂y ∂u ∂u
Se tiene que = . + . = .cos θ + .senθ . De modo que:
∂r ∂x ∂r ∂y ∂r ∂x ∂y
∂u ∂u ∂u ∂u ∂u ∂u ∂u ∂u
r = [Link] θ. + [Link]θ. =x +y . = −y +x .
∂r ∂x ∂y ∂x ∂y ∂θ ∂x ∂y
∂u ∂u ∂x ∂u ∂y ∂u ∂u
Se tiene que = . + . =− [Link]θ + [Link] θ .
∂θ ∂x ∂θ ∂y ∂θ ∂x ∂y
∂u ∂u ∂u
Entonces = −y +x . Suponga que la función u satisface las condiciones
∂θ ∂x ∂y
requeridas para que se aplique la regla de la cadena.
c. Pruebe que: (Diferenciando las ecuaciones en b)
∂u ∂2u ∂2u ∂2u ∂u ∂u
+r 2 = x +y + (cos θ) + (senθ) .
∂r ∂r ∂x∂r ∂y∂r ∂x ∂y

Prof. José Luis Quintero 92

 Demostración evidente.
∂2u ∂2u ∂2u ∂u ∂u
= − y + x −x −y .
∂θ2 ∂x ∂θ ∂y ∂θ ∂ x ∂y
Demostración evidente.
∂u ∂2u ∂2u ∂2u ∂2u ∂2u ∂2u
d. Emplee c. para obtener r + r2 2 + 2 = xr + yr −y +x .
∂r ∂r ∂θ ∂x∂r ∂y∂r ∂x∂θ ∂y∂θ
Demostración evidente.
∂u ∂u
e. Pruebe que: (Idea: sustituya y para u en b.)
∂x ∂y
∂2u ∂2u ∂2u
r =x 2 +y . Demostración evidente.
∂x∂r ∂x ∂x∂y
∂2u ∂2u ∂2u
r =x + y 2 . Demostración evidente.
∂y∂r ∂y∂x ∂y
∂2u ∂2u ∂2u
= −y 2 + x . Demostración evidente.
∂x∂θ ∂x ∂x∂y
∂2u ∂2u ∂2u
= −y + x 2 . Demostración evidente.
∂y∂θ ∂y∂x ∂y
∂u ∂2u ∂2u  ∂2u ∂2u 
f. Emplee d. y e. para establecer r + r2 2 + 2 = r2  2 + 2  .
∂r ∂r ∂θ  ∂x ∂y 

∂u ∂2u ∂2u ∂2u ∂2u ∂2u ∂2u
Usando d. se tiene r + r2 2 + 2 = xr + yr −y +x .
∂r ∂r ∂θ ∂x∂r ∂y∂r ∂x∂θ ∂y∂θ
Usando e. se tiene
∂u ∂2u ∂2u ∂2u ∂2u ∂2u ∂2u ∂2u ∂2u
r + r2 2 + 2 = x2 2 + xy + y2 2 + xy + y2 2 − xy +
∂r ∂r ∂θ ∂x ∂x∂y ∂y ∂y∂x ∂x ∂x∂y
∂2U ∂2U
x2 − xy
∂y2 ∂y∂x
En conclusión, en coordenadas polares, la ecuación de Laplace toma la forma:
∂2u 1 ∂u 1 ∂2u
+ . + . = 0.
∂r2 r ∂r r 2 ∂θ2

Ejemplo 82. Problema de Dirichlet para un círculo. Considere el problema de

resolver la ecuación de Laplace en la región circular r < a , sujeta a la condición de
frontera u(a, θ) = f(θ), en donde f es una función dada sobre 0 ≤ θ < 2π .
Solución.
Para completar el planteamiento del problema se observa que, para que u(r, θ) sea de
valores únicos, es necesario que u sea periódica en θ con período 2π . Además, se
expresa explícitamente que u(r, θ) debe estar acotada para r ≤ a , ya que esto será
importante posteriormente. Para aplicar el método se separación de variables a este
problema se supone que u(r, θ) = R(r)ϕ(θ) . Al sustituir en la ecuación diferencial se
obtiene
1 1
R ''(r)ϕ(θ) + .R '(r)ϕ(θ) + 2 .R(r)ϕ ''(θ) = 0 ,
r r
o bien,

Prof. José Luis Quintero 93

 R '' R' ϕ ''
r2 +r =− = K,
R R ϕ
en donde K es la constante de separación. Por tanto se obtienen las dos ecuaciones
diferenciales ordinarias r2R ''+ rR '− KR = 0 , ϕ ''+ Kϕ = 0 .

En este problema no hay condiciones homogéneas en la frontera; sin embargo,

recuerde que las soluciones deben ser acotadas y también periódicas en θ , con
período 2π . Se considerarán sucesivamente los casos en que K es negativa, cero y
positiva. Si K < 0, sea K = −λ2 , en donde λ > 0 ; entonces las ecuaciones ordinarias a
resolver son: r2R ''+ rR '+ λ2R = 0 , ϕ ''− λ2 ϕ = 0 . En particular la ecuación ϕ ''− λ2 ϕ = 0
tiene como solución ϕ(θ) = Aeλθ + Be−λθ , la cual es periódica solo si A = B = 0 , de modo
que K < 0, no ofrece una solución de interés. Si K = 0 , entonces las ecuaciones
ordinarias a resolver son r2R ''+ rR ' = 0 , ϕ '' = 0 . La ecuación ϕ '' = 0 tiene como
solución ϕ(θ) = A + Bθ , la cual es periódica si B = 0 , de modo que ϕ(θ) es constante.
Por otro lado la ecuación r2R ''+ rR = 0 es de Cauchy-Euler y tiene la solución
R(r) = C + D ln(r) . No es posible aceptar el término logarítmico, si U(r, θ) ha de
permanecer acotada cuando r → 0 ; de donde D = 0 . Por tanto, correspondiente a
K = 0 , se obtiene la solución U0 (r, θ) = D0 .

Por último, si K > 0 , se hace K = λ2 , en donde λ > 0 . Entonces las ecuaciones

ordinarias a resolver vienen dadas por r2R ''+ rR '− λ2R = 0 , ϕ ''+ λ2 ϕ = 0 . La ecuación
r2R ''+ rR '− λ2R = 0 es una ecuación de Cauchy-Euler y su solución es R(r) = Ar λ + Br −λ ,
mientras que la ecuación ϕ ''+ λ2 ϕ = 0 tiene solución ϕ(θ) = C cos(λθ) + Dsen(λθ) . Para
que ϕ(θ) sea periódica con período 2π es necesario que λ sea un entero positivo n.
Con λ = n se deduce que debe descartarse la solución r −λ , en virtud de que se vuelve
no acotada cuando r → 0 . Como consecuencia B = 0 y las soluciones apropiadas son
un (r, θ) = Enrn cos(nθ) + Fnrnsen(nθ) , n = 1, 2,... Estas funciones, junto con u0 (r, θ) = D0 ,
forman un conjunto de soluciones fundamentales del presente problema.

De la manera acostumbrada, suponga que u(r, θ) puede expresarse como una

combinación lineal de las soluciones fundamentales; es decir,
∞
E0
u(r, θ) =
2
+ ∑ r (E
n =1
n
n cos(nθ) + Fnsen(nθ)) .

Entonces, la condición de frontera requiere que

∞
E0
u(a, θ) = f(θ) =
2
+ ∑ a (E
n =1
n
n cos(nθ) + Fnsen(nθ))

donde,
1 2π 1 2π 1 2π
E0 =
2π ∫ 0
f(θ)dθ , En = n
aπ ∫ 0
f(θ) cos(nθ)dθ , Fn = n
aπ ∫ 0
f(θ)sen(nθ)dθ .

Prof. José Luis Quintero 94

Ejemplo 83. Problema de Dirichlet para un círculo.
∂2U 1 ∂U 1 ∂2U
+. + . =0
∂r2 r ∂r r2 ∂θ2
.
sen(θ) 0 ≤ θ ≤ π
U(1, θ) = f(θ) = 
 0 π ≤ θ ≤ 2π
Solución.
∞
E0
U(r, θ) =
2
+ ∑ r (E
n =1
n
n cos(nθ) + Fnsen(nθ))

donde,
1 2π 1 2π 1 2π
E0 =
π ∫ 0
f(θ)dθ , En =
an π ∫ 0
f(θ) cos(nθ)dθ , Fn =
anπ ∫ 0
f(θ)sen(nθ)dθ .

Para el caso de interés:

1 π 2
E0 =
π ∫ 0
sen(θ)dθ =
π
.

π π
1 1 n  1 
En =
π ∫ 0
sen(θ) cos(nθ)dθ = .
π n2 − 1 
[Link](ny) + cos y cos(ny)
n 0
1 n 1 n 2 (−1)n +1 − 1
n impar . .0 n par . .− En =
π n2 − 1 π n2 − 1 n π(n2 − 1)
π π
1 1 n  1 
Fn =
π ∫ 0
sen(θ)sen(nθ)dθ = .
π n2 − 1 
− cos [Link](ny) + cos ysen(ny) = 0
n 0
∞
(−1)n +1 − 1 n
∑
1 1
u(r, θ) = + r cos(nθ)
π π n =1 n2 − 1

Prof. José Luis Quintero 95

4.33. PROBLEMAS PROPUESTOS

1. Encuentre x • y para cada uno de los siguientes pares de vectores en R 3 .

a. x = (1 / 2,2, −1) y = (4, −2, 3)
b. x = (2 / 3,1 / 2,1) y = (−1 / 2, 4, 2)

2. Encuentre f • g para cada uno de los siguientes pares de vectores en C 0,1 .

a. f(x) = x g(x) = 1 − x2
b. f(x) = sen(πx / 2) g(x) = cos(πx / 2)
c. f(x) = x − 1
2
g(x) = 1
2
− x− 1
2

3. Encuentre f • g para cada uno de los siguientes pares de vectores en C 0,1 si el

producto interior está definido con respecto a la función de peso r(x) = ex .
a. f(x) = 1 − 2x g(x) = e−x
b. f(x) = x g(x) = 1 − x

4. Sea r cualquier función en C a,b que se anula para, cuando más, un número
finito de valores en el intervalo a,b . Demuestre que
b
f•g =
∫ a
f(x)g(x) r(x) dx

define un producto interior en C a,b .

5. Encuentre la longitud de los siguientes vectores en R 4 .

a. (3, 4, −3,1)
b. ( 12 , − 3, 2,1)
c. ( 23 , 12 , − 23 , − 12 )

6. Encuentre la distancia entre cada uno de los siguientes pares de puntos en R 4 .

a. x = (3, 0, −1,5) y = (2,2, −1, 3)
b. x = ( 23 ,1, 0, 2) y = ( 16 , 32 , − 12 , 32 )

7. Calcule f para cada uno de los siguientes vectores en C 0,1 .

a. f(x) = 1 − x2
b. f(x) = ln(x + 1)
c. f(x) = ex / 2
d. f(x) = sen( 2π x)

8. Encuentre el ángulo entre cada uno de los siguientes pares de vectores en R 3 .

a. x = (1,1,1) y = ( 12 , −1, 12 )
b. x = (1, −1, 0) y = (2, −1, 2)
c. x = (−3, −1, 0) y = (1, 2, − 5)

Prof. José Luis Quintero 96

9. Encuentre el coseno del ángulo de cada uno de los siguientes pares de vectores en
P3 si el producto interior es
1
p•q =
∫ −1
p(x)q(x)dx .

a. 1, x
b. x, x2
c. x,1 − x

10. Repita el apartado anterior usando ahora el producto interior

1
p•q =
∫ 0
p(x)q(x)dx .

11. Sean p y q polinomios en el espacio vectorial P2 (espacio de los polinomios de

coeficientes reales de grado a lo sumo dos). Discuta si la operación definida de la
forma p, q = p(−1)q(−1) + p(0)q(0) + p(1)q(1) es un producto interno en P2 .

12. Utilice la ley de los cosenos para establecer que

x•y
cos θ =
x y
para cualquier par de vectores no nulos x y y en R 2 .

13. Demuestre que x + y = x + y si, y solamente si, y = αx o x = αy para algún

número real α ≥ 0 .

14. Sean x y y vectores arbitrarios en un espacio euclidiano, y se supone x = y .

Demuestre que x + y y x − y son ortogonales. Interprete este resultado de forma
geométrica.

15. Sea (V, <, >) un espacio producto interno real y sean v, w elementos de V. Pruebe
2 2 2 2
que v+w + v−w = 2( v + w ) . Esta igualdad se conoce como ley del
paralelogramo.

16. Demuestre que las funciones dadas son ortogonales en el intervalo indicado:
a. f1(x) = x, f2 (x) = x2 , −2 ≤ x ≤ 2
3 2
b. f1(x) = x , f2 (x) = x + 1, −1 ≤ x ≤ 1
c. f1(x) = ex , f2 (x) = xe−x − e−x , 0≤x≤2
2
d. f1(x) = cos x, f2 (x) = sen x, 0≤x≤π

17. Demuestre que el conjunto de funciones dado es ortogonal en el intervalo indicado.

Encuentre la norma de cada función del conjunto.
a. senx, sen3x, sen5x, ..., 0 ≤ x ≤ 2π
b. {1, cos(npπx ), sen( mpπx )} n = 1, 2, 3,... m = 1, 2, 3,... −p ≤ x ≤ p

Prof. José Luis Quintero 97

18. Verifique mediante integración directa que las funciones son ortogonales con
respecto a la función de peso indicada, en el intervalo dado.
2
a. H0 (x) = 1, H1(x) = 2x, H2 (x) = 4x2 − 2 w(x) = e−x , −∞<x<∞
b. L 0 (x) = 1, L1(x) = −x + 1, L2 (x) = 1
2
x2 − 2x + 1 w(x) = e−x , 0 < x < ∞

19. Sea {φn(x)} un conjunto ortogonal de funciones en a ≤ x ≤ b tal que φ0 (x) = 1 y

φ1(x) = x. Demuestre que
b

∫ a
(αx + β)φn(x)dx = 0

para n = 2, 3,..., y constantes α y β cualesquiera.

20. Sea {φn (x)} un conjunto ortogonal de funciones en a ≤ x ≤ b . Demuestre que

2 2 2
φm (x) + φn (x) = φm (x) + φn (x) , m ≠ n .

21. Halle los valores de a, b y c para que f(x) = ax + 2 , g(x) = bx2 + cx + 1 y h(x) = x − 1
sean mutuamente ortogonales en el intervalo 0,1 .

22. Demuestre que las funciones f1(x) = 1 y f2 (x) = x son ortogonales, con respecto al
producto interno integral en el intervalo −1,1 , y determine las constantes A y B
de modo que la función f3 (x) = 1 + Ax + Bx2 sea ortogonal a f1 y f2 .

23. Sea C  − 12 , 12  con el producto interno

1/2
f, g =
∫ −1 / 2
r(x)f(x)g(x)dx ,

donde la función peso r está definida mediante

1
r(x) = .
1 − x2
Pruebe que la sucesión {fn }n = 0 , donde fn (x) = cos([Link](x)) es ortogonal.
∞

24. Ortogonalice cada una de las bases siguientes en R 3 .

a. (1,1, 0), (−1,1, 0), (−1,1,1)
b. (2, −1,1), (3, −2, 13 ), (1, 0,1)
c. (1, 0, 0), (0,1, 0), (−1, 0,1)

25. Ortogonalice cada uno de los conjuntos de vectores siguientes en C 0,1 .

a. ex , e − x
b. ex , e2x
c. 1, 2x, ex

26. Determine cuál de las siguientes sucesiones {ak } converge en R, y encuentre el

límite de cada sucesión convergente.

Prof. José Luis Quintero 98

 (−1)k
a. ak =
k
b. ak = 1 + (−1)k
c. ak = 1 + e−k
2k +1 + (−1)k
d. ak =
2k

27. Determine cuál de las siguientes sucesiones {xk } converge en R 3 , y encuentre el

límite de cada sucesión convergente. (En cada caso, los vectores e1 , e2 y e3 son
una base ortonormal para R 3 .
 1 k +1
a. xk =  1 +  e1 + e2 + e3
 k k +2
senk
b. xk = (1 − 2−k )e1 + e2 + (lnk)e3
k

28. Calcule los coeficientes de Fourier de la función f(x) = x2 , − π ≤ x ≤ π en C −π, π

con respecto al conjunto ortonormal
1 cos x cos 2x
, , ,...,
2π π π
y también con respecto al conjunto ortonormal
senx sen2x
, ,...
π π

29. Repita el ejercicio anterior, para la función ex , − π ≤ x ≤ π .

30. Considere que las funciones 1, senx, cos x,..., senkx, cos kx,... son una base para el
espacio euclidiano C −π, π . Encuentre los coeficientes de Fourier de f en C −π, π
en términos de esta base. ¿Cuál es el desarrollo en serie de f?

31. Encuentre todas las soluciones de cada uno de los problemas siguientes con valor
en la frontera.
a. y ''+ y = 0; y(0) = 1, y(π) = −1
b. y ''+ y = 0; y(0) = 0, y '(π) = 0
c. y ''+ 4y = sen(2x); y(0) = 0, y(π) = 0
d. y ''+ 5y = ex ; y(0) = 0, y(π) = 0
e. y ''+ 9y = 0; y '(0) = 0, y '(π) = 0
f. y ''+ 9y = cos(2x); y '(0) = 0, y '(π) = 0
g. y ''+ 9y = x; y '(0) = 0, y '(π) = 0

32. Calcule los valores propios y las funciones propias para los problemas con valores
en la frontera en los siguientes ejercicios:
a. y ''+ (1 + λ)y = 0 ; y(0) = 0, y(π) = 0

Prof. José Luis Quintero 99

 b. y ''+ 2y '+ (1 − λ)y = 0 ; y(0) = 0, y(1) = 0
c. 4y ''− 4y '+ (1 + λ)y = 0 ; y(−1) = 0, y(1) = 0
d. y ''+ 2y '+ (1 − λ)y = 0 ; y '(0) = 0, y '(π) = 0

33. La ecuación diferencial de Hermite y ''− 2xy '+ 2ny = 0, n = 0,1,2,... tiene
soluciones polinomiales Hn(x). Escriba la ecuación en su forma autoadjunta y
deduzca una relación de ortogonalidad.

34. ¿Cuál es la función propia correspondiente para el problema de valores en la

frontera y ''+ λy = 0, y '(0) = 0, y(π / 2) = 0 para λ = 25 ?

35. La ecuación diferencial de Laguerre xy ''+ (1 − x)y '+ ny = 0, n = 0,1,2,... tiene

soluciones polinomiales Ln (x) . Escriba la ecuación en su forma autoadjunta y
deduzca una relación de ortogonalidad.

36. Desarrolle la función g(x) = sen(x), 0 ≤ x ≤ π , en términos de las funciones propias

del problema de Sturm-Liouville
 y ''+ λy = 0
 .
 y(0) = 0, y(π) + y '(π) = 0

37. Determine los valores propios y las funciones propias del siguiente problema de
Sturm-Liouville
x2 y ''+ 3xy '+ λy = 0
 .
 y(1) = 0, y(e) = 0

38. Demuestre que las funciones

 πx   πx   nπx   nπx 
1, cos   , sen  a  ,..., cos  a  , sen  a  ,...
 a       
son linealmente independientes en PC −a, a .

39. Descomponga cada una de las siguientes funciones en sus partes par e impar:
x +1
a.
x −1
1
b.
x +1
c. xcox(x) − cos(2x)
∞
a0
d. f(x) =
2
+ ∑ (a
k =1
k cos(kx) + bk sen(kx)) ak ,bk cons tan tes

40. Encuentre el desarrollo en serie de Fourier de cada función dada:

a. f(x) = x, − π < x < π
b. f(x) = ex , −π<x<π

Prof. José Luis Quintero 100

 − x − 1 −π ≤ x < − 1
2 2

c. f(x) =  0 −2 ≤ x ≤ 2
1 1


 x − 2 <x≤π
1 1
2

d. f(x) = (π − x)(π + x), −π≤x≤π

x + π −π < x < 0
e. f(x) = 
x − π 0 < x < π

41. a. Encuentre el desarrollo en serie de Fourier de la función

0 −π < x ≤ 0
f(x) =  .
x 0 ≤ x < π
b. Use esta serie para demostrar que
π2 1 1 1
= 1 + 2 + 2 + 2 + ...
8 3 5 7

42. Encuentre el desarrollo en serie de Fourier de la función f(x) = cos(x) , 0 < x < π y
use este resultado para deducir que
2π 1 3 5 7
= 2 − 2 + 2
− 2
+ ... .
16 2 − 1 6 − 1 10 − 1 14 − 1

43. Encuentre el desarrollo en serie de Fourier de senos de ex , siendo 0 < x < π .

44. Sea f una función en PC 0, π  con la propiedad

π  π 
f  + x  = −f  − x  , 0 < x < π 2 .
2  2 
Demuestre que los únicos términos no cero en el desarrollo en serie de Fourier de
coseno de f son de la forma Ak cos(kx) , k impar donde
4 π2

Ak =  π ∫ 0
f(x) cos(kx)dx k = 1, 3,5, 7,...
 0 k = 0, 2, 4, 6,...


45. Encuentre el desarrollo en serie de Fourier de la función en PC −2, 2 definida por
 0 −2 < x < −1

f(x) =  x −1 < x < 1 .

0 1<x<2

46. Encuentre el desarrollo en serie de Fourier de la función

 1 8<x<9
f(x) =  .
10 − x 9 < x < 10

47. Sin recurrir a la integración, encuentre el desarrollo en serie de Fourier, x ∈ (−π, π) ,

de cada una de las siguientes funciones:
a. f(x) = sen3 (x)

Prof. José Luis Quintero 101

 x
b. f(x) = sen(x) cos2  
2
c. f(x) = sen2 (7x)sen2 (x) cos2 (x)

48. Demuestre que

 1 
n sen   n +  x 
1  2 
+ ∑
2 k =1
cos(kx) =
2sen  
x
,

2
como sigue:
Usando la fórmula de Euler eix = cos(x) + isen(x) se deduce que
n

∑ cos(kx)
1
+
2 k =1

es la parte real de la expresión

n

∑e
1 ikx
+ .
2 k =1

Evalúe esta expresión usando la fórmula para la suma de los n primeros términos
de una serie geométrica, y luego encuentre su parte real.

49. En cada uno de los ejercicios siguientes, encuentre la solución u(x,t) de la ecuación
de calor unidimensional que satisface las condiciones de frontera dadas:
 x 0<x≤L 2
a. u(0, t) = u(L, t) = 0; u(x, 0) = 
L − x L 2 ≤ x < L
b. u(0, t) = u(L, t) = 0; u(x, 0) = x(L − x)
 πx 
c. ux (0, t) = ux (L, t) = 0 u(x, 0) = ksen   , k ctte.
 L 
 πx 
d. u(0, t) = 100, u(L, t) = 0; u(x, 0) = 100 cos  
 2L 


 kUxx = Ut

e. U(0, t) = 0, U(L, t) = 0

U(x, 0) = 1 0 < x < 2
L

 
0 2 < x < L
L

uxx + x = ut , 0 < x < π , t > 0

f. u(0, t) = t , u(π, t) = 0 , t > 0
 2
 u(x, 0) = cos x , 0 < x < π

50. La temperatura en una barra aislada delgada, de longitud L, satisface las

condiciones en los puntos extremos u(0, t) = 0 y u(L, t) = 1 , y la condición inicial
 πx 
u(x, 0) = sen  .
 L 
a. Encuentre la temperatura en la barra como una función de posición y tiempo.

Prof. José Luis Quintero 102

 b. ¿Cuál es la temperatura de estado estable en la barra (esto es, la temperatura
cuando t → ∞ )?

51. Para una barra aislada de longitud 1, suponga que se mantiene la temperatura fija
con valor 0 en ambos puntos extremos y la distribución de temperatura inicial se
da según T(x) = 200sen(πx) − 100sen(2πx) . Suponga que la temperatura u(x,t)
satisface ut = uxx .
a. Encuentre u(x,t) para todo t ≥ 0 , 0 ≤ x ≤ 1 .
b. ¿Dónde se encuentra el punto más caliente de la barra cuando t = 0 ?

52. Suponga que la temperatura u(x,t) satisface ut = uxx .Suponga que para una barra
aislada de longitud 1, se mantiene la temperatura fija en 0 en ambos puntos
extremos, y la distribución inicial se da según la función T(x) = 3sen(2πx) .
a. Si la constante k en la ecuación de calor tiene valor numérico 0.05, encuentre
la distribución de temperatura u(x,1) en el tiempo t = 1 .
b. Si en el tiempo t = 1 la temperatura en cada punto ha decrecido a la mitad de
su valor inicial, ¿cuál es el valor numérico de k?

53. Resuelva el problema de valores en la frontera

 kuxx + Ae−βx = ut , β > 0, 0 < x < 1, t > 0

u(0, t) = 0, u(1, t) = 0, t > 0, u(x, 0) = f(x), 0 < x < 1
La ecuación diferencial es una forma de la ecuación de conducción de calor cuando
se genera calor dentro de una varilla delgada, por ejemplo, por decaimiento
radiactivo del material.

54. Resuelva el problema de valor de frontera

 ut = 1 + uxx , 0 < x < 1, t > 0
 .
u(0, t) = 0, u(1, t) = 0, u(x, 0) = sen(πx)

55. Conducción de calor con radiación. Una posibilidad interesante que puede
surgir es el caso donde uno de los extremos de la varilla, digamos x = 0 , se
mantiene a alguna temperatura, por ejemplo a 0 C , mientras que el otro extremo
x = L irradia a un medio circundante, el cual se asume que está también a la
temperatura 0 C . Se puede asumir como antes que la distribución de temperatura
inicial está especificada por φ(x) . Resuelva la siguiente ecuación diferencial que
modela la situación antes descrita:
 ut = a2uxx , 0 < x < L , t > 0

ux (0, t) = ux (L, t) + hu(L, t) = 0 .
 u(x, 0) = φ(x)


56. Encuentre la temperatura de estado estable en una varilla que satisface la

condición en la frontera ux (0, t) − u(0, t) = 0 en el extremo x = 0 y se mantiene a la
temperatura constante T en el extremo x = L .

Prof. José Luis Quintero 103

57. Resuelva la siguiente ecuación diferencial:
ut = 10uxx − 250u , 0 < x < 11 , t > 0

 ux (0, t) = K1 , u(11, t) = K2 .
 u(x, 0) = g(x)


58. Resuelva la ecuación diferencial parcial de la difusión de calor en una barra

delgada de longitud L, cuya difusividad térmica es α , aislada térmicamente en sus
dos extremos, si su temperatura inicial tiene la forma de un triángulo isósceles con
valor máximo de temperatura T0 > 0 en el punto medio de la barra y temperatura
nula en los dos extremos.

59. Sea una varilla de aluminio de longitud L que inicialmente está a la temperatura
uniforme de 25 C . Suponga que, en el instante t = 0 , el extremo x = 0 se enfría
hasta 0 C , mientras que el extremo x = L se calienta hasta 60 C y que, a partir
de entonces, los dos se mantienen a esas temperaturas.
a. Encuentre la distribución de temperatura en la varilla en cualquier instante t.
b. Use solamente el primer término de la serie para la temperatura u(x,t) con el
fin de hallar la temperatura aproximada en x = 5 cm , cuando t = 30 s y
cuando t = 60 s . (L = 20 cm)
c. Use los dos primeros términos de la serie para u(x,t) a fin de encontrar un
valor aproximado de u(5, 30) . ¿Cuál es la diferencia porcentual entre las
aproximaciones con uno y dos términos? ¿El tercer término tiene algún efecto
apreciable para este valor de t? (L = 20 cm)
d. Use el primer término de la serie para u(x,t) para estimar el intervalo de
tiempo que debe transcurrir antes de que la temperatura en x = 5 cm llegue a
ser menos de 1% de su valor de estado estable. (L = 20 cm)

60. Si hay flujo de calor en la superficie de la varilla en estudio hacia el medio que la
rodea cuya temperatura se mantiene constante, digamos um la ecuación de calor
se convierte en kuxx − h(u − um ) = ut . Considere la temperatura u(x,t) en una varilla
delgada con u(0, t) = u(L, t) = 0 y u(x, 0) = f(x) . En vez de estar aislado
lateralmente, la varilla pierde calor hacia el ambiente (con temperatura fija igual a
cero) a una razón proporcional a u(x,t). Demuestre que la solución de la ecuación
kuxx − hu = ut sujeta a las condiciones anteriores viene dada por la expresión
∞
 nπx  L
 nπx 
∑
2 π2t / L2 2
u(x, t) = e−ht
n =1
cne−n sen 
 L  , donde cn = L
 ∫ 0
f(x)sen 
 L 
 dx .

61. Suponga que se pierde calor a través de la superficie lateral de una varilla delgada
de longitud L, el cual pasa al medio ambiente que está a la temperatura cero. Si se
aplica la ley lineal de transferencia de calor, la ecuación adopta la forma
kuxx − hu = ut , 0 < x < L , t > 0 , donde h es una constante. Halle la temperatura
u(x,t) si la temperatura inicial en toda la varilla es f(x) y los extremos x = 0 y
x = L están aislados.

Prof. José Luis Quintero 104

62. Resuelva el problema de valor en la frontera
ut = uxx − α2u , 0 < x < L , t > 0

 ux (0, t) = K1 , u(L, t) = K2 .
 u(x, 0) = 0


63. Considere la ecuación auxx − but + cu = 0 , donde a, b y c son constantes.

a. Haga u(x, t) = eδt w(x, t) , donde δ es una constante, y encuentre la ecuación
diferencial parcial correspondiente para w.
b. Si b ≠ 0 , demuestre que puede elegirse δ de modo que la ecuación que se
encontró en la parte a) no tenga término en w.
c. Resuelva el problema
uxx − ut + u = 0 , 0 < x < π , t > 0

 ux (0, t) = u(π, t) = 0 , t > 0 .
 u(x, 0) = 0 , 0 < x < π


64. En cada uno de los ejercicios siguientes, encuentre la solución u(x,t) de la ecuación
de onda unidimensional que satisface las condiciones de frontera dadas:


 a2Uxx = Utt

a.  U(0, t) = 0, U(L, t) = 0

  2hx
 L
0 < x < 2L
U(x, 0) =  , Ut (x, 0) = 0
 ( )
2h 1 − L
x L
2
<x<L

 utt = a2uxx , 0 < x < L , t > 0


b. u(0, t) = ux (L, t) + hu(L, t) = 0 , t > 0
 u(x, 0) = f(x) , u (x, 0) = g(x)
 t

utt = a2uxx + f(x) , 0 < x < L , t > 0


c.  u(0, t) = α , u(L, t) = β
 u(x, 0) = 0 , ut (x, 0) = 0

 uxx + tsen3x = utt , 0 < x < π , t > 0

d.  u(0, t) = 0 , u(π, t) = π , t > 0
u(x, 0) = x , u (x, 0) = 0 , 0 < x < π
 t
−
 uxx + te = utt , 0 < x < π , t > 0
x

e.  u(0, t) = 0 , u(π, t) = 0 , t > 0
u(x, 0) = 1 , u (x, 0) = 0 , 0 < x < π
 t

65. Pruebe que la solución u(x,t) del problema

 utt = a2uxx , 0 < x < L , t > 0

 u(0, t) = 0 , u(L, t) = 0
u(x, 0) = f(x) , u (x, 0) = g(x)
 t

puede escribirse como u(x, t) = v(x, t) + w(x, t) , en donde v(x,t) es la solución del
mismo problema con g(x) = 0 y w(x,t) es la solución del mismo problema con

Prof. José Luis Quintero 105

 f(x) = 0 . De este modo, v(x,t) representa el movimiento de una cuerda a partir del
reposo con el desplazamiento inicial f(x) y w(x,t) representa el movimiento de una
cuerda que se pone en movimiento a partir de su posición de equilibrio con
velocidad inicial g(x). Por lo tanto, se puede obtener la solución del problema
general al resolver dos problemas subsidiarios.

66. Cuerda pulsada. Resuelva la siguiente ecuación diferencial:



 a2uxx = utt , 0 < x < L , t > 0

 u(0, t) = u(L, t) = 0 , t > 0 .

u(x, 0) =  bx 0 ≤ x ≤ L2
  , ut (x, 0) = 0
b(L − x) 2 ≤ x ≤ L
L


67. Cuerda vibrante con amortiguamiento pequeño. Resuelva la siguiente ecuación

diferencial:
a2uxx − β ut = utt , 0 < x < L , t > 0

 u(0, t) = u(L, t) = 0 , t > 0 .
 u(x, 0) = f(x) , ut (x, 0) = 0


68. Cuerda vibrante bajo la gravedad. Muestre que la solución de la ecuación

diferencial
a2uxx − g = utt , 0 < x < L , t > 0

 u(0, t) = u(L, t) = 0 , t > 0
 u(x, 0) = 0 , ut (x, 0) = 0

viene dada por
4gL2  1  πx   πat  1  3πx   3πat   gx(L − x)
u(x, t) = 3 2  3
sen   cos  L  + 3 sen  L  cos  L  + ... −
π a 1  L    3      2a2

69. Resuelva los siguientes problemas de valor en la frontera:

 a2uxx = utt , 0 < x < L , t > 0

a.  u(0, t) = 0 , ux (L, t) = K , t > 0
u(x, 0) = 0 , u (x, 0) = 0 , 0 < x < L
 t

 uxx + 1 = utt , 0 < x < 2π , t > 0


b.  u(0, t) = 1 , ux ( 2π , t) = − 2π , t > 0
 1 2
u(x, 0) = 1 + 2 x − 2 πx , ut (x, 0) = 0 , 0 < x < L
1

70. Pruebe que una solución del problema de valores en la frontera

Prof. José Luis Quintero 106

 

 uxx = utt + u , 0 < x < π , t > 0

u(0, t) = u(π, t) = 0 , t > 0 , ut (x, 0) = 0 , 0 < x < π

  x 0 < x < 2π
 u(x, 0) =  π
 π − x 2 ≤ x < π
∞
(−1)k +1
∑
4
es u(x, t) = sen((2k − 1)x) cos( (2k − 1)2 + 1 t) .
π k =1 (2k − 1)2

71. Al someter una cuerda vibratoria a una fuerza vertical externa, que varía en
función de la distancia horizontal al extremo izquierdo, la ecuación de onda tiene la
forma a2uxx + Ax = utt , donde A es una constante. Resuelva esta ecuación
diferencial parcial sujeta a u(0, t) = u(1, t) = 0, t > 0 , u(x, 0) = 0 , ut (x, 0) = 0 ,
0 < x < 1.

72. Resuelva el siguiente problema de valor en la frontera:

 a2uxx = utt , 0 < x < L , t > 0

 u(0, t) = 0 , u(L, t) = 0 , t > 0
u(x, 0) = x(L − x) , u (x, 0) = 0 , 0 < x < L
 t

73. Resuelva el siguiente problema de valor en la frontera:



 a2uxx = utt , 0 < x < L , t > 0

 u(0, t) = 0 , u(L, t) = 0 , t > 0
 ∞ ∞
 nπx   nπx 
u(x, 0) =

 n =1
∑Ansen 
 L 
 , ut (x, 0) =
n =1
∑Bnsen 
 L 
 , 0<x<L

74. Resuelva
 Uxx + Uyy = 0

U(0, y) = 0, U(1, y) = 1 − y .
 U (x, 0) = 0, U (x,1) = 0
 y y

75. Resuelva
 Uxx + Uyy = 0

U(x, 0) = x, U(x,1) = 1 − x .
 U (0, y) = 0, U (1, y) = 0
 x x

76. Resuelva
Urr + r −1Ur + r −2Uθθ = 0
 2
.
 U(1, θ) = 2 + sen (2θ)
1

Prof. José Luis Quintero 107

77. En cada uno de los siguientes ejercicios, encontrar la solución de la ecuación de
Laplace en el rectángulo 0 < x < L , 0 < y < M que satisface las condiciones de
frontera dadas:
a. u(0, y) = u(L, y) = 0, u(x, 0) = u(x,M) = 2x(x − L)

 0 0<x≤ L
2
b. u(0, y) = u(L, y) = 0, u(x, 0) = u(x,M) = 
L2
x2 − 3L
2
x+ 2
L
2
≤x<L
πy
c. u(x, 0) = u(x,M) = x(x − L), u(0, y) = u(L, y) = sen( M )

78. Deduzca la fórmula para la temperatura de estado estable en una placa cuadrada
de π × π cuyas caras están aisladas en forma tal que el calor no fluye a través de
ellas, bajo la suposición de que ux (0, y) = ux (π, y) = u(x, π) = 0, u(x, 0) = f(x) .

79. Resuelva la siguiente ecuación diferencial:

uxx + uyy = 0, 0 < x < a, 0 < y < b

 u(0, y) = 0, u(a, y) = 0, 0 < y < b .
u(x, 0) = 0, u(x,b) = f(x), 0 < x < a


80. Determine la temperatura de estado estable u(r, θ) en un anillo semicircular, si

u(a, θ) = θ(π − θ), u(b, θ) = 0, 0 < θ < π
.
u(r, θ) = 0, u(r, π) = 0, a < r < b

81. Determine la temperatura de estado estable u(r, θ) en una placa semicircular de

radio r = 2 , si
u 0 < θ < π 2
u(2, θ) =  0 ,
0 π 2 < θ < π
con u0 constante, y las orillas θ = 0 y θ = π están aisladas.

82. Determine la temperatura de estado estable u(r, θ) en una placa circular de radio
a; si se sabe que ur (a, 0) = f(θ) , 0 < θ < 2π .

RESPUESTAS A LOS PROBLEMAS PROPUESTOS

5
1. a. -5 b. 11/3 3. a. 0 b. 3-e 5. a. 35 b. 5/2 c. 6
2
7. a. 2
15
30 b. 2(1 − ln2) c. e −1 d. 1
2
2

9. a. 0 b. 0 c. − 12 10. 1
2
3, 1
4
15, 1
2

15. a. (1,1,1), (-1,1,0), (0,0,1) b. (2, −1,1), 1

18
(4, −11, −19), 2
83
(5, 7, −3)
c. (1,0,0), (0,1,0), (0,0,1)
2 2(e2 + e + 1) x
16. a. ex , e− x − 2 ex b. ex , e2x − e
e −1 3(e + 1)
c. 1, 2x − 1, ex + 6(e − 3)x + 10 − 4e
17. a. 0 b. No converge c. 1 d. 2 18. a. e1 + e2 + e3 b. No converge

Prof. José Luis Quintero 108

 eπ − e−π (−1)k
20. Respecto al primer conjunto c0 = ; ck = (eπ − e−π ), k = 1, 2,3,...
2π (1 + k2 ) π
(−1)k +1k
Respecto al segundo conjunto ck = (eπ − e−π ), k = 1, 2, 3,...
(1 + k 2 ) π
31. a. y = cos x + csenx c. No hay solución e. y = [Link](3x) g. No hay solución
32. a. λn = n2 − 1 , ϕn (x) = sen(nx), n = 1,2,... ortogonal en C 0, π
b. λn = −n2 π2 , ϕn (x) = e− x sen(nπx), n = 1, 2, 3,...; ortogonal en C 0,1 con respecto

a la función de peso e2x .

 x /2  nπx 
e sen  , n = 2, 4, 6,...
2 2   2 
c. λn = n π , ϕn (x) =  ortogonal en C −1,1 respecto a
ex / 2 cos  nπx  , n = 1,3,5,...
  2 
 
la función peso e− x .
d. λn = −n2 , ϕn (x) = e− x (n cos(nx) + sen(nx)), n = 1,2, 3,...

x2 + 1 2x 1 x
39. a. fE (x) = 2
, f0 (x) = 2
b. fE (x) = 2
, f0 (x) = −
x −1 x −1 1−x 1 − x2
c. fE (x) = − cos(2x) , f0 (x) = x cos(x)
∞ ∞
a0
d. fE (x) =
2
+ ∑
k =1
ak cos(kx) , f0 (x) = ∑ b sen(kx)
k =1
k

∞ 1 ∞ 
(−1)k +1 eπ − e−π (−1)k
40. a. 2 ∑
k =1
k
sen(kx) b.
π
 +
 2
∑1+ k
k =1
2
(cos(kx) − ksen(kx))

∞
1   n  
∑ nπ sen  2  cos(nx) + (−1) n
1 1 n n +1 2 1 2
c. − + + cos   − 2 sen    sen(nx)
4π 2 n =1
n nπ   n π
2 2
  
∞ ∞
(−1)k +1
∑ ∑ k sen(kx)
2 2 1
d. π +4 cos(kx) e. −2
3 k =1
k2 k =1
∞ k +1 ∞
π  (−1) − 1
k
(−1) 
∑ ∑
8 k
41. a. +  cos(kx) + sen(kx)  5. sen(2kx)
4 2
k π k  π k =1 4k − 1
2
k =1  
∞

∑
2 k 
43. 1 + (−1)k +1 eπ  sen(kx)
π k =1 1 + k2  

 2 4
 kπ − 2 2 k = 1,5, 9,...
 k π
∞  8
 kπx   − 2 2 k = 2, 6,10,...
∑
1
45. f(x) = + Ak cos   , donde Ak =  k π .
4 k =1  2   2 4
− − 2 2 k = 3, 7,11,...
 kπ k π
 0 k = 4, 8,12,16,...

∞
 (−1)k − 1 
∑ 
3 1
46. + cos(kπx) + sen(kπx) , 8 < x < 10
4 k =1
2 2
k π kπ 

Prof. José Luis Quintero 109

 3 1 1 1
47.a. sen(x) − sen(3x) b. sen(x) + sen(2x)
4 4 2 4
1 1 1 1 1
c. − cos(4x) + cos(10x) − cos(14x) + cos(18x)
16 16 32 16 32
∞
− (2m +1)2 π2t / L2a2   (2m + 1)πx 
∑ (2m + 1)
4L 1  
49.a. u(x, t) = 2 2
e sen  
π m =0  L 
∞
2 − (2m +1)2 π2t / L2a2   (2m + 1)πx 
∑ (2m + 1)
8L 1  
b. u(x, t) = 3 3
e sen  
π m=0  L 
∞
2k  − 4m2 π2t / L2a2  
 cos  2mπx  
∑
1 
c. u(x, t) = 1 − 2 e   
π  2
4m − 1  L  
 m =1
∞
 x  200 − n2 π2t / L2a2   nπx 
∑ n(4n
1  
d. u(x, t) = 100 1 −  + e sen  
 L  π n =1
2
− 1)  L 
∞
− cos(n2π ) +1
∑
2 2 π2 / L2 )t
e. u(x, t) = e−k(n sen(nLπx )
π n =1
n
∞
(−1)n − n
2 π2t / L2a2 
 nπx  x
∑
x 2 
50.a. u(x, t) = + e sen   b. u(x, t) =
L π n =1
n  L  L
2 2
51.a. u(x, t) = 200e−π t sen(πx) − 100e−4π t sen(2πx) b. x = 2 3
52.b. (ln 2) (4π2 )
∞
53. u(x, t) = ψ(x) + ∑A e
n =1
n
−kn2 π2t
sen(nπx) donde

A  −βx L
ψ(x) =
kβ2 
−e + (e−β − 1)x + 1 y An = 2
 ∫ 0
f(x) − ψ(x) sen(nπx)dx

∞
2(1 − cos(nπ))  −n2π2t
∑ 1 −
1
54. u(x, t) = x(1 − x) + e sen(nπx)
2 n =1
n3 π3 
∞
 h2 + λn2 L 
∑
2 2
55. u(x, t) = 2  2 2
 L(h + λn ) + h 0
n =1  
∫
φ(τ) cos(λn τ)dτ  e−a λnt cos(λnx)

donde λn son las raices positivas de λtg(λL) = h .

L L(h2 + λn2 ) + h
∫
2
Nota. Pruebe que sen(λnx) = sen2 (λnx)dx =
0 2(h2 + λn2 )
56. v(x) = T(1 + x) /(1 + L)
∞
 nπx 
∑  nπ + nπ cos(nπ) e
60x 50 70 −0.86n2 π2t / L2
59.a. u(x, t) = + sen  
L n =1  L 
b. u(5, 30) ≅ 12.6 C ; u(5, 60) ≅ 13.7 C
c. u(5, 30) ≅ 14.1 C ; 12% ; tercer tér min o ≅ −0.005 C d. t ≅ 160 s
1 ∞
L  L
 nπx    nπx  
∑∫
2 2 2 2
61. u(x, t) = e−ht 
 L ∫ 0
f(x)dx + 
L n =1  0
f(x) cos 
 L 
 dx  e−n π t / L cos 


 L 
∞
sen(n2π )
∑
8h
64.a. u(x, t) = cos(nπLat )sen(nLπx )
π2 n =1
n2

Prof. José Luis Quintero 110

 ∞
b. u(x, t) = ∑ A
n =1
n cos(aλnt) + Bnsen(aλnt) sen(aλnx) donde

1 L 1 L
An =
sen(λnx)
2 ∫ 0
f(x)sen(λnx)dx Bn =
aλn sen(λnx)
2 ∫ 0
g(x)sen(λnx)dx

λn son las raices positivas de htg(λL) = −λ .

∞ 2 (−1)n +1 + 1 
 cos(nt) + Mn t  sen(nx)
e. u(x, t) = ∑  

n =1 
nπ n2 

∞
 nπ   nπat   nπx 
∑n
4bL 1
66. u(x, t) = 2 2
sen   cos   sen  
π n =1  2   L   L 
∞
(2n − 1)πx   (2n − 1)πat 
∑ 2n − 1 sen 
2KL 1
69.a. u(x, t) = Kx +  cos  
π n =1
L   L 
∞
2  (2n − 1)π cos(nπ) + 2 
∑
x 4
b. u(x, t) = 1 − −   sen((2n − 1)x) cos((2n − 1)t)
2 π n =1  (2n − 1)3 
∞
(−1)n
∑
A 2A
71. u(x, t) = (x − x3 ) + cos(nπat)sen(nπx)
6a2 a2 π3 n =1
n3
∞
4L2  nπat   nπx 
∑n
1 
72. u(x, t) = 3 
1 + (−1)n +1  cos   sen  
π3   L   L 
n =1
∞
  nπat   nπat    nπx 
∑ A
L
73. n cos 
L  + nπa Bnsen  L   sen  L 
n =1      
∞
1 − (−1)n
∑ n senh(nπ) senh(nπx) cos(nπy)
1 2
74. U(x, y) = x+ 2 2
2 π n =1

75.
4 + 4e−nπ 4enπ + 4
∑
1
U(x, y) = + (Anenπy + Bne−nπy ) cos(nπx) An = , Bn = −
2 n impar
n2 π2 (enπ − e−nπ ) n2 π2 (enπ − e−nπ )

76. U(r, θ) = 1 − 12 r 4 cos(4θ)

∞
  (2m + 1)π(M − y)   (2m + 1)πy    (2m + 1)πx 
u(x, y) = ∑C
m =0
m senh 
  L 

+ senh 
 L 

 sen 
 L 

77.a.
16L2
donde Cm = −
(2m + 1)3 π3senh (2m + 1)πM / L 
∞
  nπ(M − y)   nπy    nπx 
u(x, y) = ∑C
n =1
n senh 
  L 

+ senh 
 L 
  sen 
 L 

b.
L2   nπ    nπ   
donde Cn = 3 3 nπsen   + 4 (−1)n − cos   
n π senh(nπM / L)   2    2   

Prof. José Luis Quintero 111

 u(x, y) =
{senh(πx / M) + senh π(L − x) / M} sen(πy / M) +
senh(πL / M)
∞
  (2M + 1)π(M − y)   (2m + 1)πy    (2m + 1)πx 
c. ∑C
m =0
m senh 
  L  + senh 
  L   sen 
  L 

8L2
donde Cm = −
(2m + 1)3 π3senh (2m + 1)πM / L 
∞ π
a0 senh(n(π − y)) cos(nx)
∑
2
78. u(x, y) =
2π
(π − y) +
n =1
an
senh(nπ)
donde an =
π ∫ 0
f(x) cos(nx)dx

∞
nπy   nπx  a
 nπx 
∑ A sen 
2
( )∫
79. u(x, y) = sen   donde An = f(x)sen   dx
a 
n nπb
n =1  a  asenh 0  a 
a
∞ n
1 − (−1)n r2n − b2n  a
∑
4
80. u(r, θ) = . 2n  r  sen(nθ)
π n =1
n 3
a − b2n  
u0 2u0
∞ sen ( n2π )  r n cos(nθ)
81. u(r, θ) =
2
+
π ∑
n =1
n 2
 

Prof. José Luis Quintero 112