FOLLETO:

CURSO DE ECUACIONES DIFERENCIALES EN

DERIVADAS PARCIALES
Ph.D. Idulfo Arrocha Rodrı́guez

Mayo de 2014.

Resumen
El presente escrito constituye un resumen del contenido del curso de Ecuaciones Dife-
renciales en Derivadas Parciales, el objetivo de mismo es brindar un panorama general de
la asignatura. Sin embargo, para una ampliación de los contenidos es necesario revisar la
bibliografı́a que aparece al final del mismo.

1
Índice
1. GENERALIDADES ACERCA DE LAS ECUACIONES DIFERENCIALES 4
1.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.2. Definición de una ecuación diferencial . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3. Clasificación de las EDPs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3.1. De acuerdo al orden . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3.2. De acuerdo al grado de linealidad . . . . . . . . . . . . . . . . . . . . . . . 7
1.4. Forma general de una EDP de segundo orden . . . . . . . . . . . . . . . . . . . . 7
1.4.1. Clasificación de una EDP de Segundo Orden . . . . . . . . . . . . . . . . 7
1.5. Solución de una EDP . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.6. Condiciones de contorno o frontera. Condiciones iniciales . . . . . . . . . . . . . . 9
1.6.1. Dominios y Fronteras . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.6.2. Condiciones de Contorno . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.6.3. Condiciones Iniciales: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
1.7. Problemas bien planteados y mal planteados . . . . . . . . . . . . . . . . . . . . . 12
1.8. Método de Separación de Variables . . . . . . . . . . . . . . . . . . . . . . . . . . 13
1.9. Generalidades de las EDPs Clásicas de Segundo Orden. . . . . . . . . . . . . . . 15
1.9.1. La Ecuación del Calor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.9.2. La ecuación de onda . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1.9.3. La Ecuación de Laplace . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

2. ECUACIONES DIFERENCIALES PARCIALES DE PRIMER ORDEN 18

2.1. Clasificación de las EDPs de primer orden . . . . . . . . . . . . . . . . . . . . . . 18
2.2. Método de las Caracterı́sticas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.3. Ecuaciones No lineales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.4. Movimiento de ondas unidireccional. . . . . . . . . . . . . . . . . . . . . . . . . . 24

3. ELEMENTOS DE LAS SERIES, INTEGRALES y TRANSFORMADA DE

FOURIER 26
3.1. Funciones Ortogonales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.1.1. Producto Interno . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.1.2. Producto interno de funciones . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.1.3. Desarrollo en series ortogonales . . . . . . . . . . . . . . . . . . . . . . . . 28
3.2. Series de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.2.1. Serie de Fourier de cosenos y senos . . . . . . . . . . . . . . . . . . . . . . 32
3.3. Integrales de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33
3.3.1. Deducción de la integral de Fourier . . . . . . . . . . . . . . . . . . . . . . 33
3.4. La Transformada de Fourier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36
3.4.1. Propiedades de la Transformada de Fourier . . . . . . . . . . . . . . . . . 37

4. PROBLEMA DE AUTOVLORES DE STURN-LIOUVILLE 41

4.1. Problema regular de autovalores de Sturn-Liouville . . . . . . . . . . . . . . . . . 41

5. SOLUCIÓN DE ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN 43

5.1. Solución de la Ecuación del Calor . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
5.2. Solución de la Ecuación de Ondas . . . . . . . . . . . . . . . . . . . . . . . . . . . 45

6. REFERENCIAS BIBLIOGRÁFICAS 47

3
1. GENERALIDADES ACERCA DE LAS ECUACIONES DI-
FERENCIALES
1.1. Introducción
Una Ecuación en Derivadas Parciales (EDP), describe una relación entre una función incógni-
ta y sus derivadas parciales. Las EDPs aparecen a menudo en todas las áreas de la fı́sica y la
Ingenierı́a. Sin embargo, recientemente se ha visto un dramático incremento en el uso de las
EDPs en áreas tales como la biologı́a, quı́mica, ciencias computacionales (particularmente con
relación al procesamiento de imágenes y gráficas) y en la economı́a (finanzas). En efecto, en
cada área donde existe una interacción entre un número de variables independientes, intentamos
definir funciones que modelen los procesos y conduzcan a ecuaciones que estableces las relaciones
entre las mismas. Cuando el valor de la función incógnita en cierto punto depende solamente
de lo que sucede en la vecindad de dicho punto, podemos en general, obtener una EDP, la cual
será complementada por condiciones adicionales tales como las condiciones iniciales (como se ve
a menudo en la teorı́a de las ecuaciones diferenciales ordinarias ( EDOs) o las condiciones de
frontera.
El análisis de una EDP tiene muchas facetas. El enfoque clásico que predominó durante el
siglo XIX fue para desarrollar métodos para encontrar soluciones explı́citas. Dada la inmensa
importancia de las EDPs en los diferentes brazos de la fı́sica, cada desarrollo matemático que
producı́a una nueva clase de EDP estuvo acompañado de un progreso significativo de la Fı́sica.
Ası́, el método de las caracterı́sticas inventado por Hamilton dió lugar a mayores avances en
la óptica y en la mecánica analı́tica. EL método de Fourier, por su parte, permitió la solución
de problemas de transferencia de calor y propagación de ondas, mientras que el método de
Green fué fundamental en el desarrollo de la teorı́a del electromagnestismo. El mayor progreso
en el desarrollo de las EDPs se ha logrado en los últimos 50 años gracias a la introducción de
métodos numéricos que permiten el uso de la computadora para resolver casi cualquier EDP,
con geometrı́as complejas y condiciones externas arbitrarias. En definitiva, no es posible estudiar
las EDPs de forma desligada de su aplicación a la ciencia. En ese sentido, en el enfoque que le
daremos al presente curso no pretendemos hacer un tratado riguroso de las EDPs, sino más bien
una descripción de los conceptos fundamentales de las EDPs ası́ como su interpretación en los
problemas que han producido su desarrollo. Creemos que el estudio es mucho más significativo
si se ve desde esta perspectiva.

1.2. Definición de una ecuación diferencial

¿Qué es una ecuación en derivadas parciales?.
Una ecuación diferencial en derivadas parciales (EDP) es una relación de la forma:

F (x, t, u, ∂x , u, ..., ∂xn , ∂t u, ..., D α u) = 0 (1)

donde u(x, t) es una función de variable espacial x = (x1 , ..., xn ) ∈ R y de la variable temporal
t.
La incógnita u representa una cantidad fı́sica como por ejemplo: la temperatura, concen-
tración de un contaminante, intensidad de una señal acústica, deformación de una estructura,
etc.
Cuando la variable espacial sólo tiene una componente, se dice que es un problema unidi-
mensional, es decir u(x, t) = (x, t). Mientras que si la variable espacial tiene dos componentes,
esto es: x = (x1 , x2 ), entonces se dice que el problema es bidimensional. De manera análoga si la
variable espacial tiene 3 componentes, entonces el problema es tridimensional. Por otro lado, si la

4
ecuación (1) no depende de la variable temporal, es decir: u = u(x), x = (x1 , ..., xn ), entonces
se dice que el problema es estacionario, mientras que si hay dependencia espacio-temporal se
dice que el problema es de evolución.
La fórmula (1) establece una relación entre la incógnita u, sus derivadas parciales hasta cierto
orden y el punto (x, t) espacio-temporal.

Ejemplos: Ecuaciones diferenciales parciales importantes.

∂u ∂u ∂2u
+u =ν 2 Ecuación unidimensional de Burges viscosa. (2)
∂t ∂x ∂x
∂2u 2
2∂ u
= c Ecuación unidimensional de ondas (3)
∂t2 ∂x2
∂u ∂2u
=k 2 Ecuación unidimensional del calor (4)
∂t ∂x
 2 
∂u 2 ∂ u ∂2u
=c + 2 Ecuación bidimensional del calor (5)
∂t ∂x2 ∂y
∂2u ∂2u
+ 2 =0 Ecuación bidimensional de Laplace (6)
∂x2 ∂y
∂2u ∂2u
+ 2 = f (x, y) Ecuación bidimensional de Poisson (7)
∂x2 ∂y
∂2u ∂2u ∂2u
+ 2 + 2 =0 Ecuación tridimensional de Laplace (8)
∂x2 ∂y ∂z
∂u
De esta manera, en la notación habitual (Legendre) para las derivadas ∂xj u = ∂xj = uxj
∂u
mientras que ∂t u = ∂t .
En esta notación se tiene que las derivadas parciales de orden superior e iteradas se escriben
como:
∂u ∂2u ∂3u
= uxj , = ux j x k = ux j x k x v (9)
∂xj ∂xj ∂xk ∂xj ∂xk ∂xv
En la fórmula (1) también se ha utilizado la notación de Schwartz llamada la notación de los
multiı́ndices.

Definición: Un vector α
~ = (α1 , ..., αn ) ∈ R tal que cada uno de los αk es un entero no ne-
gativo, se denomina mutiı́ndice de orden [~ α] = α1 + α2 + ... + αn .
De manera que dado un multiı́ndice α ~ , se suele utilizar la siguiente notación para una derivada
parcial de orden [~
α].

∂ [~α] u(x1 , ..., xn )

D [~α] u(x1 , ..., xn ) = = ∂xα11 ∂xα22 ...∂xαnn u(x1 , ..., xn ) (10)
∂xα1 ...∂xαnn
De este modo, el orden de la EDP (1) es el de la derivada de mayor orden, esto es, el máximo
de los módulos |α| de los ı́ndices α de la fórmula (1).
Cuando la incógnita de u de (1) no es una función escalar sino un vector
 
u1
u =  ...  ,
 

un

5
la expresión (1) puede también representar un sistema de ecuaciones. En ese caso F es también
una función vetorial. Si F tiene M componentes, (1) es pues un sistema de M ecuaciones con
N incógnitas.
La ecuación (1) puede ser lineal o no lineal dependiendo de que F lo sea o no en relación
a la incógnita u y sus derivadas ∂ α u. Las EDPs no-lineales, lejos de ser problemas puramente
académicos, intervienen de manera decisiva en muchos e importantes ámbitos de las ciencias y
la tecnologı́a.

Ejemplos:

1. La EDP ux + ex+y uy − u = x2 y 2 es lineal de orden 1

2. La EDP ux xu + uy + xy = 0 es no-lineal de orden 2, el término no-lineal es este caso es

ux xu.

Si F es un polinomio de grado uno en D α se dice que la EDP es lineal. En este caso la EDP es
de la forma X
aα (x)D α u − f (x) = 0, (11)
α
P
donde α significa que la suma se extiende a un conjunto finito de multi-ı́ndices α con |α|≥ 0.
Las funciones aα (x) y f (x) se suponen dadas. Normalmente cuando tratamos con una ecuación
como (11), la escribimos como: X
aα (x)D α u = f (x), (12)
α

donde f (x) es el término no homogeneo de la ecuación, también llamado término fuente. Si

f (x) = 0 diremos que la EDP lineal es homogénea. En general, las EDPs no lineales son mucho
más difı́ciles de tratar que las lineales.

1.3. Clasificación de las EDPs

. Las EDPs se pueden clasificar de acuerdo a su orden, linealidad y tipo de condición de
frontera.

1.3.1. De acuerdo al orden

Como ya mencionamos, el orden de una EDP está determinado por la derivada parcial de
mayor orden presente en la expresión, como se muestra en los siguientes ejemplos
∂u ∂u
−k =0 (13)
∂x ∂y
∂2u ∂u
2
−k =0 (14)
∂x ∂y
∂2u ∂2u ∂2u
+ 2 + 2 =0 (15)
∂x2 ∂y ∂z
 3 2
∂ u 2
∂ u ∂2u
+ + 2 =0 (16)
∂x3 ∂y 2 ∂z

como podemos ver, la ecuación (13) es de primer orden, (14 y 15) son de segundo orden, mientras
que (16) es de tercer orden.

6
1.3.2. De acuerdo al grado de linealidad
Las EDPs también se pueden clasificar según su linealidad en lineales, cuasilineales y no
lineales. Consideremos la siguiente EDP:

∂2u ∂2u ∂2u

a + b + c +d=0 (17)
∂x2 ∂x∂y ∂y 2
en la que podrı́an presentarse los siguientes casos:

1. Si los coeficientes a, b y c en (17) son constantes o funciones de la variable independiente

(x,y), entonces se trata de una ecuación diferencial parcial lineal, por ejemplo:

∂u ∂u
−k =0
∂x ∂y

2. Si los coeficientes a, b y c en (17) son funciones de la variable

 dependiente
 y/o de sus
∂u ∂u
derivadas de menor orden que el de la ecuación diferencial x, y, ∂x , ∂y , entonces se trata
de una EDP cuasilineal, por ejemplo:

∂2u ∂u
2
−u =0
∂x ∂y

3. Si los coeficientes a, b, c y d en (17) son funciones de derivadas del mismo orden que en la
EDP  
∂2u
x, y, u, v 2
∂x
entonces se trata de una ecuación diferencial no lineal, por ejemplo:
 2
∂2u ∂2u ∂2u
+ + 2 =0
∂x2 ∂y 2 ∂z

1.4. Forma general de una EDP de segundo orden

La forma general de una ecuación diferencia parcial de segundo orden se puede presentar por
la siguiente expresión:

∂2u ∂2u ∂2u ∂u ∂u

a + b + c +d +e + fu = g (18)
∂x2 ∂x∂y ∂y 2 ∂x ∂y

donde los coeficientes a, b, c, d, e, f y g de (18) son constantes o funciones de las variables inde-
pendientes.
Si el parámetro g es cero, entonces se trata de una ecuación diferencial parcial homogénea.

1.4.1. Clasificación de una EDP de Segundo Orden

Las tres formas canónicas de la EDP de segundo orden (18) son determinadas por los si-
guientes criterios:

Si b2 − 4ac < 0 es una ecuación elı́ptica

Si b2 − 4ac = 0 es una ecuación parabólica

7
 Si b2 − 4ac > 0 es una ecuación hiperbólica

Ejemplos

1. La ecuación de conducción de Calor o ecuación de difusión

∂2u ∂u
k 2
= con k>0
∂x ∂t
es una EDP parabólica pues b = 0, c = 0 y a = k y por lo tanto b2 − 4ac = 0.
En esta ecuación k es la constantede conductividad y depende del material en que se lleve
a cabo el fenómeno de difusión. La variable dependiente u(x, t) representa la temperatura
de una varilla en una posición y tiempo dado.

2. La ecuación de ondas unidimensional (3) es hiperbólica, mientras que la ecuación de calor

bidimensional (5) es parabólica. Por otro lado, las ecuaciones (6-8) son elı́pticas.

3. Clasifique las siguientes ecuaciones en parabólicas, hiperbólicas o elı́pticas.

∂2u ∂u ∂2u ∂2u ∂2u ∂2u

a) 3 = b) = c) + 2
∂x2 ∂y ∂x2 ∂y 2 ∂x2 ∂y
Solución:
a) Para este caso a = 3,e = −1, b = c = d = f = 0, por lo tanto b2 − 4ac = 0, lo que
significa que la ecuación es parabólica.
b) En este caso a = 1 y c = −1, por lo que b2 − 4ac = 4 > 0, lo cual significa que la
ecuación es hiperbólica.
c) En este caso a = 1 y c = 1, luego b2 − 4ac = −4 < 0, es decir se trata de una ecuación
elı́ptica.

