UNIVERSIDAD SAN FRANCISCO DE QUITO USFQ

Colegio de Ciencias e Ingenierías

Una aproximación a las funciones elípticas y sus

aplicaciones
Ensayo Académico
.

Eddy Santiago Achig Andrango

Matemáticas

Trabajo de titulación presentado como requisito

para la obtención del título de
Licenciado en Matemáticas

Quito, 19 de diciembre de 2016

 2

UNIVERSIDAD SAN FRANCISCO DE QUITO USFQ

COLEGIO DE CIENCIAS E INGENIERÍAS

HOJA DE CALIFICACIÓN
DE TRABAJO DE TITULACIÓN

Una aproximación a alas funciones elípticas y sus aplicaciones

Eddy Santiago Achig Andrango

Calificación:
David Fernando Hervas Ortega , PhD.
Nombre del profesor, Título académico
en Matemáticas Aplicadas

Firma del profesor

Quito, 19 de diciembre de 2016

 3

Derechos de Autor
Por medio del presente documento certifico que he leído todas las Políticas y Manuales

de la Universidad San Francisco de Quito USFQ, incluyendo la Política de Propiedad

Intelectual USFQ, y estoy de acuerdo con su contenido, por lo que los derechos de propiedad

intelectual del presente trabajo quedan sujetos a lo dispuesto en esas Políticas.

Asimismo, autorizo a la USFQ para que realice la digitalización y publicación de este

trabajo en el repositorio virtual, de conformidad a lo dispuesto en el Art. 144 de la Ley Orgánica

de Educación Superior.

Firma del estudiante: _______________________________________

Nombres y apellidos: Eddy Santiago Achig Andrango

Código: 00118621

Cédula de Identidad: 171969160-0

Lugar y fecha: Quito, diciembre de 2016

 4

Resumen

El presente trabajo realiza una aproximación a las funciones elı́pticas. En primer

lugar, se definen las funciones elı́pticas sobre C y se estudian las propiedades que
tiene. Después se plantea el problema de existencia de una función elı́ptica, para el
cual se propone una solución constructiva a partir de la ecuación de conducción de
calor. Se definen las funciones theta para poder construir las funciones elı́pticas de
Jacobi. Además, se demuestran ciertas propiedades e identidades útiles de las fun-
ciones theta y de las funciones elı́pticas de Jacobi. A partir de las funciones elı́pticas
de Jacobi, se estudian brevemente las integrales elı́pticas de primeras y de segunda
especie. Finalmente, se presentan aplicaciones de las funciones elı́pticas de Jacobi
en la fı́sica y la geometrı́a.

El enfoque de esta aproximación es seguir los primeros pasos históricos de las

funciones elı́pticas, empezando por las funciones theta que nos ayudarán a construir
las primeras funciones elı́pticas como lo hizo Jacobi. Para encontrar las funciones
theta se utiliza la ecuación de conducción de calor para mostrar que estas funciones
no se encuentran lejos de la realidad fı́sica.

El estudio de las funciones elı́pticas, ha dado soluciones analı́ticas para varios

problemas de la fı́sica y geometrı́a. Además de estos resultado, una de las grandes
aplicaciones de las funciones elı́pticas es en teorı́a de números, donde se utiliza el
estudio de las funciones elı́pticas junto con las curvas elı́pticas. Esta aplicación ha
dado grandes resultados como ser parte fundamental de la demostración del último
teorema de Fermat.

Palabras clave: Funciones elı́pticas, Funciones theta, Análisis complejo, Funciones

de Jacobi.
 5

Abstract

The present work makes an approximation to the elliptic functions. First, elliptic
functions on C are defined and their properties are studied. Then the problem of
the existence of an elliptic function is proposed, for which a constructive solution is
exposed from the heat conduction equation. Theta functions are defined to be able
to construct Jacobi’s elliptic functions. In addition, certain properties and useful
identities of the theta functions and Jacobi’s elliptic functions are demonstrated.
From the elliptic functions of Jacobi, the elliptic integrals of first and second species
are briefly studied. Finally, applications of Jacobi’s elliptic functions in physics and
geometry are presented.

The objective of this approach is to follow the first historical steps of the elliptic
functions, beginning with theta functions that will help us to construct the first
elliptic functions as Jacobi did. To find theta functions we use the heat conduction
equation to show that these functions are not far from physical reality.

The study of elliptic functions, has given analytical solutions to various problems
on physics and geometry. In addition to these results, one of the great applications
of elliptic functions is in number theory, where the study of elliptic functions toget-
her with elliptic curves is used. This application has given great results as being a
fundamental part of the demonstration of Fermat’s last theorem.

Keywords: Elliptic functions, Theta functions, Complex analysis, Jacobi’s fun-

ctions.
 6

Índice general

Introducción 8

1. Funciones elı́pticas 9
1.1. Definición . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.2. Propiedades . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

2. Funciones theta 15
2.1. La función theta como solución de la ecuación de conducción de calor 15
2.2. Propiedades . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.3. Definición de θ1 , θ2 , θ3 , θ4 . . . . . . . . . . . . . . . . . . . . . . . 20
2.4. Propiedades de las funciones θ1 , θ2 , θ3 , θ4 . . . . . . . . . . . . . . . 21

3. Funciones elı́pticas de Jacobi 29

3.1. Definición de las funciones elı́pticas de Jacobi . . . . . . . . . . . . . 29
3.2. Propiedades . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32

4. Integrales elı́pticas 36
4.1. Integrales elı́pticas de primera especie . . . . . . . . . . . . . . . . . . 36
4.2. Integrales elı́pticas de segunda especie . . . . . . . . . . . . . . . . . . 37

5. Aplicación de las funciones elı́pticas de Jacobi 39

5.1. El péndulo simple . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39
5.2. La ecuación de Duffing . . . . . . . . . . . . . . . . . . . . . . . . . . 40
5.3. Geometrı́a de la elipse . . . . . . . . . . . . . . . . . . . . . . . . . . 42

Conclusiones 43

Bibliografı́a 44
 7

Índice de figuras

1.1. La cuadrı́cula Λ generada por 1 y τ . . . . . . . . . . . . . . . . . . . 10

 8

Introducción

En el presente trabajo vamos a considerar funciones definidas sobre C y a valores

complejos. Una función f (z) se define como una función simple periódica de z si
existe una constante distinta de cero w tal que

f (z + w) = f (z)

para todo z ∈ C. Este número w es el periodo de f (z). Claramente, si n ∈ Z y n 6= 0,

entonces nw es también un periodo. Si w es un periodo tal que ningún submúltiplo
de w es periodo, entonces w es un periodo fundamental (Copson, 1935).

Por ejemplo, sabemos que ez es periódica de periodo 2πi

ez+2πi = ez = ez+2nπi

y 2πi es el periodo fundamental de ez .

Nosotros nos preguntamos si existen funciones meromorfas (funciones holo-

morfas excepto en sus polos) tales que para w1 6= w2 y ambos diferentes de cero,
tengamos:
f (z + w1 ) = f (z)
f (z + w2 ) = f (z)
donde w1 y w2 sean periodos fundamentales. Jacobi respondió a esta pregunta al
analizar el caso en el cual la razón entre w1 y w2 es real, y cuando es imaginario (Stein
y Shakarchi, 2010). En el siguiente texto utilizamos la terminologı́a de función
holomorfa para referirnos a una función analı́tica.
 9

Capı́tulo 1

Funciones elı́pticas

1.1. Definición
Nuestro estudio se centrará en funciones meromorfas f en C que tienen dos
periodos fundamentales, es decir, que existen dos números complejos diferentes de
cero w1 y w2 tales que
f (z + w1 ) = f (z) y f (z + w2 ) = f (z)
para toda z ∈ C.
Definición 1.1. Una función con dos periodos fundamentales es llamada doble
periódica.
Para tener dos periodos fundamentales es necesario que w1 y w2 sean linealmente
independientes cuando consideramos a C como un espacio vectorial definido sobre C.

Supongamos que w1 y w2 son linealmente dependientes sobre R, es decir, w2 /w1

∈ R. Si w2 /w1 es racional, digamos igual a p/q, donde p y q son relativamente
primos. Entonces por el Teorema de Bachet-Bézout (Brochero, Moreira, Saldanha,
y Tengan, 2015) existen m, n ∈ Z tal que mq + np = 1. Ahora consideremos:
f (z + mw1 + nw2 ) = f (z)
sabiendo que w2 = pq w1 , tenemos:
(mq+np)
mw1 + nw2 = mw1 + pq w1 = q
w1 = 1q w1

Concluimos que si w2 /w1 es racional, entonces f es una función periódica con un

simple periodo fundamental w0 = 1q w1 .
Si w2 /w1 es irracional igual a τ . Notemos que al fijar z0 , tenemos que
f (z0 + mw1 + nw2 ) = f (z0 )
dado que w2 = τ w1 , por lo tanto f (z0 ) = f (z0 + (m + nτ )w1 ). Dado que el conjunto
{m + nτ |m, n ∈ Z} es denso en R (Brochero y cols., 2015), y por la continuidad de
la función tenemos que f (z0 ) = f (z0 + kw1 ) para toda k ∈ R.Lo que implica que
la función f es constante. Por lo tanto es necesario que los periodos w1 y w2 sean
linealmente independientes sobre R.
 10

Describiremos a continuación una normalización que será útil. Sea τ = w2 /w1 .

Dado que τ ∈ / R, y que τ y 1/τ tienen sus partes imaginarias con signos opues-
tos, vamos a asumir (después de posibles cambios entre w1 y w2 ) que Im(τ ) > 0.
Observemos que la función f tiene periodos w1 y w2 si y solamente si la función
F (z) = f (w1 z) tiene periodos 1 y τ , y además, la función f es meromorfa si y sola-
mente si la función F es meromorfa. También, las propiedades de f son deducibles
de las propiedades de F . Por lo tanto, vamos a asumir sin perdida de generalidad
que f es una función meromorfa en C con periodos 1 y τ donde Im(τ ) > 0.

Notemos que tras aplicar varias veces las condiciones de periodicidad tenemos
f (z + n + mτ ) = f (z) (1.1)
para todos los los enteros n, m y toda z ∈ C, y por ello es natural considerar la
cuadrı́cula en C definida por:
Λ = {n + mτ |n, m ∈ Z}.
Decimos que 1 y τ generan Λ.

Figura 1.1: La cuadrı́cula Λ generada por 1 y τ

La ecuación (1.1) nos dice que f es constante bajo traslaciones de elementos de

