Dos sumas

18 views
Skip to first unread message

Dr. Wolfgang Hintze

unread,
Jan 11, 2013, 2:51:52 AM1/11/13
to
Hallar las sumas (a>-1)

s1 = Sum( 1/(n^2+a), (n,1,oo))
s2 = Sum((-1)^(n+1)/(n^2+a), (n,1,oo))

Saludos,
Wolfgang

Marko Riedel

unread,
Jan 11, 2013, 7:13:03 PM1/11/13
to
Estimados lectores,

aqu� vamos a usar el truco que consiste en hallar la integral de

g(z) = f(z) pi cot(pi z)

y de

g(z) = f(z) pi csc(pi z)

con f(z) siendo el t�rmino de la suma que queremos calcular sobre un
camino que consiste en una circunferencia C de radio R, con R -> oo. Es
un m�todo bastante cl�sico, por eso omito los detalles. Es muy facil ver
que en ambos casos la integral sobre R desaparece en el l�mite, con lo
cual los residuos a los polos deben sumar cero. (Desaparece porque
1/(n^2+a) es O(1/R^2) sobre C.)

Lo que estamos usando aqu� es esencialmente el hecho de que con m un
entero

Res_{z=m} pi cot(pi z) =
= (pi cos(pi z)/(pi cos(pi z)))_{z=m} = 1

y

Res_{z=m} pi csc(pi z) =
= (pi/(pi*cos(pi z)))_{z=m} = (-1)^m.

Primer caso:

Ponemos g(z) = 1/(z^2+a) pi cot(pi z). Luego (siendo "rq" la ra�z
cuadrada)

Res_{z=-rq(a) i} g(z) + Res_{z=+rq(a)i} g(z)
+ Res_{z=0} g(z) + 2 s1 = 0.

lo que da

(1/(2z) pi cot(pi z))_{z=-rq(a)i}
+ (1/(2z) pi cot(pi z))_{z=+rq(a)i}
+ 1/a + 2 s1 = 0

con lo cual

1/(2 rq(a) i) pi cot(pi rq(a) i)
+ 1/(2 rq(a) i) pi cot(pi rq(a) i)
+ 1/a + 2 s1 = 0

y asi finalmente

s1 = - 1/(2a) + 1/(2 rq(a)) pi coth(pi rq(a)).

Segundo caso.

Ponemos g(z) = - 1/(z^2+a) pi csc(pi z). Luego

Res_{z=-rq(a) i} g(z) + Res_{z=+rq(a)i} g(z)
+ Res_{z=0} g(z) + 2 s2 = 0.

lo que da

- 1/(2z) pi csc(pi z)_{z=-rq(a)i}
- 1/(2z) pi csc(pi z)_{z=+rq(a)i}
- 1/a + 2 s2 = 0

con lo cual

- 1/(2 rq(a) i) pi csc(pi rq(a) i)
- 1/(2 rq(a) i) pi csc(pi rq(a) i)
- 1/a + 2 s2 = 0

y asi finalmente

s2 = 1/(2a) - 1/(2 rq(a)) pi csch(pi rq(a)).


Hubo algun progreso en la migraci�n de e.c.m hacia un sitio que ofrezca
MathJAX?

Un saludo.

Marko

Luis

unread,
Jan 12, 2013, 12:02:09 PM1/12/13
to

"Marko Riedel" <markor...@yahoo.de> escribi� en el mensaje
news:87k3rjg...@yahoo.de...
Muy bonito, se�or Riedel. Yo lo hab�a comenzado a resolver
jugando con la funci�n digamma, pero no terminaba de dar
con ello. He impreso su soluci�n pues tengo en mente iniciar
un repaso exhaustivo de la Variable Compleja y entonces
quedar� registrada en mis apuntes.

Un saludo,


Marko Riedel

unread,
Jan 12, 2013, 2:44:47 PM1/12/13
to
De nada, y que tenga mucho exito en sus estudios. La Variable Compleja
es sin duda de lo m�s bonito que hay, con sus conexiones con la teor�a
de n�meros jugando un papel clave.

Un saludo.

Marko

Dr. Wolfgang Hintze