Las ecuaciones parabólicas determinan cómo una incógnita varı́a tanto en el espacio como
en el tiempo, lo cual se manifiesta por la presencia de las derivadas espacial y temporal, como
la ecuación de conducción de calor (4). Tales casos se conocen como problemas de propagación,
puesto que la solución se propaga, o cambia, con el tiempo.
La categorı́a hiperbólica, también tiene que ver con problemas de propagación. Sin embar-
go, una importante diferencia manifestada por la ecuación de onda (3), es que la incógnita se
caracteriza por una segunda derivada con respecto al tiempo. En consecuencia, la solución oscila.
Una cuerda tensa que vibra a baja amplitud es un sistema fı́sico simple que puede caracterizarse
por una EDP hiperbólica.
Comúnmente, las ecuaciones elı́pticas se utilizan para caracterizar sistemas en estado
estacionario. Como en la ecuación de Laplace (6). Esto se indica, por la ausencia de una derivada
con respecto al tiempo. Ası́, estas ecuaciones se emplean para determinar la distribución en
estado estacionario de una incógnita en dos dimensiones espaciales.

1.5. Solución de una EDP

Definición: Una función u de las variables x1 , x2 , ..., xn es una solución clásica de la EDP
(1) en una región R si se cumplen las siguientes condiciones:

1. En R la función u(x1 , ..., xn ) adminte todas las derivadas que aparecen en F .

2. La función u(x1 , ..., xn ) y todas las derivadas parciales son contı́nuas en la región R.

8
 3. Al sustituir en la EDP, u y sus derivadas, la ecuación (1) se convierte en una identidad
respecto a las variables independientes del problema.

4. En ocasiones se exige que la solución sea contı́nua hasta la frontera.

En general, la totalidad de las soluciones de una ecuación diferencial parcial es muy grande.
Por ejemplo las funciones

u = x2 + y 2 , u = ex cosy, u = ln(x2 + y 2 ) (19)

que son por completo diferentes entre sı́, son soluciones de la Ecuación Bidimensional de Laplace
(6), como podrán ustedes comprobar.
Cuando se estudiaron las Ecuaciones Diferenciales Ordinarias (EDOs) se vió que si una EDO
es lineal y homogénea, entonces, a partir de soluciones conocidas, pueden obtenerse otras solu-
ciones por superposición. Para una EDP lineal y homogénea la situación es muy parecida. De
hecho tenemos el siguiente teorema:

Definición: Principio de Superposición o linealidad

Si u1 y u2 son soluciones cualesquiera de una ecuación diferencial parcial lineal y homogénea en
alguna región R, entonces
u = c1 u1 + c2 u2
donde c1 y c2 son constantes cualesquiera, también es una solución de esa ecuación en R.
No existe una teorı́a general conocida en relación con la resolubilidad de todas las EDPs.
Tal teorı́a es extremadamente improbable que exista, dado la rica variedad de fenómenos fı́sicos,
geométricos, probabilı́sticos que puede ser modelados por las EDPs. En lugar de ello, la inves-
tigación se centra en una diversidad de casos particulares de EDPs que son importantes para
aplicaciones dentro de y fuera de las matemáticas, con la esperanza de que la percepción de los
orı́genes de estas EDPs nos den pistas sobre sus soluciones.

1.6. Condiciones de contorno o frontera. Condiciones iniciales

1.6.1. Dominios y Fronteras
En general cuando consideramos una EDP (1) la función incógnita u = u(x, t) se supone
definida sobre un conjunto dado Ω de Rn . Supondremos siempre que satisface las dos condiciones
siguientes:

Ω es un conjunto abierto. Esto es, para todo punto a ∈ Ω, existe un radio r > 0 tal que
para todo punto x ∈ Rn , tal que d(x, a) < r entonces x ∈ Ω.

Ω es conexo. Es decir, no es posible encontrar dos conjuntos abierto no vacı́os Ωi , (i =

1, 2) tales que Ω1 ∩ Ω2 = Ø y Ω1 ∪ Ω2 = Ω

En tal caso diremos que Ω es un dominio de Rn . La frontera ∂Ω de Ω es el conjunto formado

por los puntos a ∈ Rn , tales que todo radio r > 0, existen puntos x, tanto dentro (x ∈ Ω) como
fuera (x ∈
/ Ω).

Definición: Diremos que un conjunto R ∈ Rn es una región si está formado por la unión
de un conjunto abierto y conexo de Rn junto con una parte de su frontera (∂Ω). Si en R no
está incluida la frontera del abierto, se dice que es una región abierta. Por el contrario si la
frontera está incluida en R se dice que es una región cerrada.

9
1.6.2. Condiciones de Contorno
Dada una EDP sobre un dominio Ω ⊂ Rn :

F (x, D α u) = 0, x ∈ Ω, (20)

normalmente se nos pide no sólo que encontremos una función u = (x) que satisfaga la EDP en
todo punto de Ω, sino también que tal función satisfaba una serie de condiciones:

fi (x, D α u) = 0, x ∈ ∂Ωi , i = 1, ..., m,

donde las funciones fi son funciones dependientes de las variables xi y de un número finito de
derivadas D α u con |α|≥ 0. Condiciones de esta clase de denominan condiciones de contorno
(también se conocen como condiciones de frontera). Un problema consistente en resolver una
EDP sobre un dominio Ω y un conjunto de condiciones de contorno se denomina problema de
contorno (o problema de frontera).

1. Condiciones de contorno de tipo Dirichlet(de primer tipo): Este tipo de condición

de contorno se da cuando en una EDP, se le especifican los valores de la solución que
necesita la frontera del dominio. La cuestión de hallar las soluciones a esas ecuaciones con
esta condición se le conoce como Problema de Dirichlet. Esto se expresa como:

u|∂Ω = u(x) = g(x), x ∈ ∂Ω

Este tipo de condición de contorno establece un valor fijo para cada uno de los puntos
de la solución que están en la frontera. Dependiendo del problema al cual se le aplique la
condición de Dirichlet, se tendrá una interpretación distinta de este tipo de condición de
contorno.

2. Condiciones de contorno de tipo Neumann (de segundo tipo): Este tipo de con-
diciones tienen la forma: 
∂u 
= g(x), x ∈ ∂Ω,
∂n ∂Ω
∂
donde n es el vector normal a la frontera del dominio, por lo tanto ∂n es la derivada
direccional normal a la frontera.

10
 3. Condiciones de contorno de tipo Robin o mixta (de tercer tipo): Este tipo de con-
dición de contorno es una combinación lineal de las condiciones de Dirichlet y de Neumann.
Esto es:  
∂u 
au + b = g(x), x ∈ ∂Ω,
∂n ∂Ω
a y b en este caso son funciones dadas.

4. Condición de contorno de Cauchy: Un problema que combina las condiciones de

Dirichlet para una parte de la frontera y de Neumann para otra parte se dice que es un
problema de Cauchy. Esto es:

∂u 
u|∂Ω = g(x), = h(x), x ∈ ∂Ω,
∂n ∂Ω
donde g y h son funciones dadas.

5. Condiciones de contorno periódicas: Este tipo de condiciones tienen la forma:

u|∂Ω = u(x) = u(x + L), x ∈ ∂Ω

donde L es el periódo de la función u.

1.6.3. Condiciones Iniciales:

Otro tipo de condiciones que se suelen exigir a las soluciones de una EDP son las denomina-
das condiciones iniciales respecto a la variable independiente temporal t. Normalmente son un
conjunto de condiciones de la forma:
 
∂u  ∂ r−1 u 
u|t=t0 = Φ0 , = Φ1 , ... = Φr−1
∂t t=t0 ∂tr−1 t=t0

donde r > 0 y las funciones Φi dependen del resto de las variables independientes, esto es:
Φi = Φi (x).
En problemas fı́sicos en los que se analiza la evolución de un sistema suelen coexistir tanto
condiciones de contorno como iniciales.

Ejemplos:
1. Ecuación del calor unidimensional con temperatura cero en los lı́mites finitos
La ecuación en derivadas parciales serı́a:
∂u ∂2u
= k 2, 0 < x < L, t>0
∂t ∂x
Las condiciones de contorno (C.C.) son:

u(0, t) = 0, u(L, t) = 0,

mientras que las condición inicial (C.C.) tiene la forma:

u(x, 0) = f (x)

Como podemos observar, se trata de un problema de Dirichlet puesto que las condiciones
de contorno son de tipo Dirichlet, además es un problema de evolución temporal con una
condición inicial para el tiempo.

11
 2. Ecuación del calor unidimensional con fronteras aisladas
Las ecuaciónes correspodientes son las siguientes:

∂u ∂2u
EDP = k 2, 0 < x < L, t>0
∂t ∂x
∂u ∂u
C.C. (0, t) = 0, (L, t) = 0
∂x ∂x
C.I. u(x, 0) = f (x)

Como se puede observar, se trata de un problema con condiciones de contorno de tipo

Neumann y una condición inicial dependiente de x.

3. Ecuación de la cuerda vibrante con extremos fijos

Las ecuaciones correspodientes son:

∂2u 2
2∂ u
EDP = c , 0 < x < L, t > 0
∂t2 ∂x2
C.C. u(0, t) = 0, u(L, t) = f (x)
∂u
C.I. u(x, 0) = f (x) (x, 0) = g(x)
∂t
Este también es un problema de Dirichlet con condiciones iniciales para la función y su
primera derivada.

1.7. Problemas bien planteados y mal planteados

