CAPITULO III

FUNCION BESSEL
La función Bessel fue encontrada como solución de un problema de movimiento planetario por
el matemático y astrónomo Friedrich Wilhem Bessel. La aplicación de las funciones Bessel (es
más extensa que la de las funciones gamma) se extiende a los campos de ingeniería eléctrica,
acústica, hidrodinámica, termodinámica, electricidad y mecánica celeste.
Las funciones Bessel son soluciones particulares de la ecuación diferencial:
x2 y 00 + xy 0 + (x2 − n2 )y = 0 (1)
que recibe el nombre de ecuación de Bessel de orden n, aunque por supuesto,la ecuación
la ecuacion diferencial es de orden 2; n es cualquier número real o complejo. La forma de la
solución general depende del caracter de n.
El estudio de las funciones Bessel se empezará buscando una solución particular de la ecuacion
(1) para cuando n es real positivo, usando el método de Frobenius, donde se postula que la
solución es de la forma:
∞
X
y 1 = xα an x n
n=1

ya que inmediatamente se ve que x = 0 es un punto singular regular de la ecuación (1).

Entonces esta solución debe ser válida en 0 ≤ x ≤ ρ ahora sigue hacer el cuadro correspon-
0

diente:
x α
x +1 α
x +2 x + 3...α α
2 2 2 2 2
−n y −a n −a n an −a n
 
1 0 1 2 3
2
xy 1 a a 0 1
xy10 a0 c a1 (c + 1) a2 (c + 2) a3 (c + 3)
x2 y100 a0 c(c − 1) a1 c(c + 1) a2 (c + 1)(c + 2) a3 (c + 2)(c + 3)

(2)
La ecuación inicial se encuentra a partir de los términos con menor exponente (o usando
la fórmula 8 del Capitulo II):
a0 (c2 −n2 ) = 0 c = ±n (a0 6= 0 y arbitraria), como n > 0, por hipotesis, entonces :

c1 − c2 = 2n > 0
cuando 2n es igual a un entero positivo, se tiene la certeza de una solución de la forma
(1), la cual corresponde a la raiz más grande,
c1 = n

Del cuadro (2) se encuentra:

 a1 (n + 1)2 − n2 = 0
 
a1 = 0

−a0
a2 [(n + 2)2 − n2 ] + a0 = 0 a2 =
2(2n + 2)

a3 [(n + 3)2 − n2 ] + a1 = 0 a3 = 0

−a2
a4 [(n + 4)2 − n2 ] + a2 = 0 a4 =
4(2n + 4)

Por consiguiente:
−ak−2
ak = (k ≥ 2)
k(2n + k)
Nótese que todas las aes con subíndice impar valen cero. Con las aes pares se construyen
los coecientes en términos de a :
Así:
0

a 0
a4 =
2 × 4(2n + 2)(2n + 6)
−a0
a6 =
2 × 4 × 6(2n + 2)(2n + 4)(2n + 6)
↓
a0
a2k = (−1)k
2 × 4 × 6..,2k(2n + 2)(2n + 4)...(2n + 2k)
a0
= (−1)k
22k k!(n + 1)(n + 2)...(n + k)
por ende el término general de y es:
1

a0
(−1)k 2k
xn+2k
2 k!(n + 1)(n + 2)...(n + k)
Recuerden que la funcion de Bessel es una solución particular de la ecuacion diferencial
(1), por ello es conveniente seleccionar a a para que este término general que se acaba de
encontrar se simplique. Se podria tomar a = 1; pero resulta mejor:
0
0

1
a0 =
2n Γ(n + 1)
porque de esta forma tanto el exponente de Xcomo el de 2 se haceniguales; además los
factores (n + 1)...(n + k) se combinan con Γ(n + 1) para dar Γ(n + k + 1). Por consiguiente,
la solución particular, y , que se busca es:
1
 (2)
∞
X xn+2k
y1 = Jn (x) = (−1)k , (n ≥ 0)
2n+2k k!Γ(n + k + 1)
Esta es la denición de la función Besseñ de orden (n ≥ 0). La serie converge uniforme-
k=0

memnte en cualquier intervalo nito:

 n+2(k+1) 1
 
a2(k+1)  xn+2(k+1)
2 )k+1)!Γ(n+k+2) 
| |= 1
a2k xn+2k

 2n+2kk!Γ(n+k+1) 
n+2k
 
 2 k!Γ(n + k + 1)   2
=  n+2k 2  x 
2 2 (k + 1)k!(n + k + 1)Γ(n + k + 1) 
 
 1 
 −→k→∞ 0
=  2
2 (k + 1)(n + k + 1) 
También hay convergencia absoluta para cualquier X nita, puesto que:
 
 1  1
 2n+2k k!Γ(n + k + 1)  < k!
 

para toda K y toda X. Esta es la prueba M de Weierstrass. Por consiguiente la serie que
dene a J (x) es uniformememnte convergente.
Ejemplo: Calcular J (x) Se ha obtenido:
n
1/2

(3)
∞ n+2k
X x k
J (x) =
n (−1) n+2k
k=0
2 k!Γ(n + k + 1)

Por consiguiente: ∞
X x2k+1/2
J1/2 (x) = (−1)k
22k+1/2 k!Γ( 2k+2
2
)
Dentro de la sumatoria multiplíque por (x/2) /2 y fuera de ella por (2/x) /2 Ahora use la
k=0
1 1

propiedad recursiva de la función gamma para obtener:

(4)
∞ 2k+1
X
1/2 x
k
J (x) = (2/x)
1/2 (−1) k+1 2k+1 1
2 k!( )...( )Γ(1/2)2 2

En el denominador de ecuación (4), después del factorial y antes de Γ el número 2 aparece

