PONTIFICIA UNIVERSIDAD CATÓLICA DE CHILE

FACULTAD DE FÍSICA
FIZ0313 - Métodos de la Fı́sica Matemática II

Profesor: Roberto Rodrı́guez

Alumno: Néstor Espinoza

Solución Tarea 2
Fecha de entrega: Lunes 20 de Abril.

1. Problema 1
(a) La fórmula de Rodrigues correspondiente al polinomio de Hermite n-ésimo es:
2 dn  −x2 
Hn (x) = (−1)n ex e
dxn
Ası́, tendrı́amos que:
∞ ∞ 
X tn X 2 d
n 
−x2
 tn
Hn (x) = (−1)n ex e (1)
n=0
n! n=0 dxn n!
Pero, de acuerdo con el Teorema de Couchy, es posible expresar la derivada n-ésima de una función de
acuerdo a la igualdad:
dn n! f (z)
I
(f (z))|z = dz, z, z0 ∈ C (2)
(z − z0 )n+1
n 0
dx 2πi
Ası́, reemplazando ésto en (1):
∞  2
!
n! e−z tn
 I
n x2
X
(−1) e n+1 dz
n=0
2πi (z − x) n!

La expresión puede ser ordenada “sacando” la integral de la sumatoria, notando que los n! se cancelan.
Ası́:
2 2 2 2 2 2 n
ex e−z tn (−1)n ex e−z (−1)n tn
∞ ∞ ∞
ex e−z

−t
I X I X I X
n+1 dz = 2πi (z − x) (z − x)
n dz =
2πi (z − x) z − x
dz
n=0 2πi (z − x) n=0 n=0

De ésta manera, podemos ordenar mas compactamente:

2 2 ∞  n
ex e−z X −t
I
dz
2πi (z − x) n=0 z − x
 La suma correspondiente es, de hecho, una suma geométrica1. Ası́, imponiendo la condición |t/(z−x)| <
1: 2 2
ex e−z
 
z−x
I
dz
2πi (z − x) z − x + t
O mejor. simplificando los términos (z − x):
2 2 2 2
ex e−z ex e−z
I I
dz = dz
2πi (z − x + t) 2πi (z − (x − t))
Ocupando la ecuación (2) de forma inversa (ésto es, pasando de la forma integral a la forma diferencial)
encontramos, dada la forma de la integral, que debemos ocupar el teorema para n = 0, ésto es, no hay
derivada, por ende solo debemos evaluar la función en z0 = x − t. Ası́:
2 2
ex e−z
I  
2 2 2 2 2 2
dz = ex e−(x−t) = ex e−x +2xt−t = e2xt−t
2πi (z − (x − t))
O sea, hemos demostrado que:
∞
X tn 2
Hn (x) = e2xt−t = ψ(x, t) (3)
n=0
n!

(b) Usando la función generatriz deducida en el ejercicio anterior, tomemos el producto interno (ψ(x, t), ψ(x, t)):
Z ∞
2
ψ(x, t)2 e−x dx
−∞

Por un lado, tendrı́amos, según la ecuación (3):

Z ∞ Z ∞ ∞ ∞
2 2
 Z Z
2 2 2 2 2
+2t2
ψ(x, t)2 e−x dx = e2xt−t e−x dx = e4xt−2t e−x dx = e−(x−2t) dx
−∞ −∞ −∞ −∞

O sea: Z ∞ Z ∞
2t2 −(x−2t)2 2t2 2
e e dx → µ = x − 2t → e e−µ du
−∞ −∞
R∞ −µ2
√
Y como −∞
e du = π:
∞ √
Z
2 2 2
e2t e−µ du = πe2t (4)
−∞
Y por otro lado también:
∞
! ∞
! ∞
∞ ∞
tn tk ∞
t2n −x2
Z Z Z
2 −x2 2
X X X
ψ(x, t) e dx = Hn (x) Hk (x) e−x dx = Hn (x)2 e dx
−∞ −∞ n=0
n! k! −∞ n=0 (n!)2
k=0