¿Cuál deberı́a ser la solució ideal de una ecuación diferencial?. Idealmente, deberı́amos ob-
tener soluciones explı́citas (analı́ticas) en términos de funciones elementales. En la práctica esto
solamente es posible para EDPs muy simples o soluciones muy simples para problemas muy
generales. En general es imposible encontrar soluciones explı́citas para todas las soluciones de
todas las EDPs. Con mucha frecuencia las EDPs surgen al confeccionar modelos de un deter-
minado fenómeno, usando diferentes hipótesis. Evidentemente nos interesan aquellos esquemas
que describen adecuadamente la realidad y permiten además efectuar predicciones. Si nuestro
modelo representa adecuadamente el fenómenos fı́sico que estamos estudiando, cabe esperar que
posea una solución única. A priori no es evidente si un problema de ecuaciones en derivadas
parciales está bien planteado aunque fı́sicamente parezca razonable que la solución haya de ser
única, puede ocurrir que no tenga solución o que, de tener, no sea única. Aunque estemos segu-
ros de que la solución existe y es única, esto no basta para afirmar que el problema está bien
planteado.
En efecto, en un problema realista los datos se obtienen a partir de experimentos, ası́ que
forzosamente no serán exactos, de manera que, para que el problema tenga sentido, ha de ocurrir
también que un pequeño cambio en esos datos no influya de manera notable en la solución. Pero
aún hay más: como en la mayor parte de la ocasiones estas ecuaciones tienen que resolverse
numéricamente, lo que siempre implica aproximaciones y errores de redondeo, es preciso que
estas aproximaciones no produzcan grandes desviaciones del resultado exacto, ya que en caso
contrario el modelo usado serı́a de muy poca utilidad. El estudio de los problemas bien planteados
es lo que se suele denominar teorı́a clásica de las ecuaciones en derivadas parciales, (esto no
significa que los problemas mal planteado sean irrelevantes, ya que, por ejemplo, juegan un
papel destacado en mecánica de fluidos). En ese sentido no son suficientes las soluciones clásicas
y deberemos usar las soluciones débiles y generalizadas, que fueron introducidas por Sobolev
alrededor de 1930; se trata de soluciones que pueden no ser suficientemente regulares o incluso

12
tratarse de distribuciones (o funciones generalizadas). Algunas ecuaciones, por ejemplo la de
Laplace, admiten estas soluciones clásicas. Otras ecuaciones, sin embargo, no poseen en general
este tipo de soluciones clásicas. Este es el caso la ley de conservación escalar:
ut + ∂x F (u) = 0,
que rige diversos fenómenos unidimensionales en dinámica de fluidos, y que en particular sirve
para elaborar modelos sobre la formación y propagación de ondas de choque (una onda de choque
es precisamente una curva en la cual la solución u(x, t) presenta una discontinuidad).
Concretando un poco más, diremos que un problema para una cierta EDP está bien planteado
si se cumplen las tres condiciones siguientes:
1. El problema posee una solución (existencia).
2. La solución es única (unicidad).
3. La solución es estable, es decir, depende de forma continua de los datos del problema (de
las condiciones iniciales y de las de contorno) (estabilidad).
Las dos primeras exigencias son razonables, pero la tercera pudiera parecer caprichosa. Sin
embargo es especialmente importante para problemas que aparecen en aplicaciones fı́sicas, ya
que, es deseable que la solución que buscamos, además de ser única, cambie muy poco cuando se
modifican ligeramente las condiciones que especifican el problema. Algunos tipos de ecuaciones,
en particular las EDP de segundo orden, han sido estudiadas exhaustivamente y se conocen
resultados para saber qué tipos de problemas están bien planteados. Si una o más de las 3
condiciones listadas anteriormente no se cumple, decimos que el problema está mal planteado.
Por ejemplo, consideremos la ecuación no lineal:
u2x + u2y = −1,
dado que el lado izquierdo de la ecuación es no negativo,claramente esta ecuación no se satisface
para ninguna función real u(x, y). Podemos decir, con certeza que los problemas fundamentales
de la fı́sica matemática son todos bien planteado. Sin embargo, en ciertas aplicaciones de in-
genierı́a pueden aparecer problemas mal planteados. En la práctica, este tipo de problemas no
tienen solución. Por lo tanto, cuando nos enfrentamos a un problema mal planteado, el primer
paso deberı́a ser modificarlo adecuadamente con el fin de hacerla bien planteado.

1.8. Método de Separación de Variables

A pesar de que existen varios métodos que pueden utilizarse para encontrar soluciones par-
ticulares de una EDP lineal, con el método de separación de variables, el objetivo es encontrar
soluciones particulares en forma del producto de una función de x y de una función de y.
u(x, y) = F (x)G(y)
Luego, mediante esta suposición con frecuencia es factible reducir una EDP lineal de dos variables
en dos EDOs (Ecuaciones Diferenciales Ordinarias). Ası́ tenemos por ejemplo:
∂u du ∂u du ∂2u d2 u ∂2u d2 u
=G , =F , = G , = G , (21)
∂x dx ∂y dy ∂x2 dx2 ∂y 2 dy 2
Ejemplo:
Encontrar las soluciones producto de:
∂2u ∂u
2
=4
∂x ∂y

13
Solución:
Suponemos una solución de la forma: u(x, y) = F (x)G(y), calculamos sus derivadas según las
fórmulas (21) y sustituimos:
dF dG
G 2 =4
dx dy
separamos las variables y obtenemos
d2 F dG
2
=
4F dx Gdy
La expresión anterior nos indica que ambos miembros de la ecuación son independientes de x e y.
En otras palabras, cada miembro de la ecuación debe ser constante. Para fines prácticos, resulta
conveniente escribir dicha constante de separación como −λ. A partir de las dos igualdades:
d2 F dG
= = −λ
4F dx2 Gdy
obtenemos las dos ecuaciones diferenciales lineales ordinarias
d2 F dG
+ 4λF = 0 y + λG = 0
dx2 dy
que también podemos escribir como:

F ′′ + 4λF = 0, G′ + λG = 0

Haciendo un análisis, se obtiene que la ecuación auxiliar de la ecuación para F y la solución

para G son respectivamente:

m2 + 4λ = 0 y G = Ce−λy ,

luego las raı́ces de la ecuación auxiliar serán:

√
m = ±2 −λ

En ambos casos se observa que la solución depende del valor de λ. Para evitar tener que lidiar
con raı́ces inexactas, muchos autores, sobre todos los de libros de ingenierı́a, establecen que
λ = −α2 . Ası́ analizaremos tres casos: λ = 0, λ = −α2 < 0 y λ = −α2 > 0, veamos:
1. Caso 1: (λ = 0) Se tiene que m = 0. por lo que F = c1 + c2 x y G = c3 . Luego la solución
general en este caso será u = A1 + B1 x.

2. Caso 2: (λ = −α2 < 0). En este caso, las soluciones son reales y distintas puesto que
−λ > 0, ası́ que m = ±2α. Luego las soluciones de cada ecuación serán:
2y
F = c4 e−2αx + c5 e2αx , y G = c6 eα

Trabajando un poco en la solución para F se obtiene una combinación de funciones hi-

perbólicas, esto es:
F = c4 cosh 2αx + c5 sinh 2αx
Por lo tanto la solución general en este caso será:
2y 2y
u = A2 eα cosh 2αx + B2 eα sinh 2αx

donde A2 = c4 c6 y B2 = c5 c6 .

14
 Figura 1:

3. Caso 3: (λ = α2 > 0). En este caso las soluciones son imaginarias puesto que −λ < 0 y
m = ±2αi. Ası́ que:
F = c7 e−2αix + c8 e2αix ,
lo cual se puede escribir también como: F = c7 cos 2αx + c8 sin 2αx, si usamos la fórmula
2
de euler. Mientras que G = C9 e−α y .
Por lo tanto la solución general en este caso será:
2 2
u = A3 e−α y cos 2αx + B3 e−α y sin 2αx,

donde A3 = c7 c9 y B3 = c8 c9 .

1.9. Generalidades de las EDPs Clásicas de Segundo Orden.

Para finalizar este módulo, vamos a exponer brevemente, las carácterı́sticas más generales
de las EDPs clásicas, como lo son la ecuación del calor, la ecuación de ondas y la ecuación de
Laplace, mismas que abordaremos en forma más profunda en los siguientes módulos.

1.9.1. La Ecuación del Calor

La EDP
∂u ∂2u
=k 2 (22)
∂t ∂x
se conoce como Ecuación unidimensional del calor. Esta ecuación se presenta en la teorı́a
del flujo de calor, esto es, la transferencia de calor por conducción en una varilla o un alambre
delgado. La función u(x, t) es la temperatura.
Supongamos que una varilla circular delgada del longitud L tiene área transversal A y coin-
dice con el eje x en el intervalo [0, L](figura 1). Supongamos también que:
Dentro de la varilla, el flujo de calor tiene lugar sólo en la dirección x.

La superficie lateral, o curva, de la varilla está aislada; esto es no escapa calor de su

superficie.

No se está generando calor dentro de la varilla.

La varilla es homogénea, esto es, su masa por unidad de volumen ρ (densidad) es constante.

El calor especı́fico γ y la conductividad térmica K del material de la varilla son constantes.

K
Es muy comunmente establecer k = γρ . A este constante positiva se le llama difusividad
térmica.

15
 Figura 2:

1.9.2. La ecuación de onda

La EDP
∂2u 2
2∂ u
= a (23)
∂t2 ∂x2
se conoce como ecuación de onda. Para describir el fenómeno que describe la ecuación de onda
consideremos una cuerda de longitud L, como la cuerda de una guitarra, tensada entre dos
puntos localizados en el eje x, digamos, x = 0 y x = L. Cuando la cuerda comienza a vibrar,
supongamos que el movimiento tiene lugar en el plano xy de tal manera que cada punto de
la cuerda se mueve en dirección perpendicular al eje x (vibraciones trasversales) (figura 2). De
esta manera u(x, t) expresa el desplazamiento vertical de cualquier punto de la cuerda medido
a partir del eje x para t > 0. Además suponemos que:

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 en comparación con la longitud de la cuerda.

La pendiente de la curva es pequeña en todo momento.

La tensión T actúa en dirección tangente a la cuerda y su magnitud T es igual en todos

los puntos.

La tensión es grande en comparación con la fuerza de gravedad.

No actúan otras fuerzas externas sobre la cuerda.

T
El término a2 se define como a2 = ρ

1.9.3. La Ecuación de Laplace

La ecuación de Laplace tiene la forma:

∂2u ∂2u
+ 2 =0
∂x2 ∂y

la misma se puede abreviar como ∇2 u = 0, donde

∂2u ∂2u
∇2 u = + 2
∂x2 ∂y

16
 Figura 3:

se llama laplaciano bidimensional de la función u. En tres dimensiones, el laplaciano de u es

∂2u ∂2u ∂2u

∇2 u = + 2 + 2
∂x2 ∂y ∂z
La ecuación de Laplace en dos y tres dimensiones se presenta en problemas independientes
del tiempo que involucran potenciales como el electrostático, el gravitacional y la velocidad en
mecánica de fluidos. Además, la solución de la ecuación de Laplace también puede interpretarse
como la distribución de temperatura de estado estable (figura 3).

17
2. ECUACIONES DIFERENCIALES PARCIALES DE PRI-
MER ORDEN
Definición: Sea Ω un dominio sobre Rn . Una ecuación diferencial de primer orden en Ω
está dada por:
F (x, u, ∇u(x)) = 0, x ∈ Ω, (24)
∂u ∂u ∂u
donde F es continua en (x, u, p) ∈ Ω×R×Rn , u es una función C 1 (Ω) y ∇u = ∂x 1
+ ∂x 2
+...+ ∂x n
.
Las ecuaciones de primer orden aparecen en una variedad de procesos fı́sicos y de ingenierı́a,
tales como el transporte de material en el flujo de un fluido y la propagación de frentes de
ondas en óptica. Sin embargo, son menos frecuentes que las ecuaciones de segundo orden. Por
simplicidad limitaremos el estudio a las ecuaciones en dos variables, sobre la cuales se puede
establecer una clasificación atendiendo al grado de linealidad.
De acuerdo con (24), la ecuación general de primer orden para la función u(x, y) tiene la
forma:
F (x, y, u, ux , uy ) = 0 (25)

2.1. Clasificación de las EDPs de primer orden

El grado de linealidad que tiene una EDP de primer orden se puede usar para establecer una
clasificación de las mismas.
Ecuaciones lineales: una EDP lineal de primer orden para la función incógnita u(x, y) se
puede escribir como:

p(x, y)ux + q(x, y)uy = r(x, y)u + s(x, y) (26)

donde p(x, y), q(x, y), r(x, y) y s(x, y) son funciones arbitrarias de x e y y s(x, y) es llamado
el término no homogéneo. Por ejemplo, una ecuación lineal tı́pica de primer orden es:

ux + xuy = u + 2

EDPs de primer orden semilineal: una EDP semilineal de primer orden tiene la forma:

p(x, y)ux + q(x, y)uy = f (x, y, u) (27)

donde las funciones p y q sólo depende de x e y y f es una función no lineal arbitraria de

u. Un tı́pico ejemplo de una EDP semilineal de primer orden es:

ux + (1 + x)uy = (1 + x + y)u2

donde la función f (x, y, u) = (1 + x + y)u2 es no lineal debido al término u2 .

EDPs de primer orden cuasilineal: Una EDP cuasilineal de primer orden tiene la siguiente
ecuación:
p(x, y, u)ux + q(x, y, u)uy = f (x, y, u), (28)
donde p y q pueden depender o no depender de x e y y por lo menos uno de ellos depende
de u. Cuando f está presente, pueda que dependa o no de x, y y u, aunque esto no altera
la naturaleza de la ecuación. Una EDP tı́pica cuasilineal es:

ux + uuy = u,

en este caso la cuasilinearidad se debe a la presencia del término uuy .

18
 Tanto las EDPs lineales como las no lineales de primer orden aparecen a menudo en sistemas
que involucran varias variables dependientes, y en muchas ocasiones es posible, no para todas,
pero si para algunas, eliminar alguna de las variables dependientes, obteniendose una única
ecuación de orden superior dependiente del resto de las variables dependientes.
Aunque en principio se puede encontrar una solución general de una EDP de primer orden,
a diferencia de la solución general de una ODE lineal que contiene una constante arbitraria, la
solución general de una EDP lineal contiene una función arbitraria. Esto lo podemos ver en el
siguiente ejemplo:
ux + xuy = u + 2, (29)
para la que podemos probar que tiene la solución general
u(x, y) = C exp(x + φ(ξ)) − 2, (30)
donde ξ 2 = x2 − 2y, φ es una función diferenciable con argumento ξ y C una constante.
Para encontrar soluciones particulares supongamos, por ejemplo, que una solución de la
ecuación anterior satisface la condición de contorno u(x, 0) = −1. Haciendo y = 0 en la solución
general, y reemplazando el valor de ξ 2 se sigue que sobre el eje x − ξ = x, podemos encontrar
de la condición u(x, 0) = −1 que la función arbitraria φ se debe escoger de tal modo que:
−1 = C exp[x + φ(x)] − 2 y 1 = C exp x + φ(x)

2.2. Método de las Caracterı́sticas

Usando el método de las caracterı́sticas, podemos resolver problemas de valor inicial (PVI)
para EDPs de primer orden lineales, cuasilineales y no lineales en general, discutiendo el com-
portamiento de la solución en algunos ejemplos concretos. Este método está basado en una
interpretación geométrica de la EDP que conduce a un sistema de ecuaciones diferenciales ordi-
narias (EDO).
La forma general de la EDP lineal con coeficientes constantes es:
a(x, y)ux + b(x, y)uy = c(x, y)u + d(x, y) (31)
Supongamos que estamos localizados en el punto P (x, y, u) de la superficie solución u = u(x, y)
y nos movemos en dirección del vector (a, b, c). Para el caso en que la ecuación es homogénea
se tiene que sobre cualquier punto de la superficie, la dirección de la normal está dada por el
vector (ux , uy , −1). Entonces se tiene que:
(a, b, c) · (ux .uy , −1) = 0 por ser perpendiculares
Ası́, (a, b, c) es perpendicular a la normal y debe encontrarse en el plano tangente a la
superficie u = u(x, y). Ası́ que cualquier superficie solución de la EDP de primer orden pasa a
través del punto P (x, y, u) debe ser tangente al vector (a, b, c). Además, como (a, b, c) es siempre
tangente a superficie, nunca nos salimos de la superficie.
Si a, b, c, y d son funciones suaves (o regulares), la solución u de (31) es una superficie en el
espacio (x, y, u) que podemos parametrizar utilizando dos parámetros, digamos s y ξ. Para hallar
la expresión de u = u(x, y) suponemos dada una curva inicial C de parámetro ξ , que será el dato
inicial o condición inicial (CI) que acompaña a (31). Si por cada punto de C podemos trazar
otra curva de parámetro s, entonces la superficie ası́ determinada será la solución buscada. El
lado izquierdo de (31) tiene la forma de una derivada direccional de u(x, t) en la dirección del
vector [a, b] en los puntos (x, y) donde [a, b] está definido y no es nulo. Luego las ecuaciones
dx dy
= a(x, y); = b(x, y), s = parámetro, (32)
ds ds

19
 Figura 4:

determinan una familia uniparamétrica de curvas en el plano (x, y):

x = x(s); y = y(s) (33)

cuyos vectores tangentes [x′ (s), y ′ (s)] tienen la misma dirección de [a, b] (figura 4). Por lo tanto
la derivada total de u sobre las curvas (33) será:
du ∂u dx ∂u dy
= + (34)
ds ∂x ds ∂y ds
Comparando esta derivada direccional con la ecuaciones (31 y 32) se tiene que:
du ∂u ∂u
=a +b = cu + d
ds ∂x ∂y
obteniéndose el siguiente sistema de EDOs
dx dy du
= a; = b; = cu + d (35)
ds ds ds
Despejando el difencial del parámetro (ds) e igualando las tres ecuaciones anteriores se obtiene:
dx dy du
= = (36)
a b cu + d
que se conocen como las Ecuaciones caracterı́sticas auxiliares. Todo esto conduce al siguiente
teorema:

Teorema:
Dada la ecuación general de la EDP lineal de primer orden

aux + buy = cu + d (37)

Las superficies ortogonales a la solución, son las generadas por las curvas de la ecuación auxiliar
dx dy du
= = (38)
a b cu + d

20
La fórmula anterior se conoce como ecuación caracterı́stica auxiliar de la EDP lineal de primer
orden y es la que utilizaremos para resolver una EDP de primer orden con coeficientes constan-
tes. Es decir, resolviendo el par de ecuaciones resultantes se obtiene la solución general de la
EDP. Sin embargo, no estamos interesados en la solución general, sino en hallar una solución
particular u = u(x, y) que será una unica superficien en el espacio y que contenga una curva
dada C. Esta curva será la CI para (31). La EDP y la CI constituyen un problema de valor
inicial (PVI) o más especificamente un Problema de Cauchy.

Ejemplos:

1. Consideremos la EDP lineal

aux + buy = f (x, y)
De acuerdo con (38), la ecuación auxiliar correspondiente será:

dx dy du
= =
a b f (x, y)

Resolviendo el primer sistema obtenemos que la solución es: bx−ay = ξ. Luego despejando
x se obtiene x = ay+ξ
b :
Luego la otra ecuación nos queda:
dy du
=  ,
b f ay+ξ
, y
b

Resolviendo se obtiene: du = F (y, c)dy. Cuya solución es u = G(y, ξ) + c1 . Ası́, la solución

general se obtiene reemplazando c1 por φ(ξ) y ξ por bx−cy, por lo tanto la solución general
tiene la forma:
u(x, y) = G(y, bx − ay) + φ(bx − ay)

2. La EDP lineal
aux + buy + ku = f (x, y) (39)
tiene como ecuación auxiliar
dx dy du
= = ,
a b f − ku
Luego la solución del par izquierdo es bx − ay = ξ. Mientras que el otro par
dx du
=
a f − ku
puede reducirse a una EDO lineal con u como variable dependiente y x como variable
independiente. Esta ecuación esta dada por

du ku f (x, y)
+ =
dx a a

que se puede resolver mediante el método del factor integrante ekx/a o reducirla a una
ecuación homogénea como la del caso anterior.

21
3. Resolver el problema de Cauchy

ux + 3uy = 2u sujeto a u(x, 0) = ex

Solución:
Esta es una ecuación lineal y como el problema de Cauchy es una curva en el eje x usaremos
las ecuaciones caracterı́sticas. Esto es:
dx dy du
= =
1 3 2u
La solución del par izquierdo es: x = 13 y + ξ ó y = 3x + ξ.

Mientras que el otro par es:

du
= dx
2u
y su solución es: ln u = 2x + f (ξ), o también u = e2x+f (ξ) . Reemplanzando el valor del
parámetro ξ se obtiene:
u = e2x+f (y−3x)

Para encontrar el valor de la función f (y − 3x) reemplazamos la condición de Cauchy

u(x, 0) = ex , esto es:
u(x, 0) = e2x+f (−3x) = ex ,
de donde se obtiene que f (−3x) = −x ó lo que es lo mismo f (x) = 13 x. Luego f (y − 3x) =
1
3 (y − 3x). Por lo tanto la solución del problema de Cauchy será:

u = ex e1/3x

Le queda al lector comprobar que efectivamente esta es la solución del EDP propuesta.

4. Resolver el problema de Cauchy ux + uy = eu , sujeta a u(0, y) = y

Solución: Note que esta es una ecuación semilineal debido al término no homogéneo eu ,
pero en este caso la condición de Cauchy se especifica sobre el eje y. Sin embargo, como
veremos, también podemos usar el método de las caracterı́sticas.
De esta manera las ecuaciones caracterı́sticas serán:
du
dx = dy =
eu

Al igual que en los casos anteriores, la solución del par izquierdo es x + ξ = y. Mientras
que el otro par tiene como solución e−u = −x + f (ξ) o también
 
1
u = ln
f (ξ) − x

Despejando el valor del parámetro ξ y reemplazandolo en la ecuación para u se obtiene:

 
1
u = ln
f (y − x) − x

22
 Para encontrar el valor de la función f aplicamos la condición de Cauchy u(0, y) = y, esto
es:  
1
u(0, y) = ln =y
f (y)
de donde se deduce que f (y) = e−y . Por lo tanto la solución del problema de Cauchy será:
 
1
u = ln x−y para ex−y > x.
e −x
Como podrán comprobar, esta expresión satisface la condición de Cauchy especificada y la
ecuación de la EDP original, ası́ que esta es una solución clásica. La restricción ex−y > 0,
es para asegurar que u(x, y) es real, lo cual significa que la solución no está definida sobre
todo el plano (x, y) y por lo tanto no es una solución global.

2.3. Ecuaciones No lineales

Consideremos el siguiente problema de Cauchy:
ux + uuy + u = 0, sujeta au(0, y) = 1 + y
Como podemos observar esta ecuación es cuasilineal debido a la presencia del término uuy . Si
embargo aún es posible aplicar el método de las caracterı́sticas como sigue:
dx dy du
= =
1 u −u
El primer par de igualdades produce la EDO
dy
=u
dx
La cual no podemos resolver hasta encontrar una expresión para u. De manera que resolvemos
primero el segundo par de igualdades, esto es:
du
dx =
−u
La solución en este caso es: ln u = −x + ln g(ξ), ası́ que u = g(ξ)e−x .
Se sigue de la condición de Cauchy que u = 1 + ξ en el punto (0, ξ), ası́ haciendo x = 0 y
reemplazando u por 1 + ξ en el último resultado, encontramos que
g(ξ) = 1 + ξ
y ası́ se tiene que
u = (1 + ξ)e−x
Ahora utilizamos este resultado para resolver el primer par de ecuaciones.
dy
= (1 + ξ)e−x
dx
Integramos esta ecuación usando el hecho de que la caracterı́stica pasa a través del punto (0, ξ),
y se obtiene: Z x
y = ξ + (1 + ξ) e−η dη
0
que da como resultado
y = (1 − e−x ) + ξ(2 − e−x )
Cuando elemininamos ξ la solución se transforma en
 
1+y
u= e−x , con la condición x 6= − ln 2
2 − e−x

23
2.4. Movimiento de ondas unidireccional.
La EDP de primer orden
∂u ∂u
+c = 0, c = cte, x ∈ R, t > 0 (40)
∂t ∂x
representa un movimiento de ondas unidireccional en R (conocida también como ecuación de
convección). Resolver un Problema de Valor Inicial (PVI) para (40) con la condición inicial

u(x, 0) = F (x), x ∈ R, (41)

siendo F una función dada.

Solución:
Aplicando el método de las caracterı́sticas, se tiene que las ecuaciones auxiliares son:

dx
dt = , du = 0
c
que, junto a las CI conforman los tres PVI para estas EDOs, cuyas soluciones son

dx
dt = , du = 0
c
Cuyas soluciones son x − ct = ξ y u = f (ξ). Reemplazando el valor del parámetro ξ en la
ecuación para u, se obtiene la solución general

u(x, t) = f (x − ct)

Finalmente, aplicando la condición inicial u(x, 0) = F (x), se obtiene la solución paticular

u(x, t) = F (x − ct) (42)

Si F ∈ C 1 (R) entonces (42) es la solución del PVI. Nótese que la solución está determinada sólo
en términos del dato inicial y .

Interpretación fı́sica de la solución

Si F representa la forma inicial de una onda, entonces u(x, t) = F (x − ct) muestra que la
onda no cambia de forma cuando t crece (e.d. a medida que el tiempo transcurre). Además, si
c > 0, entonces la onda viaja a la derecha con velocidad dx/dt = c; y si c < 0 la onda viaja
a la izquierda. Las CC son x − ct =. De este modo, para cada valor de , u es constante sobre
las curvas caracterı́sticas. Esto significa que el dato inicial se transmite a lo largo de las curva
caracterı́stica, en el sentido que cualquiera sea el valor de u a t = 0 en algún x, ese valor es el
mismo en todos los puntos x pertenecientes a las curvas caracterı́sticas (ver figura 5).

24
 Figura 5: Propagación unidireccional de ondas

Tarea No 2
1. Clasifique las siguientes EDPs en lineales, semilineales, cuasilineales o no lineales.
a) ux + u2 uy = x + 2y f) ux + (1 + ux )uy = u2
b) 3ux + 4uy = sin x g) ux sin y + uy cos x = 1 + x2 + y
2
c) ux + xuy = u + 1 h) ux + (1 + u)uy = 2xy
d) ux + 2uy = cos u i) (x2 + 1)ux + u2y = 2x + 3
e) (x + 1)ux + yuy = 2u + e x j) (xy + 2)ux + (1 + y + u)uy = u
2. En los siguientes ejercicios use la solución general de la EDP (29) dada por (30) para
encontrar la solución que satisfaga la condición dada, imponiendo cualquier restricción
que se requiera para que la solución sea válida.
a) u(x, 0) = 2, y > 0
b) u(x, 1) = −1, y > 0
c) u(x, 2) = x − 2, y > 2
3. En los siguientes ejercicios resuelva el problema de Cauchy. Verifique que el resultado
obtenido es una solución y explique si existe alguna restricción que sea necesario imponerle
a esta.
a) ux + 2uy = 2, u(x, 0) = x
b) 3ux + 2uy = x, u(x, 0) = 1
c) yux + 3uy = y(1 + u), u(0, y) = y 2
d) 2ux + uy = cos x, u(x, 0) = 2x
e) ux + xuy = u + 2, u(0, y) = 3y
f) ux + uuy + u = 0, u(0, y) = sin y
4. Resuelva los siguientes problemas de valor inicial.
a) 2ut + 4ux = 3u, dado que u(x, 0) = sin 2x
b) ut − 2ux = x, dado que u(x, 0) = x2
c) ut − 3ux = 2u + 1, dado que u(x, 0) = 21 cos x
d) ut − 4ux = 3x, dado que u(x, 0) = ex

25
3. ELEMENTOS DE LAS SERIES, INTEGRALES y TRANS-
FORMADA DE FOURIER
3.1. Funciones Ortogonales
3.1.1. Producto Interno
Si ~u y ~v son dos vectores en el espacio tridimensional, entonces el producto interno < ~u, ~v >
o ~u · ~v posee las siguientes propiedades

1. < ~u, ~v > = < ~v , ~u >

2. < k~u, ~v >= k < ~u, ~v >, k un escalar

3. < ~u, ~v >= 0 si u = 0 y < ~u, ~v > > 0 si ~u 6= 0

4. < ~u + ~v , w
~ > = < ~u, w
~ > + < ~v , w
~>

3.1.2. Producto interno de funciones

Sean f1 y f2 dos funciones definidas sobre un intervalo [a, b]. El producto interno de dos
funciones f1 y f2 se define como:
Z b
< f1 , f2 >= f1 (x)f2 (x)dx (43)
a

Definición: (Vectores Ortogonales:) Dos vectores ~u y ~v son ortogonales si:

< ~u, ~v >= 0 (44)

Definición: (Funciones Ortogonales:)

Dos funciones f1 y f2 son ortogonales sobre un intervalo [a, b] si:
Z b
< f1 , f2 >= f1 (x)f2 (x)dx = 0 (45)
a

Ejemplo:
Las funciones f1 (x) = x2 y f2 (x) = x3 son ortogonales en el intervalo [−1, 1] ya que:
Z 1
1 Z 1
2 3 x6  5
< f1 , f2 >= x · x dx = x dx = =0
−1 −1 6 −1

A diferencia del cálculo vectorial, donde la palabra ortogonal es sinónimo de perpendicular, en

este contexto el término ortogonal no tiene significado geométrico.

Definición: Conjunto ortogonal

Un conjunto de funciones de valor real {φ0 (x), φ1 (x), φ2 (x), ...} se dice que es ortogonal en un
intervalo [a, b] si
Z b
< φm , φn >= φm (x)φn (x)dx = 0, m 6= n (46)
a

26
 Dado un vector ~u, la norma o módulo k ~u k se puede expresar en términos del producto
interno, esto es:
< ~u, ~v >=k ~u k2 (norma cuadrada) (47)
Por lo tanto p
k ~u k= < ~u, ~v > (Norma) (48)
De igual modo, para un conjunto de funciones reales {φn (x)}n∈N , las fórmulas correspon-
dientes serán:
k φn (x) k2 =< φn , φn > (49)
p
k φn (x) k= < φn , φn > (50)
Por lo tanto, de la definición de función ortogonal se tiene que:
Z b
k φn (x) k2 = φ2n (x)dx (51)
sa
Z b
k φn (x) k= φ2n (x)dx (52)
a

Si {φn (x)} es un conjunto ortogonal de funciones en el intervalo [a, b] tal que:

k φn (x) k= 1, n = 0, 1, 2, ...
entonces, se dice que {φn (x)} es un conjunto ortonormal.
Ejemplo 1: Demuestre que el conjunto {1, cos x, ..., cos nx} es ortogonal en el intervalo [−π, π].
Solución:
Tenemos que probar que
Z π
< φn , φm >= φn (x)φm (x)dx = 0, m 6= n
−π
Caso 1: Z π Z π
π
1 
< φ0 , φn >= φ0 (x)φn (x)dx = cos xdx = sen nx = 0
−π −π n −π

Caso 2:
Z π Z π
< φm , φn >= φm (x)φn (x)dx = cos nx cos mxdx
−π −π
Z
1 π
= [cos(m + n)x + cos(m − n)x] dx
2 −π
 
1 sen(m + n)x sen(m − n)x π
= + =0
2 m+n m−n π

Ejemplo 2: Encuentre la norma de cada una de las funciones del ejemplo 1.

Solución:
Para φ0 (x) = 1 se tiene:
p
k φ0 (x) k = < φ0 , φ0 >
sZ
π
= φ20 (x)dx
−π
sZ
π √
= dx = 2π
−π

27
Para φn (x) = cos nx
p
k φn (x) k = φ2
sZ n
π
= (cos nx)2 dx
−π
r
1 1
= [x + sen 2x]π−π dx
2 2n
√
= π

Conjunto de funciones ortonormales

Cualquier conjunto ortogonal de funciones diferente de cero {φn (x)}, n = 0, 1, 2, ..., se puede
normalizar, dividiendo cada función entre su norma. Esto es, del ejemplo anterior, se tiene que:
1 cos x cos nx
{ √ , √ , √ , ...}
2π π π
es ortogonal en [−π, π].

Supongamos que ~v1 , ~v2 y ~v3 son tres vectores distintos de cero, ortogonales entre si en el
espacio vectorial R3 . Entonces, dichos vectores pueden usarse como una base para dicho espacio,
esto es:
~u = c1~v1 + c2~v2 + c3~v3 (53)
donde c1 , c2 y c3 son escalares. Cada ci se puede expresar en terminos de ~u y del vector ~vi
correspondiente. Es decir, por ejemplo para c1 se tiene:

< ~u, ~v1 > = < c1~v1 , ~v1 > + < c2~v2 , ~v1 > + < c3~v3 , ~v1 >
= c1 < ~v1 , ~v1 > +c2 < ~v2 , ~v1 > +c3 < ~v3 , ~v1 >
= c1 k~v1 k2

Luego se tiene que

< ~u, ~v1 >
c1 =
k~v1 k2
De igual forma, podemos hacer lo mismo para c2 y c3 , esto es:
< ~u, ~v2 > < ~u, ~v3 >
c2 = , c3 =
k~v2 k2 k~v3 k2
Por lo tanto de (53) se tiene que:
3
< ~u, ~v1 > < ~u, ~v2 > < ~u, ~v3 > X < ~u, ~vn >
~u = 2
~v1 + 2
~v2 + 2
~v3 = ~vn (54)
k~v1 k k~v2 k k~v3 k k~vn k2
n=1

3.1.3. Desarrollo en series ortogonales

Siguiendo un proceso analógo al que usamos en el caso anterior,dado un conjunto ortogonal
de funciones {φn (x)}, pero esta vez infinito sobre un intervalo [a, b], si f (x) es una función
definida en [a, b], es posible encontrar coeficientes cn , n = 0, 1, 2, ... tal que:

f (x) = c0 φ0 (x) + c1 φ1 (x) + ... + cn φn (x) + ... (55)

28
Esto es, multiplicando esta ecuación por φm (x) e integrando en el intervalo [a, b] se tiene:
Z b Z b Z b Z b
φm (x)f (x)dx = c0 φm (x)φ0 (x)dx + c1 φm (x)φ1 (x)dx + ..., cn φm (x)φn (x)dx + ...
a a a a
= c0 < φm , φ0 > +c1 < φm , φ1 > +..., cn < φm , φn > +...
Cada término del lado derecho es cero, excepto cuando m = n debido a la ortogonalidad de las
funciones φn , por lo tanto se tiene:
Z b Z b
φm (x)f (x)dx = cn φ2n (x)dx
a a
Por lo tanto Rb
a φm (x)f (x)dx
cn = Rb
2
a φn (x)dx
Es decir de (55) se tiene
∞ ∞
X X < f (x), φn (x) >
f (x) = cn φn (x) = φn (x) (56)
n=0 n=0
kφn (x)k

Definición: Conjunto ortogonal con respecto a una función de peso

Se dice que un conjunto de funciones de valor real {φ0 (x), φ1 (x), φ2 (x), ...} es ortogonal respecto
a una función de peso w(x) en un intervalo [a, b] si
Z b
w(x)φm (x)φn (x)dx = 0, m 6= n. (57)
a
Ejemplo: El conjunto {1, cos x, cos 2x, ...., cos nx} es ortogonal respecto a la función de peso
w(x) = 1 en [−π, π].
Si {φn (x)} es ortogonal respecto a una función de peso w(x) en [a, b], entonces desarrollando una
función f (x) como en (55) y multiplicando por w(x)φn (x) e integrando sobre [a, b] se obtiene
Z b Z b Z b
w(x)φm (x)f (x)dx = c0 w(x)φm (x)φ0 (x)dx + c1 w(x)φm (x)φ1 (x)dx + ...
a a a
Z b
+ cn w(x)φm (x)φn (x)dx + ...
a
= c0 < w(x)φm , φ0 > +c1 < w(x)φm , φ1 > +..., cn < w(x)φm , φn > +...
Luego como {φn (x)} es ortogonal respecto del peso w(x) se tiene que < w(x)φn , φm >= 0
excepto si m = n.
Luego Z b Z b
w(x)φn (x)f (x)dx = cn w(x)φ2n (x)dx
a a
Ası́ que
Rb
w(x)φn (x)f (x)dx
a
cn =
Rb
2
a w(x)φn (x)dx
Conjuntos completos El procedimiento bosquejado para determinar los coeficientes cn
fue formal; esto es, las preguntas básicas acerca de que si un desarrollo ortogonal de una serie
como la (56) es en realidad factible o pudiera ser ignorada. Asimismo, para desarrollar f en una
serie de funciones ortogonales, desde luego es necesario que f no sea ortogonal a cada φn del
conjunto ortogonal {φn (x)}. Para evitar este último problema debemos suponer, en lo que resta
del análisis, que un conjunto ortogonal es completo. Esto significa que la única función continua
ortogonal a cada miembro del conjunto es la función cero.