k=0

2k+1
veces, por tanto:
∞
1/2
X x2k+1
J1/2 (x) = (2/πx) (−1)k
k=0
2k k!(2k + 1)..,1
Ahora se combinan 2 y k! para obtener: 2k(2k − 2)..,2. Este resultado se intercala con los
k

factores: (2k + 1)(2k − 1)..,1 para producir: (2k + 1)1!, tal que:
r ∞
2 X x2k+1
J1/2 (x) = (−1)k
πx k=0 (2k + 1)!
Esta fórmula es la serie del seno:
J1/2 (x) = (
2 1/2
πx
) senx (5)
 Fórmulas diferenciales con J0 y J1
De J (x) para n = 0, 1, se tiene:
n

x2 x4 k x2k
J0 (x) = 1 − + + ... + (−1) + ...
22 24 (2!)2 22k (k!)2

x x3 x2k+1
J1 (x) = − 3 + ... + (−1)k 2k+1 + ...
2 2 2! 2 k!(k + 1)!
Nótese que la derivada del segundo término de J (x) es el primer término de J (x) con
signo contrario, etc. En general la derivada del término que contiene a:
0 1

d k+1 x2k+2 k (2k + 2)x

2k+1
(−1) = −(−1) 2k+2
dx 22k+2 [(k + 1)!]2 2 [(k + 1)!]2
2(k + 1)x2k+1
= −(−1)k
2x 22k+1 (k + 1)k!(k + 1)!
x2k+1
= −(−1)k
22k+1 k!(k + 1)!
que es el simétrico (negativo) del término que contiene a x 2k+1
en J (x) observe que la
derivada del primer término de J (x) es cero. Por consiguiente:
1
0

J00 (x) = −J1 (x) (6)

Ahora multiplique cada término de J (x) por x y derive:
1

d x3 k x
2k+1
[xJ1 (x)] = x − 2 + ... + (−1) 2k + ...
dx 2 2 (k!)2
es decir:
d
dx
[xJ1 (x)] = xJ0 (x) (7)
Las raíces de J (X) y J (x)
0 1

Un teorema muy importante del algebra establece que una ecuación de grado n:
a0 xn + a1 xn−1 + ... + an−1 x + an = 0 (a 6= 0) (8)
tiene n-raices, reales o complejas. Suponga ahora que n → ∞, entonces se podría pensar que
el número de raíces tiende a innito, de acuerdo al teorema mencionado. ¾es esto siempre
cierto? De hecho puede resultar que ecuación (8) no tenga ninguna. Por ejemplo, la ecuación:
x2 x3 xn
1+x+ + + ... + + ... = 0
2! 3! n!
no tiene ninguna raiz, ya que esta ecuación es:
ex = 0

y no existe ningún número, real o complejo que la satisfaga.

Cuando se iguala a cero una serie innita para formar una ecuación, resulta en general muy
dicil determinar la naturaleza de sus raices, o encontrarlas aun sabiendo su naturaleza.
La naturaleza de las raices de J (x) = 0 y J (x) = 0 se conoce y existen tablas de estas raices,
las cuales son muy útiles en varios problemas de la física: Estas dos ecuaciones no tienen
0 1

raices complejas y cada una posee un número innito de raices reales diferentes.
Entre dos raices consecutivas de J (x) = 0 cae una y sólo una raiz de J (x) = 0, y entre dos
raices consecutivas de J (x) = 0 sólo hay una raiz de J (x) = 0, es decir, las raices de las dos
0 1

ecuaciones separan a una de la otra. En cada ecuación la diferencia entre dos raices conse-
1 0

cutivas tiende a π, cuando x → ∞. La interpretación de este hecho es decirque las funciones

J (x) y J (x) son "quasi-periódicas"de periodo casi 2π . En otras palabras estas dos funciones
se comportan aproximadamente como las funciones cosx y senx.
0 1

Las funciones J (x) se pueden calcular truncando la serie que las dene, o a través de aproxi-
maciones por medio de polinomios, como se ilustró durante el estudio de la función Gamma.
n

En particular las funciones J (x) y J (x) se pueden aproximar (−3 ≤ x ≤ 3) por:

0 1

x x x
J0 (x) = 1 − 2,2499997( )2 + 1,2656208( )4 − 0,3163866( )6 +
3 3 3
x 8 x 10 x 12
... + 0,0444479( ) − 0,0039444( ) + 0,0002100( ) +
3 3 3
−8
... + ε...|ε| < (5)10