En donde se usó el hecho de la ortogonalidad entre los polinomios de Hermite. Ası́, podemos escribir
la ecuación anterior ocupando el hecho que la integral de una suma es la suma de las integrales:
∞ Z ∞ ∞
t2n −x2 t2n
Z ∞
X X 2
Hn (x)2 2
e dx = 2
Hn (x)2 e−x dx
n=0 −∞
(n!) n=0
(n!) −∞

1 Recordemos que una suma geométrica es del tipo:

∞
X 1
xn =
n=0
1−x

La suma converge al valor dado siempre y cuando x sea menor que uno.
 Ası́, como esta última integral es (Hn (x), Hn (x)) = ||Hn (x)||2 :
∞
X t2n
2
||Hn (x)||2 (5)
n=0
(n!)

Ası́, igualando las ecuaciones (4) y (5) tendrı́amos:

∞
X t2n √ 2

2
||Hn (x)||2 = πe2t
n=0
(n!)

Obsérvese que la exponencial de la derecha puede ser expandida en Serie de Taylor2, de modo que
quedarı́a:
∞ ∞
X t2n 2
√ X (2t2 )n
||H n (x)|| = π
n=0
(n!)2 n=0
n!

Ası́ igualando los coeficientes de las potencias t2n :

√
||Hn (x)||2 (2)n π √
= → ||Hn (x)||2 = 2n πn!
(n!)2 n!

(c) (i) Tenemos:

∞
X tn 2
Hn (x) = e2xt−t (6)
n=0
n!
Si derivamos con respecto a t, el lado derecho queda:
∂ψ 2
= (2x − 2t)e2xt−t = (2x − 2t)ψ(x, t)
∂t
Pero, ésto es equivalente a:
∞ ∞
X tn X tn+1
(2x − 2t)ψ(x, t) = 2x Hn (x) −2 Hn (x)
n=0
n! n=0
n!

A su vez, el lado izquierdo de la ecuación (6) queda:

∞ ∞
∂ψ X ntn−1 X tn−1
= Hn (x) = Hn (x)
∂t n=1
n! n=1
(n − 1)!

Ası́, uniendo las dos expresiones obtenidas para ∂ψ(x, t)/∂t e igualando potencias de t:
∞ ∞ ∞
X tn X tn X tn
Hn+1 (x) = 2x Hn (x) − 2 Hn−1 (x)
n=0
n! n=0
n! n=1
(n − 1)!

O mejor:
∞ ∞ ∞
X tn X tn X tn
Hn+1 (x) − 2x Hn (x) + 2 Hn−1 (x) =0
n=0
n! n=0
n! n=1
(n − 1)!
2 Recordar que la expansión en Serie de Taylor de la exponencial es:
∞
X xn
ex =
n=0
n!
 Pero el lado izquierdo puede ordenarse mas aún:
∞
X tn tn tn
H1 (x) − 2xH0 (x) + Hn+1 (x) − 2xHn (x) + 2Hn−1 (x) =0
n=1
n! n! (n − 1)!

Obsérvese que como H1 (x) = 2x y H0 (x) = 1, H1 (x) − 2xH0 = 0. Ası́, quedamos con:
∞
tn
 
X Hn+1 (x) 2xHn (x)
− + 2Hn−1 (x) = 0
n=1
n − 1! n n

Pero como lo que está fuera del paréntesis no es siempre cero, lo de adentro debe serlo. Ası́:

Hn+1 (x) = 2xHn (x) − 2nHn−1 (x), n ≥ 1 (7)

Que era lo que se querı́a demostrar.

(ii) Si derivamos con respecto a x la ecuación (6), el lado derecho queda:
∂ψ 2
= 2te2xt−t = 2tψ(x, t)
∂x
Pero: ∞ ∞ ∞
X tn X tn+1 X tn
2tψ(x, t) = 2t Hn (x) =2 Hn (x) =2 Hn−1 (x)
n=0
n! n=0
n! n=1
(n − 1)!
Pero a la vez, el derivar con respecto a x también significa que el lado izquierdo de la ecuación
(6):
∞
∂ψ X tn
= Hn (x)′
∂x n=0
n!
O sea:
∞ ∞
X tn X tn
Hn (x)′ =2 Hn−1 (x)
n=0
n! n=1
(n − 1)!
Obsérvese que el lado izquierdo de ésta última expresión puede ser escrito como:
∞ ∞
X tn X tn
H0 (x)′ + Hn (x)′ =2 Hn−1 (x)
n=1
n! n=1
(n − 1)!