29
3.2. Series de Fourier
En la sección anterior estudiamos que si φ0 (x), φ1 (x), φ2 (x), ... es un conjunto de funcio-
nes con valores reales que son ortogonales en el intervalo [a, b] y si f es una función defini-
da en el mismo intervalo, entonces podemos desarrollar formalmente f en una serie ortogonal
c0 φ0 (x) + c1 φ1 (x) + c2 φ2 (x) + .... En esta sección desarrollaremos las funciones en términos de
un conjunto ortogonal especial de funciones trigonométricas.

Series trigonométricas El conjunto de funciones trigonométricas

 
π 2π 3π π 2π 3π
1, cos x, cos x, cos x, ..., sin x, sin x, sin x, ... (58)
p p p p p p

es ortogonal en el intervalo [−p, p]. Este conjunto será de especial importancia posteriormente
en la solución de ciertos tipos de problemas con valores en el lı́mite que involucran ecuaciones
diferenciales lineales parciales. En esas aplicaciones, necesitaremos desarrollar una función f
definida sobre [−p, p] en una serie ortogonal que consista en las funciones trigonométricas dadas
en (58), es decir,
∞  
a0 X nπ nπ
f (x) = + an cos x + bn sin x (59)
2 n=1
p p
Los coeficientes a0 , a1 , a2 , ..., b1 , b2 , ..., pueden determinarse exactamente en la misma forma que
en el análisis general de las expansiones de series ortogonales de la sección anterior.

Definición: Series de Fourier

La serie de Fourier de una función f definida en el intervalo (−p, p) está dada por
∞  
a0 X nπ nπ
f (x) = + an cos x + bn sin x (60)
2 p p
n=1

donde
Z
1 p
a0 = f (x)dx (61)
p −p
Z
1 p nπ
an = f (x) cos dx (62)
p −p p
Z
1 p nπ
bn = f (x) sin dx (63)
p −p p

Ejemplo: Expanda la función


0, −π < x < 0
f (x) =
π − x, 0 ≤ x ≤ π

en una serie de Fourier

Solución:

30
Calculamos los coeficientes de acuerdo con la definición anterior.
Z Z π 
1 π 1 π
a0 = f (x)dx = (π − x)dx =
π −π π 0 2
Z π Z π 
1 1
an = f (x) cos nxdx = (π − x) cos nxdx
π −π π 0
− cos nπ + 1
=
n2 π
1 − (−1)n
=
n2 π
De manera similar, calculamor el valor de bn
Z
1 π 1
bn = f (x) sin nxdx =
π −π n
Por lo tanto,
∞  
π X 1 − (−1)n 1
f (x) = + cos nx + sin nx
4 n=1 n2 π n

El teorema siguiente proporciona condiciones suficientes para la convergencia de una serie de

Fourier en un punto.

Teorema: Condiciones de convergencia para una serie de Fourier

Sean f y f ′ funciones continuas en el intervalo (−p, p); esto es, establezcamos f y f ′ continuas
excepto en un número finito de puntos en el intervalo y con discontinuidades finitas sólo en
estos puntos. Entonces, la serie de Fourier de f en el intervalo converge a f (x) en un punto de
continuidad. En un punto de discontinuidad, la serie de Fourier converge al promedio
f (x+) + f (x−)
2
donde f (x+) y f (x−) denotan el lı́mite de f en x por la derecha a izquierda respectivamente.

Ejemplo: Convergencia de un punto de discontinuidad

La función del ejemplo anterior satisface las condiciones de convergencia del teorema de conver-
gencia. En consecuencia, por cada x en el intervalo (−π, π), excepto en x = 0, la serie de Fourier
covergerá f (x). En x = 0 la función es discotnı́nua, entonces la serie converge a
f (0+) + f (0−) π+0 π
= =
2 2 2
Extensión periódica
Una serie de Fourier no sólo representa la función en el intervalo (−p, p), sino que también
proporciona la extensión periódica de f fuera de este intervalo. Ahora podemos aplicar el teo-
rema anterior a la extensión periódica de f , o suponer desde el principio que la función dada es
periódica con periodo T = 2p; esto es, f (x + T ) = f (x). Cuando f es una función continua y
existen las derivadas derecha e izquierda en x = −p y x = p, respectivamente, entonces la serie
converge al promedio [f (p−) + f (p+)]/2

Secuencia de sumas parciales

Es interesante observar cómo la secuencia de las sumas parciales {SN (x)} de una serie de Fourier
se aproxima a una función. Por ejemplo, las primeras tres sumas parciales de (3.2) son
π π 2 π 2 1
S1 (x) = , S2 (x) = + cos x + sin x, S3 (x) = + cos x + sin x + sin 2x
4 4 π 4 π 2

31
3.2.1. Serie de Fourier de cosenos y senos
El esfuerzo que se lleva a cabo en la evaluación de los coeficientes a0 , an y bn al desarrollar
una función f en una serie de Fourier se reduce de manera significativa cuando f es una función
par o impar. Se dice que una función f es:

Par si f (−x) = f (x) e impar si f (−x) = −f (x)

En un intervalo simétrico tal como (−p, p), la gráfica de una función par tiene simetrı́a respecto
al eje y, mientras que la gráfica de una función impar tiene simetrı́a en relación con el origen.

Ejemplos
Las funciones polinomiales que consisten en todas las potencias pares de x son pares,
mientras que las gráficas de los polinomios constituidos por todas las potencias impares
de x son simétricas con relación al orı́gen.

Las funciones trigonométricas coseno y seno son funciones pares e impares respectivamente
ya que:
cos(−x) = cos x y sen(−x) = − sen x

Las funciones exponenciales f (x) = ex y f (x) = e−x no son pares ni impares.

Teorema: Propiedades de las funciones pares e impares

1. El producto de dos funciones pares es par.

2. El producto de dos funciones impares es par.

3. El producto de una función par y una impar es impar.

4. Al suma (resta) de dos funciones pares es par.

Ra Ra
5. Si f es par, entonces −a f (x) = 2 0 dx
Ra
6. Si f es impar, entonces −a f (x) = 0

Si f es una función par de (−p, p), entonces, en vista de las propiedades anteriores, los coeficientes
(61), (62) y (63) se convierten en
Z Z
1 p 2 p
a0 = f (x)dx = f (x)dx
p −p p 0
Z Z
1 p nπ 2 p nπ
an = f (x) cos x dx = f (x) cos xdx
p −p p p 0 p
| {z }
par
Z p
1 nπ
an = f (x) sen x dx = 0
p −p | p
{z }
impar

De manera similar, cuando f es impar en el intervalo (−p, p), se tiene:

Z
2 p nπ
an = 0, n = 0, 1, 2, ..., bn = f (x) sen xdx
p 0 p

32
Ejemplo:
Expanda f (x) = x, −2 < x < 2 en una serie de Fourier.
Solución:
Dado que la función f (x) = x es impar en el intervalo (−2, 2), desarrollamos f en una serie de
senos con p = 2. Por lo tanto a0 = an = 0 y
Z 2
nπ 4(−1)n+1
bn = sen xdx =
0 2 nπ

Por lo tanto,
∞
4 X 4(−1)n+1 nπ
f (x) = sen x
π nπ 2
n=1

Ejemplo: La función 
−1, −π < x < 0
f (x) =
1, 0 ≤ x ≤ π
es impar en el intervalo (−π, π). Con el valor de p = π tenemos que:
Z
2 π 2 1 − (−1)n
bn = (1) sen nxdx = ,
π 0 π n

y ası́
∞
2 X 1 − (−1)n
f (x) = sen nx (64)
π n
n=1

Fenómeno de Gibbs Con ayuda de un sistema asistido por computadora, en la figura (6)
se han trazado las gráficas S1 (x), S2 (x), S3 (x), S15 (x) de las sumas parciales de los términos
diferentes a cero de (64). Como se puede observar en la figura (6 d), la gráfica de S1 5(x) tiene
picos pronunciados cerca de las discontinuidades en x = 0, x = π, x = −π, , etc. Este ”disparo”
por las sumas parciales SN de los valores funcionales cerca de un punto de discontinuidad no la
empareja sino que permanece constante, aun cuando el valor de N se considera elevado. Este
comportamiento de una serie de Fourier cerca de un punto en el cual f es discontinua se conoce
como Fenómeno de Gibbs.

3.3. Integrales de Fourier

En las secciones anteriores las series de Fourier se usaron para representar una función f
definida en un intervalo finito (−p, p). Cuando f y f ′ son continuas en dicho intervalo finito,
una serie de Fourier representa la función en el intervalo y converge hacia la extesión periódica
de f fuera del intervalo. De esta manera podrı́amos suponer que las series de Fourier solo
se asocian con funciones periódicas. Ahora procederemos a deducir, una forma de representar
ciertas funciones no periódicas que estén definidas en un intervalo infinito (−∞, ∞) o semiinfinito
(0, ∞).

3.3.1. Deducción de la integral de Fourier

Supongamos que una función f está definid en (−p, p). Si utilizamos las definiciones integrales
de los coeficientes (61),(62) y (63) en la expresión (60), entonces la serie de Fourier de f en el

33
 Figura 6: Sumas parciales de la serie de senos de (64) en el intervalo (−π, π)

intervalo es
Z p ∞ Z p  Z p  
1 1X nπ nπ nπ nπ
f (x) = f (t)dt + f (t) cos tdt cos x+ f (t) sen tdt sen x
2p −p p n=1 −p p p −p p p
(65)
Si establecemos αn = nπ/p, ∆α = αn+1 − αn = π/p, entonces (65) se convierte en
Z p  ∞ Z p  Z p  
1 1X
f (x) = f (t)dt ∆α+ f (t) cos αn tdt cos αn x + f (t) sen αn dt sen αn x ∆α
2π −p π n=1 −p −p
(66)
Expandiendo ahora el intervalo (−p, p), haciendo que p → P∞∞. Como p → ∞, entonces ∆α → 0,
y el lı́mite del segundo término tiene la forma: lı́m∆α→o R n=1 F (αn )∆α, la cual no es más que
∞
la definición de la integral de Riemann. Por lo tanto, si −∞ f (t)dt existe, el lı́mite del primer
término de (66) es cero y el lı́mite del a suma se conviert en
Z Z ∞  Z ∞  
1 ∞
f (x) = f (t) cos αtdt cos αx + f (t) sen αtdt sen αx dα (67)
π 0 −∞ −∞

Esta expresión se llama Integral de Fourier de f en (−∞, ∞)

Definición: Integral de Fourier

La Integral de Fourier de una función f definida en el intervalo (−∞.∞) está dada por
Z
1 ∞
f (x) = [A(α) cos αx + B(α) sen αx] dα (68)
π 0

34
donde
Z ∞
A(α) = f (x) cos αxdx (69)
Z−∞
∞
B(α) = f (x) sen αxdx (70)
−∞

Las condiciones suficientes en las que la integral de Fourier converge hacia f (x) son similares a,
pero ligeramente más estrictas que, las condiciones de a series de Fourier, pues se exige además
que f y f ′ sean
R ∞continuas en cada intervalo finito y f sea absolutamente integrable en (−∞, ∞),
es decir que −∞ |f (x)|dx converge.

Ejemplo:
Encuentre la representación de la integral de Fourier de la función

 0, x < 0
f (x) = 1, 0 ≤ x ≤ 2

0, x > 2.
Como podemos ver, esta función satisface las condición suficiente de convergencia de la integral
de Fourier, entonces se tiene que:
Z ∞
A(α) = f (x) cos αxdx
−∞
Z 0 Z 2 Z ∞
= f (x) cos αxdx + f (x) cos αxdx + f (x) cos αxdx
−∞ 0 2
Z 2
sen 2α
= cos αxdx =
0 α
Z ∞ Z 2
1 − cos 2α
B(α) = f (x) sen αxdx = sen αxdx =
−∞ 0 α

Sustituyendo los coeficientes en la fórmula de la integral de Fourier tenemos

Z     
1 ∞ sen 2α 1 − cos 2α
f (x) = cos αx + sen αx dα
π 0 α α

Integrales de seno y coseno Cuando f es una función par en el intervalo (−∞, ∞), enton-
ces f (x) cos αx también es par, mientras que f (x) sen αx es impar. Como consecuencia de las
propiedades de las funciones pares B(α) = 0, por lo que la integral de Fourier nos queda
Z
2 ∞
f (x) = A(α) cos αxdα
π
Z ∞0
donde A(α) = f (x) cos αxdx
0