x x
J1 (x) = x{1/2 − 0,56249985( )2 + 0,21093573( )4 +
3 3
x 6 x 8 x
... − 0,03954289( ) + 0,0044319( ) − 0,00031761( )10 +
3 3 3
x 12
... + 0,00001109( ) + ε ... |ε| < (1,3)10−8
3
Mayor precisión se obtiene con la aproximación dada por los programas en FORTRAN:
BESSEL0 [calcula J (x)] y BESSEL1 [calcula J (x)], para cualquier x real. Para el cálculo
numérico de otras funciones Bessel (n, entero >1) basta utilizar las propiedades recursivas
0 1

(ver mas adelante)

 PROGRAMA BESSEL0

WRITE(*,*) ' DAME EL VALOR DE X'

READ(*,*) X
XTEMP=BESSJ0(X)
WRITE(*,*)X,XTEMP
END
FUNCTION BESSJ0(X)
REAL*8 Y,P1,P2,P3,P4,P5,Q1,Q2,Q3,Q4,Q5,R1,R2,R3,R4,R5,R6
* S1,S2,S3,S4,S5,S6
DATA P1,P2,P3,P4,P5/1.DO,-.1098628627D-2,.2734510407D-4,
* -.2073370639D-5,.2093887211D-6/, Q1,Q2,Q3,Q4,Q5/-.1562499995D-
*1,
* .1430488765D-3,-.6911147651D-5,.7621095161D-6,-.934945152D-7/
DATA R1,R2,R3,R4,R5,R6/57568490574.D0,-13362590354.D0,651619640.7D
*0,
* -11214424.18D0,77392.33017D0,-184.9052456D0/,
* S1,S2,S3,S4,S5,S6/57568490411.D0,1029532985.D0,
* 9494680.718D0,59272.64853D0,267.8532712D0,1.D0/
IF(ABS(X).LT.8.)THEN
Y=X**2
BESSJ0=(R1+Y*(R2+Y*(R3+Y*(R4+Y*(R5+Y*R6)))))
* /(S1+Y*(S2+Y*(S3+Y*(S4+Y*(S5+Y*S6)))))
ELSE
AX=ABS(X)
Z=8./AX
Y=Z**2
XX=AX-.785398164
BESSJ0=SQRT(.636619772/AX)*(COS(XX)*(P1+Y*(P2+Y*(P3+Y*(P4+Y
* *P5)))))-Z*SIN(XX)*(Q1+Y*(Q2+Y*(Q3+Y*(Q4+Y*Q5)))))
ENDIF
RETURN
END
 PROGRAMA BESSEL1

WRITE(*.*) ' DAME EL VALOR DE X'

READ(*.*) X
XTEMP=BESSJ1(X)
WRITE(*,*)X,XTEMP
END
FUNCTION BESSJ1(X)
REAL*8 Y,P1,P2,P3,P4,P5,Q1,Q2,Q3,Q4,Q5,R1,R2,R3,R4,R5,R6,
* S1,S2,S3,S4,S5,S6
DATA R1,R2,R3,R4,R5,R6/72362614232.D0,-789059235.D0,242396853.1D0
*,
* -2972611.439D0,15704.48260D0,-30.16036606D0/,
* S1,S2,S3,S4,S5,S6/144735228442.D0,2300535178.D0,
* 18583304.74D0,99447043394D0,376.9991397D0,1.D0/
DATA P1,P2,P3,P4,P5/1.D0,.183105D-2,-.3516396496D-4,.247520174D-5
*,
* -.240337019D-6/, Q1,Q2,Q3,Q4,Q5/.04687499995D0,-.2002690873D-3
*,
* .8449199096D-5,-.88228987D-6,.105787412D-6/
IF(ABS(X).LT.8. )THEN
Y=X**2
BESSJ1=X*(R1+Y*(R2+Y*(R3+Y*(R4+Y*(R5+Y*R6)))))
* /(S1+Y*(S2+Y*(S3+Y*(S4+Y*(S5+Y*S6)))))
ELSE
AX=ABS(X)
Z=8./AX
Y=Z**2
XX=AX-2.356194491
BESSJ1=SQRT(.636619772/AX)*(COS(XX)*(P1+Y*(P2+Y*(P3+Y*(P4+Y
* *P5))))-Z*SIN(XX)*(Q1+Y*(Q2+Y*(Q3+Y*(Q4+Y*Q5)))))
* *SIGN(1.,X)
ENDIF
RETURN
END

Un ejercicio interesante es la traducción de estos dos programas al APL; también se pue-

de practicar con el primer polinomio dado arriba.
 Cuadro (1)
Valores de las funciones Bessel de orden 0, 1 y 2
(Fuente: Handbook of Mathematical Functions - ver capitulo 1)
x J0 (x) J1 (x) J2 (x)
0.0 1.00000 00000 00000 0.00000 00000 0.00000 00000
0.1 0.99750 15620 66040 0.04993 75260 0.00124 89587
0.2 0.99002 49722 39576 0.09950 08326 0.00498 33542
0.3 0.97762 62465 38296 0.14831 88163 0.01116 58619
0.4 0.96039 82266 59563 0.19602 65780 0.01973 46631
0.5 0.93846 98072 40813 0.24226 84571 0.03060 40235
0.6 0.91200 48634 97211 0.28670 09881 0.04366 50967
0.7 0.88120 08886 07405 0.32899 57415 0.05878 69444
0.8 0.84628 73527 50480 0.36884 20461 0.07581 77625
0.9 0.80752 37981 22545 0.40594 95461 0.09458 63043
1.0 0.76519 76865 57967 0.44005 05857 0.11490 34849
1.1 0.71962 20185 27511 0.47090 23949 0.13656 41540
1.2 0.67113 27442 64363 0.49828 90576 0.15934 90183
1.3 0.62008 59895 61509 0.52202 32474 0.18302 66988
1.4 0.56685 51203 74289 0.54194 77139 0.20735 58995
1.5 0.51182 76717 35918 0.55793 65079 0.23208 76721
1.6 0.45540 21676 39381 0.56989 59353 0.25696 77514
1.7 0.39798 48594 46109 0.57776 52315 0.28173 89424
1.8 0.33998 64110 42558 0.58151 69517 0.30614 35353
1.9 0.28181 85593 74385 0.58115 70727 0.32992 57277
2.0 0.22389 07791 41236 0.57672 48078 0.35283 40286
2.1 0.16660 69803 31990 0.56829 21358 0.37462 36252
2.2 0.11036 22669 22174 0.55596 30498 0.39505 86875
2.3 0.05553 97844 45602 0.53987 25326 0.41391 45917
2.4 +0.00250 76832 97244 0.52018 52682 0.43098 00402
2.5 -0.04838 37764 68198 0.49709 41025 0.44605 90584
2.6 -0.09680 49543 97038 0.47081 82665 0.45897 28517
2.7 -0.14244 93700 46012 0.44160 13791 0.46956 15027
2.8 -0.18503 60333 64387 0.40970 92469 0.47768 54954
2.9 -0.22431 15457 91968 0.37542 74818 0.48322 70505
3.0 -0.26005 19549 01933 0.33905 89585 0.48609 12606
3.1 -0.29206 43476 50698 0.30092 11331 0.48620 70142
3.2 -0.32018 81696 57123 0.26134 32488 0.48352 77001
3.3 -0.34429 62603 98885 0.22066 34530 0.47803 16865
3.4 -0.36429 55967 62000 0.17922 58517 0.46972 25683
3.5 -0.38012 77399 87263 0.13737 75274 0.45862 91842
3.6 -0.39176 89837 00798 0.09546 55472 0.44480 53988
3.7 -0.39923 02033 71191 0.05383 39877 0.42832 96562
3.8 -0.40255 64101 78564 +0.01282 10029 0.40930 43065
3.9 -0.40182 60148 87640 -0.02724 40396 0.38785 47125
4.0 -0.39714 98098 63847 -0.06604 33280 0.36412 81459
4.1 -0.38866 96798 35854 -0.10327 32577 0.33829 24809
4.2 -0.37655 70543 67568 -0.13864 69421 0.31053 47010
4.3 -0.36101 11172 36535 -0.17189 65602 0.28105 92288
4.4 -0.34225 67900 03886 -0.20277 55219 0.25008 60982
4.5 -0.32054 25089 85121 -0.23106 04319 0.21784 89837
4.6 -0.29613 78165 74141 -0.25655 28361 0.18459 31052
4.7 -0.26933 07894 19753 -0.27908 07358 0.15057 30295
4.8 -0.24042 53272 91183 -0.29849 98581 0.11605 03864
4.9 -0.20973 83275 85326 -0.31469 46710 0.08129 15231
5.0 -0.17759 67713 14338 -0.32757 i91376 0.04656 i51163
(−4)6 (−4)5 (−4)3
h i h h
11 8 7

