Tarea 2

Matemáticas Avanzadas. Semestre 2024-2

Profesores: Adán G. Andrade y Mauricio Garcı́a Vergara
Ayudante: Karen Escutia Solı́s

Fecha de entrega: 15 de abril a las 12h. Nota: En las secciones donde se te pide graficar te
sugiero hacerlo con algún software de tu elección. Puedes también tomar la referencia de las gráficas
de algún libro o lugar confiable (en ese caso no olvides añadir la referencia).

Ejercicio 1.
Responde a lo que se te indica y justifica tus razonamientos.
a) Explica por qué es natural que surjan dos familias de soluciones linealmente independientes
entre sı́, al resolver la ecuación de calor en una simetrı́a polar (r, θ).
b) Escribe la ecuación de Bessel paramétrica, y explica a qué está relacionado el parámetro λ de
dicha ecuación.
c) Grafica las funciones Bessel de primer y segundo tipo de orden cero J0 (r), Y0 (r).
d) Explica, en términos fı́sicos, por qué la función K0 no puede ser una solución al resolver la
ecuación de Laplace en un cilindro con condiciones de Dirichlet en la frontera.
e) Grafica Jp para n = 0, 12 , 1, 2, 5.
f) Explica qué es el método de Frobenius y especifica qué condiciones se necesitan para poder
aplicarlo.
g) Explica por qué (y bajo qué condición) la ecuación diferencial de Legendre es un caso particular
de la ecuación asociada de Legendre.
h) Grafica las polinomios de Legendre Pn (cos θ) para n = 0, 1, 2, 3.
i) Grafica las funciones asociadas de Legendre Pnm (cos θ) para m = 1 y n = 1, 2, 3.
j) Explica por qué es usual encontrar en los problemas de coordenadas esféricas se use el término
simetrı́a angular para referirse a la simetrı́a respecto al ángulo azimutal φ, sin generar confusión
con la coordenada angular restante θ.
k) En cada caso, da un ejemplo de bajo qué simetrı́a y qué condiciones de contorno puedes esperar
las siguientes soluciones:

a) Funciones Bessel de primer tipo de orden cero: J0 .

b) Funciones Bessel de primer tipo de orden p: Jp .
c) Funciones Bessel modificadas de primer tipo y orden 0: I0 .
d) Funciones asociadas de Legendre: Pnm (cos θ).
e) Polinomios de Legendre: Pn (cos θ).
f) Armónicos esféricos: Yn,m (θ, φ).

1
Ejercicio 2.
La finalidad de este ejercicio es utilizar la teorı́a de las funciones Bessel para obtener una expansión
en series Bessel para la función 1 − r2 .

a) Escribe la relación de recurrencia para las funciones de Bessel de orden p ≥ 0.

b) A partir de la expresión del inciso a) muestra que

d p
[x Jp (x)] = xp Jp−1 (x) . (1)
dx

(Nota: cuidado al introducir xp en la suma. Necesitarás multiplicar por un 1 adecuado para

mantener la forma de la factorización del término x/2. Además, recuerda la propiedad:
Γ(k + p + 1) = (k + p)Γ(p + k)).

c) Muestra que

d  −p
x Jp (x) = −x−p Jp+1 (x) .

(2)
dx
P∞ (−1)k 2k
Pista: Muestra que la derivada es: k=1 k!Γ(k+p+1)22k+p
x2k−1 . Luego propón el cambio k → k+1,
para comenzar la suma desde k = 0.
d) Muestra que

Z
xp+1 Jp (x)dx = xp+1 Jp+1 (x) + C . (3)

e) Muestra que:

a
ap+2
Z
α
xp+1 Jp ( x)dx = Jp+1 (α) . (4)
0 a α

Nota: encontrarás conveniente el cambio de variable: u = αx/a.

f) Considera una membrana circular (de radio a) que vibra con simetrı́a radial y que se encuentra
anclada a lo largo de su circunferencia. Escribe la ecuación que gobierna las oscilaciones y que
tiene por soluciones a u(r, t). Plantea además, las condiciones de contorno e iniciales, suponiendo
que la rapidez inicial es g(r) y la forma inicial es f (r).
g) Muestra que la solución al problema está dada por:

∞
X
u(r, t) = J0 (λn r) [An cos (cλn t) + Bn sin (cλn t)] . (5)
n=1

Escribe la expresión para los coeficientes An , Bn y especifica qué es λn .

h) Se estudia el caso particular, donde la forma inicial de la membrana está dada por la función
f (r) = 1 − r2 . Además, se considera que se tiene velocidad inicial cero. Y, por simplicidad, se
toma el caso a = c = 1. Muestra que los coeficientes de la solución en este caso toman la forma:

2
 Bn = 0 (6)
4
An = 2 2 J2 (αn ) . (7)
αn J1 (αn )

Pista: encontrarás útil lo calculado en los incisos d) y e). Además, el cambio de variable u =
αn2 − s2 te será de utilidad para realizar una integración por partes.
i) Usando la expresión (estarı́a genial si te tomas el tiempo de demostrarla)

2p
Jp+1 (x) = Jp (x) − Jp−1 (x) , (8)
x
muestra que el coeficiente puede simplificarse hasta tomar el valor

8
An = . (9)
αn3 J1 (αn )

j) Escribe la solución u(r, t) al problema y muestra que para el caso t = 0 se tiene que

∞
X 8
1 − r2 = J (αn r) ,
3 J (α ) 0
(10)
α
n=1 n 1 n

con 0 < r < 1. Esta es la expansión de la función 1 − r2 en términos de la serie Bessel.