De manera similar, cuando f es impar en (−∞, ∞), f (x) cos αx y f (x) sen αx son funciones
impares y pares, respectivamente. Por lo tanto A(α) = 0 y
Z
2 ∞
f (x) = B(α) sen αxdα
π 0
Z ∞
donde B(α) = f (x) sen αxdx
0

35
 Ejemplo: Encuentre la representación integral de Fourier de la función

1, |x| < a
f (x) =
0, |x| > a.
Solución: Dado que f es una función par, la podemos representar mediante la integral de cosenos
de Fourier, esto es:

Z ∞ Z a Z ∞
A(α) = f (x) cos αxdx = f (x) cos αxdx + f (x) cos αxdx
Z0 a 0 0
sen aα
= cos αxdx =
0 α
por lo que Z ∞
2 sen aα cos αx
f (x) = dα
π 0 α

3.4. La Transformada de Fourier

Una transformada integral es una transformación que a partir de funciones dadas produce
nuevas funciones que dependen de una variable diferente y aparecen en la forma de una in-
tegral. Estas transformaciones son de interés principalmente como herramientas para resolver
ecuaciones diferenciales ordinarias, ecuaciones diferenciales parciales y ecuaciones integrales. La
Transformada de Fourier es una de las más importantes en la ingenierı́a ya que las mismas son
de vital importancia en la resolución de ecuaciones diferenciales parciales.
El punto de partida para la deducción de la Transformada de Fourier es la forma compleja de
la representación de la integral de Fourier definida en (68) para una función f (x). Para derivar
esta representación en la que f (x) esta definida sobre el intervalo (−∞, ∞), sustituimos en (68),
las expresiones de A(α) y B(α) dadas por las fórmulas (69) y (70). De esta manera, se obtiene:
Z Z ∞ 
1 ∞
f (x) = f (u) (cos αu cos αx + sen αu sen αx) du dα
π 0 −∞
Z Z ∞ 
1 ∞
= f (u) cos α(u − x)du dα
π 0 −∞
Z Z ∞ 
1 ∞
= f (u) cos α(x − u)du dα
π 0 −∞

Como el integrando de la última integral es una función para para α, el intervalo de integración
con respecto a α puede doblarse y para compensar mediante la introducción de un factor de 21 ,
de manera que obtenemos:
Z ∞ Z ∞ 
1
f (x) = f (u) cos α(x − u)du dα
2π −∞ −∞
Para completar la expresión compleja, sumamos cero al miembro derecho en la forma
Z ∞ Z ∞ 
i
0= f (u) sen α(x − u)du dα
2π −∞ −∞
lo cual es cierto puesto que sen α(x− u) es una función impar para α. De manera que obtenemos:
Z ∞ Z ∞ 
1
f (x) = f (u) (cos α(x − u) + i sen α(x − u)) du dα
2π −∞ −∞

36
Usando ahora la fórmula de Euler eiθ = cos θ + i sen θ obtenemos la forma compleja de la integral
de Fourier.
Z ∞ Z ∞ 
1
f (x) = f (u) exp{iα(x − u)}du dα
2π −∞ −∞
Z ∞  Z ∞ 
1 1
= √ exp(iαx) √ f (u) exp(−iαu)du dα (71)
2π −∞ 2π −∞
Si ahora definimos a F(α) como la integral de la derecha de la expresión anterior y dado que u
es una variable muda la podemos cambiar nuevamente por x se tiene que:
Z ∞
1
F(α) = √ f (x) exp(−iαx)dx
2π −∞
Ası́ que (71) se convierte en:
Z ∞
1
f (x) = √ F(α) exp(iαx)dα (72)
2π −∞

Por lo tanto a F(α) se le llama Transformada de Fourier de f (x) y dado que la integral
(72) nos permite recuperar f (x) mediante F(α), a esta se le llama Transformada inversa de
Fourier. Al par de integrales F(α) y f (x) se les llama Par de Tansformadas de Fourier.

Ejemplo:
Encuentre la transformada de Fourier de:
 
1, |x| < a 1, 0<x<a 1
a) f (x) = b) g(x) = c) p(x) =
0, |x| > a. 0, en cualquier otro punto x2 + a2
Solución:

Z a  iα 
1 −iαx 1 e a − e−iαa
a) F(α) = √ e dx = √
2π −a α 2π i
r
2 sen αa
=
π α
Z a  
1 −iαx 1 1 − e−iαa
b) G(α) = √ e dx = √
2π 0 2π iα
Z ∞ −iαx Z ∞ Z ∞
1 e 1 cos αx i sen αx
c) P(α) = √ 2 + a2
dx = √ 2 + a2
dx − √ dx (73)
2π −∞ x 2π −∞ x 2π −∞ 2 + a2
x
Como el integrando de la segunda integral es impar, esta es igual a cero, por lo tanto:
Z ∞
1 cos αx
P(α) = √ dx = (74)
2π −∞ x2 + a2

3.4.1. Propiedades de la Transformada de Fourier

1. Si F1 (ω) = F[f1 (x)] y Si F2 (ω) = F[f2 (x)] son las transformadas de Fourier de dos
funciones f1 y f2 , respectivamente, entonces:

F[a1 f1 (x) + a2 f2 (x)] = a1 F[f1 (x)] + a2 F[f2 (x)], (75)

donde a1 y a2 son constantes arbitrarias

37
 1 ω


2. Sea a una constante real y F(ω) = F[f (x)], entonces F[af (x)] = |a| F a .

3. Si F[f (x)] = F(ω), entonces F[f (−x)] = F(−ω)

4. Si F(ω) = F[f (x)], entonces F[f (x − x0 )] = F(ω)e−iωx0

5. Si ω0 es una constante real y F(ω) = F[f (x)], entonces F[f (x)e−iω0 x ] = F(ω − ω0 )

6. Si F[f (x)] = F(ω) y f (x) → ±∞, entonces F[ dfdx

(x)
] = iωF(ω) = iωF[f (x)]

7. Si F1 (ω) = F[f1 (x)] y F2 (ω) = F[f2 (x)], entonces F[f1 (x) ∗ f2 (x)] = F1 (ω)F2 (ω), donde
Z ∞
f1 (x) ∗ f2 (x) = f1 (x)f2 (t − x)dx, (76)
∞

es la convolución de dos funciones f1 y f2 .

38
Tarea No 3

1. En los siguientes problemas, demuestre que las funciones dadas son ortogonales.

a) f1 (x) = ex , f2 (x) = xe−x − e−x ; [0, 2]

b) f1 (x) = cos x, f2 (x) = sen2 x; [0, π]
c) f1 (x) = x, f2 (x) = cos 2x; [−π/2, π/2]
d) f1 (x) = ex , f2 (x) = sen x; [π/4, 5π/4]

2. En los siguientes problemas, demuestre que cada conjunto de funciones es ortogonal en el

intervalo indicado. Encuentre la norma de cada función del conjunto.
n o
a) sen nπ
p x , n = 1, 2, 3, ...; [0, p]
n o
b) 1, cos nπ
p x , n = 1, 2, 3, ...; [0, p]
n o
c) 1, cos nπ
p x, sen mπ
p x , n = 1, 2, 3, ..., m = 1, 2, 3...; [0, p]

3. Sea {φn (x)} un conjunto ortogonal de funciones en [a, b] tal que φ0 (x) = 1 y φ1 (x) = x.
Rb
Demuestre que a (αx + β)φn (x)dx = 0 para n = 2, 3, ..., y para cualquier constante α y β

4. Sea {φn } un conjunto ortogonal de funciones en [a, b]. Demuestre que

kφm (x) + φn (x)k2 = kφm (x)k2 + kφn (x)k2 , m 6= n

5. Se dice que una función con valores reales es periódica con periodo T si f (x + T ) = f (x).
Por ejemplo, 4π es un periodo de sen x ya que sen(x+4π) = sen x. El valor más pequeño de
T para el que f (x + T ) = f (x) es válida se llama periodo fundamental de f . Por ejemplo,
el periodo fundamental de f (x) = sen x es T = 2π. ¿Cuál es el periodo fundamental de
cada una de las funciones siguientes?

a) f (x) = cos 2πx

b) f (x) = sen L4 x
P  
c) f (x) = A0 + ∞ n=1 An cos nπ
p x + B n sen nπ
p x , An y Bn dependen solamente de n.

6. En los siguientes problemas, encuentre la serie de Fourier de f en el intervalo dado.


1, −1 < x < 0
a) f (x) =
x, 0 ≤ x < π
 2
π , −π < x < 0
b) f (x) = 2 2
π −x , 0 ≤ x< π

0, −π < x < 0
c) f (x) =
sen x, 0 ≤ x < π
7. Compruebe las 6 propiedades de las funciones pares e impares

8. En los siguientes problemas, desarrolle la función dada en una serie apropiada de cosenos
y senos.

39
 
−1, −π < x < 0
a) f (x) =
1, 0≤x<π
b) f (x) = x2 , −1 < x < 1
c) f (x) = π 2 − x2 , −π < x < π

x − 1, −π < x < 0
d) f (x) =
x + 1, 0 ≤ x < π
9. En los siguientes problemas encuentre la representación de la integral de Fourier para la
función dada.


 0, x < −1

−1, −1 < x < 0
a) f (x) =

 2, 0<x<1

0, x>1

 0, x < π
b) f (x) = 4, π < x < 2π

0, x > 2π

0, x<0
c) f (x) = −x
e , x>0
10. Encuentre la transformada de Fourier de las siguientes funciones
 −x
e , x>0
a) f (x) =
0, x<0

 −1, −1 < x < 0
b) f (x) = 1, 0 < x < 1π

0, en caso contrario
11. Demostrar cada una de las propiedades de la transformada de Fourier.

40
4. PROBLEMA DE AUTOVLORES DE STURN-LIOUVILLE
Cuando estudiamos el método de separación de variables, se observó que dada una EDP de
segundo orden, en algunos casos es posible separar las variables y obtener ecuaciones ordinarias,
una de las cuales tiene la forma de una EDO de segundo orden lineal:
y ′′ + λy = 0 (77)
Vimos también que, haciendo un análisis de los posibles valores de los autovalores, existen
tres casos y sus respectivas soluciones:

 c1 + c2 x, λ=0
y(x) = c1 cosh αx + c2 senh αx, λ = α2 < 0

c1 cos αx + c2 sen αx, λ = α2 > 0
Si a este problema le añadimos condiciones de valor de frontera tales como:
y(0) = 0, y(L) = 0, (78)
Podemos comprobar que dicho problema tiene soluciones no triviales solamente cuando λ =
α2 > 0. Luego para este caso, si reemplaamos las condiciones de frontera (78), se obtiene:
y(0) = c1 cos 0 + c2 sen 0 = 0 ⇒ c1 = 0
y(L) = c2 sen αL = nπ n = 1, 2, 3, ...
nπ
Es decir αn = L y por lo tanto, puesto que λ = α2 , se tiene que:

n2 π 2
λn = , n = 1, 2, 3, ...
L2
Los valores de λn son llamados Valores Propios, mientras que las correspondientes soluciones:
nπ
yn = c2 sen x, n = 1, 2, 3, ...
L
son llamadas funciones propias del problema.
Nótese además que este conjunto de autofunciones o funciones propias es ortogonal en el
intervalo [0, L]. Cabe señalar que este problema constituye un caso particular de un importante
problema general de valores de frontera llamado Problema de autovalores de Sturn Liouville.

4.1. Problema regular de autovalores de Sturn-Liouville

. Un problema regular de autovalores de Sturn-Liouville consiste en una ecuación diferencial
de Sturn-Liouville de la forma:
 
d dφ
p(x) + q(x)φ + λp(x)φ = 0 (79)
dx dx
sujeta a las condiciones de frontera
A1 φ(a) + B1 φ′ (a) = 0 (80)
′
A2 φ(b) + B2 φ (b) = 0 (81)
donde p, q, r y r ′ son funciones contı́nuas con valores reales en el intervalo [a, b] , r(x) > 0 y
p(x) > 0 para todo x en el intervalo. Por otro lado, se supone también que los coeficientes de
las condiciones de frontera son reales e independientes de λ y distintos de cero.
Algunos casos particulares de las ecuaciones de Sturn-Liouville son:

41
 2
1. Caso simple: ddt2φ + λφ = 0,
en este caso r = 0, q = 0 y p = 1.