Los programas anteriores dan el resultado con siete cifras decimales para J (x) y J (x)
Las primeras cindo raices de J (x) = 0, J (x) = 0 son (con 4 cifras decimales):
0 1
0 1
 J0 (x) = 0 J1 (x) = 0
Raíz Diferencia Raíz Diferencia
2.4048 3.8317
3.1153 3.1839
5.5201 7.0156
3.1336 3.1579
8.6537 10.1735
3.1378 3.1502
11.7915 13.3237
3.1394 3.1469
14.9309 16.4706
Note que las diferencias de las raices de J (x) = 0 se acercan a π por debajo, y las de
J (x) = 0, por arriba. Las raices de J (x) = 0 están entre las raices de J (x) = 0 para formar
0

una serie creciente. Un resultado interesante es que la suma de los reciprocos al cuadrado de
1 1 0

las raices positivas de J (x) = 0 vale 1/4:

0

1
(2,4048)−2 + (5,5201)−2 + (8,6537)−2 + ... =
4
positivas de cosx = 0:  π −2  −2  −2
3π 5π 1
+ + + ... =
2 2 2 2
2
 
−2 −2 −2 −2 π
= 4π (1 + 3 + 5 + ...) = (4π )
8
 Dos integrales que involucran a J0 (x) y J1 (x)

La ecuación de Bessel de orden n es:

x2 y 00 + xy 0 + (x2 − n2 )y = 0 ...(1)
y su solución es y = J (x) Sea ax = x ∴ y = y /a , y
n
0 0 00
= y 00 /a2 ∴

x2 y 00 + xy 0 + (a2 x2 − n2 )y = 0 (9)
cuya solución, por tanto, es: y = J (ax) , para n = 0, se tiene:
n

xy 00 + y 0 + a2 xy = 0 → y = J0 (ax)

Ahora se pueden obtener resultados interesantes y útiles al considerar dos de estas ecuaciones
de Bessel de orden cero y sus soluciones:
00 0
xu + u + α xu = 0 → u = J (αx) 2
o (10)
00 0
xv + v + β xv = 0 → v = J (βx) 2
0 (11)
Ahora multiplique a la ecuación (11) por u y a la ecuacion (10) por v, enseguida reste los
resultados:
u(xv 00 + v 0 + β 2 xv) − v(xu00 + u0 + α2 xu) = 0 ∴
(β 2 − α2 )xuv = x(u00 v + uv 00 ) + (u0 v − uv 0 )
el lado derecho de esta ecuación es una derivada exacta:
d
[x(u0 v − uv 0 )] = x(u0 v 0 + u00 v − uv 00 − u0 v 0 ) + u0 v − uv 0 =
dx
= x(u00 v − uv 00 ) + (u0 v − uv 0 )
tal que la integración entre x = 0 y x = 1 da:
Z 1
1
2
(β − α ) 2
xuv dx = |x(u0 v − uv 0 )|0 = |u0 v − uv 0 |x=1
0

d
yaqueu = J0 (αx) u0 = J0 (αx)
dx
d d(αx)
u0 = J0 (αx) = −αJ1 (αx)
d(αx) dx
análogamente:
v = Jo (βx) → v 0 = −βJ1 (βx) ∴

(12)
1
−αJ1 (α)J0 (β) + βJ0 (α)J1 (β)
Z
xJ0 (αx)J0 (βx)dx =
β 2 − α2
Si α y β son raices distintas de J (x) = 0, es decir: J (α) y J (β) = 0, entonces la última
0

ecuación se reduce a:
0 o 0

Z 1
xJ0 (αx)J0 (βx)dx = 0 (α = β)
0
 Cuando β = α la expresión (12) queda indeterminada (0/0). Se puede calcular su límite