Pero como H0 (x) = 1, entonces H0 (x)′ = 0. Por lo tanto reescribiendo la expresión dejando un
lado igual a 0:
∞
tn
 
X Hn (x)′
− 2Hn−1 (x) = 0
n=1
(n − 1)! n
Como lo de afuera del paréntesis no puede ser 0, lo de adentro debe serlo. Ası́:

Hn (x)′ = 2nHn−1 (x), n ≥ 1 (8)

Que era lo que se querı́a demostrar.

(d) Usando las ecuaciones (7) y (8) tendrı́amos que:

Hn (x)′′ = 2nHn−1 (x)′ = 2Hn (x) + 2xHn (x)′ − Hn+1 (x)′ (9)

Pero ahora podemos expresar el último término con la ecuación (8), quedándonos con:

Hn+1 (x)′ = 2(n + 1)Hn (x)

 Ası́, reemplazando ésto en (9):

Hn (x)′′ = 2Hn (x) + 2xHn (x)′ − 2(n + 1)Hn (x) = 2xHn (x)′ − 2nHn (x)

O sea:
Hn (x)′′ − 2xHn (x)′ + 2nHn (x) = 0
Que era lo que se querı́a demostrar.
(e) Para resolver la integral la estrategia será formársela ocupando las relaciones de recurrencia para poder
obtener un valor que conozcamos. Obsérvese que queremos formarnos el término x2 , ası́, multiplicando
la ecuación (7) por x:
xHn+1 (x) = 2x2 Hn (x) − 2nxHn−1 (x)
O mejor:
xHn+1 (x)
x2 Hn (x) = + nxHn−1 (x)
2
Usando nuevamente la ecuación (7) en el lado derecho de esta última expresión, reemplazando los
términos que tienen a x como coeficiente, con Hn+2 y Hn :
 
Hn+2 (x) + 2(n + 1)Hn (x) Hn (x) + 2(n − 1)Hn−2 (x)
x2 Hn (x) = +n
4 2
Hn+2 (x) (n + 1)Hn (x) nHn (x)
= + + + n(n − 1)Hn−2 (x)
4 2 2
Hn+2 (x) Hn (x)
= + nHn (x) + + n(n − 1)Hn−2 (x)
4 2
2
Ası́, multiplicando esta expresión por Hn e−x e integrando en el intervalo [−∞, ∞]:
Z ∞ Z ∞ 
2 −x2 Hn+2 (x) Hn (x) 2
x Hn (x)Hn (x)e dx = + nHn (x) + + n(n − 1)Hn−2 (x) Hn (x)e−x dx
−∞ −∞ 4 2
Z ∞ 
Hn (x) 2
= nHn (x) + Hn (x)e−x dx
−∞ 2
||Hn (x)||2
= n||Hn (x)||2 +
2

Ası́, el valor de la integral depende de ||Hn (x)||2 , cuyo valor encontramos en la parte (b). Ası́:
√ √ √
Z ∞
2
x2 Hn (x)Hn (x)e−x dx = n2n πn! + 2n−1 πn! = 2n−1 (2n + 1) πn!
−∞

2. Problema 4
(a) Para evaluar la integral nos fijamos que el denominador es muy parecido a la función generatriz ψ(x, t)
de los polinomios de Legendre Pn (x):
∞
1 X
ψ(x, t) = √ = Pn (x)tn (10)
1 − 2xt + t2 n=0

Ası́, hacemos el cambio de variables −x = cos(θ), t = b/a, de modo que la integral queda:
Z 1 Z 1 Z 1 Z 1
t−x
√ dx = ψ(x, t)(t − x)dx = tψ(x, t)dx − xψ(x, t)dx
−1 1 − 2xt + t2 −1 −1 −1
 Escribiendo las funciones generatriz como sumatorias, según la ecuación (10):