2. Flujo de calor en una varilla uniforme:

 
d dφ
k0 + λcρφ = 0,
dx dx

en este caso r = k0 , q = α y p = cρ

3. Vibración en una cuerda no uniforme:

d2 φ
T0 + αφ + λρ0 φ = 0
dt2
En este caso r = T0 , q = α y p = cρ

4. Flujo de calor circular  

d dφ
r + λrφ = 0,
dr dr
en este caso x = r, q = 0 y p = x

Existe una cantidad importante de teoremas relativos al problema de Sturn-Liouville

El siguiente teorema es una lista de algunas de las propiedades del problema regular de
Sturn-Liouville.

Teorema: (Propiedad del problema regular de Sturn-Liouville)

1. Existe un número infinito de valores propios reales que pueden disponerse en orden ascen-
dente λ1 < λ2 < λ3 , ... < λn < ... de tal menera que λn → ∞ a medida que n → ∞.

2. Para cada valor propio existe solamente una función propia (excepto para múltiples cons-
tantes diferentes de cero).

3. Las funciones propias correspondientes a los diferentes valores propios son linealmente
independientes.

4. El conjunto de funciones propias correspondientes al conjunto de valores propios es orto-

gonal respecto a la función de peso p(x) en el intervalo [a, b].

5. La autofunción φn (x) forman un conjunto completo de modo que cualquier función contı́nua
a trozos f(x) se puede repreentar (expandir) mediante una serie generalizada de Fourier de
dicha autofunción.
∞
X
f (x) ≈ an φn (x),
n=1

más aún, dicha serie infinita converge a:

f (x+ ) + f (x− )
para a<x<b
2

42
5. SOLUCIÓN DE ECUACIONES DIFERENCIALES DE SE-
GUNDO ORDEN
5.1. Solución de la Ecuación del Calor
Consideremos una varilla delgada de longitud L con temperatura inicial f (x) en toda ella y
cuyos extremos se mantienen a una temperatura de cero en todo el tiempo t > 0, tal y como
se muestra en la figura (7). Si la varilla satisface los supuestos de la sección 1.9.1, entonces la

Figura 7: Econtrar la temperatura u en una varilla finita.

temperatura u(x, t) se puede determinar mediante el problema de valores en la frontera.

∂2u ∂u
k 2
= , 0 < x < L, t > 0
∂x ∂t
C.C. u(0, t) = 0, u(L, t) = 0 t > 0
C.I. u(x, 0) = f (x), 0 < x < L
Utilizando ahota el método de separación de variables, suponiendo una solución de la forma:
u(x, t) = X(x)T (t) y λ una constante de separación, se obtiene la siguiente igualdad.
d2 X dT
kT (t) 2
= X(x) = −λ
dx dt
separando variables se tiene:
1 d2 X 1 dT
2
= = −λ
X dx kT dt
de lo anterior se obtienen dos ecuaciones diferenciales ordinarias de primer y segundo orden
d2 X dT
2
+ λX = 0 y + kλT = 0 (82)
dx dt
Por otro lado, las condiciones de frontera se convierten en:
u(0, t) = X(0)T (t) = 0yu(L, t) = X(L)T (t) = 0
Como esto tiene que ser válido para todo tiempo t, se tiene X(0) = 0 y X(L) = 0.
Estas condiciones de frontera homogéneas, junto con la EDO de segundo orden para X(x)
constituye un problema de autovalores de Sturn-Liouville. Esto es:
d2 X
+ λX = 0, X(0) = 0, X(L) = 0 (83)
dx2
De modo que las soluciones generales correspondientes a esta ecuación son:
X(x) = C1 + C2 x si λ = 0,
X(x) = C1 cosh αx + C2 senh αx, si λ = α2 < 0,
X(x) = C1 cos αx + C2 sen αx, si λ = α2 > 0,

43
al aplicar las condiciones de frontera para los casos λ = 0 y λ < 0, se obtiene la solución trivial
X(x) = 0, la cual no nos brinda ninguna información, de modo que nos quedamos sólo con la
solución para λ > 0. Luego aplicando la primera condición de frontera X(0) = 0 a la solución
dada obtenemos C1 = 0. Por lo tanto X(x) = C2 sen αx. Luego la segunda condición de frontera
X(L) = 0, implica que:
X(L) = C2 sen αL = 0
Si C2 = 0 ⇒ X = 0 por lo que u = 0. Cuando C2 6= 0 ⇒ sen αL = 0. Esto es cierto
si αL = nπ,
n = 1, 2, 3, ... En consecuencia esta ecuación tiene solución no trivial cuando
2
λn = α2n = nπ
L , n = 1, 2, 3, .... Luego se tiene que:
nπ
X(x) = C2 sen αx = C2 sen x, n = 1, 2, 3, ...
L
De esta forma λn y Xn (x) son, respectivamente, los valores propios y las funciones propias del
problema planteado.
Por otra parte, la solución general para la ecuación dependiente del tiempo es:

T (t) = C3 e−λkt
n2 π 2
Sustituyendo el valor de λn = L2
se tiene:

n2 π 2
 
−k t
T (t) = C3 e L2 (84)

Por lo tanto, la solución general de la ecuación del calor será:

n2 π 2 nπ
 
−k t
un (x, t) = X(x)T (t) = An e L2 sen x (85)
L
donde las constantes C1 y C2 se han reemplazado por una única contante An .
Esta solución satisface la ecuación diferencial parcial del calor con condiciones de frontera
para cada valor positivo de n. Sin embargo, para que satisfaga la condición inicial hacemos:
n2 π 2 nπ nπ
 
−k (0)
un (x, 0) = f (x) = An e L2 sen x = An sen x,
L L
de manera que un (x, 0) = An sen nπ
L x. Pero un (x, t) no será solución de la EDP del calor a menos
que calculemos los coeficientes An . Con ese objetivo, aplicamos el principio de superposición
generalizado y se obtiene:
∞ ∞ n2 π 2 nπ
 
X X −k t
U (x, t) = un = An e L2 sen x (86)
n=1 n=1
L

Lo que se obtiene es una serie de Fourier de senos. Luego si sustituimos en esta última expresión
la condición inicial se tiene:
∞
X nπ
U (x, 0) = f (x) = An sen x
L
n=1

Por lo tanto esta última expresión es la expansión de medio intervalo de f en una serie de senos.
De modo que: Z
2 L nπ
An = f (x) sen xdx
L 0 L

44
Concluimos que una solución al problema de valor de frontera está dada por la serie infinita:
∞ Z L   2 2
2X nπ −k n π2 t nπ
U (x, t) = f (x) sen xdx e L sen x (87)
L 0 L L
n=1

Por ejemplo, si la temperatura inicial es u(x, 0) = 100, L = π y k = 1, entonces:

Z
2 π
An = 100 sen nxdx
π 0
Z π
200
= sen nxdx
π 0
 
200 1 − (−1)n
=
π n

Por lo tanto, la serie resultante será:

∞  
200 X 1 − (−1)n −n2 t
U (x, t) = e sen nx
π n
n=1

5.2. Solución de la Ecuación de Ondas

El desplazamiento vertical u(x, t) de una cuerda de longitud L que se encuentra vibrando
libremente en el plano vertical esta determinado por la ecuación:

∂2u ∂2u
a2 = , 0 < x < L, t > 0
∂x2 ∂x2
C.C. u(0, t) = 0, u(L, t) = 0 t > 0
∂u
C.I. u(x, 0) = f (x), |t=0 = g(x), 0 < x < L
∂t
Solución: Supongamos que u(x, t) = X(x)T (t), mediante separación de variables obtenemos:

d2 X d2 T
a2 T = X
dx2 dt2
1 d2 X 1 d2 T
= X = −λ
x dx2 a2 T dt2
por lo que se obtienen dos ecuaciones de segundo orden.

d2 X d2 T
= −λX y = −λa2 T
dx2 dt2
Por otro lado, de las condiciones de frontera se obitnene que:

X(0) = 0 y X(L) = 0

La ecuación odinaria para X junto con estas condiciones de frontera forman un problema regular
de Sturn-Liouville

x′′ + λX = 0 (88)
X(0) = 0 y X(L) = 0 (89)

45
De los tres posibles valores del parámetro λ, (λ = 0; λ = −α2 < 0, λ = α2 > 0), solamente la
última nos lleva a soluciones no triviales. La solución general de esta ecuación correspondiente
a λ = α2 > 0, α > 0, es:
x = c1 cos αx + c2 sen αx
Reemplazando las condiciones de frontera X(0) = 0 y X(L) = 0 se obtiene que: c1 = 0 y
c2 sen αL = 0. Esto implica que: α = nπ
L . Por lo tanto los valores propios y las funciones propias
son respectivamente:
 nπ 2 nπ
λn = y X(x) = c2 sen x, n = 1, 2, 3, ...
L L
Por otro lado, la solución general de la EDO para la variable temporal t es:
nπ nπ
T (t) = C3 cos x + c4 sen x
L L
Renombrando a c3 ∆c2 como An y a c2 c4 como Bn , la solución que satisfacen tanto la ecuación
de onda como las condiciones de frontera es:
 nπa nπa  nπ
un (x, t) = An cos t + Bn sen t sen x (90)
L L L
por lo tanto, la solución general será:
∞ 
X nπa nπa  nπ
U (x, t) = An cos t + Bn sen t sen x (91)
L L L
n=1

Aplicando ahora la condición incial u(x, 0) = f (x), obtenemos:

∞
X nπ
u(x, 0) = f (x) = An sen x
n=1
L

Dado que esta última expresión es una serie de Fourier de medio intervalo, se tiene que:
Z L
2 nπ
An = f (x) sen x (92)
L 0 L

Para determinar el coeficiente Bn aplicamos la segunda condición incial de tipo Newmann. Para
ello derivamos con respecto al tiempo la solución general (91) y reemplazamos:
∞ 
∂u X nπa nπa nπa nπa  nπ
= −An sen t + Bn cos t sen x
∂t L L L L L
n=1
 ∞ 
∂u  X nπa  nπ
= g(x) = B n sen x
∂t t=0
 L L
n=1

Despeando el coeficiente, de acuerdo con los coeficientes de Fourier

Z L
2 nπ
Bn = g(x) sen x (93)
nπa 0 L

Por lo tanto los coeficientes correspondientes a la solución general (91) son (92) y (93.)

46
6. REFERENCIAS BIBLIOGRÁFICAS

[1] E. Zachmanoglon and [Link]. Introduction to Partial Differential Equations with Applica-
tions. Dover, 1986.

[2] L. C. Evans. Partial Differential Equations. Graduate Studies in Mathematics, Vol 19, AMS
1997.

[3] R. Haberman Elementary Applied Partial Differential Equations with Fourier Series and
Boundary Value Problems. Prentice Hall,1987.

[4] A. Jeffrey Advanced Engineering Mathematics. Harcourt Academy Press, 2002.

[5] R.K. Nagel, E. B. Staff and A. D. Snider Ecuaciones Diferenciales y Problemas con Valores
en la Frontera. Pearson Educations. Mexico 2005.

[6] J. B. Espinoza y D. Elizarraz M. Ecuaciones Diferenciales, Técnicas de Solución y Aplica-

ciones. Universidad Autónoma Metropolitana. Azcapotzalco, México 2004.

[7] Y. Pinchover and J. Rubinstein. Introduction to Partial Differential Equations. Cambridge

University Press, 2005.

[8] I. P. Alonso. Primer Curso de Ecuaciones en Derivadas Parciales. Universidad Autónoma

de Madrid. Departamento de Matemáticas.

[9] Manuel Mañas B. y Luis Martı́nez A. Ecuaciones Diferenciales II. Universidad Complutense
de Madrid. Departamento de Fı́sica Teórica.

[10] Erwin Kreyszig Matemáticas Avanzadas Para Ingenierı́a Vol II. Tercera Edición. Limusa
Wiley, 2003.

[11] Qing Han. An Introduction to Linear Partial Differential Equations.

[12] Denis G. Zill., Jacqueline M. Dewar. Matemáticas Avanzadas Para Ingenierı́a Vol. 2.
McGraw-Hill. Tercera Ed. 2008.

[13] Premk Kithe. Partial differential equations and Mathematica. CRC Press, 1997.

[14] A.D. Polyanin, V.F Zaitev, A. Moussiaux. Handbook of First Order Partial Differential
Equations. Taylor and Francis, 2002.

[15] L. Lapidus, G. Pinder. Numerical Solutions of Partial Differential Equations in Science and
Engineering, John Wiley and Sons Inc, 1999.

47