para ello considere que α es una raiz de J (x) = 0 y β una variable que tiende a α, entonces
el cociente del lado derecho de ecuación (12) es:
0

αJ1 (α)J0 (β)

−
β 2 − α2
Ahora se aplica la regla de Hôpital, diferenciando numerador y denominador con respecto
a β:
αJ1 (α)J0 (β) αJ1 (α)J1 (β)
lı́m − 2 2
= lı́m −
βtoα β −α βtoα 2β

= 1/2J12 (α) ∴
Z 1
xJ02 (αx)dx = 1/2J12 (α)

Desarrollo de f (x) en términos de la función Bessel de orden cero, J (x).

0

Sean α , α , ..., α , ... las raices positivas de J (x) = 0. Un gran número de funciones pueden
0

desarrollarse como:
1 2 n 0

f (x) = A J (α x) + ... + A J (α x) + ... (0 < x < 1)

1 0 1 n 0 n (13)
En esta clase de problemas hayq ue encontrar los coecientes de .A".
Para lograrlo se multiplica la ecuación (13) por xJ (α x) y enseguida se integra de o a 1, con
la hipótesis de que la integral es válida. Por el resultado anterior, todas las integrales del lado
0 n

derecho de ecuación (13) se hacen cero, excepto la única que contiene al factor de A , es decir: n
Z 1 Z 1
An 2
xJ0 (αn x)f (x)dx = An xJ02 (αn x)dx = J (αn )
0 0 2 1
por ende:
(14)
Z 1
2
An = 2 xJ0 (αn x)f (x)dx
J1 (αn )
La ecuación (13) da el desarrollo de f (x) y la fórmula (14) los coecientes.
0

Ejemplo:
Desarrollar el número 1 en una serie de funciones de Bessel de orden cero.
En este ejemplo:
Z 1
2
f (x) = 1, An = 2 xJ0 (αn x)dx
J1 (αn )
La integral se puede calcular, usando una d elas fórmulas vistas, cambiando a x por α x:
0

d [αn xJ1 (αn x)] = αn xJ0 (αn x) d(αn x)dx

 
xJ1 (αn x)
d = zJ0 (αn x)dx , porconsiguiente :
αn
 xJ1 (αn x) 1
 
2 2
An =   =
J1 r(αn )  αn 
0 αn J1 (αn )
La serie correspondiente para el número 1 es:
1=
2J0 (α1 x) 2J0 (α2 x)
+
α1 J1 (α1 ) α2 J1 (α2 )
+ ... +
2J0 (αn x)
αn J1 (αn )
+ ... (15)
Esta serie es válida para 0 < x < 1; sin embargo, no es válida para x = 1, ya que se reduce
a cero, porque todos los numeradores son raices de la ecuación J (x) = 0 (Todas las α son
sus raices positivas), es decir, existe una contradiccioón:
0

1 = 0 + 0 + ... + 0 + ... = 0 :
pero para x = 0 la serie obtenida sí es válida. Como J (0) = 1, entonces:
0

∞
X 1
1=2
n=1
αn J1 (αn )
para encontrar el resultado:
1 1
= displaystylesum∞
n=1
2 αn J1 (αn )
 La función Bessel J−n (x)

Si cambiamos n por -n en la denición que se dió de la función J (x) [ecuación (3)] se

tiene:
n

∞
X x−n+2k
J−n (x) = (−1)k
2−n+2k k!Γ(−n + k + 1)
Cuando n no es un entero positivo la serie converge para todos los valores reales de x,
k=0

excepto cero y dene la funcion J (x).

Cuando ntoN > 0, entonces Γ(−n + k + 1) se vuelve innita para k = 0, 1, 2, ..., (N − 1); por
−n

consiguiente estos términos tienden a cero. La parte diferente de cero principia cuando k = N
y toma la forma:
∞
x−N +2k xN
X 
k
(−1) −N +2k = (−1)N +
k=N
2 k!Γ(−N + k + 1) 2N N !

xN +4

x
− N +2 + − ... = (−1)N JN (x)
2 (N + 1)! 2N +4 (N + 2)!2!
la cual dene el valor de J(x), donde N es un entero positivo; por ende, la denición
completa de la función de Bessel de orden negativo es:
∞
X x−n+2k
J−n (x) = (−1)k
k=0
2−n+2k k!Γ(−n + k + 1)
(0 ≤ n 6= 1, 2, 3, ...)

J−N (x) = (−1)N JN (x) (N = 1, 2, 3, ...)

J−1/2 (x)al igual que J (x) se puede expresar en forma nita.
Ponga en la función J (x):
1/2
n = 1/2 −n

∞ −1/2+2k
−1/2+2k k! Γ 2k+1
( )
X
J−1/2 (x) = (−1)k x 2 2