Λ.
Definición 1.2. El paralelogramo fundamental (P0 ) es el conjunto
{z ∈ C|z = a + bτ donde 0 ≤ a < 1 y 0 ≤ b < 1}
El paralelogramo fundamental es de gran importancia ya que f está totalmente
determinada por su comportamiento en P0 . Como consecuencia de la definición
tenemos que: [
C= (n + mτ + P0 )
n,m∈Z

es decir, la cuadrı́cula Λ es una cobertura disjunta del plano complejo.

 11

Definición 1.3. Dos números complejos son congruentes módulo Λ si

z = w + n + mτ para algún n, m ∈ Z
y escribimos z ∼ w.
En otras palabras, z y w difieren por un punto en la cuadrı́cula Λ, es decir
z − w ∈ Λ. Por (1.1) concluimos que f (z) = f (w) cuando z ∼ w.
Lema 1.4. La relación de congruencia módulo es una relación de equivalencia
Demostración.
1. Reflexividad. Tenemos que z = z + 0 + 0τ , y 0 + 0τ ∈ Λ. Tenemos que z ∼ z.

2. Simetrı́a. Si tenemos z ∼ w, entonces z = w + h, donde h ∈ Λ. Notemos que

w = z + (−h), donde −h ∈ Λ, por lo tanto w ∼ z.

3. Transitividad. Sea z ∼ w y w ∼ u, entonces existen h, k ∈ Λ tal que z = w + h

y w = u + k. Por lo tanto z = (u + k) + h = u + (k + h), donde k + h ∈ Λ.
Concluimos que z ∼ u
Proposición 1.5. Sea f una función meromorfa con dos periodos 1 y τ que genera
la cuadrı́cula Λ. Entonces:
1. Todo punto en C es congruente a un único punto en el paralelogramo funda-
mental

2. Todo punto en C es congruente a un único punto en cualquier paralelogramo

del tipo P0 + α, donde P0 es el paralelogramo fundamental y α es cualquier
número complejo.

3. La función f está completamente determinada por sus valores en cualquier

paralelogramo del tipo P0 + α, donde P0 es el paralelogramo fundamental y α
es cualquier número complejo.
Demostración. (Stein y Shakarchi, 2010)
1. Sea z = x+iy con x, y ∈ R. Entonces podemos escribir z = a+bτ con a, b ∈ R.
Esto es posible ya que 1 y τ son linealmente independientes sobre los reales, es
decir, generan el plano complejo. Después, escogemos n, m como los mayores
enteros ≤ a y ≤ b respectivamente. Consideremos ahora w = z − n − mτ
donde w ∼ z y w ∈ P0 . Para demostrar unicidad, supongamos que existen w
, w0 ∈ P0 con w ∼ z y w0 ∼ z. Tenemos entonces w ∼ w0 (Por la transitividad
de ∼ Lema 1.4). Si escribimos w = a + bτ y w0 = c + dτ , con a, b, c, d ∈ R,
sigue que w − w0 = (a − c) + (b − d)τ ∈ Λ. Y por lo tanto, a − c, b − d ∈ Z. Sin
embrago, ya que 0 ≤ a, b, c, d < 1, entonces −1 < a − c < 1 y −1 < b − d < 1.
Por lo que concluimos que a = c, b = d, y por lo tanto w = w0 .

2. Sea un punto z0 en el plano complejo y sea P un paralelogramo de la forma

P0 + α, donde α ∈ C. Consideremos z0 − α, por la parte anterior tenemos que
existe un único w0 ∈ P0 tal que z0 − α ∼ w0 , es decir, existe h ∈ Λ tal que
z0 − α = w0 + h. Por lo tanto, z0 = w0 + α + h, donde w0 + α ∈ P . Por lo
tanto podemos concluir que w0 + α es el único punto congruente a z0 en P .
 12

3. De la parte 1. podemos concluir que la función f está totalmente determinada

por sus valores en P0 . Y por la parte 2. aplicada a todos los puntos en P0 ,
concluimos que la función f está completamente determinada por sus valores
en cualquier paralelogramo del tipo P0 + α, donde P0 es el paralelogramo
fundamental y α es cualquier número complejo.
Teorema 1.6. Una función analı́tica en todo C y doble periódica es constante.
Demostración. (Stein y Shakarchi, 2010) Por la Proposición 1.5 tenemos que la
función está completamente determinada por sus valores en P0 , dado que la clausura
de P0 es un compacto, podemos concluir que la función está acotada. Por el teorema
de Liouville, sabemos que toda función entera y acotada es constante.
Definición 1.7. Una función meromorfa doble periódica no constante es llamada
una función elı́ptica.

1.2. Propiedades
Teorema 1.8. Las funciones elı́pticas tienen finitos ceros y polos en el paralelogramo
fundamental P o
Demostración. (Stein y Shakarchi, 2010) Sabemos que una función meromorfa
sólo puede tener finitos ceros y polos en cualquier disco, concluimos que las funcio-
nes elı́pticas van a tener finitos ceros y polos en el paralelogramo fundamental. En
especial, podemos concluir que tendrá finitos ceros y polos en todos los paralelogra-
mos P del tipo P0 + α, donde P0 es el paralelogramo fundamental y α es cualquier
número complejo. Como es usual, contamos los ceros y polos con su multiplicidad.
Teorema 1.9. El número total de polos de una función elı́ptica en P0 es siempre
≥ 2.
Demostración. (Stein y Shakarchi, 2010) Supongamos que f no tiene polos en
la frontera ∂P0 del paralelogramo fundamental. Por el teorema del residuo tenemos:
Z X
f (z)dz = 2πi resf
∂P0
Considerando la integral y rompiendo en sus cuatro caminos tenemos que:
Z Z 1 Z 1+τ Z τ Z 0
f (z)dz = f (z)dz + f (z)dz + f (z)dz + f (z)dz
∂P0 0 1 1+τ τ

por la periodicidad de f , tenemos que

Z 1+τ Z τ Z τ Z 0
f (z)dz = f (z)dz y también tenemos, f (z)dz = f (z)dz
1 0 1+τ 1
R P
y por tanto ∂P0 f (z)dz = 2πi resf = 0. Por lo tanto concluimos que f tiene
por lo menos dos polos en P0 .
Ahora supongamos que f tiene un polo en ∂P0 , escojamos un h ∈ C, con ||h||
pequeño, tal que si P = h + P0 , entonces f no tiene polos en ∂P . Con el argumento
anterior, sabemos que P tiene al menos dos polos en P , y por lo tanto f tiene al
menos dos polos en P0 .
 13

Lema 1.10. Si f es una función holomorfa periódica con periodo ξ, entonces f 0

también es una función periódica con periodo ξ.

Demostración. (Stein y Shakarchi, 2010)Por definición tenemos que para todo

z0 ∈ C
f (z0 + h) − f (z0 )
f 0 (z0 ) = lı́m
h→0 h
por la periodicidad de f , tenemos:

f (z0 + h + ξ) − f (z0 + ξ) f (z0 + h) − f (z0 )

f 0 (z0 + ξ) = lı́m = lı́m = f 0 (z0 )
h→0 h h→0 h

En especial concluimos que si f es una función elı́ptica con periodos 1 y τ , entonces

f 0 /f es una función doble periódica con periodos 1 y τ .

El número total de polo (contados acorde a su multiplicidad) se una función

elı́ptica es llamado su orden. En el siguiente demostramos que el número de ceros
contados acorde a su multiplicidad es igual al orden de la función elı́ptica.

Teorema 1.11. Toda función elı́ptica de orden m tiene m ceros en P0

Demostración. (Stein y Shakarchi, 2010) Vamos a asumir en primer lugar que

la función f no tiene ceros o polos en la frontera de P0 , sabemos por el principio del
argumento de Cauchy que

f 0 (z)
Z
dz = 2πi(N − P )
∂P0 f (z)

donde N y P son el número de ceros y polos de f en P0 respectivamente.

De la integral de la izquierda veamos que:
Z 1 0 Z 1+τ 0 Z τ 0 Z 0 0
f 0 (z)
Z
f (z) f (z) f (z) f (z)
dz = dz + dz + dz + dz
∂P0 f (z) 0 f (z) 1 f (z) 1+τ f (z) τ f (z)

por la periodicidad de f 0 /f (Lema 1.10), tenemos que

1+τ Z τ 0 Z τ 0 Z 0 0
f 0 (z)
Z
f (z) f (z) f (z)
dz = dz y también tenemos, dz = dz
1 f (z) 0 f (z) 1+τ f (z) 1 f (z)

Tenemos: ∂P0 f 0 (z)/f (z)dz = 2πi(N − P ) = 0. Por lo tanto concluimos que

R

N = P.
Ahora supongamos que f tiene un polo en ∂P0 , escojamos un h ∈ C, con ||h|| pe-
queño, tal que si P = h + P0 , entonces f no tiene polos en ∂P . Con el argumento
anterior, sabemos que P tiene el mismo número de ceros y polos en P , y por lo tanto
f tiene el mismo número de polos y ceros en P0 .

Como consecuencia, si f es una función elı́ptica entonces la ecuación f (z) = c

tiene tantas soluciones como el orden de f para cada c ∈ C, simplemente porque
f − c es elı́ptica y tiene tantos polos como f .
 14

Teorema 1.12. Si f (z) y g(z) son funciones elı́pticas con polos en los mismos
puntos y con la misma parte principal en estos puntos, entonces f (z) = g(z) + c
para una constante c

Demostración. (Jones y Singerman, 1987) Notemos que la función f − g es una

función entera doble periódica, y por el Teorema 1.6 tenemos que f − g es constante.

Teorema 1.13. Si f (z) y g(z) son funciones elı́pticas con ceros y polos del mismo
orden en los mismos puntos, entonces f (z) = cg(z) para una constante c

Demostración. (Jones y Singerman, 1987) Notemos que la función f /g es una

función entera doble periódica, por lo tanto por el Teorema 1.6, f /g es constante.

A pesar de la naturaleza simple de las propiedades de las funciones elı́pticas, aún

queda abierta la pregunta acerca de la existencia de las funciones elı́pticas. Ahora el
objetivo es dar una demostración constructiva a esta pregunta a partir del estudio
de las funciones theta que se estudiarán en el siguiente capı́tulo.
 15

Capı́tulo 2

Funciones theta

2.1. La función theta como solución de la ecuación

de conducción de calor
Introduciremos la función theta considerando un problema especı́fico de conduc-
ción de calor. Sea θ la temperatura en el tiempo t en cualquier punto de un material
sólido isotrópico cuya conductividad térmica es uniforme. Entonces, si ρ es la densi-
dad del material, s es su calor especı́fico, y k su conductividad térmica, θ satisface
la ecuación diferencial

κ∇2 θ = ∂θ/∂t
donde κ = k/sρ es la difusividad. Centraremos el problema en el caso particular
en el cual no hay variación de temperatura en las direcciones x , y en el marco
Cartesiano. El calor fluye en todos los lugares paralelo a eje z y la ecuación de
conducción de calor se reduce a la forma

κ∂ 2 θ/∂z 2 = ∂θ/∂t (2.1)

por lo tanto, θ queda en función de z y t, es decir, θ = θ(z, t) (Lawden, 2013).

El problema especı́fico que estudiaremos es la conducción del calor en un bloque

infinito de material limitado por los planos z = 0, π, cuando las condiciones sobre
cada plano en el lı́mite se mantienen uniformes para cada tiempo t. La conducción
de temperatura es entonces enteramente en la dirección del eje z y la ecuación (2.1)
se aplica.

Primero, supongamos que como condiciones de frontera tenemos que las caras
del bloque son mantenidas a temperatura cero, es decir, θ = 0 para z = 0, π y todo
t. Inicialmente en t = 0, supongamos θ = f (z) para 0 < z < π. Podemos observar
que la ecuación de calor y las condiciones de frontera son homogéneas y lineales. En
primer lugar, encontramos todas las soluciones a las condiciones lineales homogéneas
de la forma
θ(z, t) = h(z)g(t) (2.2)
Por el principio de superposición, una superposición lineal arbitraria de estas solu-
ciones va a seguir siendo una solución.
Sustituimos θ(z, t) = h(z)g(t) en la ecuación (2.1) y obtenemos
 16

h(z)g 0 (t) = κh00 (z)g(t)

Separamos los términos que involucran t a la izquierda y los términos que invo-
lucran z a la derecha:

g 0 (t) h00 (z)

=κ
g(t) h(z)
Como el lado izquierdo no depende de z y el lado derecho no depende de t,
tenemos que ambos lados son iguales a una constante λ.

g 0 (t) h00 (z)

= =λ
κg(t) h(z)

Y obtenemos dos ecuaciones diferenciales ordinarias

g 0 (t) = λκg(t) (2.3)

y
h00 (x) = λh(x) (2.4)
Las condiciones de frontera nos dicen que θ(0, t) = θ(π, t) = 0. Reemplazando en la
ecuación (2.2) tenemos
h(0)g(t) = h(π)g(t)
Para todos los valores de t. La solución g(t) = 0 para todos los valores de t, nos
lleva a la solución trivial θ(z, t) = 0 para todo z y t. Por lo tanto consideraremos el
caso en el que g(t) no es siempre 0, tenemos:

h(0) = h(π) = 0

El caso en el cual g(t) es siempre cero nos lleva a la respuesta trivial θ = 0. Por lo
tanto, para encontrar soluciones no triviales, debemos resolver la siguiente ecuación
diferencial ordinaria con condiciones de frontera:

h00 (x) = λh(x) , h(0) = 0 = h(π) (2.5)

Si λ = 0, entonces (2.5) se convierte en:

h00 (x) = 0 , h(0) = 0 = h(π)

por lo tanto esto nos lleva a la solución trivial h = 0.

Si λ > 0, entonces (2.5) se convierte en:

h00 (x) − λh(x) = 0

La solución general es:

√ √
h(x) = a cosh( λx) + b sinh( λx)

Dada la condición de frontera, tenemos:

√
h(0) = a = 0 entonces, h(x) = b sinh( λx)
 17

y también √
h(π) = b sinh( λπ) = 0 entonces, h(x) = 0
por lo tanto esto nos lleva a la solución
√ trivial h = 0.
Si λ < 0, en este caso definimos w = −λ, y (2.5) se convierte en:

h00 (x) + w2 h(x) = 0 , h(0) = 0 = h(π)

, con solución general:

h(x) = a cos(wx) + b sin(wx)
Dada la condición de frontera, tenemos:

h(0) = a = 0 entonces, h(x) = b sin(wx)

y también
h(π) = b sin(wπ) = 0
entonces b = 0 o sin(wπ) = 0. La posibilidad b = 0, nos lleva a la solución trivial.
La segunda posibilidad se cumple en el caso w = n, donde n ∈ Z. Por lo tanto las
únicas soluciones no triviales para (2.5) son múltiplos constantes de:

h(x) = sin(nx) , n ∈ Z

con λ = −n2 .
Para cada una de estas soluciones, se debe encontrar la función g(t) correspondiente
con respecto a la ecuación (2.3).

g 0 (t) = λκg(t)

donde λ = −n2 , esta ecuación diferencial tiene como solución

2 κt
g(t) = ae−n

donde a es una constante de integración por lo tanto encontramos que las soluciones
no triviales al problema de la ecuación de calor junto con las condiciones de frontera
dadas son múltiplos constantes de:
2 κt
θn (z, t) = sin(nz)e−n

Por el principio de superposición, sigue que:

θ(z, t) = b1 sin(z)e−κt + b2 sin(2z)e−2κt + ... (2.6)

es también una solución al problema planteado.

Ahora determinaremos los coeficientes bn en (2.6) de tal manera que la condición

inicial f (z) = θ(z, 0) se satisfaga. Fijando t = 0 en (2.6), tenemos:
∞
X
f (z) = bn sin(nz)
n=1

Aplicando la teorı́a de series de Fourier, podemos encontrar bn :

2 π
Z
bn = f (z) sin(nz)dz
π 0
 18

En el caso especial donde f (z) = πδ(z − 21 )π (donde δ(z) es la función de impulso

unitario de Dirac), el bloque está inicialmente a temperatura cero en todas partes,
excepto en el plano z = 21 π, donde la temperatura en muy alta. Calculamos que:
Z π    
1 1
bn = 2 δ z − π sin(nz)dz = 2 sin nπ
0 2 2
Por lo tanto, la difusión de calor sobre el bloque está gobernada por la ecuación:
∞
2 κt
X
θ(z, t) = 2 (−1)n e−(2n+1) sin((2n + 1)z)
n=0

Escribiendo e−4κt = q, θ asume la forma:

∞
2
X
θ1 (z, q) = 2 (−1)n q (n+1/2) sin((2n + 1)z) (2.7)
n=0

La primera función theta θ1 (z, q) está definida por la ecuación (2.7) para todos
los números complejos de z y q tal que |q| < 1. Reemplazando la función seno por
su representación de Euler, obtenemos:
∞
2
X
iθ1 (z, q) = (−1)n q (n+1/2) ei(2n+1)z (2.8)
n=−∞

Para establecer la convergencia, notamos un como el enésimo término en la serie.

Entonces:  
 un+1 
 un  = |q
  e | = |q|2n+2 e−2y
2n+2 2iz

escribiendo z = x + iy. Siendo |q| < 1, la razón tiende a cero cuando n → +∞,
por el test DÁlambert sabemos que la serie converge en +∞. Cuando n → −∞
consideramos la razón |un /un+1 |, claramente esta razón converge a cero, y por lo
tanto la serie converge en −∞ también.

Como notación alternativa, escribimos:

q = eiπτ
donde la parte imaginaria de τ debe ser positiva para obtener |q| < 1. Entonces
tenemos: ∞
2
X
iθ1 (z, τ ) = (−1)n ei(n+1/2) πτ ei(2n+1)z
n=−∞

2.2. Propiedades
De manera general tenemos que la función theta está dada por las series:
∞
2
X
Θ(z, τ ) = eiπn τ e2πinz
n=−∞

que converge para toda z ∈ C, y τ con parte imaginaria mayor que cero (H).
 19

Proposición 2.1. La función Θ satisface las siguiente propiedades:

1. Θ es entera en z ∈ C y holomorfa en τ ∈ H
2. Θ(z + 1|τ ) = Θ(z|τ )
3. Θ(z + τ |τ ) = Θ(z|τ )e−πiτ e−2πiz
4. Θ(z|τ ) = 0 cuando z = 1/2 + τ /2 + n + mτ donde n, m ∈ Z
Demostración. (Stein y Shakarchi, 2010)
1. Supongamos que Im (τ ) = t ≥ t0 > 0 y z = x + iy pertenece a un con-
junto limitado de C, digamos |z| ≤ M . Entonces la serie definida por Θ es
absolutamente y uniformemente convergente, porque
∞
2 2
X X
|eπin τ e2πinz | ≤ C e−πn t0 e2πnM < ∞
−∞ n≥0

Por lo tanto, para un τ ∈ H fijo la función Θ(.|τ ) es entera, y para cada z ∈ C

fija, la función Θ(z|.) es holomorfa en el medio plano superior (H).
2. Tenemos:
∞ ∞
iπn2 τ 2πi(n+1)z 2
X X
Θ(z + 1|τ ) = e e = eiπn τ e2πinz e2πi
n=−∞ n=−∞

Por la fórmula de Euler sabemos que e2πi = 1. Por lo tanto Θ(z +1|τ ) = Θ(z|τ )
3. Tenemos:
∞
2
X
Θ(z + τ |τ ) = eπin τ e2πin(z+τ )
n=−∞
∞
2
X
= eπi(n+1) τ e2πinz e−πiτ
n=−∞
∞
2
X
= eπi(n+1) τ e2πi(n+1)z e−πiτ e−2πinz
n=−∞

= Θ(z|τ )e−πiτ e−2πiz

4. Por los literales (2) y (3) anteriores, es suficiente mostrar que Θ(1/2+τ /2|τ ) =
0. Tenemos que:
∞
2
X
Θ(1/2 + τ /2|τ ) = eπin τ e2πin(1/2+τ /2)
n=−∞
∞
2 +n)τ
X
= (−1)n eπi(n
n=−∞

Notemos que para n ≥ 0, n y −n − 1 tienen diferente pariedad y además:

(−n − 1)2 + (−n − 1) = n2 + n. Por lo tanto:
∞ ∞
n πi(n2 +n)τ 2 +n)τ
X X
Θ(1/2 + τ /2|τ ) = (−1) e − (−1)n eπi(n =0
n=0 n=0
 20

Un hecho sobresaliente de la función theta es su naturaleza dual. Vista como

una función de z, se encuentra en el campo de las funciones elı́pticas , dado que
Θ es periódica con periodo 1 y çuasiperiodo”τ como se mostró arriba. Cuando la
consideramos como una función de τ , Θ revela su naturaleza modular y su conexión
cercana con la función de partición y el problema de representar enteros como la
suma de cuadrados.

2.3. Definición de θ1, θ2, θ3 , θ4

En la sección 2.1, definimos θ1 como la solución para un problema de conduc-
ción de calor especı́fico. Otra función theta aparece como solución si cambiamos las
condiciones de frontera. Supongamos que las caras del bloque son isoladas, de tal
manera que el calor no puede escapar a través de ellas. Esta condición de frontera
se escribe como ∂θ/∂z = 0 cuando z = 0, π, y el método de separación de variables
nos lleva a la solución:
∞
1 X 2
θ(z, t) = a0 + an e−n κt cos nz
2 n=1

donde
2 π
Z
an = f (z) cos nzdz
π 0
En el caso f (z) = πδ z − 12 π tenemos que


 
1
an = 2 cos nπ
2
y la solución al problema es:
∞
2
X
θ4 (z, q) = 1 + 2 (−1)n q n cos(2nz)
n=1
−4κt
donde q = e .

Si cambiamos la función coseno por su expresión por exponenciales según la

fórmula de Euler, tenemos una otra representación de θ4
∞
2
X
θ4 (z) = (−1)n q n e2inz (2.9)
n=−∞

Definimos entonces:
  ∞
1 X 2
θ2 (z) = θ1 z + π = q (n+1/2) ei(2n+1)z (2.10)
2 n=−∞
  ∞
1 X 2
θ3 (z) = θ4 z + π = q n e2inz (2.11)
2 n=−∞

Si definimos q = eiπτ , con la parte imaginaria de τ es positiva. Donde tenemos que:

∞ z 
2
X
θ3 (z, τ ) = eiπn τ e2inz = Θ ,τ (2.12)
−∞
π
donde Θ es la serie definida en la sección 2.2.
 21

2.4. Propiedades de las funciones θ1, θ2, θ3 , θ4

Estas cuatro funciones theta, son importantes ya que nos permiten definir las
funciones elı́pticas de Jacobi. A continuación vamos a demostrar algunas relaciones
útiles entre θ1 , θ2 , θ3 y θ4 . Definiendo q = eiπτ :

Lema 2.2. Las funciones θ1 , θ2 , θ3 , θ4 , satisfacen:

     
1 1 1 1
θ1 (z) = −θ2 z + π = −iµθ4 z + πτ = −iµθ3 z + π + πτ
2 2 2 2

donde µ = q 1/4 eiz .

Demostración. (Lawden, 2013) Tenemos que:

  ∞
1 X 2
θ2 z + π = eiπτ (n+1/2) ei(2n+1)(z+π/2)
2 n=−∞
∞
2
X
= eiπτ (n+1/2) ei(2n+1)z einπ eiπ/2
n=−∞
∞
2
X
= eiπτ (n+1/2) ei(2n+1)z einπ eiπ/2
n=−∞
∞
2
X
=i (−1)n eiπτ (n+1/2) ei(2n+1)z
n=−∞

= −θ1 (z)

  ∞
1 X 2
θ4 z + πτ = (−1)n eiπτ n e2in(z+πτ /2)
2 n=−∞
∞
2
X
−iz −iπτ /4
=e e (−1)n eiπτ (n+1/2) ei(2n+1)z
n=−∞
−1
= µ iθ1 (z)

Para terminar la demostración, basta probar que: θ3 (z + π/2) = θ4 (z):

  ∞
1 X 2
θ3 z+ π = eiπτ n e2in(z+π/2)
2 n=−∞
∞
2
X
= (−1)n eiπτ n e2inz
n=−∞

= θ4 (z)

El Lema 2.2 junto con la ecuación 2.12 nos muestra una manera natural de
escribir las cuatro funciones theta definidas en esta sección, por medio de la función
Θ definida por una serie en la sección 2.2. A continuación, demostraremos una
propiedad importante de las funciones theta que juntamente con el Lema 2.2, será
de gran importancia en el siguiente capı́tulo.
 22

Lema 2.3. Se cumple:

θ42 (x)θ42 (y) − θ12 (x)θ12 (y) = θ4 (x + y)θ4 (x − y)θ42 (0)

Demostración. Vamos a demostrar esta identidad multiplicando series absoluta-

mente convergentes por lo cual el reordenar los términos de las series queda justifi-
cado. Según la ecuación (2.8) tenemos:
2 2
XX
θ1 (x, q)θ1 (y, q) = − (−1)m+n q (m+1/2) +(n+1/2) ei(2m+1)x+i(2n+1)y (2.13)
m n

Cambiando las variables enteras (m, n) a (r, s) con la relación:

m + n = r, m − n = s
Además, podemos escribir:
 2  2
1 1 1 1
m+ + n+ = (r + 1)2 + s2
2 2 2 2

(2m + 1)x + (2n + 1)y = (r + 1)(x + y) + s(x − y)

Entonces,
2 +(1/2)s2
X
θ1 (x, q)θ1 (y, q) = − (−1)r q (1/2)(r+1) ei(r+1)(x+y)+is(x−y)

Si (m, n) son ambos pares o ambos impares, (r, s) serán ambos pares. Si por el
contrario, (m, n) tienen pariedad opuesta, (r, s) serán ambos impares. Entonces per-
mitiendo a (r, s) recorrer todas las parejas de números pares y todas las parejas de
números impares, cada pareja de enteros (m, n) se obtendrá una sola vez. Sigue que:
2 2
XX
θ1 (x, q)θ1 (y, q) = − q 2(r+1/2) +2s ei(2r+1)(x+y)+2is(x−y)
r s
2 +2(s+1/2)2
XX
+ q 2r e2ir(x+y)+i(2s+1)(x−y)
r s

Cada una de las dobles series de arriba, puede ser expresada como el producto de
dos series simples, como sigue:
2
X X 2
θ1 (x, q)θ1 (y, q) = − q 2(r+1/2) ei(2r+1)(x+y) q 2s e2is(x−y)
r s
2r2 2ir(x+y) 2(s+1/2)2 i(2s+1)(x−y)
X X
+ q e q e
r s

Sustituyendo con las ecuaciones (2.10) y (2.11), obtenemos:

θ1 (x, q)θ1 (y, q) = θ3 (x + y, q 2 )θ2 (x − y, q 2 ) − θ2 (x + y, q 2 )θ3 (x − y, q 2 ) (2.14)

Según la ecuación (2.11) tenemos:

2 +m2
XX
θ3 (x, q)θ3 (y, q) = qn e2inx+2imy (2.15)
m n

Cambiando las variables enteras (m, n) a (r, s) con la relación:

 23

m + n = r, m − n = s
Además, podemos escribir:
r 2 + s2
m2 + n2 =
2
2nx + 2my = r(x + y) + s(x − y)
Entonces,
2 +s2 )/2
X
θ3 (x, q)θ3 (y, q) = q (r eir(x+y)+is(x−y)
Por el mismo argumento realizado arriba, notamos que debemos sumar sobre todas
las parejas de números pares e impares (r, s). Sigue que:
XX 2 2
θ3 (x, q)θ3 (y, q) = q 2r +2s e2ir(x+y)+2is(x−y)
r s
2 +2(s+1/2)2
XX
+ q 2(r+1/2) ei(2r+1)(x+y)+i(2s+1)(x−y)
r s

Cada una de las dobles series de arriba, puede ser expresada como el producto de
dos series simples, como sigue:
X 2 X 2
θ1 (x, q)θ1 (y, q) = − q 2r ei2r(x+y) q 2s e2is(x−y)
r s
2(r+1/2)2 i(2r+1)(x+y) 2
X X
+ q e q 2(s+1/2) ei(2s+1)(x−y)
r s

Sustituyendo con las ecuaciones (2.10) y (2.11), obtenemos:

θ3 (x, q)θ3 (y, q) = θ3 (x + y, q 2 )θ3 (x − y, q 2 ) + θ2 (x + y, q 2 )θ2 (x − y, q 2 ) (2.16)

Realizando un incremento de 12 π en x, y, tenemos:

θ4 (x, q)θ4 (y, q) = θ3 (x + y, q 2 )θ3 (x − y, q 2 ) − θ2 (x + y, q 2 )θ2 (x − y, q 2 ) (2.17)

Elevando al cuadrado y restando las ecuaciones (2.14) y (2.17), obtenemos:

θ42 (x, q)θ42 (y, q)−θ12 (x, q)θ12 (y, q) = [θ32 (x+y, q 2 )−θ22 (x+y, q 2 )][θ32 (x−y, q 2 )−θ22 (x−y, q 2 )]
(2.18)
Reemplazando y = 0 en la ecuación (2.17) tenemos:

θ32 (x, q 2 ) − θ22 (x, q 2 ) = θ4 (x, q)θ4 (0, q) (2.19)

Utilizando las ecuaciones (2.18) y (2.19), tenemos:

θ42 (x)θ42 (y) − θ12 (x)θ12 (y) = θ4 (x + y)θ4 (x − y)θ42 (0) (2.20)

Las funciones definidas en la sección anterior tienen propiedades de periodicidad que

se pueden escribir como en la siguiente proposición:
Proposición 2.4. La funciones θ1 , θ2 , θ3 , θ4 , satisfacen:
1. θ1 (z) = −θ1 (z + π) = −λθ1 (z + πτ ) = λθ1 (z + π + πτ )

2. θ2 (z) = −θ2 (z + π) = λθ2 (z + πτ ) = −λθ2 (z + π + πτ )

 24

3. θ3 (z) = θ3 (z + π) = λθ3 (z + πτ ) = λθ3 (z + π + πτ )

4. θ4 (z) = θ4 (z + π) = −λθ4 (z + πτ ) = −λθ4 (z + π + πτ )

donde λ = ei(2z+πτ ) .

Demostración. (Lawden, 2013)

1. Usando la ecuación (2.8), tenemos:

∞
2
X
θ1 (z + π) = −i (−1)n eiπτ (n+1/2) ei(2n+1)(z+π)
n=−∞
∞
2
X
= −i (−1)n eiπτ (n+1/2) ei(2n+1)z eiπ(2n+1)
n=−∞
∞
2
X
=i (−1)n eiπτ (n+1/2) ei(2n+1)z
n=−∞

= −θ1 (z)

∞
2
X
θ1 (z + πτ ) = −i (−1)n eiπτ (n+1/2) ei(2n+1)(z+πτ )
n=−∞
∞
2
X
= −i (−1)n eiπτ (n+1/2) ei(2n+1)z eiπτ (2n+1)
n=−∞
∞
2
X
=i (−1)n+1 eiπτ ((n+1)+1/2) ei(2(n+1)+1)z e−iπτ e−2iz
n=−∞
−1
= −λ θ1 (z)

Utilizando lo demostrado arriba, resulta:

θ1 (z + π + πτ ) = −θ1 (z + πτ )
= λ−1 θ1 (z)

2. Usando la ecuación (2.10), tenemos:

∞
2
X
θ2 (z + π) = eiπτ (n+1/2) ei(2n+1)(z+π)
n=−∞
∞
2
X
= eiπτ (n+1/2) ei(2n+1)z eiπ(2n+1)
n=−∞
∞
2
X
=− eiπτ (n+1/2) ei(2n+1)z
n=−∞

= −θ2 (z)
 25

∞
2
X
θ2 (z + πτ ) = eiπτ (n+1/2) ei(2n+1)(z+πτ )
n=−∞
∞
2
X
= eiπτ (n+1/2) ei(2n+1)z eiπτ (2n+1)
n=−∞
∞
2
X
=− eiπτ ((n+1)+1/2) ei(2(n+1)+1)z e−iπτ e−2iz
n=−∞
−1
= λ θ2 (z)
Utilizando lo demostrado arriba, resulta:
θ2 (z + π + πτ ) = −θ2 (z + πτ )
= −λ−1 θ1 (z)

3. Utilizando la ecuación (2.11), tenemos:

∞
2
X
θ3 (z + π) = eiπτ n e2inz e2inπ
n=−∞

= θ3 (z)

∞
2
X
θ3 (z + πτ ) = eiπτ n e2inz e2inπτ
n=−∞
∞
2
X
= eiπτ (n+1) e2i(n+1)z e−iπτ e−2iz
n=−∞

= λ−1 θ3 (z)
Utilizando lo demostrado arriba, resulta:
θ3 (z + π + πτ ) = θ3 (z + πτ )
= λ−1 θ1 (z)

4. Utilizando la ecuación (2.9) tenemos:

∞
2
X
θ4 (z + π) = (−1)n eiπτ n e2inz e2inπ
n=−∞

= θ4 (z)

∞
2
X
θ4 (z + πτ ) = (−1)n eiπτ n e2inz e2inπτ
n=−∞
∞
2
X
=− (−1)n+1 eiπτ (n+1) e2i(n+1)z e−iπτ e−2iz
n=−∞

= −λ−1 θ4 (z)
Utilizando lo demostrado arriba, resulta:
θ4 (z + π + πτ ) = θ4 (z + πτ )
= −λ−1 θ4 (z)
 26

A continuación estudiaremos las derivadas de razones con las funciones θ1 , θ2 , θ3 ,

θ4 . Para lo cual primero vamos a demostrar una identidad importante:
Lema 2.5. Se cumple la identidad
θ10 (0) = θ2 (0)θ3 (0)θ4 (0) (2.21)
Demostración. (Lawden, 2013) En la ecuación (2.14) incrementamos el valor de
y en 21 π, en la misma ecuación incrementamos el valor de x, y en 12 π, en la ecuación
(2.16) incrementamos el valor de y en 12 π, y utilizando el Lema 2.2, obtenemos las
siguientes ecuaciones respectivamente:
θ1 (x, q)θ2 (y, q) = θ1 (x + y, q 2 )θ4 (x − y, q 2 ) + θ4 (x + y, q 2 )θ1 (x − y, q 2 ) (2.22)
θ2 (x, q)θ2 (y, q) = θ2 (x + y, q 2 )θ3 (x − y, q 2 ) + θ3 (x + y, q 2 )θ2 (x − y, q 2 ) (2.23)
2 2 2 2
θ3 (x, q)θ4 (y, q) = θ4 (x + y, q )θ4 (x − y, q ) − θ1 (x + y, q )θ1 (x − y, q ) (2.24)
Si derivamos parcialmente la ecuación (2.22) con respecto a x, y utilizamos x = y =
0, obtenemos:
θ10 (0, q)θ2 (0, q) = 2θ10 (0, q 2 )θ4 (0, q 2 ) (2.25)
Después, al utilizar x = y = 0 en las ecuaciones (2.23) y (2.24), y dado que θ1 (0) = 0,
tenemos:
θ22 (0, q) = 2θ2 (0, q 2 )θ3 (0, q 2 )
θ3 (0, q)θ4 (0, q) = θ42 (0, q 2 )
Dividimos la ecuación (2.25) para las identidades de arriba y obtenemos:
θ10 (0, q) θ10 (0, q 2 )
=
θ2 (0, q)θ3 (0, q)θ4 (0, q) θ2 (0, q 2 )θ3 (0, q 2 )θ4 (0, q 2 )
Repitiendo el resultado, podemos demostrar que:
n
θ10 (0, q) θ10 (0, q 2 )
=
θ2 (0, q)θ3 (0, q)θ4 (0, q) θ2 (0, q 2 )θ3 (0, q 2n )θ4 (0, q 2n )
n

n
para todo entero positivo n. Dado que |q| < 1, al tener n → ∞ tenemos q 2 → 0,
por lo tanto:
θ10 (0, q) θ10 (0, q)
= lı́m (2.26)
θ2 (0, q)θ3 (0, q)θ4 (0, q) q→0 θ2 (0, q)θ3 (0, q)θ4 (0, q)
Con referencia a las ecuaciones (2.8), (2.9), (2.10) y (2.11), sabemos que:
∞
2
X
θ10 (0) = (−1)n (2n + 1)q (n+1/2) = 2q 1/4 + O(q 9/4 )
n=−∞
∞
2
X
θ2 (0) = q (n+1/2) = 2q 1/4 + O(q 9/4 )
n=−∞
∞
2
X
θ3 (0) = q n = 1 + O(q)
n=−∞
∞
2
X
θ4 (0) = (−1)n q n = 1 + O(q)
n=−∞

Por lo que el lı́mite de la ecuación (2.26) es 1. De donde obtenemos:

θ10 (0) = θ2 (0)θ3 (0)θ4 (0)
 27

Teorema 2.6. Se cumplen las siguientes derivadas:

d
(θ1 /θ4 ) = θ42 (0)θ2 (x)θ3 (x)/θ42 (x) (2.27)
dx
d
(θ2 /θ4 ) = −θ32 (0)θ1 (x)θ3 (x)/θ42 (x) (2.28)
dx
d
(θ3 /θ4 ) = −θ22 (0)θ1 (x)θ2 (x)/θ42 (x) (2.29)
dx
Demostración. (Lawden, 2013) Intercambiando x, y en la ecuación en la ecuación
(2.24) y multiplicándola por la ecuación (2.22) obtenemos:

θ1 (x, q)θ2 (y, q)θ3 (y, q)θ4 (x, q) = [θ1 (x + y, q 2 )θ4 (x − y, q 2 ) + θ4 (x + y, q 2 )θ1 (x − y, q 2 )]
× [θ4 (x + y, q 2 )θ4 (x − y, q 2 ) + θ1 (x + y, q 2 )θ1 (x − y, q 2 )]
(2.30)

Intercambiando x, y en la ecuación (2.30) y sumando con la misma ecuación, tene-

mos:
θ1 (x, q)θ2 (y, q)θ3 (y, q)θ4 (x, q) + θ1 (y, q)θ2 (x, q)θ3 (x, q)θ4 (y, q)
(2.31)
= 2θ1 (x + y, q 2 )θ4 (x + y, q 2 )[θ42 (x − y, q 2 ) + θ12 (x − y, q 2 )]

Al reemplazar y = 0 en la ecuación (2.22) y x = 0 en la ecuación (2.24), tenemos

las siguientes identidades:

θ1 (x, q)θ2 (0, q) = 2θ1 (x, q 2 )θ4 (x, q 2 )

θ3 (0, q)θ4 (y, q) = θ42 (y, q 2 ) + θ12 (y, q 2 )

Usando estas identidades en la ecuación (2.31), se reduce a la forma:

θ1 (x + y)θ4 (x − y)θ2 (0)θ3 (0) = θ1 (x)θ4 (x)θ2 (y)θ3 (y) + θ2 (x)θ3 (x)θ1 (y)θ4 (y) (2.32)

Derivamos la ecuación (2.32) parcialmente con respecto a y, y obtenemos:

[θ10 (x + y)θ4 (x − y) − θ1 (x + y)θ40 (x − y)]θ2 (0)θ3 (0)

(2.33)
= θ1 (x)θ4 (x)[θ20 (y)θ3 (y) + θ2 (y)θ30 (y)] + θ2 (x)θ3 (x)[θ10 (y)θ4 (y) + θ1 (y)θ40 (y)]

Si ponemos y = 0 en la ecuación (2.33), tenemos que:

θ10 (x)θ4 (x) − θ1 (x)θ40 (x)

θ1 (x)θ4 (x)[θ20 (0)θ3 (0) + θ2 (0)θ30 (0)] + θ2 (x)θ3 (x)[θ10 (0)θ4 (0) + θ1 (0)θ40 (0)] (2.34)
=
θ2 (0)θ3 (0)

Utilizando que θ1 (0) = θ20 (0) = θ30 (0) = θ40 (0) = 0 en la ecuación (2.34), obtenemos:

θ2 (x)θ3 (x)θ10 (0)θ4 (0)

θ10 (x)θ4 (x) − θ1 (x)θ40 (x) = (2.35)
θ2 (0)θ3 (0)

Utilizando la identidad (2.21), tenemos:

d
(θ1 /θ4 ) = θ42 (0)θ2 (x)θ3 (x)/θ42 (x)
dx
 28

Con lo que queda demostrada la identidad (2.27).

Incrementamos el valor de y en πτ /2 + π/2, en la ecuación (2.20). Y obtenemos:

θ42 (x)θ22 (y) = θ2 (x + y)θ2 (x − y)θ42 (0) + θ12 (x)θ32 (y) (2.36)

Con y = 0 en la última ecuación tenemos:

θ42 (x)θ22 (0) = θ22 (x)θ42 (0) + θ12 (x)θ32 (0)

Podemos escribir esta última ecuación como:

 2  2
θ1 (x) 2 θ2 (x)
θ3 (0) + θ42 (0) = θ22 (0)
θ4 (x) θ4 (x)

Si diferenciamos parcialmente con respecto a x, tenemos:

   
θ1 (x) d 2 θ2 (x) d
(θ1 /θ4 )θ3 (0) + (θ2 /θ4 )θ42 (0) = 0
θ4 (x) dx θ4 (x) dx

Utilizando la identidad (2.27), obtenemos:

d
(θ2 /θ4 ) = −θ32 (0)θ1 (x)θ3 (x)/θ42 (x)
dx
Al incrementar los valores de x, y en π/2 en la ecuación (2.20), tenemos:

θ4 (x + y)θ4 (x − y)θ42 (0) = θ32 (x)θ32 (y) + θ22 (x)θ22 (y) (2.37)

Al poner y = 0 en la anterior igualdad, tenemos:

θ42 (x)θ42 (0) = θ32 (x)θ32 (0) + θ22 (x)θ22 (0)

Podemos escribir esta última ecuación como:

 2  2
θ2 (x) 2 θ3 (x)
θ2 (0) + θ32 (0) = θ42 (0)
θ4 (x) θ4 (x)

Si diferenciamos parcialmente con respecto a x, tenemos:

   
θ2 (x) d 2 θ3 (x) d
(θ2 /θ4 )θ2 (0) + (θ3 /θ4 )θ32 (0) = 0
θ4 (x) dx θ4 (x) dx

Utilizando la identidad (2.28), tenemos:

d
(θ3 /θ4 ) = −θ22 (0)θ1 (x)θ2 (x)/θ42 (x)
dx
 29

Capı́tulo 3

Funciones elı́pticas de Jacobi

3.1. Definición de las funciones elı́pticas de Jacobi

Las funciones elı́pticas sn(u), cn(u) y dn(u) quedan definidas como razones de
funciones theta.
θ3 (0) θ1 (z)
sn(u) = .
θ2 (0) θ4 (z)
θ4 (0) θ2 (z)
cn(u) = .
θ2 (0) θ4 (z)
θ4 (0) θ3 (z)
dn(u) = .
θ3 (0) θ4 (z)
donde z = u/θ32 (0).

Proposición 3.1. sn(u) y cn(u) satisfacen: sn2 (u) + cn2 (u) = 1

Demostración. (Lawden, 2013) Por el Lema 2.2, después de aumentar π/2+τ π/2
al valor de y en (2.20), tras simplificar obtenemos:

θ42 (x)θ22 (y) − θ12 (x)θ32 (y) = θ2 (x + y)θ2 (x − y)θ42 (0) (3.1)

Reemplazando y = 0,

θ42 (x)θ22 (0) − θ12 (x)θ32 (0) = θ22 (x)θ42 (0)

θ32 (0)θ12 (x) θ42 (0)θ22 (x)

+ =1
θ22 (0)θ42 (x) θ22 (0)θ42 (x)
Utilizando x = u/θ32 (0), tenemos:

sn2 (u) + cn2 (u) = 1

Definición 3.2. k = θ22 (0)/θ32 (0) es el módulo de las funciones elı́pticas y k 0 =

θ42 (0)/θ32 (0) es el módulo complementario.

Proposición 3.3. Se cumple que:

1. dn2 (u) + k 2 sn2 (u) = 1

2. dn2 (u) − k 2 cn2 (u) = k 02

 30

3. k 2 + k 02 = 1

Demostración. (Lawden, 2013)

1. En la ecuación (2.20) aumentamos π/2 al valor de y. Por el Lema 2.2, obtene-

mos:
θ3 (x + y)θ3 (x − y)θ42 (0) = θ42 (x)θ32 (y) − θ12 (x)θ22 (y)
Utilizando y = 0, obtenemos:

θ42 (x)θ32 (0) = θ12 (x)θ22 (0) + θ32 (x)θ42 (0)

De donde:
θ12 (x)θ22 (0) θ32 (x)θ42 (0)
+ =1
θ42 (x)θ32 (0) θ42 (x)θ32 (0)
θ24 (0) θ12 (x)θ32 (0) θ32 (x)θ42 (0)
+ =1
θ34 (0) θ42 (x)θ22 (0) θ42 (x)θ32 (0)
Con x = u/θ32 (0), concluimos:

dn2 (u) + k 2 sn2 (u) = 1

2. En la ecuación (2.20) aumentamos π/2 al valor de x y de y. Por el Lema 2.2,

obtenemos:

θ4 (x + y)θ4 (x − y)θ42 (0) = θ32 (x)θ32 (y) − θ22 (x)θ22 (y)

Utilizando y = 0, obtenemos:

θ32 (x)θ32 (0) − θ22 (x)θ22 (0) = θ42 (x)θ42 (0) (3.2)

θ32 (x)θ32 (0) θ22 (x)θ22 (0)

− =1
θ42 (x)θ42 (0) θ42 (x)θ42 (0)
θ42 (0)θ32 (x) θ24 (0)θ42 (0)θ22 (x) θ44 (0)
− =
θ32 (0)θ42 (x) θ34 (0)θ24 (0)θ42 (x) θ34 (0)
Con x = u/θ32 (0), concluimos:

dn2 (u) − k 2 cn2 (u) = k 02

3. Utilizando x = 0 en (3.2), tenemos:

θ24 (0) + θ44 (0) = θ34 (0)

θ24 (0) θ44 (0)

+ 4 =1
θ34 (0 θ3 (0)
De donde sigue:
k 2 + k 02 = 1

Concluiremos esta sección comprobando la doble periodicidad de las funciones defi-

nidas.
 31

Proposición 3.4. Se cumple que:

1. sn(u) tiene dos periodos 2πθ32 (0) y πτ θ32 (0)

2. cn(u) tiene dos periodos 2πθ32 (0) y πθ32 (0) + πτ θ32 (0)

3. dn(u) tiene dos periodos πθ32 (0) y 2πτ θ32 (0)

Demostración. (Lawden, 2013) Usando la proposición 2.4 y las definiciones dadas

al principio del capı́tulo, tenemos:

1.
θ3 (0) θ1 (µ/θ32 (0) + 2π)
sn(u + 2πθ32 (0)) =
θ2 (0) θ4 (µ/θ32 (0) + 2π)
θ3 (0) θ1 (µ/θ32 (0))
=
θ2 (0) θ4 (µ/θ32 (0))
= sn(u)

θ3 (0) θ1 (µ/θ32 (0) + πτ )

sn(u + πτ θ32 (0)) =
θ2 (0) θ4 (µ/θ32 (0) + τ π)
θ3 (0) −λ−1 θ1 (µ/θ32 (0))
=
θ2 (0) −λ−1 θ4 (µ/θ32 (0))
= sn(u)

2.
θ3 (0) θ2 (µ/θ32 (0) + 2π)
cn(u + 2πθ32 (0)) =
θ2 (0) θ4 (µ/θ32 (0) + 2π)
θ3 (0) θ2 (µ/θ32 (0))
=
θ2 (0) θ4 (µ/θ32 (0))
= cn(u)

θ3 (0) θ2 (µ/θ32 (0) + π + πτ )

cn(u + πθ32 (0) + πτ θ32 (0)) =
θ2 (0) θ4 (µ/θ32 (0) + π + τ π)
θ3 (0) −λ−1 θ2 (µ/θ32 (0))
=
θ2 (0) −λ−1 θ4 (µ/θ32 (0))
= cn(u)

3.
θ3 (0) θ3 (µ/θ32 (0) + π)
dn(u + πθ32 (0)) =
θ2 (0) θ4 (µ/θ32 (0) + π)
θ3 (0) θ3 (µ/θ32 (0))
=
θ2 (0) θ4 (µ/θ32 (0))
= dn(u)
 32

θ3 (0) θ1 (µ/θ32 (0) + 2πτ )

dn(u + 2πτ θ32 (0)) =
θ2 (0) θ4 (µ/θ32 (0) + 2τ π)
θ3 (0) λ−2 θ1 (µ/θ32 (0))
=
θ2 (0) λ−2 θ4 (µ/θ32 (0))
= dn(u)

La proposición 3.4, nos garantiza que las funciones sn, cn y dn son funciones elı́pti-
cas ya que son funciones meromorfas con doble periodo. Para facilitar la notación,
definiremos:
1
K = πθ32 (0)
2
1
iK 0 = πτ θ32 (0)
2
Por lo que la proposición 3.4 resulta:

sn(u) = sn(u + 4K) = sn(u + 2iK 0 )

cn(u) = cn(u + 4K) = cn(u + 2K + 2K 0 )
dn(u) = dn(u + 2K) = dn(u + 4iK 0 )

3.2. Propiedades
Teorema 3.5. Las derivadas de las funciones elı́pticas de Jacobi son:
d
sn(u) = cn(u)dn(u)
du
d
cn(u) = −sn(u)dn(u)
du
d
dn(u) = −k 2 sn(u)cn(u)
du
Demostración. (Lawden, 2013) Utilizando la ecuación (2.27) al diferenciar la
definición de sn, tenemos:

d θ3 (0) θ42 (0)θ2 (z)θ3 (z)

sn(u) = = cn(u)dn(u)
du θ2 (0) θ32 (0)θ42 (z)

Utilizando la ecuación (2.28) al diferenciar la definición de cn, tenemos:

d θ4 (0) θ32 (0)θ1 (z)θ3 (z)

cn(u) = − = −sn(u)dn(u)
du θ2 (0) θ32 (0)θ42 (z)

Utilizando la ecuación (2.29) y la Definición 3.2, al diferenciar la definición de cn,

tenemos:
d θ4 (0) θ22 (0)θ1 (z)θ2 (z)
dn(u) = − 2 2
= −k 2 sn(u)cn(u)
du θ3 (0) θ3 (0)θ4 (z)
Otra propiedad importante son las correspondientes a la adición, para lo cuál primero
demostraremos las siguientes identidades:
 33

Lema 3.6. Se cumplen las siguiente identidades:

θ2 (x + y)θ4 (x − y)θ2 (0)θ4 (0) = θ2 (x)θ4 (x)θ2 (y)θ4 (y) − θ1 (x)θ3 (x)θ1 (y)θ3 (y) (3.3)
θ3 (x + y)θ4 (x − y)θ3 (0)θ4 (0) = θ3 (x)θ4 (x)θ3 (y)θ4 (y) − θ1 (x)θ2 (x)θ1 (y)θ2 (y) (3.4)
Demostración. (Lawden, 2013) Al multiplicar la ecuación (2.22) por la ecuación
(2.24), obtenemos:
θ1 (x, q)θ2 (y, q)θ3 (x, q)θ4 (y, q) = [θ1 (x + y, q 2 )θ4 (x − y, q 2 ) + θ4 (x + y, q 2 )θ1 (x − y, q 2 )]
× [θ4 (x + y, q 2 )θ4 (x − y, q 2 ) − θ1 (x + y, q 2 )θ1 (x − y, q 2 )]
(3.5)
Intercambiando x, y en la ecuación en la ecuación (3.5) y sumando con la misma
ecuación, tenemos:
θ1 (x, q)θ2 (y, q)θ3 (x, q)θ4 (y, q) + θ1 (y, q)θ2 (x, q)θ3 (y, q)θ4 (x, q)
(3.6)
= 2θ1 (x + y, q 2 )θ4 (x + y, q 2 )[θ42 (x − y, q 2 ) − θ12 (x − y, q 2 )]
Al reemplazar y = 0 en la ecuación (2.22) y en la ecuación (2.24), tenemos las
siguientes identidades:
θ1 (x, q)θ2 (0, q) = 2θ1 (x, q 2 )θ4 (x, q 2 )
θ3 (x, q)θ4 (0, q) = θ42 (x, q 2 ) − θ12 (x, q 2 )
Usando estas identidades en la ecuación (3.36), se reduce a la forma:
θ1 (x + y)θ3 (x − y)θ2 (0)θ4 (0) = θ1 (x)θ3 (x)θ2 (y)θ4 (y) + θ2 (x)θ4 (x)θ1 (y)θ3 (y)
Incrementando x por π/2, la identidad anterior se transforma en:
θ2 (x + y)θ4 (x − y)θ2 (0)θ4 (0) = θ2 (x)θ4 (x)θ2 (y)θ4 (y) − θ1 (x)θ3 (x)θ1 (y)θ3 (y)
Con lo que queda demostrada la primera identidad.
Multiplicando la ecuación (2.16) y la ecuación (2.17), encontramos que:
θ3 (x, q)θ3 (y, q)θ4 (x, q)θ4 (y, q) = θ32 (x + y, q 2 )θ32 (x − y, q 2 ) − θ22 (x + y, q 2 )θ22 (x − y, q 2 )
(3.7)
Reemplazando y = 0, en la ecuación (2.16) y (2.17), tenemos
2θ32 (x, q 2 ) = θ3 (x, q)θ3 (0, q) + θ4 (x, q)θ4 (0, q)
2θ22 (x, q 2 ) = θ3 (x, q)θ3 (0, q) − θ4 (x, q)θ4 (0, q)
Sustituyendo estas dos identidades en la ecuación (3.7), tenemos:
2θ3 (x)θ3 (y)θ4 (x)θ4 (y) = [θ3 (x + y)θ4 (x − y) + θ4 (x + y)θ3 (x − y)]θ3 (0)θ4 (0) (3.8)
Incrementando el valor de x, y por τ π/2 por el Lema 2.2, tenemos:
2θ2 (x)θ2 (y)θ1 (x)θ1 (y) = [θ4 (x + y)θ3 (x − y) − θ3 (x + y)θ4 (x − y)]θ3 (0)θ4 (0) (3.9)
Restando la ecuación (3.9) de la ecuación (3.8), encontramos:
θ3 (x + y)θ4 (x − y)θ3 (0)θ4 (0) = θ3 (x)θ4 (x)θ3 (y)θ4 (y) − θ1 (x)θ2 (x)θ1 (y)θ2 (y)
Con lo que se demuestra la segunda identidad.
 34

Teorema 3.7. Las fórmulas de adición de las funciones elı́pticas de Jacobi son:

sn(u)cn(v)dn(v) + sn(v)cn(u)dn(u)
sn(u + v) =
1 − k 2 sn2 (u)sn2 (v)
cn(u)cn(v) − sn(u)sn(v)dn(u)dn(v)
cn(u + v) =
1 − k 2 sn2 (u)sn2 (v)
dn(u)dn(v) − k 2 sn(u)sn(v)cn(u)cn(v)
dn(u + v) =
1 − k 2 sn2 (u)sn2 (v)

Demostración. (Lawden, 2013) Utilizando la definición de sn(u), tenemos que:

θ3 (0) θ1 (x + y)
sn(u + v) = .
θ2 (0) θ4 (x + y)

donde x = u/θ32 (0), y = v/θ32 (0).

Si dividimos la ecuación (2.32) para la ecuación (2.20), tenemos que:

θ1 (x + y) θ2 (0)[θ1 (x)θ4 (x)θ2 (y)θ3 (y) + θ2 (x)θ3 (x)θ1 (y)θ4 (y)]

= 4
θ4 (x + y) θ2 (0)θ3 (0)[θ42 (x)θ42 (y) − θ12 (x)θ12 (y)]

Usando las definiciones al inicio de este capı́tulo, tenemos que:

sn(u)cn(v)dn(v) + sn(v)cn(u)dn(u)
sn(u + v) =
1 − k 2 sn2 (u)sn2 (v)

Utilizando la definición de cn(u), tenemos que:

θ4 (0) θ2 (x + y)
cn(u + v) = .
θ2 (0) θ4 (x + y)

donde x = u/θ32 (0), y = v/θ32 (0).

Si dividimos la ecuación (3.3) para la ecuación (2.20), obtenemos:

θ2 (x + y) θ4 (0)[θ2 (x)θ4 (x)θ2 (y)θ4 (y) − θ1 (x)θ3 (x)θ1 (y)θ3 (y)]

=
θ4 (x + y) θ2 (0)[θ42 (x)θ42 (y) − θ12 (x)θ12 (y)]

Usando las definiciones al inicio de este capı́tulo, tenemos:

cn(u)cn(v) − sn(u)sn(v)dn(u)dn(v)
cn(u + v) =
1 − k 2 sn2 (u)sn2 (v)

Utilizando la definición de dn(u), tenemos que:

θ4 (0) θ3 (x + y)
dn(u + v) = .
θ3 (0) θ4 (x + y)

donde x = u/θ32 (0), y = v/θ32 (0).

Si dividimos la ecuación (3.4) para la ecuación (2.20), obtenemos:

θ3 (x + y) θ4 (0)[θ3 (x)θ4 (x)θ3 (y)θ4 (y) − θ1 (x)θ2 (x)θ1 (y)θ2 (y)]

=
θ4 (x + y) θ3 (0)[θ42 (x)θ42 (y) − θ12 (x)θ12 (y)]
 35

Usando las definiciones al inicio de este capı́tulo, tenemos:

dn(u)dn(v) − k 2 sn(u)sn(v)cn(u)cn(v)
dn(u + v) =
1 − k 2 sn2 (u)sn2 (v)

Con lo que se demuestra el teorema.

Utilizando el Teorema 3.5 se demuestran las siguientes identidades:

Z
1
sn(u)du = ln(dn(u) − kcn(u)) (3.10)
k
Z
1
cn(u)du = arcsin(ksn(u)) (3.11)
k
Z
dn(u)du = arcsin(sn(u)) (3.12)

Definimos entonces am(u) como:

Z u
am(u) = dn(v)dv
0

. Utilizando la ecuación (3.12), tenemos:

sn(u) = sin(am(u)) (3.13)

 36

Capı́tulo 4

Integrales elı́pticas

4.1. Integrales elı́pticas de primera especie

A lo largo de este capı́tulo, el argumento u y el módulo k de todas las funciones
elı́pticas será asumido real y, también, supondremos que 0 < k < 1, si no se establece
otra cosa.
Empezaremos por calcular la derivada de la función inversa sn−1 (x, k). En este
capı́tulo todas las funciones inversas serán tomadas en el rango (0, K), esto se hace
para que los valores que tome estén bien definidos evitando la doble periodicidad.
Con u = sn−1 (x, k), entonces x = snu, y por lo tanto usando el Teorema 3.5, la
Proposición 3.1 y la Proposición 3.3:
dx p
= cn(u)dn(u) = (1 − x2 )(1 − k 2 x2 )
du
Integrando esta ecuación en el rango (0, x) (0 ≤ x ≤ 1) para x y el rango correspon-
diente (0, u) (0 ≤ u ≤ K) para u, da como resultado:
Z x
−1 dt
sn (x, k) = u = p (4.1)
0 (1 − t )(1 − k 2 t2 )
2

para (0 ≤ x ≤ 1). Esta integral es llamada una integral elı́ptica de primera

especie.
Dado que sn−1 = K, un caso especial de la ecuación (4.1) es:
Z 1
dt
K= p (4.2)
0 (1 − t2 )(1 − k 2 t2 )

Esta es la integral elı́ptica completa de primera especie.

Ya que K 0 (k) = K(k 0 ), también tenemos
Z 1
0 dt
K (k) = p
0 (1 − t2 )(1 − k 2 t2 )

Si reemplazamos t = sin α en la ecuación (4.1), esta se reduce a la forma:

Z φ
−1 dα
F (φ, k) = sn (sin φ, k) = p (4.3)
0 1 − k 2 sin2 α
 37

donde x = sin φ. Si reemplazamos φ = π/2, tenemos:

Z π/2
dα
K(k) = p (4.4)
0 1 − k 2 sin2 α
Esta ecuación se puede expandir en potencias ascendentes de k 2 , y puede ser inte-
grado término a término:
"  2  2  2 #
π 1 1,3 1,3,5
K= 1+ k2 + k4 + k 6 − ... (4.5)
2 2 2,4 2,4,6

Si tenemos u = cn−1 (x), tal que x = cn(u). Entonces

d p
x = −sn(u)dn(u) = − (1 − x2 )(k 02 + k 2 x2 )
du
Integrando en el intervalo (0, u) para u y (1, x) para x nos lleva al resultado:
Z 1
−1 dt
cn (x, k) = 2 02 2 2
(4.6)
x (1 − t )(k + k t )

4.2. Integrales elı́pticas de segunda especie

La integral elı́ptica de segunda especie está definida por:
Z xr
1 − k 2 t2
dt
0 1 − t2
Si tenemos t = sin φ, entonces:
Z φ q
E(k, φ) = 1 − k 2 sin2 φdφ
0

Además, si utilizamos sn(u) = t = sin φ tenemos que cualquier integral elı́ptica de

la segunda especie puede ser reducida a evaluar la función elı́ptica de Jacobi
E(u, k) definida como:
s
Z u Z τ
1 − k 2 t2
E(u, k) = dn2 vdv = dt
0 0 1 − t2

donde t = sn(v) y τ = sn(u)

La forma de Legendre de esta integral se deriva al hacer la transformación

sn(v) = sin θ, sn(u) = sin φ. Viendo como una función de φ y k, la integral será
denotada por D(φ, k). Entonces
Z φq
D(φ, k) = (1 − k 2 sin2 θ)dθ
0

Por la ecuación (3.13), φ = am(u), la relación entre las dos notaciones es expresada
por:
E(u, k) = D(am(u), k) (4.7)
 38

La integral elı́ptica completa de segunda especie denotada por E es definida

por la ecuación:
Z K Z π/2 p
2
E = E(k) = dn (u)du = E(K, k) = 1 − k 2 sin2 θdθ (4.8)
0 0

Esta ecuación puede ser expandida en potencias de k 2 , y después ser integrada

término a término y tenemos como resultado:
"  2  2  2 #
π 1 1 1,3 1 1,3,5
E= 1− k2 − k4 − k 6 − ... (4.9)
2 2 3 2,4 5 2,4,6
Teorema 4.1. La fórmula de adición de la función elı́ptica de segunda especie es:
E(u1 + u2 ) = E(u1 ) + E(u2 ) − k 2 sn(u1 )sn(u2 )sn(u1 + u2 )
Demostración. (Bowman, 1953)Por la fórmula de adición de dn en el Teorema
3.7, tenemos que:
2dn(x)dn(y)
dn(x + y) + dn(x − y) =
1 − k 2 sn2 (x)sn2 (y)
−2k 2 sn(x)sn(y)cn(x)cn(y)
dn(x + y) − dn(x − y) =
1 − k 2 sn2 (x)sn2 (y)
Multiplicando ambas identidades obtenemos:
2 2 −4k 2 sn(x)sn(y)cn(x)cn(y)dn(x)dn(y)
dn (x + y) − dn (x − y) =
(1 − k 2 sn2 (x)sn2 (y))2
y, integrando ambos lados con respecto a y junto con el Teorema 3.5, resulta:
2sn(x)cn(x)dn(x)
E(x + y) + E(x − y) = − +C (4.10)
sn2 (x)(1− k 2 sn2 (x)sn2 (y))
donde C es una constante que depende solo de x. Si y = x, tenemos:
2sn(x)cn(x)dn(x)
E(2x) = − +C
sn2 (x)(1 − k 2 sn4 (x))
y restando esta ecuación de la ecuación (4.10) tenemos:
 2
sn2 (x) − sn2 (y)
 
2k sn(x)cn(x)dn(x)
E(x + y) + E(x − y) − E(2x) =
1 − k 2 sn4 (x) 1 − k 2 sn2 (x)sn2 (y)
Usando nuevamente el Teorema 3.7 tenemos que:
sn2 (x) − sn2 (y)
sn(x + y)sn(x − y) = (4.11)
1 − k 2 sn2 (x)sn2 (y)
De igual manera,
2sn(x)cn(x)dn(x)
sn(2x) = (4.12)
1 − k 2 sn4 (x)
Por lo tanto, utilizando las ecuaciones (4.11) y (4.12), tenemos:
E(x + y) + E(x − y) − E(2x) = k 2 sn(2x)sn(x + y)sn(x − y)
Finalmente, si utilizamos u1 = x + y, u2 = x − y, entonces la ecuación anterior
resulta:
E(u1 + u2 ) = E(u1 ) + E(u2 ) − k 2 sn(u1 )sn(u2 )sn(u1 + u2 )
 39

Capı́tulo 5

Aplicación de las funciones

elı́pticas de Jacobi

5.1. El péndulo simple

Definición 5.1. Un péndulo simple consiste en una partı́cula suspendida de un
hilo inextensible y sin peso, moviéndose libremente en un plano vertical bajo la
fuerza de la gravedad.

Sea l la longitud del hilo de suspensión, g la aceleración gravitacional, y m la

masa de la partı́cula. Entonces, si β es el ángulo que forma el hilo con el eje vertical
y v es la velocidad de la masa en cualquier tiempo t, por la conservación de energı́a
tenemos:
1 2
mv − mgl cos β = constante
2
Dado que v = lβ̇, con ω 2 = g/l, esta ecuación puede ser escrita como:

β̇ 2 − 2ω 2 cos β = constante

donde β̇ es la derivada de β con respecto al tiempo. Sea α la amplitud de la oscilación,

entonces β̇ = 0 cuando β = α, por lo tanto
    
2 2 2 2 1 2 1
β̇ = 2ω (cos β − cos α) = 4ω sin α − sin β
2 2

Integrando, tenemos el siguiente resultado:

1 β
Z
dβ
ωt = (5.1)
2 0 sin 2 α − sin2
2 1 1




2
β

donde t es el tiempo que toma la partı́cula en ir desde la posición más baja hasta
llegar a la posición en la cual el hilo está inclinado a un ángulo β con respecto al eje
vertical.
Haciendo un cambio de variable a φ, donde
1 1
sin β = sin α sin φ
2 2
 40

Usando la ecuación (4.3) de la integral elı́ptica de primera especie, tenemos:

1 φ
Z  
dφ −1 1
ωt = q
= sn sin φ, sin α
2 0 2 1 2
1 − sin 2 α sin φ 2

Invirtiendo obtenemos:  
1
sin φ = sn ωt, sin α
2
Entonces:
1 1
sin β = sin αsn(ωt, k) (5.2)
2 2
donde k = sin 12 α mostrando que sin 21 β oscila con amplitud sin 12 α y periodo

T = 4K/ω

donde Z π/2
dφ
K= q
1 − sin2 12 α sin2 φ


0

notemos que:
s 
1 2 1 2
p
cos β = 1 − sin α sin φ = (1 − k 2 sn2 ωt) = dnωt (5.3)
2 2

Diferenciando la ecuación (5.2), tenemos:

1 1 1
β̇ cos β = ω sin αcnωtdnωt
2 2 2
Y utilizando la ecuación (5.3):
1
β̇ = 2ω sin αcnωt
2
(Lawden, 2013).

5.2. La ecuación de Duffing

Definición 5.2. La ecuación de Duffing gobierna las oscilaciones de una masa
adherida al final de un resorte cuya tensión (o compresión) T está relacionada con
su extensión x por una ecuación de la forma.

T = αx + βx3

donde α es siempre positivo.

Si β = 0, el resorte obedece la ley de Hooke y las oscilaciones son armónicas

simples. Si β > 0, la tensión incrementa con la extensión más rápido que en la ley
de Hooke y el resorte se dice duro. Si β < 0, la tensión incrementa más lento que en
 41

la ley de Hooke y el resorte se dice suave.

La ecuación del movimiento de la masa puede ser puesta en la forma:

ẍ + x + x3 = 0

que es la forma canónica de la ecuación de Duffing.

Consideremos, primero el caso del resorte duro, para el cual  > 0. Supongamos que
d 1 2
en el momento inicial tenemos: t = 0, x = a, ẋ = 0. Dado que ẍ = dx 2
ẋ , podemos
integrar con respecto a x para obtener:
1 2 1 2 1 4 1 2 1 4
ẋ + x + x = a + a
2 2 4 2 4
Integrando con respecto a t, obtenemos:
r Z a
2 dx
t= p
 x (a2 − x2 )(2/ + a2 + x2 )
" s #
1 x a2
=p cn−1 ,
2
(1 + a ) a 2 + 2a2

Haciendo referencia a la ecuación (4.6). Invirtiendo tenemos que:

p
x = acn( (1 + a2 )t) (5.4)

con módulo:
a2
k2 =
2 + 2a2
Si trabajamos con el primer orden de . Entonces, k 2 = 21 a2 y la ecuación (4.5) nos
dice que K = 12 π(1 + 18 a2 ). Y en ese caso tenemos:

3
T = 2π(1 − a2 )
8
Para un resorte√suave,  < 0 y escribiremos  = −η. Ya que x − ηx3 tiene el
máximo en x = 1/ 3η, teóricamente la tensión deberı́a disminuir mientras la ex-
tensión aumenta para valores suficientemente
√ grandes de x, sin embargo este no es
un caso real, por lo que asumiremos a < 1/ 3η.

Teniendo, x = 0 en t = 0, la ecuación para el tiempo es:

r Z x
2 dx
t= p
η 0 (a2 − x2 )(2/η − a2 − x2 )
r " s #
2 x ηa2
= sn−1 ,
2 − ηa2 a 2 − ηa2

haciendo referencia la ecuación (4.1). Invirtiendo la última ecuación obtenemos:

r
1
x = asn( (1 − ηa2 )t)
2
 42

donde el módulo es determinado por

ηa2
k2 =
2 − ηa2
Por lo tanto el periodo de oscilación está dado por:
4K
T =q
(1 − 12 ηa2 )

y, para pequeños η consideraremos solo el primer grado de η, esto junto a la ecuación

(4.5) nos lleva a
3
T = 2π(1 + ηa2 )
8
(Lawden, 2013).

5.3. Geometrı́a de la elipse

Tomando una elipse parametrizada como:

x = a sin θ, y = b cos θ

con a > b y el ángulo excéntrico θ es medido desde el eje menor. Si s es la longitud de

arco medido en sentido
p horario alrededor
p de la curva empezando desde
p el punto 2(0, b),
2 2 2
entonces: ds
p= (dx + dy ) = (a cos θ + b sin θ)dθ = a (1 − e2 sin θ)dθ
2 2 2

donde e = (1 − b2 /a2 ) es la excentricidad. Entonces, la longitud de arco desde el

punto (0, b) a cualquier punto P donde θ = φ está dado por:
Z φ q
s=a (1 − e2 sin2 θ)dθ = aE(u, e)
0

donde φ = am(u, e). Esto es por la ecuación (4.7)

Si ponemos φ = 12 π, u = K tenemos aE para la longitud de un cuadrante de la
elipse.
Si introducimos un parámetro alternativo u al poner θ = am(u, e), entonces por la
ecuación (3.13), la forma paramétrica puede ser expresada como:

x = asnθ, y = bcsθ

Esto nos muestra como las funciones elı́pticas nos ayudan a describir objetos geométri-
cos (Lawden, 2013).
 43

Conclusiones

Se ha presentado una aproximación histórica a las funciones elı́pticas por medio

de su construcción usando funciones theta. Estas funciones fueron encontradas co-
mo lo hizo Jacobi y ası́ llegamos al origen del estudio sobre las funciones elı́pticas.
Además, gracias a estas funciones se han podido resolver de manera analı́tica pro-
blemas en geometrı́a y fı́sica. Es importante además que las derivadas e integrales
de las funciones elı́pticas de Jacobi sn, cn y dn pueden ser expresadas en base a sı́
mismas.

En el primer capı́tulo hemos trabajado analı́ticamente mostrando las propiedades

que tienen las funciones elı́pticas. La relevancia de esta parte está dada por su aplica-
ción sobre curvas elı́pticas. Los capı́tulos siguientes estudian de manera constructiva
cómo hallar funciones elı́pticas y demuestran sus propiedades. Esta perspectiva, per-
mite que el tema de funciones elı́pticas pueda ser entendido con facilidad aún sin
mucho conocimiento de análisis complejo, y se pueda apreciar completamente las
aplicaciones en el capı́tulo final.

Uno de los objetivos de este trabajo fue presentar aplicaciones de las funciones
elı́pticas por medio del trabajo de Jacobi. Sin embargo para aplicaciones recientes
es necesario un paso más profundo dentro del estudio de estas funciones. El estudio
de las funciones elı́pticas de Weistrass es el siguiente paso para poder comprender
mejor las aplicaciones que tienen estas funciones en teorı́a de números que se ha
convertido en una de las más importantes aplicaciones de las funciones elı́pticas.
 44

Bibliografı́a

Bowman, F. (1953). Introduction to Elliptic Functions: With Applications. English

Universities Press.
Brochero, F., Moreira, C., Saldanha, N., y Tengan, E. (2015). Teoria dos números
- um passeio pelo mundo inteiro com primos e outros números familiares.
Projeto Euclides, IMPA.
Copson, E. (1935). An Introduction to the Theory of Functions of a Complex
Variable. Clarendon Press.
Jones, G., y Singerman, D. (1987). Complex Functions: An Algebraic and Geometric
Viewpoint. Cambridge University Press.
Lawden, D. (2013). Elliptic Functions and Applications. Springer New York.
Stein, E., y Shakarchi, R. (2010). Complex Analysis (n.o v. 13). Princeton University
Press.