Ejercicio 3.
La finalidad de este ejercicio es introducir las funciones Bessel esféricas de primer tipo.
a) La ecuación diferencial para la distribución de temperatura u en una esfera sólida, bajo condi-
ciones de simetrı́a radial y estacionarias, está dada por:
d2 u 2 du [kr2 − n(n + 1)]
+ + u = 0, 0 < r < a, u(a) = 0 , (11)
dr2 r dr r2
donde u representa la temperatura, r es la distancia radial desde el centro de la esfera, n =
0, 1, 2, ... es un parámetro que depende de las condiciones de contorno y de la simetrı́a del pro-
blema y k es algún número real no negativo. Se buscan soluciones en el intervalo 0 < r < a. Haz
el cambio de variable

u = r−1/2 w , (12)

y sustituye en la ecuación (11) para mostrar que toma la forma

1
r2 w′′ + rw′ + [kr2 − (n + )2 ]w = 0 , w(a) = 0 . (13)
2
b) Identifica qué valor debe tomar p para conectar la expresión (13) con la la forma paramétrica
de la ecuación de Bessel de orden p.
c) Explica qué condición (del problema fı́sico) se debe satisfacer para que sea necesario elegir
k = λ2 , donde

αn+1/2,j
λ = λn,j = , (14)
a
considerando que αn+1/2,j denota el j-ésimo cero positivo de la función Jn+1/2 .

3
 d) Las correspondientes soluciones resultan de la forma

wn,j (r) = An Jn+ 12 (λn,j r) , n = 0, 1, 2, ... , j = 1, 2, ... , (15)

con An alguna constante por determinar.
Si se revierte el cambio de variable en (12) uno es capaz de identificar que las soluciones al
problema original deben ser múltiplos escalares de:

r−1/2 Jn+ 21 (λn,j r) .

Ası́, se definen las funciones Bessel esféricas de primero tipo como
 π 1/2
jn (r) = Jn+ 12 (λn,j r) , n = 0, 1, 2, ... . (16)
2r
Estas funciones especiales permiten escribir la solución al problema como

un,j (r) = An jn (λn,j r) , n = 0, 1, 2, ... , j = 1, 2, ... . (17)

Las jn codifican la información del orden fraccionario p = n+1/2 de las funciones Bessel originales
Jp . Grafica j0 , j1 , j2 , j3 .
e) Explica a partir de las gráficas, por qué las funciones Bessel esféricas jn perfectamente pueden
representar (en alguna combinación lineal) una solución para el problema planteado. (Pista:
piensa en la idea de soluciones acotadas).
f) A partir de las relaciones de ortogonalidad de las funciones Bessel Jp muestra que las funciones
Bessel esféricas también son ortogonales en el intervalo [0, a] con respecto a una función de peso
x2 , es decir,

Z a
jn (λn,j x)jn (λn,j ′ x)x2 dx = 0 , j ̸= j ′ , (18)
0
a
a3 2
Z
jn2 (λn,j x)x2 dx = j (αn+1/2,j ) . (19)
0 2 n+1

Aquı́ αn+1/2,j es el j-ésimo cero positivo de jn (y de forma equivalente el de Jn+1/2 .)

Ejercicio 4.
El tejido cardı́aco es un material complejo con propiedades eléctricas que determinan el ritmo
cardı́aco de una persona por medio de procesos de sistole y diasole. La ecuación de Laplace se puede
utilizar para modelar el potencial eléctrico dentro del tejido cardı́aco, lo que permite comprender mejor
la propagación de los impulsos eléctricos y la función cardı́aca. Consideremos un tejido cardı́aco de
forma cilı́ndrica con radio R y altura 2a. La ecuación de Laplace que describe el potencial eléctrico
u(r, θ, z) dentro del tejido es:

1 ∂2u ∂2u
 
1 ∂ ∂u
∆u = r + 2 2 + 2 =0
r ∂r ∂r r ∂θ ∂z
con las siguientes condiciones de frontera ”laterales”:
1. u|z=a = 0, u|z=−a = 0 (Potencial eléctrico nulo en la parte superior e inferior del tejido cardı́aco).
2. u|r=R = f (z, θ) − a, −a < z < a, f (z, θ + 2π) = f (z, θ) (Distribución de estimulación
eléctrica en el borde del tejido cardı́aco).
Calcula el potencial eléctrico en el tejido cardı́aco.

4
Ejercicio 5.
Supongamos que estamos estudiando la difusión de una sustancia quı́mica a través de un anillo de
tejido biológico que representa una estructura anular en el cuerpo humano, como un vaso sanguı́neo o
una sección de un órgano tubular.
La concentración de la sustancia quı́mica en el tejido se representa por la función u(r, θ), donde
r es la distancia radial desde el centro del anillo y θ es el ángulo azimutal que describe la posición
alrededor del anillo. La ecuación de Laplace describe la difusión de la sustancia quı́mica en el tejido,
y se expresa como:

∇2 u = 0
Las condiciones de contorno para este problema son las siguientes:

En el cı́rculo interno del anillo, donde el radio es r1 = 1, la concentración de la sustancia quı́mica

es constante y igual a 1 en todas las posiciones alrededor del anillo: u(1, θ) = 1 para todo θ.
En el cı́rculo externo del anillo, donde el radio es r2 = 3, la concentración de la sustancia quı́mica
varı́a sinusoidalmente con el ángulo azimutal θ, con una amplitud decreciente a medida que
aumenta el número de onda n. Además, u(3, θ) = 1 + θ.

Este modelo puede ayudar a comprender cómo se distribuye y difunde una sustancia quı́mica en
estructuras anulares biológicas y cómo las condiciones en los bordes del anillo afectan esta distribución.

a) Realiza un esquema del dominio de la ecuación de Laplace.

b) Resuelve la ecuación de Laplace y da la expresión de la distribución de nutrientes.