k=0
∞
X x2k
= (2/x)1/2 (−1)k  
2k−1 2k−3 Γ)(1/2)
k=0 (22k k! 2 2
... 2

∞
X x2k
= (2/πx)1 /2 (−1)k
k=0
2k k! (2k − 1)(2k − 3)..,1
∞
X
1 x2k
= (2/πx) /2 (−1)k = (2/πx)1 /2cosx
(2k)!
La función J cuando se substituye por y en la ecuación diferencial de Bessel la satisface,
k=0

obviamente, porque J (x) diere de la solución J (x) sólo en el signo de n y ésta aparece
−n

solamente al cuadrado en la ecuación diferencial de Bessel. Así, cuando n diferente de entero

−n n

y de cero, se tienen dos soluciones particulares independientes:

y1 = Jn (x) , y2 = J−n (x)
 de la ecuación diferencial de Bessel de orden n. Por tanto, la solución general es:
y = AJn (x) + bJ−n (x) (16)
donde A, B son constantes arbitrarias. Así, por ejemplo, la solución general de la ecuación
diferencial de Bessel de orden 1/2:
x2 y 00 + xy 0 + (x2 − 0,25)y = 0 es :

y = AJ1/2 (x) + BJ−1/2 (x)

= A(2/πx)1/2 senx + B(2/πx)1/2 cosx

senx cosx
= C1 √ + C2 √
x x
Cuando n = 0, entonces:
y = (A + B)J (x) = CJ (x), C constante arvitraria;
0 0

cuando n es un entero positivo (N > 0), entonces:

y = [A + (−1) B]J (x):
N

por consiguiente, cuando n es cero o un entero la ecuación (16) no d ala solución general. No
N

son sucientes las funciones J (x), J (x) para obtenerla. Las funciones que se han estudia-
do hasta ahora se llaman funciones Bessel de la primera clase : El problema en cuentión se
n −n

resuelve con las funciones de Bessel de la segunda clase, las cuales se verán más adelante.
 Las funciones ber y bei

Considere la ecuación diferencial que frecuentemente aparece en problemas de la física:

00 0
p
xy + y − ixy = 0 (donde i = (−1) (17)
que tiene la forma de la ecuación que ya se estudió:
xy 00 + y 0 + a2 xy = 0 → y = J0 (ax)
por tanto a 2
=i , tal que una solución particular es:
y = J0 (ax) = J0 [(−1)1/2 x] = J0 (i3/2 x)
(17)
Substituyendo (i 3/2
x) por x en la serie que se obtuvo para J (x),
Se tiene:
0

i 3 x2 i6 x4 i9 x6
J0 (i3/2 x) = 1 − + − + ...
22 24 (2!)2 26 (3!)2
ix2 x4 ix6
=1+ 2 − 4 − + ...
2 2 (2!)2 26 (3!)2
x4 x8
 
1 − 2 2 + 2 2 2 2 − ... +
24 2 4 6 8
x2 x6 x1 0
 
+i 2 − 2 2 2 + 2 2 2 2 2 + ...
2 2 4 6 2 4 6 8 10
" ∞
#
4k
X x
= 1+ (−1)k 2 2 2 +
k=1
2 4 6 ... (4k)2
" ∞ #
4k−2
X x
+ i (−1)k 2 2 2
k=1
2 4 6 ... (4k − 2)2

Por consiguiente para x real J (i x) es compleja. Este resultado se usa para denir las
3/2

funciones ber (Bessel - real) y bei (Bessel - imaginaria) como las series:
0

" ∞
#
X x4k
1+ (−1)k 2 2 2
2 4 6 ...(4k)2
y
k=1

" ∞
#
X x4k−2
i (−1)k 2 2 2
k=1
2 4 6 ... (4k − 2)2

respectivamente, tal que:

J0 (i3/2 x) = ber x + i bei x
Ahora se derivarán dos fórmulas que conectan a las funciones
 ber y bei:
Sí al segundo término de la serie ber, a saber: − primero se diferencia y despuésx4
22 42
se multiplica por x, el resultado es: , resultado que también se obtiene, pero con signo
x4

opuesto, sí al primer término de la serie bei, a saber , primero se multiplica por x y enseguida
21 42
x2

se integra de 0 a x. Relaciones similares se obtienen entre el tercer término de ber y el segundo

22

de bei; para así llegar a:

∞
0
X x4k
x ber x = (−1)k =
k=1
22 42 62 ...(4k − 2)2 (4k)
Z x
=− x(bei x)dx

en otras palabras:
0

d
(x ber0 x) = −x bei x
dx
Note que resultados idénticos se obtienen multiplicando cualquier término de la serie ber
por x y después integrando de 0 a x, y al diferenciar el término correspondiente de la serie
bei y después multiplicar por x. El resultado es:
x ∞
x2 X x4k+2
Z
x(ber x)dx = + (−1)k 2 2 2
0 2 k=1
2 4 6 ...(4k)2 (4k + 2)
∞
X x4k−2
=− (−1)k = x bei0 x
22 42 62 ...(4k − 4)2 (4k − 2)
por tanto:
k=1

d
(x bei0 x) = x ber x
dx
{ Al multiplicar por x se obtiene:
∞
X x4k+1
x(ber x) = x + (−1)k
22 42 62 ...(4k)2
integrando:
k=1

x ∞
x2 X x4k+2 (4k + 2)
Z
x(ber x)dx = + (−1)k 2 2 2
0 2 2 4 6 ...(4k)2 (4k + 2)2
El último se consiguió trayendo dentro de la sumatoria al primer término de la serie,
k=1

2
(x /2) }
Ejemplo:
Demostrar que:
Z a
x (ber2 x) + (bei2 x) dx =
 
0

a [(ber a)(bei0 a) − (bei a)(ber0 a)]

Se ha visto que:
d
x (ber x) = (x[bei0 x])
dx
 d
(bei x) = − (x[ber0 x])
dx
Integrando por partes:
Z a Z a
2
x(ber x)dx = [x (ber x)] (ber x)dx
0 0
Z a
d
= [x (bei0 x x)](ber x)dx
0 dx
Z a
0
= a(bei a)(ber a) − x(bei0 x)(ber0 x)dx
0
Z a Z a
x(bei2 x)dx = [x(bei x)] (bei x)dx
0 0
Z a
d
=− [x(ber0 x)](bei x)dx
0 dx
Z a
0
= −a(ber a)(bei a) + x(ber0 x)(bei0 x)dx

por consiguiente ya que las integrales se cancelan, se tiene:

0

Z a
x (ber2 x) + (bei2 x) dx =
 
0

a [(ber a)(bei0 a) − (bei a)(ber0 a)]

 Funciones Bessel de segunda clase

La función de Neumann se dene como:

Nn (x) =
cos(nπ)Jn (x) − J−n (x)
sen(nπ)
(18)
Para n no entero satisface la ecuación diferencial de Bessel, ya que es una combinación
lineal de dos soluciones conocidas: J (x) y J (x), Sin embargo, para n entero, la ecuación
(18) queda indeterminada. La denición de N (x) se eligió de forma deliberada para que po-
n −n

seyera estas dos propiedades.

Evaluando a N (x) mediante la regla de L'Hòpital se obtiene:
d
dx
[cos(nπ)Jn (x) − J.n (x)]
N (x) = d
dx
sen(nπ)

n (x) (x) 
−πsen(nπ)Jn (x) + cos(nπ) zJzn − zJ−n
zn
=

πcos(nπ)


n=ent.
 
1 zJn (x) zJ−n (x) 
= − (1)n
π zn zn n=entero
Lo anterior conduce a una expresión que depende de la función digamma (z) y el loga-
ritmo de x, El cálculo es relativamente sencillo; pero muy laborioso (ponerlo aqui).
En esta forma la solución general de la ecuación diferencial de Bessel es:
y = AJn (x) + BNn (x)
A y B son constantes arbitrarias.
Las funciones Neumann satisfacen las relaciones de recursión que tienen las funciones Bessel:
 Relaciones Recursivas

d n
dx
[x Jn (x)] = xn Jn−1 (x) (19)
d  −n
dx
x Jn (x) = −x−n Jn+1 (x)

(20)
Jn−1 (x) + Jn+1 (x) =
2n
x
Jn (x) (21)
Jn−1 (x) − Jn+1 (x) = 2Jn0 (x) (22)
n n
Jn0 (x) = − Jn (x) + Jn−1 (x) = Jn (x) − Jn+1 (x)
x x
(23)
0
Jn (x)J−n (x) − J−n (x)Jn0 (x) = −
2πsen(nπ)
nx
(24)
Jn (x)Nn+1 (x) − Jn+1 (x)Nn (x) = −
2
πx
(25)
Demostración de la propiedad (19):
Primero se multiplica por x y después se deriva:
n

∞
d n d X x2n+2k
[x Jn (x)] = (−1)k n+2k
dx dx k=0 2 k!Γ(n + k + 1)
∞
X (2n + 2k)x2n+2k−1
= (−1)k n+2k
k=0
2 k!Γ(n + k + 1)
como se sabe:
Γ(n + k + 1) = (n + k)Γ(n + k) ∴

∞
d n X x2n+2k−1
[x Jn (x)] = (−1)k n+2k−1
dx k=0
2 k!Γ(n + k)
∞
n
X xn+2k−1
=x (−1)k = xn Jn−1 (x)
k=0
2n+2k−1 k!Γ(n + k)
 Ecuación diferencial general con soluciones Funciones Bessel

a2 − n 2 c 2
 
001 − 2a 0 α−1 2
y + y + (bcx ) + y=0
x x2
tiene la solución:
y = xZn (bxc )
donde Z es J ó N ó cualquier combinación de ellas, a, b, c, n son constantes.
Ejemplo:
y 00 + 9xy = 0
Por tanto:
1 − 2a = 0, (bc)2 = 0, 2(c − 1) = 1, a2 − n2 c2 = 0 ∴

a = 1/2, b2 = 9/c2 , c = 3/2 to b = 2, n = a/c = 1/3

Por consiguiente:
y = x1/2 Z1/3 (2x3/2 )
 

Por ende, la solución general es:

y = x1/2 AJ1/3 (2x3/2 ) + BN1/3 (2x3/2 )
 

donde A, B son constantes arnitrarias.

Como una aplicación de las funciones Bessel considere de nuevo el problema del péndulo sim-
ple, con la siguiente caracteristica:
su longitud 1 ahora crece con velocidad constante.
Encontrar la ecuación de movimiento y la solución para osciladores pequeñas.
La energa cintica es : T = 1/2 m (l˙2 + l2 Θ2 )
La energa potencial es : V = −mgl cosΘ,
entonces:
L = T − V = 1/2m (l˙2 + l2 Θ2 ) + mgl cosΘ
La ecuación Lagrangiana de movimiento para la variable Θ es:
d
(ml2 Θ) + mgl senΘ = 0
dt
Por consiguiente:
 dΘ dl
ml2 + 2ml Θ + mgl senΘ = 0 ∴
dt dt

l Θ + 2lΘ + gsenΘ = 0
Sea l la longitud de la cuerda, en el tiempo t : l = l + vt, donde v es la constante de
crecimiento, entonces: l˙ = v Para el caso de oscilaciones pequeñas senΘ puede sustiruirse por
0

Θ, por tanto:

(l0 + vt)Θ + 2vΘ + gΘ = 0

sea: l0 + vt

l0

x= =
v v
entonces: dθ dθ
dx = dt, Θ = = = Θ0 ∴ Θ = Θ00
entonces:
dt dx

vxΘ + 2vΘ + gΘ = 0
ó
2 g θ
Θ00 + Θ0 + =0
x v x
Esta última ecuación tiene la forma de la ecuación diferencial general con soluciones Fun-
ciones de Bessel.
Por tanto:
g
1 − 2a = 2, a2 − n2 c2 = 0, (bc)2 = , 2c − 2 = −1
v
a2 1p p
c = 1/2, a = −1/2, n2 = = 1, b = g/v = 2 g/v
c2 c
Por ende:
 p 
−1/2 1/2
Θ=x Z1 2 g/vx
Para simplicar la notación ponga:
p
u = 2 g/vx1/2
entonces, por lo visto anteriormente, se tiene:
Θ = A u−1 J1 (u) + Bu−1 N1 (u)

usando la propiedad diferencial que se demostró se encuentra dΘ/du:

dΘ
= − A u−1 J2 (u) + B u−1 N2 (u)
 
du
Las constantes A y B se encuentran a partir de las condiciones iniciales; por ejemplo,
cuando se estudió el péndulo simple clásico con oscilaciones pequeñas se puso:
Θ = Θ0 , Θ = 0 , para t = 0
y se vió que la solución general es:
Θ = A coswt +; B senwt

Θ0 = a cos0 + B sen0 → A = Θ0
Θ = −A wsenwt + B coswt = 0 → B = 0, por tanto :
Θ = Θ0 coswt
En el problema bajo discusión se pueden tomar estas mismas condiciones iniciales:
l0
Θ = Θ̇0 , Θ = 0 , para t = 0; x =
v
p p p 2p
u0 = 2 g/v x1/2 = 2 g/v l0 /v = gl0
v
 1/2
du p −1/2 g/v p
= (1/2)(2) g/v x = = g/l0
dx l0 v
dΘ
dΘ p p
= dx
du
= Θ0 l0 /g = Θ̇ l0 /g = 0
du
entonces para t = 0:
dx

Θ0 = A u−1 −1
0 J1 (u0 ) + B u0 N1 (u0 )

0 = A J2 (u0 ) + B N2 (u0 )
por la propiedad recursiva:
d
u−1 J1 (u) = −u−1 J2 (u)


du
Se puede probar que: 2
J1 (u)N2 (u) − J2 (u)N1 (u) = −
πu
Ahora use las relaciones de recursión, la denición de Nn (x) y la propiedad:
0 2
Jn (x)N−n (x) − Nn (x)Jn0 (x) = − sen(nπ)
πx
Jn0 (x)J−n (x) − Jn (x)J−n
0
(x)
=
sen(nπ)
para obtener:
u0 Θ0 N2 (u0 ) = A[J1 (u0 )N2 (u0 ) − N1 (u0 )J2 (u0 )]
 
2
=A −
πu0
u0 Θ0 J2 (u0 ) = B[N1 (u0 )J2 (u0 ) − N2 (u0 )J1 (u0 )]
 
2
=B
πu0
Por consiguiente:
πu20 πu20
A=− Θ0 N2 (u0 ), B= Θ0 J2 (u0 )
2 2
Entonces la solución general es:
π 2
Θ= u Θ0 u−1 [−N2 (y0 )J1 (u) + J2 (u0 )N1 (u)]
2 0
esta solución se puede simplicar si se ajustan las constantes v, l , tal que:0

2p
u0 = gl0 es una raiz de J2 (u) = 0, entonces
v
B = 0 , por tanto :
 p 
−1 −1/2
Θ=Au J1 (u) = C x J1 2 gx/v , donde
Ap l0
C= v/g, x = + t
2 v
para este caso especial Θ̇ es un múltiplo de J (u). Así:
2

Θ = 0, corresponde a las raices de J1 (u) = 0

Θ̇ = 0, corresponde a las raices de J2 (u) = 0
Un cuarto de periodo correspondiente al tiempo desde Θ = 0 hasta Θ̇ = 0, las cuales se
pueden encontrar a partir de las raices de J (u) = 0 y J (u) = 0
1 2
 Función Generadora

Proceda a introducir una función de dos variables:

g(x, t) = e1/2x (t − 1/t)

Esta función se puede desarrollar en una serie de Laurent, para obtener:

∞
X
1/2x (t − 1/t)
e = Jn (x) tn
n=0

, que constituye el coeciente de t , es la función de Bessel de primera clase de orden n.

Jn (x) n

Desarrollando las exponenciales se tiene el producto de series de MacLaurin en xt/2, −x/2t,

respectivamente:
∞   ∞  x k t−k
1/2xt −1/2xt−1
X x r tr X
e e = (−1)k
r=0
2 r! k=0 2 k!
Para cierta k se tiene tn
(≥ 0) , apartir de:
x tn+k  x k t−k xn
(−1)k = n
2 (n + k)! 2 k! 2 n!
Por consiguiente, el coeciente de t es: n

∞
X xn+2k
Jn (x) = (−1)k
k=0
2n+2k k!Γ(n + k + 1)
De acuerdo a este procedimiento, el cual origina esta serie a partir de sus propiedades de
convergencia, es evidente que estas serie da inmediatamente J (x). Además si se substituye
un número real por uno complejo, Z , entonces se extiende la denición de la función Bessel.
−n
 Ejercicios para resolver
1) Derive:
 n
n n 1 d
Jn (x) = (−1) x J0 (x)
x dx
2) Demuestre:
 n
1 d
(xv Jv (x)) = xv−n Jv−n (x)
x dx
3) Demuestre:
 n  
1 d Jv (x) Jv+n (x)
= (−1)n
x dx xv xv+n
4) Demuestre:
 n
1 d sex(x)
q
n+1/2n
Jn+1/2 (x) = (−1) 2/πx
x dx x
(n = 0, 1, 2, ...)
5) Demuestre:
∞
X
1=2 J2n (x) + J0 (x)
n=1

6) Derive: Z ∞
1
e−ax J0 (bx)dx = √
0 a2 + b 2
7) Demuestre:
Z ∞
1
Jn (x) = cos(nΘ − xsenΘ)dΘ ... (n, entero)
π 0

8) Demuestre que entre dos raices reales consecutivas de J (x) = 0 existe una, y sólo una
raiz de J (x) = 0.
0
1