unread,
Jan 14, 2013, 4:36:35 PM1/14/13
to
On 12 Jan., 20:44, Marko Riedel <markoriede...@yahoo.de> wrote:
> "Luis" <wile...@hotmail.com> writes:
> > "Marko Riedel" <markoriede...@yahoo.de> escribió en el mensaje
> >news:87k3rjg...@yahoo.de...
> >> "Dr. Wolfgang Hintze" <w...@snafu.de> writes:
>
> >>> Hallar las sumas (a>-1)
>
> >>> s1 = Sum( 1/(n^2+a), (n,1,oo))
> >>> s2 = Sum((-1)^(n+1)/(n^2+a), (n,1,oo))
>
> >>> Saludos,
> >>> Wolfgang
>
> >> Estimados lectores,
>
> >> aquí vamos a usar el truco que consiste en hallar la integral de
>
> >>     g(z) = f(z) pi cot(pi z)
>
> >> y de
>
> >>     g(z) = f(z) pi csc(pi z)
>
> >> con f(z)  siendo el término  de la suma  que queremos calcular  sobre un
> >> camino que consiste en una circunferencia C  de radio R, con R -> oo. Es
> >> un método bastante clásico, por eso omito los detalles. Es muy facil ver
> >> que en ambos  casos la integral sobre R desaparece en  el límite, con lo
> >> cual  los residuos  a los  polos  deben sumar  cero. (Desaparece  porque
> >> 1/(n^2+a) es O(1/R^2) sobre C.)
>
> >> Lo que  estamos usando aquí  es esencialmente el  hecho de que con  m un
> >> entero
>
> >>    Res_{z=m} pi cot(pi z) =
> >>  = (pi cos(pi z)/(pi cos(pi z)))_{z=m} = 1
>
> >> y
>
> >>    Res_{z=m} pi csc(pi z) =
> >>  = (pi/(pi*cos(pi z)))_{z=m} = (-1)^m.
>
> >> Primer caso:
>
> >> Ponemos  g(z) =  1/(z^2+a)  pi cot(pi  z).  Luego (siendo  "rq" la  raíz
> >> Hubo algun progreso en la migración de e.c.m hacia un sitio que ofrezca
> >> MathJAX?
>
> >> Un saludo.
>
> >> Marko
>
> > Muy bonito, señor Riedel. Yo lo había comenzado a resolver
> > jugando con la función digamma, pero no terminaba de dar
> > con ello. He impreso su solución pues tengo en mente iniciar
> > un repaso exhaustivo de la Variable Compleja y entonces
> > quedará registrada en mis apuntes.
>
> > Un saludo,
>
> De nada, y  que tenga mucho exito en sus  estudios. La Variable Compleja
> es sin duda de  lo más bonito que hay, con sus  conexiones con la teoría
> de números jugando un papel clave.
>
> Un saludo.
>
> Marko
>
>

Muy bien, Marko.

Conoczco tu "truco" como transformación de Watson-Sommerfeld, pero no
me di cuenta en ese caso.

¿ Que tal con los dos productos ?

p1 = Prod(1 + 1/(n^2+a), (n,1,oo) )
p2 = Prod(1+(-1)^(n+1)/n^2, (n,1,oo) )

Saludos,
Wolfgang

Marko Riedel

unread,
Jan 14, 2013, 11:38:48 PM1/14/13
to
Estimados lectores,

por el momento opto por un método simple, y Vds. ya me dirán si se
considera admisible o no.

Para la primera parte, introducimos

f(z) = prod_{n>=1} (1 + 1/(n^2+z^2))

asi que p1 = f(rq(a)).

Vemos que las singularidades de f(z) se hallan en los puntos

z = +/- i n

y son polos simples, entonces f(z) debe ser multiplo de

pi i z / sin(i pi z) = prod_{n>=1} 1/(1 + z^2/n^2).

Luego hay que hallar su cociente, que es

prod_{n>=1} (1 + 1/(n^2+z^2)) (1 + z^2/n^2)
= prod_{n>=1} (1 + (1+z^2)/n^2)

y este último es obviamente

sin(pi i rq(1+z^2)) / (pi i rq(1+z^2))


con lo cual

p1 = pi i rq(a) / sin(i pi rq(a))
* sin(pi i rq(1+a)) / (pi i rq(1+a))

= rq(a)/rq(1+a) sinh(pi rq(1+a))/sinh(pi rq(a)).

Esto funciona porque los polos no tienen punto de acumulación en C y el
producto original converge para todo z salvo dichos polos.

Para la segunda parte, introducimos

f(z) = prod_{m>=1} (1 - z^2/(2m)^2) prod_{m>=1} (1 + z^2/(2m-1)^2).

Todo está bien en lo que se refiere a la convergencia. El valor de p2
será entonces f(1).

Vemos de inmediato que el primer factor es

sin(1/2 pi z)/(1/2 pi z)

y el segundo,

cos(1/2 pi i z),

asi que

p2 = 1/(1/2 pi) cosh(pi/2) = 2/pi cosh(pi/2).

Un saludo.

Marko

Dr. Wolfgang Hintze

unread,
Jan 15, 2013, 7:46:36 AM1/15/13
to
Yo sí considero el método admisible, pero prefiero usar la
representación sugiente de la función Gamma(z)

(1) 1/Gamma(z) = z Exp(EulerGamma z) Product( (1 + z/n) Exp(- z/n) ),
(n,1,oo) )

en que aparecen los factores (1+z/n) inmediatamente.

Lo explicaré si encuentro el tiempo.

Saludos,
Wolfgang
Reply all
Reply to author
Forward
0 new messages