Z 1 X∞ Z 1 X∞
1· Pn (x)tn+1 dx − x Pn (x)tn dx
−1 n=0 −1 n=0

Pero, como P0 (x) = 1 y P1 (x) = x:

Z 1 1
4b
Z
P0 (x)2 tdx − P1 (x)2 tdx = t(||P0 (x)||2 − ||P1 (x)||2 ) =
−1 −1 3a

(b) Similarmente al caso anterior, hacemos el mismo cambio de variables, lo que nos queda, como x2 =
cos(θ)2 = 1 − sin(θ)2 :
Z 1 2 Z 1 Z 1 Z 1
x (1 − t) + t
√ dx = ψ(x, t)(x2 (1 − t) + t)dx = x2 ψ(x, t)(1 − t)dx + ψ(x, t)tdx
−1 1 − 2xt + t2 −1 −1 −1

Ahora bien, como P2 (x) = 21 (3x2 − 1) implica que x2 = 2P2 (x)+1

3 . Ası́, ocupando la forma de sumatoria
de la función generatriz:
Z 1 Z 1 Z 1 ∞
2P2 (x) + 1 X
x2 ψ(x, t)(1 − t)dx + ψ(x, t)tdx = Pn (x)tn (1 − t)dx + t||P0 (x)||2
−1 −1 −1 3 n=0

Ası́: ∞ ∞
1 Z 1
2P2 (x) + 1 X 2P2 (x) + 1 X
Z
n
Pn (x)t dx − Pn (x)tn+1 dx + t||P0 (x)||2
−1 3 n=0 −1 3 n=0
Lo que se reduce a:
2 2
1
t ||P2 (x)||2 − ||P2 (x)||2 t + ||P0 (x)||2 − ||P0 (x)||2 t + t||P0 (x)||2


3 3
O mejor:  
4 − 4t 2 − 2t b
t2 + + 2t, t =
15 3 a

3. Problema 5
(a) La ecuación que define el potencial pedido viene definida por:
∞  
X bn n
Φ(r, θ) = + a n r Pn (cos(θ)) (11)
n=0
rn+1

Primero, observamos que queremos evitar la discontinuidad cuando r → ∞, por lo que tomamos an = 0.
Ası́, vemos que la función potencial en la superficie de la esfera, i.e. en (R, θ) puede escribirse como:
Φ(R, θ) = ksin(θ/2)
 
1 − cos(θ)
= k
2
k
= (P0 (cos(θ)) − P1 (cos(θ)))
2
Ası́, evaluando en (R, θ) el potencial:
∞
X bn k
Φ(R, θ) = n+1
Pn (cos(θ)) = (P0 (cos(θ)) − P1 (cos(θ)))
n=0
R 2
Multiplicando a ámbos lados por Pk′ (cos(θ))sin(θ) y luego integrando en el intervalo [−π, π] (para usar
ortogonalidad):
∞
X bn k
δk′ ,n ||Pk′ (cos(θ))||2 = δk′ ,0 ||P0 (cos(θ))||2 − δk′ ,1 ||P1 (cos(θ))||2


R n+1 2
n=0

Ası́, para cada valor correspondiente a bk′ (ésto es, cuando k ′ = n):
′
Rk +1 k
δk′ ,0 ||P0 (cos(θ))||2 − δk′ ,1 ||P1 (cos(θ))||2


bk ′ =
2||Pk′ (cos(θ))||2

Solamente podemos obtener b0 y b1 (pues para otros valores de k ′ no podrı́amos haber despejado el
término ||Pk′ (cos(θ))||2 ). Ası́ reemplazando para k ′ = 0 y k ′ = 1 los valores de éstos coeficientes son:

KR
b0 =
2
KR2
b1 = −
2
Por tanto, reemplazando ésto en la función potencial original dada por la ecuación (11):

b0 P0 (cos(θ)) b1 P1 (cos(θ))
Φ(r, θ) = +
r r2
Lo que reemplazando los valores encontrados para b0 y b1 queda:
 
KR Rcos(θ)
Φ(r, θ) = 1−
2r r