Vol IV Num 1 Completo PDF
Vol IV Num 1 Completo PDF
Vol IV Num 1 Completo PDF
Se permite la reproducción parcial o total de los contenidos de la publicación para fines educativos, dándose
el debido crédito a sus autores y a la propia revista. Se prohibe, sin embargo, la reproducción parcial o total
de este texto por cualquier medio o formato incluyendo el electrónico, con fines lucrativos.
Revista Pensamiento Matemático
Grupo de Innovación Educativa Pensamiento Matemático
y
Grupo de Investigación Matemática Aplicada a la Ingeniería Civil
Universidad Politécnica de Madrid
Comité Científico
Mariló López González, Adela Salvador Alcaide, Sagrario Lantarón Sánchez, Ascensión Moratalla de la Hoz,
Javier Rodrigo Hitos, José Manuel Sánchez Muñoz, Raquel Caro Carretero, Fernando Chamizo Lorente,
Luis Garmendia Salvador, José Juan de Sanjosé Blasco, Arthur Pewsey, Alfonso Garmendia Salvador,
Fernanda Ramos Rodríguez, Milagros Latasa Asso, Nieves Zuasti Soravilla
1 de abril de 2014
I
Índice de Artículos
Editorial del Número 1 (Vol. IV) . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
Experiencias Docentes
El abandono académico: análisis y propuestas paliativas. Dos proyectos de la Universidad
Politécnica de Madrid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
Ana Casaravilla
Historias de Matemáticas
Los curiosos problemas de mezclas de Alcuino de York . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
Antonio M. Oller Marcén
Aritmética de los valores de la función zeta de Riemann en argumentos enteros . . . . . . . . . . . . 33
Anier Soria Lorente
Cuentos Matemáticos
Matemáticas para todos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
José Miguel Bel Martínez
Investigación
Arquitectura paramétrica discreta: exploración en el ámbito de la geometría ortogonal . . . . . 49
Oscar del Castillo Sánchez
Cálculo Fraccionario y dinámica newtoniana . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77
Antón Lombardero Ozores
Programación en paralelo GPGPU del método en diferencias finitas FTCS para la ecuación
del calor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez
Estimación de sucesos raros mediante Redes Bayesianas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 127
Francisco Soler Flores y José Ángel Olivas Varela
Juegos Matemáticos
Números triangulares cuadrados (la cuadratura del triángulo) . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
Dionisio Pérez Esteban
Una Sucesión Matemática Curiosa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 141
Marco Vinicio Vásquez Bernal
Críticas
Informe sobre el libro: “El enigma de Fermat. Tres siglos de desafío a la matemática”,
Albert Violant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 149
Javier Rodrigo Hitos
Entrevistas
Pablo Amster: La Música de la Matemática . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 157
Rosa María Herrera
III
Editorial del Número 1 (Vol. IV)
Equipo Editorial
Revista de Investigación
1 de abril de 2014
Resumen
Presentamos el primer número del cuarto volumen de la Revista “Pensamiento Matemá-
tico” que nuevamente incluye diversos artículos encuadrados en las diferentes secciones ha-
bituales de la Revista, cubriendo un amplio espectro de la Matemática y sus aplicaciones.
Abstract
We present the first number of the fourth volume of the “Mathematical Thinking” Journal.
It includes some papers framed in the different sections of the Journal. The articles cover a
wide range of the Mathematics topics and its applications.
Introducción
Comenzamos nuestro cuarto año de publicación con el primer número del año 2014. Cuatro
años de esfuerzo, trabajo, dedicación y perseverancia, características éstas con las que todo el
equipo que hace PM se siente identificado.
Queremos aprovechar este número para agradecer el esfuerzo de todos aquellos compañe-
ros que, de una manera u otra, colaboran desinteresadamente con esa Revista. Nos reconforta
en gran medida saber que nuestro proyecto resulta cada vez de mayor interés para la comu-
nidad docente con cada nuevo número publicado, creciendo cada día en número de lectores y
colaboradores.
Seguimos con el espíritu de que Pensamiento Matemático esté dirigida a un público lo más
amplio posible pretendiendo introducir a sus lectores en la belleza y aplicaciones de las matemá-
ticas. Por ello en este número se encuentran artículos que nuevamente tratan diversos campos
de la Matemática, sus aplicaciones y curiosidades.
A continuación presentamos brevemente los artículos que componen el presente número en
sus diferentes secciones.
Experiencias Docentes
El artículo “El abandono académico: análisis y propuestas paliativas. Dos proyectos de la Universi-
dad Politécnica de Madrid” presenta el abandono de los estudios superiores, el cual preocupa en
las universidades de todo el mundo y tiene una incidencia relevante en las titulaciones de in-
geniería y arquitectura en España. Son múltiples los análisis y estudios sobre el abandono que,
1
Equipo Editorial Editorial
con diversos enfoques, pueden encontrarse en la literatura. Sin embargo, no hay unidad de cri-
terio para conceptualizar el fenómeno, establecer medidas que permitan comparar la dimensión
real del abandono en instituciones, países y situaciones distintas y, consecuentemente, propo-
ner medidas eficientes para mejorar las tasas de permanencia y minimizar las consecuencias
negativas –tanto de carácter personal como social– del abandono de estudios superiores a nivel
general. Se presentan aquí dos proyectos promovidos por la Universidad Politécnica de Madrid
(UPM) cuyo objetivo es profundizar en el conocimiento del abandono y determinar propuestas
paliativas con aplicabilidad en diferentes contextos. Uno de ellos es un Proyecto Transversal de
Innovación Educativa, de carácter interno a la UPM, que involucra a un gran número de profe-
sores y centros de esta universidad. El otro es un Proyecto Internacional financiado por la Unión
Europea, en el que participan veinte universidades de dieciséis países, doce de América Latina
y cuatro europeos, que lidera la UPM. Ambos coinciden en el tiempo (2011-2014) y comparten,
entre otras características, una visión integradora del abandono como fenómeno que debe ser
abordado desde la perspectiva de responsabilidad conjunta del estudiante y de la institución.
Historias de Matemáticas
En “Los curiosos problemas de mezclas de Alcuino de York” se presentan los problemas “de mez-
clas” contenidos en la colección medieval Propositiones ad acuendos juvenes de Alcuino de York
y se analizan posibles métodos de resolución, puramente aritméticos, que pudo haber segui-
do el autor del texto original. Pensamos que un análisis de este tipo puede tener interesantes
implicaciones didácticas a la hora de afrontar el paso de la Aritmética al Álgebra.
Rabano Mauro apoyado por Alcuino ofrece una obra a Otgar de Maguncia.
“Aritmética de los valores de la función zeta de Riemann en argumentos enteros” trata sobre el
carácter aritmético de los valores de la función zeta de Riemann en argumentos enteros y en
particular en impares, el cual hoy día aún continúa siendo un problema abierto dentro de la
comunidad matemática. Este artículo se dedica a presentar los principales resultados alcanza-
dos por varios matemáticos desde el siglo XVII hasta la actualidad, correspondientes al carácter
aritmético que siguen los valores de la función zeta de Riemann en argumentos enteros.
1 1 1
ζ ( s) = ∑ = 1+ s + s +··· ,
n ≥1
n s 2 3
Cuentos Matemáticos
“Matemáticas para todos” es un relato que trata de poner de manifiesto cómo tras la “ignoran-
cia” matemática puede haber atrevimiento por el que se puede llegar a cualquier resultado por
extraño que éste resulte, todo ello llevado en un tono narrativo ameno y divertido.
Investigación
“Arquitectura paramétrica discreta: exploración en el ámbito de la geometría ortogonal” trata sobre
la idea de arquitectura paramétrica, la cual tiene una vertiente instrumental en el uso de softwa-
re orientado a la gestión de la información del proyecto arquitectónico, en particular cuando se
emplean geometrías complejas. Sin embargo, la idea de que la forma, sin importar su geometría,
esté matemáticamente definida mediante un número limitado de parámetros puede ser de gran
alcance en lo que se refiere a las posibilidades de generación de variaciones por modificación
de dichos parámetros, manipulación de la forma en procesos de optimización y por la relación
entre las estructuras matemático-formales producidas y la idea clásica de unidad formal de la
obra de arte.
“Ciudad – Alea”; imagen obtenida mediante generación de una malla tridimensional consistente en una alternancia a la que se
aplica una eliminación aleatoria de planos.
“Cálculo Fraccionario y dinámica newtoniana” presenta el Cálculo Fraccionario como una útil
herramienta que abre la posibilidad de extender los conceptos de derivada e integral a órdenes
no enteros. En este contexto, cabe plantearse una generalización de la segunda Ley de Newton
en la que se sustituya la derivada segunda interviniente por otra de orden α (1, 5 < α ≤ 2). En
este trabajo se estudian algunos sistemas clásicos (péndulo, proyectil y resorte) bajo este nuevo
enfoque.
O
L
T
θ
mg sen θ
mg cos θ
mg
Fuerzas que intervienen en el péndulo
En “Programación en paralelo GPGPU del método en diferencias finitas FTCS para la ecuación del
calor” se muestran las ventajas de la programación en paralelo resolviendo numéricamente la
ecuación del calor en dos dimensiones a través del método de diferencias finitas explícito cen-
trado en el espacio FTCS. De las conclusiones de este trabajo se pone de manifiesto la impor-
tancia de la programación en paralelo para tratar problemas grandes, en los que se requiere un
elevado número de cálculos, para los cuales la programación secuencial resulta impracticable
por el elevado tiempo de ejecución. En la primera sección se describe brevemente los conceptos
básicos de programación en paralelo. Seguidamente se resume el método de diferencias finitas
explícito centrado en el espacio FTCS aplicado a la ecuación parabólica del calor. Seguidamente
se describe el problema de condiciones de contorno y valores iniciales específico al que se va
a aplicar el método de diferencias finitas FTCS, proporcionando pseudocódigos de una imple-
mentación secuencial y dos implementaciones en paralelo. Finalmente tras la discusión de los
resultados se presentan algunas conclusiones.
δΩ
Ω
z B
δB
x
Esfera arbitraria B contenida en una región Ω. δB denota el límite de B.
En “Estimación de sucesos raros mediante Redes Bayesianas”podremos ver cómo la teoría y apli-
caciones de sucesos de baja probabilidad de ocurrencia o sucesos raros está en boga en los
últimos años debido a su importancia práctica en muy diferentes campos, como seguros, fi-
nanzas, ingeniería o ciencias del medio ambiente. Este trabajo presenta una metodología para
predecir sucesos raros basada en Redes Bayesianas que permite a su vez estudiar situaciones
alternativas.
0.25
caso 1
caso 2
caso 3
caso 4
caso 5
0.2
0.15
Probabilidad
0.1
0.05
sucesos raros
0
0 2 4 6 8 10 12 14 16
Sucesos raros y sus distribuciones según los valores de las variables predictoras
Juegos Matemáticos
En “Número triangulares cuadrados (la cuadratura del triángulo)” se pone de manifiesto cómo
entre los números triangulares, se encuentran camuflados algunos cuadrados. En este artículo
se les desenmascara, y en el proceso nos topamos con algunos personajes interesantes: fraccio-
nes continuas, ecuación de Pell, aproximación de soluciones, etc.
En “Una Sucesión Matemática Curiosa” se pone de manifiesto que existe una sucesión de nú-
meros naturales, que inician en cero y va aumentando cada vez el entero correspondiente al
orden del elemento de la misma, que tiene un comportamiento a primera vista curioso, mas al
realizar un estudio de la misma se entiende a cabalidad este comportamiento. La sucesión gene-
ra además ciertos resultados interesantes, como el poder determinar resultados para la potencia
de cualquier número entero, en función de sus elementos.
Críticas
En “Informe sobre el libro: ‘El enigma de Fermat. Tres siglos de desafío a la matemática´, Albert
Violant” se presenta un informe sobre dicho libro, perteneciente a la colección “El mundo es
matemático”, de la que ya se comentó otro libro en un artículo anterior del autor.
Entrevistas
“Soy matemático. En una época era músico, o al menos eso creía. Después empecé a decir que era buen
matemático entre los músicos y buen músico entre los matemáticos. Pero ahora me asumo como matemá-
tico, aunque siga haciendo algo de música. Quizás la diferencia es que ahora que no me siento músico la
disfruto más”. Este es el autorretrato arrancado a Pablo Amster al empezar este trabajo, cono-
ciéndole un poco esto es muy escaso, pero es una herramienta válida para poner en contacto al
lector y al matemático.
Pablo Amster.
No queremos cerrar esta breve introducción a nuestro nuevo número sin agradecer a la
Universidad Politécnica de Madrid el haber otorgado al GIE Pensamiento Matemático el Premio
al Mejor Grupo de Innovación Educativa 2013 de la Universidad. Son estos reconocimientos
los que animan a los equipos a seguir trabajando en la difusión, investigación y enseñanza de
cualquier disciplina, por ello muchas gracias.
“La perseverancia es el trabajo duro que hay que hacer cuando uno se cansa del trabajo duro
que ya hizo.”
Newt Gingrich
“Hay una fuerza motriz más poderosa que el vapor, la electricidad y la energía atómica: la
voluntad.”
Albert Einstein
H.W. Arnold
El Comité Editorial
1 de abril de 2014
Resumen
El abandono de los estudios superiores preocupa en las universidades de todo el
mundo y tiene una incidencia relevante en las titulaciones de ingeniería y arquitectura en
España. Son múltiples los análisis y estudios sobre el abandono que, con diversos enfoques,
pueden encontrarse en la literatura. Sin embargo, no hay unidad de criterio para
conceptualizar el fenómeno, establecer medidas que permitan comparar la dimensión real
del abandono en instituciones, países y situaciones distintas y, consecuentemente, proponer
medidas eficientes para mejorar las tasas de permanencia y minimizar las consecuencias
negativas –tanto de carácter personal como social– del abandono de estudios superiores a
nivel general. Se presentan aquí dos proyectos promovidos por la Universidad Politécnica
de Madrid (UPM) cuyo objetivo es profundizar en el conocimiento del abandono y
determinar propuestas paliativas con aplicabilidad en diferentes contextos. Uno de ellos es
un Proyecto Transversal de Innovación Educativa, de carácter interno a la UPM, que
involucra a un gran número de profesores y centros de esta universidad. El otro es un
Proyecto Internacional financiado por la Unión Europea, en el que participan veinte
universidades de dieciséis países, doce de América Latina y cuatro europeos, que lidera la
UPM. Ambos coinciden en el tiempo (2011-2014) y comparten, entre otras características,
una visión integradora del abandono como fenómeno que debe ser abordado desde la
perspectiva de responsabilidad conjunta del estudiante y de la institución.
7
Ana Casaravilla Experiencias Docentes
Abstract
Dropout in higher education is a concern in universities worldwide and has a
significant impact in our country, especially in engineering and architecture degrees. There
are multiple and different perspectives to analyze this phenomenon; however, there is not a
unified approach to conceptualize it. Adequate ways to measure and compare the dropout
dimension in different institutions or countries have not been found so far. Consequently, it
is difficult to find effective procedures to improve retention rates and to minimize negative
consequences -both personal and social- of the abandonment of higher education. This
article is about two projects which have been promoted by the Technical University of
Madrid (UPM) whose main objective is to deepen the knowledge of academic attrition and
to determine palliative approaches to be applied in different contexts. One of them is a
UPM Inner Transverse Educational Innovation Project, involving a large number of
teachers and schools of the UPM. The second one is an International Project lead by the
UPM and funded by the European Union. This project involves twenty universities and
sixteen countries, twelve in Latin America and four in Europe. Both projects have been
developed simultaneously (2011-2014) and share, among other characteristics, the
integrated viewpoint that the dropout phenomenon should be addressed by considering
both the student and the institution responsibility.
1. Introducción
El abandono de los estudios superiores preocupa en las universidades de todo el mundo y
tiene una incidencia relevante en nuestro país, y muy especialmente, en las titulaciones de
ingeniería y arquitectura. Son múltiples los análisis sobre el abandono que, con diferentes
propósitos, pueden encontrarse en la literatura especializada. En ellos se estudian las causas
que influyen en el abandono de estudios, se elaboran modelos de predicción con el objetivo de
medir el riesgo de abandono, o se presentan diversas teorías que, a lo largo del tiempo y con
diferentes enfoques, han intentado explicar los motivos y circunstancias por los que un
estudiante decide dejar los estudios que comenzó, en un cierto momento. Estos análisis suelen
presentar el abandono como un hecho puntual, que se produce al desvincularse el alumno de
los estudios iniciado. Son, por ello, en exceso simplistas ya que no toman en consideración lo
que sucede después de que el alumno abandona la carrera, ni las consecuencias de este hecho
en su futuro académico, profesional y personal.
Por otro lado, no hay unidad de criterio a la hora de conceptualizar el fenómeno (ni tan
siquiera una definición de “abandono” generalmente asumida por la comunidad
universitaria). Tampoco se han podido establecer –hasta ahora– formas de medición que
permitan comparar la dimensión real del abandono en instituciones, países y situaciones
distintas. Consecuentemente, ante esta falta de criterios homogéneos, no es posible proponer
medidas eficientes de carácter general para mejorar las tasas de permanencia y minimizar las
consecuencias negativas –tanto de carácter personal como social– del abandono de estudios
superiores.
Se presentan aquí dos proyectos promovidos por la Universidad Politécnica de Madrid
(UPM) cuyo objetivo común es doble. En primer lugar, profundizar en el conocimiento del
abandono en toda su complejidad, presentando un marco conceptual, ciertos instrumentos y
modelos de predicción, y algunas pautas para el estudio del fenómeno, con aplicabilidad en
diferentes contextos. Y, en segundo término, ofrecer propuestas, suficientemente contrastadas,
para la reducción de los índices de abandono en nuestras universidades.
Uno de ellos es un Proyecto Transversal de Innovación Educativa, de carácter interno a la
UPM, que involucra a un gran número de profesores y centros de esta universidad. El
segundo es un Proyecto Internacional, financiado por la Unión Europea, en el que participan
veinte universidades de dieciséis países, doce de América Latina y cuatro europeos, que lidera
asimismo la UPM. Ambos proyectos coinciden en el tiempo (2011-2014) y comparten, entre
otras características, una visión integradora del abandono como fenómeno que debe ser
abordado desde la perspectiva de la responsabilidad conjunta del estudiante y de la
institución.
Figura 1. Portal GUIA, referencia para el análisis del abandono en la educación superior en Latinoamérica.
1Se considera que un estudiante matriculado en una titulación, en una universidad, en el curso n,
“abandona”, cuando no se matricula en la misma titulación en el curso n+1.
Entre ellos, un caso particular es el de los alumnos que se matriculan el curso n+1 en otra titulación de la
misma universidad. A éstos se les considera “reubicados” dentro de la propia institución.
Figura 2. “Análisis del Absentismo y Abandono en las titulaciones de grado de la UPM y propuestas para la mejora de los
índices de permanencia”, libro digital publicado en CDP.
del primer año como indicador de calidad, y se define un protocolo para la implantación,
seguimiento y evaluación de las actuaciones dirigidas a su mejora.
Figura 3. Presentación de la ponencia ”Protocolo de calidad para la tasa de permanencia a un año en la UPM”, por Ana
Casaravilla en III-CLABES. Noviembre 2013.
Referencias
[1] ARRIAGA, J. et al. La acción académica posterior tomada por el estudiante que abandona.
Elemento necesario para un “Sistema nacional de registro de Información sobre abandono”. Actas
de la Segunda Conferencia Latinoamericana sobre Abandono en la Educación Superior
(II-CLABES), pp. 83-90. Madrid. 2012.
[3] COROMINAS, E. La transición a los estudios universitarios. Abandono o cambio en el primer año
de Universidad, pp.127-151. Revista de Investigación Educativa, 19,1. Murcia. 2001.
[4] GARCÍA, A. et al. Protocolo de calidad para la tasa de permanencia a un año en la Universidad
Politécnica de Madrid. Actas de la Tercera Conferencia Latinoamericana sobre Abandono
en la Educación Superior (III-CLABES), pp. 985-995. Madrid. 2013.
[5] GLENN, P. A. & RYAN, M. P. Increasing one-year retention rates by focusing on academic
competence: an empirical odyssey. Journal of College Student Retention, v.4, n.3, pp. 297-324.
Minnesota. 2003.
[7] ROHR S. L. How well does the SAT and GPA predict the retention of science, technology,
engineering, mathematics, and business students. Journal of College Student Retention, v.14, n.2,
pp.195-208. Minnesota. 2012.
[8] VILLAR, A. Del abandono de los estudios a la reubicación universitaria. Revista de la Asociación de
Sociología de la Educación (RASE), 3(2), pp.267-283. Madrid. 2011.
[9] VRIES, W. et al. ¿Desertores o decepcionados? Distintas causas para abandonar los estudios
universitarios. Revista de la educación superior. Asociación Nacional de Universidades e
Instituciones de Educación Superior de México (ANUIES) , XL(4), n.160, pp.31-49. México.
2011.
[10] COMISIÓN EUROPEA. ESTRATEGIA EUROPA 2020. 2010. Accesible el 27 de septiembre de 2013 en
http://www.mecd.gob.es/inee/publicaciones/indicadores-educativos/Indicadores-
Internacionales/UE.html#Objetivos_ET2020_2011
Sobre la autora:
Nombre: Ana Casaravilla.
Directora adjunta del proyecto GUIA. Catedrática de Bachillerato y Profesora de la ETSEM-
UPM (jubilada). Premio de Innovación Educativa UPM 2007 y Mención Especial en los
premios UPM de Innovación Educativa 2012.
Correo Electrónico: [email protected]
Institución: Universidad Politécnica de Madrid, España.
1 de abril de 2014
Resumen
En este trabajo se presentan los problemas “de mezclas” contenidos en la colección
medieval Propositiones ad acuendos juvenes de Alcuino de York y se analizan posibles
métodos de resolución, puramente aritméticos, que pudo haber seguido el autor del texto
original. Pensamos que un análisis de este tipo puede tener interesantes implicaciones
didácticas a la hora de afrontar el paso de la Aritmética al Álgebra.
Abstract
In this work we present the “mixing” problems contained in the medieval collection
Propositiones ad acuendos juvenes by Alcuin of York and possible, purely arithmetic, solving
methods that could have been used by the author are analyzed. We think that such an
analysis can have interesting didactical implications when facing the transition from
Arithmetic to Algebra.
1. Introducción
El Álgebra, entendida como el arte de resolver ecuaciones, es posiblemente uno de los
mayores avances en la historia de las Matemáticas. Problemas de enunciado complicado y
difíciles de manejar se convierten en ecuaciones para algunas de las cuales, además, se posee
un procedimiento algorítmico de resolución. Sin embargo el Álgebra también tiene sus
peligros si no se usa adecuadamente. Uno de los más importantes es que es muy fácil obviar
el significado de las manipulaciones que uno está efectuando. Se olvida a veces que los
números que aparecen en el problema no son realmente entes abstractos sino que
generalmente simbolizan cantidades de magnitud y al operar con ellos no deberíamos perder
17
Antonio M. Oller Marcén Historias de Matemáticas
de vista lo que significa cada una de las operaciones y resultados parciales (Livneh y
Linchevski, 2007).
Pensamos que el uso de textos antiguos en el aula puede ser una buena forma de mejorar
los procesos de aprendizaje de las Matemáticas (Richard, Meavilla y Fortuny, 2010) y, en
particular, de facilitar el proceso de transición de la Aritmética al Álgebra. En concreto el uso
de textos de épocas en las que el Álgebra no era conocida o no estaba aún completamente
generalizada (Férnandez, 1997). Así, puede ser útil realizar actividades de diversos tipos, por
ejemplo:
2. Comparar una solución algebraica con una aritmética del mismo problema,
comentando las ventajas e inconvenientes de cada una de ellas.
Evidentemente, para poder realizar este tipo de actividades es necesaria por parte del
docente una labor previa de investigación y de análisis de textos que se adecuen al propósito
perseguido (Meavilla, 2004). En este trabajo hemos analizado una serie de problemas que
consideramos interesantes desde un punto de vista aritmético dentro de la colección medieval
de problemas titulada Propositiones ad acuendos juvenes de Alcuino de York. Hemos elegido
una serie de problemas que, bajo una clasificación moderna, podrían clasificarse como “de
mezclas” y que admiten un tratamiento muy rico sin necesidad de recurrir al Álgebra.
Alcuino de York nació en la actual Inglaterra en la primera mitad del siglo VIII. Estudió
en la escuela benedictina de York, por aquel entonces un importante centro cultural, donde
también enseñó y que llegó a dirigir. Hacia el año 782 se trasladó a la Escuela Palatina en
Aquisgrán por invitación de Carlomagno; allí pasó dos periodos: de 782 a 790 y de 793 a 796.
Tras este último periodo se retiró como abad a la abadía de San Martín en Tours, donde
murió.
Figura 1. Rabano Mauro apoyado por Alcuino ofrece una obra a Otgar de Maguncia1
La obra en la que nos vamos a detener se titula en el latín original Propositiones ad acuendos
juvenes o, lo que es lo mismo, Problemas para aguzar a los jóvenes. Se trata de una colección de 53
(o 56, según las fuentes) problemas de diverso carácter: los hay aritméticos, geométricos e
incluso aparecen ejemplos de lo que hoy llamaríamos “Matemática recreativa”. El interés de
la obra radica en que, en muchos casos, se trata de problemas que aparecen por primera vez
en la tradición occidental y que recogen tradiciones orientales (donde los textos recopilatorios
de problemas y soluciones tenían una tradición muy anterior) de origen incierto. El texto
presenta los enunciados de los problemas junto con sucintas soluciones que, en el caso de
1
Alcuino de York Carolingian Manuscript manuscriptum Fuldense ca. 831/40, Österreichische Nationalbibliothek
Wien.
respuestas numéricas, incluyen tan sólo la respuesta junto con una comprobación de su
corrección.
Para estudiar el texto hemos recurrido a dos fuentes que, además, están disponibles en la
red. Son las siguientes:
Contenida en Beati Flacci Albini seu Alcuini, Abbatis et Caroli Magni Imperatoris
Magistri. Opera Omnia: Operum pars octava: Opera Dubia. Ed. D. Frobenius [Froster],
Ratisbona, 1777. Revisada y republicada por J.P. Migne como: Patrologiae Cursus
Completus: Patrologiae Latinae, Tomus 101, Paris, 1863, columnas 1143-1160.
Como hemos dicho, ambos textos (las reediciones de J.P. Migne) se hayan disponibles en
la página web [1]. En el primero de ellos aparecen 53 problemas, en el segundo son 56; los
problemas que aparecen en ambos textos tienen enunciados y títulos idénticos
En concreto son ocho los problemas de este tipo que se encuentran en las Propositiones.
Sus enunciados son los siguientes, se ha seguido la numeración de Hadley y Singmaster
(1992):
“Dijo un mercader: quiero comprar 100 cerdos por cien denarios; sin embargo
pagaría 10 denarios por un macho, 5 por una hembra y uno por dos lechones.
Diga, aquel que sepa, ¿cuántos machos, cuántas hembras y cuántos lechones debe
haber de forma que no sobre ni falte ninguno?”
“Cierto paterfamilias disponía de 100 sirvientes, entre los que pretendió repartir
100 modios de maíz, de manera que los hombres recibieran tres modios, dos las
mujeres y los niños medio modio. Diga pues, aquél que vale, ¿cuántos hombres,
cuántas mujeres y cuántos niños había?”
“Quiso cierto hombre comprar 100 animales de distintas clases por 100 sueldos,
de forma que los caballos se compran por tres sueldos, los bueyes por 1 y 24
ovejas por 1 sueldo. Diga, aquél que vale, ¿cuántos los caballos, cuántos los
bueyes y cuántas fueron las ovejas?”
“Cierto hombre quiso comprar en oriente 100 animales de diversas clases por 100
sueldos, para lo que ordenó a un sirviente que por un camello se recibieran 5
sueldos, 1 por un asno y 20 ovejas se compraran por 1 sueldo. Diga, quién quiera,
¿cuántos camellos, asnos y ovejas se obtuvieron a cambio de los 100 sueldos?”
“Cierto obispo ordenó repartir 12 panes entre el clero. Previó así que cada
sacerdote recibiera dos panes; un diácono medio y un lector la cuarta parte.
Obrando así el número de clérigos y de panes resulta el mismo. Diga, quién
quiera, ¿cuántos sacerdotes, diáconos y lectores debe haber?”
Si analizamos con detenimiento los enunciados podemos ver que en todo ellos hay algo a
repartir (modios de maíz, panes o dinero) y sobre todo que, una vez hecho el reparto todos los
receptores han recibido en promedio lo mismo: una unidad. Si compramos 100 animales por
100 denarios el precio medio es de un denario por animal, si repartimos 12 panes entre 12
personas cada persona recibe un pan por término medio y lo mismo si repartimos 20 modios
de maíz entre 20 personas. Así pues, todos estos problemas se reducen a tener que calcular la
cantidad correspondiente de cada componente de una mezcla conocidos los precios de cada
componente y el precio medio final (que en este caso, para redondear las cosas, es 1).
Que, de hecho, estos problemas pueden atacarse empleando técnicas similares a las
utilizadas en los problemas clásicos de mezclas es algo que ya observó Leonardo de Pisa
(Fibonacci). En su importante obra Liber Abaci podemos encontrar, por ejemplo, el siguiente
enunciado que es esencialmente idéntico a los que acabamos de presentar (Sigler, 2002, p.
249):
“Un hombre compra carne de cerdo a tres denarios la libra, de ternera a dos denarios la
libra y de cabra a medio denario la libra. ¿Cómo ha de hacer para comprar siete libras de
carne por siete denarios en total?”
sueldos o 100 denarios) y se hace variar el precio de cada una de las tres clases de
animales. En este grupo se incluyen los problemas 5, 38 y 39.
También en este caso vamos a definir una familia de problemas R(n) cuyo
enunciado sería: “Se han de repartir n modios de maíz entre n sirvientes de tal
forma que los hombres reciben 3 modios, las mujeres 2 y los niños medio.
¿Cuántos hombres, mujeres y niños ha de haber?”
Además de dar estos valores, como ya hemos dicho, Alcuino efectúa los cálculos
oportunos para comprobar que se trata efectivamente de soluciones del problema propuesto.
Sin embargo no hay ningún tipo de indicación o pista sobre cómo ha obtenido dicha solución.
Puesto que, dada la época en la que fue escrita la obra, es poco probable que Alcuino
estuviera familiarizado con técnicas algebraicas hemos de suponer que – de existir – el
método utilizado hubo de ser puramente aritmético (cuando no de mero ensayo y error). En
las secciones que siguen vamos a tratar de presentar, aunque sea de modo hipotético, algunas
posibles formas de pensar que pudo haber seguido Alcuino a la hora de resolver los
problemas y que justificarán por qué las respuestas dadas son las que son y no otras de entre
las posibles. Los problemas se analizarán por bloques, teniendo en cuenta la clasificación
anterior. Tan sólo el Problema 47 recibirá, dado su carácter único, un tratamiento
individualizado.
6. El problema 47
Recordemos el enunciado del problema: “Cierto obispo ordenó repartir 12 panes entre el
clero. Previó así que cada sacerdote recibiera dos panes; un diácono medio y un lector la
cuarta parte. Obrando así el número de clérigos y de panes resulta el mismo. Diga, quién
quiera, ¿cuántos sacerdotes, diáconos y lectores debe haber?”. La solución de Alcuino: 5
sacerdotes, 1 diácono y 6 lectores.
Evidentemente, este sistema tiene infinitas soluciones que vienen dadas por:
Ahora bien, no hay que olvidar que, por las condiciones del problema, x, y y z han de ser
enteros no negativos por lo que ha de ser un entero no negativo que cumpla ⁄ .
Es decir, . Así pues, las posibles soluciones son las ternas (4,8,0) y (5,1,6). La primera,
fácil de encontrar por tanteo, posiblemente no fuera aceptable por parte del autor ya que
implica que no hay lectores, lo que en cierto modo contraviene las condiciones del problema.
Observamos en primer lugar que lo que queremos hacer es repartir 12 panes entre 12
personas de tal forma que cada uno reciba en promedio un pan. Como cada sacerdote recibe 2
panes, cada vez que “cogemos” un sacerdote la media empeora en una unidad. Por su parte,
cada diácono mejora la media en media unidad (pues recibe medio pan cada uno) y, por
último, al recibir un cuarto de pan, cada lector mejora la media en tres cuartos.
Ahora, para que la media sea de un pan por persona, las desviaciones de la media que
acabamos de señalar deben compensarse; es decir, el número de sacerdotes debe ser igual a la
mitad de los diáconos más las tres cuartas partes de los lectores. Como el número de
sacerdotes es entero positivo se sigue que la mitad de los diáconos más las tres cuartas partes
de los lectores tiene que ser entero positivo. Esto último sólo puede pasar si el doble de los
diáconos más el triple de los lectores es múltiplo de 4.
1. Cada uno de ellos es mayor o igual que 1 y menor o igual que 10.
De esto se deduce que el múltiplo de 4 buscado es mayor que 5 y menor o igual que 32; es
decir, uno de entre la lista: 8, 12, 16, 20, 24, 28, 32. Para cada uno de estos valores tenemos las
siguientes posibilidades (teniendo en cuenta las dos indicaciones anteriores):
Al igual que sucedía más arriba el sistema obtenido tiene infinitas soluciones que vienen
dadas, por ejemplo, por:
{
La forma de presentar estas soluciones (hay arbitrariedad a la hora de elegir el parámetro)
está motivada por el hecho de que, en los casos considerados por Alcuino, n es un múltiplo de
5 y así para que las soluciones dadas sean enteras es suficiente con que lo sea . Además, para
que las soluciones sean no negativas habrá que exigir que . Así, se tienen las
siguientes soluciones para los casos estudiados por Alcuino:
En primer lugar hay que observar que conocida una solución para un cierto problema
R(n) es fácil construir soluciones para otros R(m). Por ejemplo es muy fácil observar que si
(x,y,z) es una solución de R(n), entonces (kx,ky,kz) es solución de R(kn). Esta afirmación no
requiere de demostración puesto que se sigue de forma casi directa de las condiciones del
problema. Pese a todo, Alcuino no parece haber observado este hecho puesto que la solución
que ofrece para R(90) no es el triple de la que ofrece para R(30) ni la de R(100) es cinco veces la
de R(20).
Otra relación que salta a la vista entre los problemas que propone Alcuino es que hay dos
parejas de problemas cuyo parámetro difiere en lo mismo: 10 unidades. En este sentido, si
somos capaces de construir una solución de R(n+k) a partir de una solución de R(n) de un
modo concreto podremos construir del mismo modo una solución de R(m+k) a partir de una
de R(m). La propiedad en cuestión es la siguiente: Si (x,y,z) es solución de R(n) y (x+a,y+b,z+c)
es solución de R(n+k), entonces conociendo una solución (x’,y’,z’) de R(m) se tiene que
(x’+a,y’+b,z’+c) será solución de R(m+k). Además se debe cumplir que
.
Recapitulemos las soluciones dadas por Alcuino a los problemas de reparto de maíz que
plantea en su obra (como antes x es el número de hombres, y el de mujeres y z el de niños):
En todo el relato anterior lo único que no encaja con las soluciones dadas por Alcuino es
la de R(90). Sin embargo estamos convencidos de que es poco probable que a la hora de
resolver los problemas la idea de Alcuino fuera resolver uno tras otro repitiendo algún
método similar al del Problema 47; sobre todo teniendo en cuenta que los cuatro problemas
aparecen sucesivamente en el texto. Más aún, el Problema 33a, que es el que corresponde a
R(90) no aparece en algunas ediciones del texto por lo que no sería descabellado pensar que
en un proceso de copia o trascripción se perdiera la solución o se sustituyera la original por
otra. En cualquier caso, todo esto queda en el terreno de la especulación.
Donde x es el número de animales del primer tipo (que cuestan a sueldos cada uno), y la
cantidad de animales del segundo tipo (cuyo precio unitario es de b sueldos) y z es la cantidad
de animales del tercer tipo (que cuesta 1 sueldo cada grupo de c de ellos). Podemos
transformar el sistema en otro equivalente:
{
( ) ( ) ( )
( ) ( )
{
Discutir el carácter entero y no negativo de las soluciones es en este caso bastante más
complicado que en los ejercicios anteriores – y lo omitiremos – debido al mayor número de
parámetros. Sin embargo es interesante indicar que cuando b=1, como sucede en los
Problemas 38 y 39 (que además son consecutivos en el texto), la discusión es algo más sencilla.
A la hora de intentar indagar el proceso que pudo haber seguido Alcuino para dar sus
soluciones a estos problemas, hemos de tratar por separado el problema C(10,5,2) de los otros.
En los otros dos casos, C(3,1,24) y C(5,1,20) se tiene que el precio de los animales del segundo
tipo coincide con el precio medio que se quiere pagar. Esto hace que el problema admita un
análisis mucho más sencillo.
El problema C(3,1,24) admite un análisis similar, sólo que en este caso cada animal del
primer tipo debe compensarse con del tercero y de ahí se sigue el resultado. De hecho, el
problema C(a,1,c) admite en general la solución ( ) ( ( )) que se
obtiene razonando de forma idéntica a la anterior y que, en general, no será única.
En el caso del problema C(10,5,2) el análisis en más complejo puesto que no podemos
despreocuparnos de ninguno de los tipos de animales. Lo que podemos hacer es estudiar
separadamente la relación entre cada uno de los dos primeros tipos de animales (los que
tienen un precio superior a la media) y los animales del tercer tipo (que tienen un precio
inferior a la media). Haciendo esto se obtiene que por cada macho debe haber 18 lechones y
que por cada hembra ha de haber 8 lechones. Como en total ha de haber 100 animales se sigue
que, a lo sumo, hay 9 hembras (si hubiera 10, habría 80 lechones y si hubiera al menos 1
macho con los 18 lechones que conlleva excederíamos los 100 animales). Así pues,
procediendo por tanteo, la primera opción es comprar 9 hembras, que se han de compensar
con 72 lechones; quedan por comprar 19 animales que, casualmente (si es que existen las
casualidades), se pueden completar con 1 macho y los 18 lechones que conlleva. Así pues
podemos afirmar que la solución dada por Alcuino es la primera que se obtiene por tanteo
después de un sencillo análisis inicial.
9. Conclusiones
Hemos presentado un análisis pormenorizado de los problemas “de mezclas” que
aparecen en las Propositiones ad acuendos juvenes de Alcuino de York. Además de resolverlos
mediante técnicas algebraicas modernas, hemos llevado a cabo una discusión sobre los
procedimientos que pudo haber seguido el autor a la hora de presentar sus soluciones,
aunque se ha de tener en cuenta que las ideas presentadas no pueden ser más que meras
hipótesis dado que Alcuino tan sólo indica las soluciones numéricas a sus problemas sin dar
la más mínima pista sobre cómo llegó hasta ellas.
Referencias
[1] COOPERATORUM VERITATIS SOCIETAS. Documenta Catholica Omnia,
http://www.documentacatholicaomnia.eu
[2] FERNÁNDEZ, F. (1997) Aspectos históricos del paso de la aritmética al Álgebra, Uno: Revista de
didáctica de las matemáticas, Nº 14, pp. 75-98.
[5] GÓMEZ, B. (2006). Los ritos en la enseñanza de la regla de tres. En José Mariano Vallejo, el
matemático ilustrado. Una mirada desde la Educación Matemática, pp. 49-69. Publicaciones
Universidad de Córdoba.
[6] HADLEY, J. Y SINGMASTER, D. (1992). Problems to sharpen the young. An annotated translation
of ‘Propositiones ad acuendos juvenes’ the oldest mathematical problem collection in Latin
attributed to Alcuin of York, Math. Gaz. Nº 76, pp. 102-126.
[9] MEAVILLA, V. (2011). El lobo, la cabra y la col: Antología de problemas matemáticos recreativos.
Editorial Almuzara, Córdoba.
[10] RICHARD, P.R; MEAVILLA, V. Y FORTUNY, J.M. (2010). Textos clásicos y geometría dinámica:
estudio de un aporte mutuo para el aprendizaje de la geometría. Enseñanza de las Ciencias, Nº
28(1), pp. 95-112.
[11] SIGLER, L.E. (2002). Fibonacci’s Liber Abaci. A translation into modern English of Leonardo
Pisano’s Book of Calculation. Springer Verlag, New York.
Sobre el autor:
Nombre: Antonio M. Oller Marcén
Correo Electrónico: [email protected]
Institución: Centro Universitario de la Defensa de Zaragoza, España.
1 de abril de 2014
Resumen
Hoy en día, el carácter aritmético de los valores de la función zeta de Riemann en argu-
mentos enteros y en particular en impares, continúa siendo un problema abierto dentro de
la comunidad matemática. Este artículo, se dedica a presentar los principales resultados al-
canzados por varios matemáticos desde el siglo XVII hasta la actualidad, correspondientes
al carácter aritmético que siguen los valores de la función zeta de Riemann en argumentos
enteros.
Abstract
At present, arithmetical character of the values of the Riemann’s zeta function in integer
arguments and in particular in odd arguments, keep being an open problem into the mathe-
matical community. This paper, dedicates to show the principal results obtained for several
mathematicians from the century XVII until the present time.
1. Introducción
El problema de calcular la suma de la serie
1 1 1
ζ ( s) = ∑ s
= 1+ s + s +··· , (1)
n ≥1
n 2 3
para s ≥ 2 entero, había atraído la atención de varios matemáticos desde el siglo XVII, en
particular para s = 2. Por lo que a principios del siglo XVIII las series infinitas pasan a ser
33
Anier Soria Lorente Historias de Matemáticas
uno de los temas estrellas dentro del universo matemático de la época. El estudio de que si
convergían hacia un número o si se hacían cada vez más grandes sería uno de los retos de
cualquier matemático de aquel entonces. En el siglo XVIII se interesaron en este problema varios
matemáticos como Jacob Bernoulli, Daniel Bernoulli y Christian Goldbach, quienes obtuvieron
algunos resultados preliminares sobre la suma de esta serie en el caso s = 2, los cuales pronto
serían superados por Leonhard Euler que, en este marco conceptual, hizo su primer contacto
con dicha serie en el caso s = 2 y pronto mejoraría los cálculos de sus predecesores. El problema
de calcular la suma de la serie en el caso s = 2 no era fácil debido a su lenta convergencia, por
ejemplo, para calcular el número al que converge, con una precisión de seis decimales, hay que
sumar al menos un millón de términos de la serie. En efecto, como
1 1 1 1 1 1 1
− = < 2 < = − ,
k k+1 k ( k + 1) k k ( k − 1) k−1 k
de tal forma que aproximar la serie con n lugares decimales requiere calcular la suma de alme-
nos 10n términos.
Se conoce que, Leibniz a propuesta de su mentor Huygens, calculó la suma de los recíprocos
de los números triangulares. Descubrir que
1 1 1
∑ n ( n + 1) = ∑ − = 1,
n ≥1 n ≥1
n n+1
satisfizo tanto al joven Leibniz que impulsó su afición por las matemáticas, que luego le llevaría
a co-descubrir el Análisis Matemático. Si la suma de los inversos de los números triangulares
constituyó un problema fácil para Leibniz, no ocurrió lo mismo con la suma de los inversos de
los números cuadrados ζ (2). El problema le fue planteado a Leibniz por Oldenburg, secretario
de la Royal Society en 1673, aunque ya había sido abordado veinte años antes por Pietro Mengoli
y por el mismo Walis (que dió el valor de 1, 645 como aproximación de la suma de la serie).
Leibniz comunicó a sus corresponsales Jacob y Johann Bernoulli el problema, y les dijo que en
apariencia debía tener una solución tan simple como la de los números triangulares inversos. Y
así lo pensaron Jacob y Johann Bernoulli, pero pronto se dieron cuenta de que algo no marchaba
bien. No fue difícil demostrar por comparación que su suma estaba acotada superiormente por
la suma de la serie de los inversos de los números triangulares
1 2
2
< , n > 1.
n n ( n + 1)
Pero el resultado preciso de la suma se les negaba, hasta tal punto que lanzaron públicamente
este grito de socorro, “Grande será nuestra gratitud si alguien encuentra y nos comunica lo que
hasta ahora ha escapado a nuestros esfuerzos”. Desde entonces al problema se le conoce como
Problema de Basilea y fue Johann Bernoulli, hermano menor de Jacob y en ese entonces mentor
de Euler, quien seguramente le sugirió a este último investigar la evaluación de esa suma.
En 1729 Euler recibió una carta de su amigo Christian Goldbach donde le señala un método
de aproximación que lo lleva a estimar el valor entre 1, 64 y 1, 66. Goldbach reta a Euler para que
lo mejore. En el momento de recibir este desafío, Euler, que contaba sólo con 22 años de edad,
se encontraba en la Academia de San Petersburgo enfrascado en varios problemas concretos de
mecánica. Sin embargo, no se olvidó del reto y dos años más tarde, hizo pública una asombro-
sa aproximación de seis cifras decimales exactas, 1, 643934, transformando habilidosamente la
Varios años más tarde, un día que leía con interés una de las obras de Newton, la genialidad de
Euler se desbordó al encontrar la idea de que el desarrollo en serie de la función seno estaba re-
lacionado con la solución exacta del problema. Lo ingenioso sería utilizar el desarrollo del seno
no solo como sumas sino también como producto de infinitos factores. Newton, precisamente
en esta obra, utilizaba con mucha eficacia la relación entre los coeficientes de las potencias y las
raíces de los polinomios. Esto mismo intentó Euler con la serie de los inversos de los cuadrados.
Basándose en el desarrollo en series de potencias del seno
z2n+1
sen z = ∑ (−1)n .
n ≥0 (2n + 1)!
Luego, utilizando el hecho de que los ceros de la función P (z) se producen para los valores en
que el numerador se anula (con excepción de z = 0, donde P (0) = 1, es decir, para todo z = πn,
donde n = ±1, ±2, ±3, . . . Entonces factoriza como si P (z) fuese un polinomio
!
z2 1 1 z2
P ( z ) = ∏ 1 − 2 2 = 1 − 2 ∑ 2 z2 + · · · = 1 − +···
n ≥1
π n π n ≥1
n 3!
De esta manera, comparando los términos de segundo grado, le lleva a Euler a obtener el ma-
ravilloso resultado ζ (2) = π 2 /6. Así le comunicaba Euler su extraordinario hallazgo a Daniel
Bernoulli, el hijo de Johann, en una carta fechada en 1735, dando solución así al llamado Proble-
ma de Basilea [26], que sin lugar a dudas le abrió las puertas para ingresar a la élite matemática
de su época. La deducción de Euler es una joya de las Matemáticas, la misma expresa muy
bien el estilo de esa época prodigiosa. Cuando Euler hizo este cálculo hacía veintinueve años
que Jacob había fallecido. Johann Bernoulli, reconciliado ya con la figura de su hermano mayor,
comentó a Euler: “De este modo el deseo más ferviente de mi hermano se ha cumplido . . . ¡Si
estuviera aquí!”.
Entre 1740 y 1744, utilizando las mismas herramientas Euler encontró la suma de las series
de los inversos de las potencias pares de los números naturales hasta el orden 26
1315862
ζ (26) = π 26 .
11094481976030578125
Todos estos triunfos, estimularon formidablemente a Euler para continuar extendiendo estos
resultados. El espíritu inquieto y perspicaz de Euler no podía sentirse satisfecho con lo encon-
trado hasta el momento. Además, aún faltaban las sumas en el caso de n impar. De modo que
en 1750 publica otro artículo donde señala los valores aproximados de las series armónicas de
orden impar n = 2k + 1, para k = 1, . . . , 7. Y en su famoso Tratado de Cálculo Diferencial de
1755, al fin consigue exponer la elegante fórmula [14]
muy importante en las matemáticas, pues generalizó un problema que se había resistido mucho
tiempo, haciéndose muy famoso por el camino. Pero Euler no se quedó aquí, entusiasmado
expuso la siguiente conjetura sobre el caso impar
1 p
¿ ∑ 2k +1
= π 2k+1 ?
n ≥1 n q
donde p y q son números enteros. Los esfuerzos de Euler para probar la validez de esta conjetura
fueron vanos. No obstante, puede servirle de consuelo que aún hoy, a más 255 años después,
no se ha conseguido ni validar ni refutar su conjetura, la misma sigue siendo todo un misterio.
Actualmente la serie (1) es conocida bajo el nombre de función zeta de Riemann [17, 18, 25,
46], quien la extendió al campo complejo, mostrando, y previendo, muchas de sus interesantes
propiedades. Desafortunadamente Euler no obtuvo nada acerca de los casos impares, todos sus
intentos por evaluar la función zeta de Riemann en argumentos impares ζ (2k + 1), k ∈ N,
fueron fallidos. Hasta la fecha, el carácter aritmético de estos números, exceptuando ζ (3); es
decir, si son racionales o irracionales, sigue siendo un problema abierto dentro de la comunidad
matemática.
Después de los estudios iniciados por Euler, nada se supo sobre la naturaleza aritmética de
ζ (2k + 1) para k ∈ N, hasta que Roger Apéry, en el año 1979, sorprendió a la comunidad ma-
temática con una demostración de la irracionalidad de ζ (3) [2]. De ahí el conocido teorema de
Apéry, ζ (3) ∈/ Q. Él sólo brindó una breve descripción de su impresionante demostración, cu-
yos resultados se pueden encontrar de forma detallada en [2, 15, 47, 65]. Como reconocimiento
a este resultado, la constante ζ (3), se denomina actualmente, constante de Apéry.
El método de obtención de los aproximantes, aunque ingenioso, no hacía, sin embargo, uso
alguno de resultados que no hubieran sido conocidos por los matemáticos del siglo XVIII. ¡Una
demostración que se le había escapado al gran Euler! Una excelente exposición puede encon-
trarse en Vander Poorten [65], quien dio una conferencia sobre la demostración de Apéry en el
congreso internacional celebrado ese mismo año en Helsinki.
A partir de la demostración de la irracionalidad de ζ (3) dada por Apéry, se organizaron
múltiples seminarios en los que se pretendía entender dicha demostración, con la finalidad
de responder a las interrogantes acerca de las propiedades aritméticas de la función zeta de
Riemann en enteros impares, lo cual tuvo lugar en el Institute for Advanced Study, dirigido por
E. Bombieri. No obstante, hasta la fecha no se sabe siquiera si ζ (5) es irracional o no; algunos de
los pocos resultados relacionados con dicho número se pueden encontrar en [51, 53, 61, 67, 68,
70]. Sin embargo, el resultado de Apéry inspiró a varios autores [8, 10, 12, 36, 43, 44, 58, 59, 71] a
construir diferentes métodos para explicar la irracionalidad de ζ (3). Sorprendentemente, estos
métodos conducen a la misma sucesión de aproximantes racionales de la clase de Apéry. En [8]
se muestra que todos estos métodos coinciden, teniendo como origen común, un problema de
aproximación simultánea. A este hecho se le denomina “fenómeno de Apéry”.
Uno de ellos, se basa en una integral doble que involucra los polinomios de Legendre, la
cual fue considerada por Beukers [8, 10, 11, 37]
Z 1Z 1
log xy
rn = q n ζ (3 ) − p n = − Ln ( x ) Ln (y) dx dy (2)
0 0 1 − xy
( xyz (1 − x ) (1 − y) (1 − z))n
Z 1Z 1Z 1
= dx dy dz,
0 0 0 (1 − (1 − xy) z)n+1
36 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Aritmética de los valores de la función zeta de Riemann en argumentos enteros Anier Soria Lorente
donde
n 2 2 n 2 2
n+k n n+k n
qn = ∑ , y pn = ∑ γn,k , (3)
k =0
k k k =0
k k
n
1 k
(−1) j−1 n + j −1 n −1
∑ ,
j3 j∑
γn,k = +
j =1 =1 2j3 j j
denotan los polinomios de Legendre, ortogonales con respecto a la medida de Lebesgue en
√ 4n
(0, 1). Por otra parte, Beukers demostró de una manera sencilla que rn = O 2−1 , lo
cual le permitió dar una nueva demostración de la irracionalidad de ζ (3). Es de vital impor-
tancia destacar, que para ciertas modificaciones de la integral de Beukers (2), se ha mejorado la
medida de irracionalidad [30, 48] dada por Apéry [2, 65, 66].
De una manera similar a [23], Beukers [12] consideró un problema de aproximación racio-
nal en un intento por formular de un modo más natural la demostración de Apéry. Para ello
introdujo la función racional
(n − z + 1)2n
Rn ( z ) = , (4)
(−z)2n+1
donde (t)n = t (t + 1) · · · (t + n − 1) denota el símbolo de Pochhammer [24], a partir del cual
mediante el cálculo de su desarrollo en fracciones simples, dedujo de forma directa los aproxi-
mantes racionales de Apéry (3).
Sorokin en [58, 63], obtuvo los aproximantes racionales de Apéry (3) del mismo modo que
Beukers, considerando el problema de aproximación racional
An (z) f 1 (z) + Bn (z) f 2 (z) − Cn (z) = O z−n−1 ,
An (z) f 2 (z) + 2Bn (z) f 3 (z) − Dn (z) = O z−n−1 ,
log2 x
Z 1 Z 1 Z 1
dx log x 1
f 1 (z ) = , f 2 (z) = − dx, f 3 (z ) = dx.
0 z−x 0 z−x 2 0 z−x
De este modo, demostró que la solución de este problema viene dada por las relaciones de
ortogonalidad
Z 1
( An (x ) − Bn (x ) log x ) x k dx = 0,
0
k = 0, . . . , n − 1,
Z 1
k
(( An ( x ) − Bn ( x ) log x ) log x ) x dx = 0,
0
junto con la condición An (1) = 0. Luego, usando la convolución de Mellin [62, 63, 64] como un
ingrediente crucial, consiguió lo siguiente
Z 1
( An ( x ) − Bn ( x ) log x ) log x
rn = dx
0 1−x
Z 1Z 1
log xy
= − Ln ( x ) Ln (y) dx dy,
0 0 1 − xy
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 37
Anier Soria Lorente Historias de Matemáticas
lo cual implica la irracionalidad de ζ (3) de acuerdo con la estimación de Beukers dada en (2),
véase [10].
Otra de las principales aportaciones que han sido desarrolladas para explicar la irracionali-
dad de ζ (3), fue propuesta por Nesterenko, quien en 1996 inspirado en la obra de Gutnik [23]
consideró la siguiente modificación de la función racional (4) de Beukers
(1 − z)2n
Rn ( z ) = , (5)
(z)2n+1
y de esta manera demostró que la sucesión residuo dada en (2) se podía escribir de la siguiente
forma
2
(1 − ν)2n π 2
!
∂ (1 − k ) n 1
Z
rn = − ∑ = dν,
k ≥1
∂k ( k )2 2πi L (ν)2
n +1
sin πν
n +1
donde L es la línea vertical Re z = C, 0 < C < n + 1, orientada de arriba hacia abajo. De hecho,
él aplicó a esta integral, el conocido método de Laplace [34, 35], lo cual le permitió llegar al
mismo comportamiento de la sucesión residuo (2) encontrado por Beukers.
Más tarde, en el 2002, Zudillin basado en los resultados de Nesterenko, utilizó la función
racional (5) y haciendo uso del algoritmo de Zeilberger [3, 4, 5, 6, 41, 42], obtuvo la relación de
recurrencia de Apéry
(n + 1)3 yn+1 − (2n + 1)(17n2 + 17n + 5)yn + n3 yn−1 = 0, n ≥ 1, (6)
la cual satisfacen la sucesión de los numeradores p n y denominadores qn de los aproximantes
racionales (3) a ζ (3) con condiciones iniciales
p0 = 0, p1 = 6, q0 = 1, q1 = 5,
donde a partir de la misma se prueba evidentemente la irracionalidad de ζ (3). Además, a partir
de la relación de recurrencia anterior y de los aproximantes racionales (3) a ζ (3), Apéry [65]
dedujo el siguiente desarrollo en fracciones continuas
6| 1 | 64 | n6 |
ζ (3 ) = − − −···− 2
−···
| 5 | 117 | 535 | (2n + 1) (17n + 17n + 5)
También, de forma similar, en 1996 Nesterenko [36] propuso el siguiente desarrollo en fracciones
continuas
1| 2| 1| 4| 2| 6| 4|
2ζ (3) = 2 + + + + + + + +··· ,
|2 |4 |3 |2 |4 |6 |5
donde los numeradores an , n ≥ 1, y los denominadores bn , n ≥ 2, están definidos mediante
b4k+1 = 2k + 2, a4k+1 = k (k + 1) , b4k+2 = 2k + 4,
a4k+2 = (k + 1) (k + 2) , b4k+3 = 2k + 3,
a4k+3 = ( k + 1 )2 , b4k = 2k, a4k = (k + 1)2 .
Luego, en el año 2009, Nesterenko publicó una nueva demostración de la irracionalidad de ζ (3)
[38], demostrando para ln3 Bn , ln3 Dn ∈ Z, lo siguiente
∂
(−1)n ln3 ∑ R (k) = (−1)n−1 ln3 (2Bn ζ (3) − Dn ) < (4/5)n ,
k ≥1
∂k
donde
[(n−1)/2] [n/2]
k−j k−j
R (k) = k −2 ∏ ∏ ,
j =1
k+j j =1
k+j
[(n−1)/2] [n/2]
(3) (2)
Bn = ∑ bk y Dn = ∑ 2bk Hk + ak Hk ,
k =0 k =1
con
h i h i
n −1 n n −1 n
+k +k
n −1 2 2 2 2
bk = (−1) ,
k k k k
[(n−1)/2] [n/2] !
2 1 1 1 1
ak = bk + ∑ − + ∑ − ,
k j =0
k−j k+j j =0
k−j k+j
lo cual prueba evidentemente el teorema de Apéry. Aquí, ln denota el mínimo común múltiplo
(r ) (1)
de {1, 2, . . . , n} y Hk el k-ésimo número armónico de orden r (Hk = Hk y H0 = 0) [16].
Se conoce además, que de los pocos resultados desarrollados para ζ (4), se encuentran los
aportados por Zudilin en [74, 75], donde basado en la serie de tipo hipergeométrica
2 2
(−1)n+1
!
∂ (1 − k)n (k + n + 1)n (2k + n)
rn,z = ∑ ∂k = qn,z ζ (4) − pn,z (7)
6 k ≥1 (k)4n+1
√ 3n
= O 3−2 3 ,
− 3n3 9n2 − 1 yn−1 = 0, n ≥ 1, (8)
donde los aproximantes racionales involucrados en (7), la satisfacen con condiciones iniciales
q0,z = 1, q1,z = 12, p0,z = 0, p1,z = 13,
y a partir de la misma obtuvo el desarrollo en fracciones continuas,
13 | 17 · 2 · 3 · 4 | 27 · 5 · 6 · 7 |
ζ (4) = + + +···
| P (0 ) | P (1) | P (2)
n7 (3n − 1) (3n) (3n + 1) |
+ +··· ,
| P ( n)
siendo P (n) = 3 (2n + 1) 3n2 + 3n + 1 15n2 + 15n + 4 . Cabe destacar, que la relación de
recurrencia de segundo orden (8), no brinda aproximantes diofánticos que prueban la irracio-
nalidad de ζ (4), sin embargo, presentan un eficiente y rápido algoritmo para el cálculo de esta
constante. También, existen otros resultados, relacionados con el número ζ (4), los cuales han
sido desarrollados en su gran mayoría por Zudilin (véase [14, 27, 60, 72, 74, 75, 77]).
En vista a extender el resultado anterior para ζ (5), Zudilin en [70] auxiliándose de las series
hipergeométricas
n (1 − k ) n ( k + n + 1)n
Fn = n!4 (−1)n ∑ k+ ,
k ≥1
2 (k)6n+1
n (−k)n+1 (k + n)n+1
n!4 (−1)n+2 ∑ k +
F̃n = ,
k ≥1
2 (k)6n+1
dedujo la siguiente relación de recurrencia de tercer orden
donde
la cual satisfacen la sucesión de los numeradores pn,5 y denominadores qn,5 de los aproximantes
racionales a ζ (5) con condiciones iniciales
87 1190161
p0,5 = 0, p1,5 = , p2,5 = − ,
2 64
q0,5 = −1, q1,5 = 42, q2,5 = −17934.
Además, comprobó que la sucesión rn,5 = qn,5 ζ (5) − pn,5 > 0 también satisface la relación de
recurrencia (9) y verificó que la misma y la sucesión de los denominadores qn,5 , satisfacen los
siguientes límites
log |rn,5 |
lı́m = log |µ2 | = −1.08607936 . . . ,
n→∞ n
log |qn,5 |
lı́m = log |µ3 | ,
n→∞ n
donde
son las raíces del polinomio característico µ3 + 2368µ2 − 752µ − 16 de la relación de recurrencia
(9). Con estos resultados, Zudilin presentó un eficiente y rápido algoritmo para el cálculo de
esta constante ζ (5), puesto que la sucesión de números racionales pn,5 /qn,5 converge a ζ (5) con
una velocidad |µ2 /µ3 | < 1.42521964 ·10−4.
En el 2001, Ball y Rivoal [9, 49, 50, 54] probaron que la sucesión {ζ (2k + 1)} k≥1 contiene una
infinidad de irracionales. En otra dirección, para el 2002, el ruso Vladimir Zudilin, un investi-
gador de la Universidad Lomonosov de Moscú, probó que al menos uno de los números ζ (5),
ζ (7), ζ (9) y ζ (11) es irracional [68]. Lo cierto es, que son escasos los resultados que existen
acerca de las propiedades aritméticas de los números ζ (2k + 1), con k ∈ N, los cuales pueden
consultarse en [19, 20, 21, 22, 33, 45, 52, 55, 56, 57, 69, 73, 76]
Referencias
[1] A BRAMOWITZ, M. and S TEGUN, I. A., Handbook of mathematical functions with formulas,
graphs, and mathematical tables, Dover, New York, 1972.
[2] A PÉRY, R., Irrationalité de ζ (2) et ζ (3), Astérisque, Vol. 61, pp. 11–13, 1979.
[4] A BRAMOV, S.A. and L E, H. Q., A criterion for the applicability of Zeilberger’s algorithm to rational
functions, Discrete Math., Vol. 259, pp. 1–17, 2002.
[5] A BRAMOV, S.A., When does Zeilberger’s algorithm succeed?, Appl. Math., Vol. 30, pp. 424–441.
[6] A BRAMOV, S.A. and L E, H. Q., Telescoping in the context of symbolic summation in Maple, J. of
Symb. Comput., Vol. 38, pp. 1303–1326, 2004.
[7] A LMKVIST, G., VAN S TRATEN, D. and Z UDILIN, W., Apéry limits of differential equations of
order 4 and 5, Modular Forms and String Duality, Fields Inst. Commun. Ser., Amer. Math.
Soc. & Fields Inst., Providence, Vol. 54, pp. 105–123, 2008.
[8] A RVESÚ, J., Orthogonal forms: A key tool for deducing Apéry’s recurrence relation, J. Approx.
Theory, 2012.
[9] B ALL, K. and R IVOAL, T., Irrationalité d’une infinité de valeurs de la fonction zeta aux entiers
impairs, Inv. Math., Vol. 146, pp. 193–207, 2006.
[10] B EUKER, F., A note on the irrationality of ζ (2) and ζ (3), Bull. London Math. Soc., Vol. 11, pp.
268–272, 1979.
[11] B EUKER, F., Legendre polynomials in irrationality proofs, Bull. Austral. Math. Soc., Vol. 22, pp.
431–438, 1980.
[12] B EUKER, F., Padé approximations in number theory, Padé approximation and its applications,
(Amsterdam, 1980), Lecture Notes in Math., Vol. 888, pp. 90–99, Springer, Berlin-New York,
1981.
[13] B EUKER, F., Consequences of Apéry’s work on ζ (3), preprint of talk presented at the Rencon-
tres Arithmétiques de Caen, ζ (3) Irrationnel: Les Retombées, 1995.
[14] B ALANZARIO, E. P., Método elemental para la evaluación de la función zeta de Riemann en los
enteros pares, Miscelánea Matemática, Vol. 33, pp. 31–41, 2001.
[15] C OHEN, H., Démonstration de I’irrationalité de ζ (3) (d’aprés R. Apéry), Séminaire de théorie
des nombres, Grenoble, 1978.
[16] C OHEN, J. and G UY, R., The Book of Numbers, Copernicus, Springer-Verlag New York, Inc.,
1997.
[18] E DWARDS, H. M., Riemann’s Zeta Function, Academic Press, Dover, New York, 1974.
[19] F ISCHLER, S., Irrationalité de valeurs de zêta [d’après Apéry, Rivoal, ...], Séminaire Bourbaki,
Astérisque, Vol. 294, Nº 990, pp. 27–62, 2004.
[20] F ISCHLER, S., Restricted rational approximation and Apéry-type constructions, Indagationes
Mathem, Vol. 20, pp. 201–215, 2009.
[21] F ISCHLER, S., W. Zudilin, A refinement of Nesterenko’s linear independence criterion with appli-
cations to zeta values, Math. Ann., Vol. 347, pp. 739–763, 2010.
[22] F ISCHLER, S., Nesterenko’s criterion when the small linear forms oscillate, Arch. Math., Vol. 98,
pp. 143–151, 2012.
[23] G UTNIK, L. A., On the irrationality of certain quantities involving ζ (3), Acta Arith., Vol. 42,
pp. 255–264, 1983.
[24] G ASPER, L. and R AHMAN, M., Basic Hypergeometric Series, Encyclopedia of Mathematics
and its Applications, Cambridge University Press, Cambridge, 2004.
[25] G OURDON, X. and S EBAH, P., The Riemann Zeta-function ζ (s): generalities, Numbers, cons-
tants and computation, 2004.
[26] G RANERO, F., El Problema de Basilea: historia y algunas demostraciones, La Gaceta de la RSME,
Vol. 12, Nº 4, pp. 721–737, 2009.
[27] H ANCL, J., A simple proof of the irrationality of ß4 , Amer. Math. Monthly, Vol. 93, pp. 374–375,
1986.
[28] H ATA, M., Legendre type polynomials and irrationality measures, J. Reine Angew. Math., pp.
99–125, 1990.
[29] H ATA, M., On the linear independence of the values of polylogarithmic functions, J. Math. Pures
et Appl., Vol. 69, pp. 133–173, 1990.
[30] H ATA, M., A new irrationality measure for ζ (3), Acta Arith., Vol. 92, pp. 47–57, 2000.
[31] J IN, Y. and D ICKINSON, H., Apéry sequences and Legendre transforms, J. Austral. Math. Soc.,
Vol. 68, pp. 349–356, 2000.
[32] J AIN, L. and T ZERMIAS, P., Beukers’ integrals and Apéry’s recurrences, J. Integ. Seq., Article
05.1.1, Vol. 8, 2000.
[33] K RATTENTHALER, C. and R IVOAL, T., Hypergéométrie et fonction zêta de Riemann, Mem.
Amer. Math. Soc., Vol. 186, Nº 875, 2007.
[34] L ÓPEZ, J. L. and PAGOLA, P. J., The Laplace’s and Steepest Descents Methods Revisited, Inter-
national Mathematical Forum, Vol. 2, Nº 7, pp. 297–314, 2007.
[35] L ÓPEZ, J. L. and PAGOLA, P. J., Fórmulas explícitas para los coeficientes de los métodos de Laplace
y Saddle Point, XI Congreso de Matemática Aplicada, Ciudad Real, pp. 1–9, septiembre 2009.
[36] N ESTERENKO, Yu. V., A few remarks on ζ (3), Math. Notes, Vol. 59, Nº 6, pp. 625–636, 1996.
[37] N ESTERENKO, Yu. V., Integral identities and constructions of approximations to zeta values, J.
Théor. Nombres Bordx, Vol. 15, pp. 535–550, 2003.
[38] N ESTERENKO, Yu. V., An Elementary Proof of the Irrationality of ζ (3), Moscow University
Mathematics Bulletin, Vol. 64, Nº 4, pp. 165–171, 2009.
[39] N IKISHIN, E. M. and S OROKIN, V. N., Rational Approximations and Orthogonality, Transla-
tions of Mathematical Monographs, Amer. Math. Soc. Providence, RI, Vol. 92, 1991.
[40] O SLER, T. J., A Computer Hunt for Apery’s, Math. Spectrum, Vol. 35, Nº 1, pp. 5–8, 2002.
[41] PAULE, P., A Mathematica Version of Zeilberger’s Algorithm for Proving Binomial Coefficient
Identities, J. Symb. Comput., Vol. 20, pp. 673–698, 1995.
[42] PAULE, P., Symbolic summation–some recent developments, In Computer Algebra in Science
and Engeneering–Algorithms, Systems, and Applications, pp. 138–162, 1995.
[43] P RÉVOST, M., A new proof of the irrationality of ζ (2) and ζ (3) using Padé approximants, J.
Comput. Appl. Math., Vol. 67, pp. 219–235, 1996.
[44] P RÉVOST, M. and S OROKIN, V. N., Hermite-Padé Approximants and Apery Theorem, Technical
Report L.M.P.A., 2007.
[45] P ILEHROOD, Kh. H., Simultaneous generation for zeta values by the Markov-WZ method, D.
Math. Theor. Comput. Sci. DMTCS, pp. 115–124, 2008.
[46] R IEMANN, B., Über die Anzahl der Primzahlen unter einer gegebenen Grösse, Monastsber. Akad.
Berlin, pp. 671–680, 1859.
[47] R EYSSAT, E., Irrationalité de ζ (3), Selon Apéry in: 20 année. Seminaire Delange–Pisot–
Poitou, 1978.
[48] R HIN, G. and V IOLA, C., The Group Structure for ζ (3), Acta Arith., Vol. 97, pp. 269–293, 2001.
[49] R IVOAL, T., La fonction zêta de Riemann prend una infinité de valeurs irrationnelles aux entiers
impairs, C.R. Acad. Sci. Paris Sér. I Math., pp. 267–270, 2000.
[50] R IVOAL, T., Propriétés diophantiennes des valeurs de la fonction zêta de Riemann aux entiers
impairs, Thèse de doctorat, Laboratoire SDAD, Université de Caen, 2001.
[51] R IVOAL, T., Irrationalité d’au moins un des neuf nombres ζ (5) , ζ (7) , . . . , ζ (21), Acta Arith.,
pp. 157–167, 2002.
[52] R IVOAL, T. and Z UDILIN, W., Diophantine properties of numbers related to Catalan’s constant,
Math. Ann., Vol. 326, pp. 705–721, 2003.
[53] R IVOAL, T., Séries hypergéométriques et irrationalité des valeurs de la fonction zêta, Actes des
journées arithmétiques de Lille (July, 2001), J. Théorie Nombres Bordeaux, pp. 351–365, 2003.
[54] R IVOAL, T., Valeurs aux entiers de la fonction zêta de Riemann, Quadrature, Vol. 49, 2003.
[55] R IVOAL, T., Simultaneous generation of Koecher and Almkvist-Granville’s Apéry-like formulae,
Experiment. Math., Vol. 13, pp. 503–508, 2004.
[56] R IVOAL, T., Hypergeometric constructions of rational approximations of (multiple) zeta values,
in Algebraic and Analytic Aspects of Zeta Functions and L-Functions, G. Bhowmik et. al.
(eds.), MSJ Memoirs, Mathematical Society of Japan, Tokyo, Vol. 21, pp. 167–183, 2010.
[57] S ÁNCHEZ, P. A., Valores irracionales de la función zeta de Riemann en los impares, Morfismos,
Vol. 10, Nº 2, pp. 15–43, 2006.
[58] S OROKIN, V. N., Hermite-Padé approximations for Nikishin systems and the irrationality of ζ (3),
Communications of the Moscow Math. Soc., Vol. 49, pp. 176–177, 1993.
[59] S OROKIN, V. N., Apéry’s theorem, Vestnik Moskov. Univ. Ser. I Mat. Mekh. [Moscow Univ.
Math. Bull.], Nº 3, pp. 48–52, 1998.
[60] S OROKIN, V. N., One algorithm for fast calculation of ß4 , Preprint (Russian Academy of Scien-
ces, M.V. Keldysh Institute for Applied Mathematics, Moscow), 2002.
[61] S ALIKHOV, V., On multiple integrals represented as a linear form in 1, ζ (3), ζ (5), . . . , ζ (2k − 1),
J. Math. Sci., Vol. 146, Nº 2, 2007.
[62] S MET, Ch. and VAN A SSCHE, W., Mellin transforms for multiple Jacobi-Piñeiro polynomials and
a q-analogue, J. Approx. Theory, pp. 782–806, 2010.
[63] VAN A SSCHE, W., Multiple orthogonal polynomials, irrationality and transcendence, Contemp.
Math., Vol 236, pp. 325–342, 1999.
[64] VAN A SSCHE, W., Hermite-Padé Rational Approximation to Irrational Numbers, Comput. Meth.
Funct. Theory, Vol. 10, Nº 2, pp. 585–602, 2010.
[65] VAN DER P OORTEN, A. J., A proof that Euler missed ... Apéry’s proof of the irrationality of ζ (3),
An informal report, Math. Intelligencer, Vol. 1, pp. 195–203, 1979.
[66] Z UDILIN, W., Difference equations and the irrationality measure of numbers, Proceedings of the
Steklov Institute of Mathematics, Vol. 218, pp. 160–174, 1997.
[67] Z UDILIN, W., One of the eight numbers ζ (5) , ζ (7) , . . . , ζ (17) , ζ (19) is irrational, Mat. Za-
metki [Math. Notes], pp. 472–476, 2001.
[68] Z UDILIN, W., One of the eight numbers ζ (5) , ζ (7) , ζ (9) , ζ (11) is irrational, Uspekhi Mat.
Nauk [Russian Math. Surveys], pp. 149–150, 2001.
[69] Z UDILIN, W., Irrationality of values of zeta function at odd integers, Uspekhi Mat. Nauk [Rus-
sian Math. Surveys], pp. 215–216, 2001.
[70] Z UDILIN, W., A third-order Apéry-like recursion for ζ (5), Mat. Zametki [Math. Notes], pp.
796–800, 2002.
[71] Z UDILIN, W., An elementary proof of Apéry’s theorem, E-print math. NT/0202159, Moscow
Lomonosov State University, pp. 1–8, 2002.
[72] Z UDILIN, W., Difference equation and permutation group for ζ (4), Actes des 12‘emes rencon-
tres arithmétiques de Caen (June 29–30, 2001), J. Théorie Nombres Bordeaux (2002).
[73] Z UDILIN, W., Irrationality of values of the Riemann zeta function, Izv. Ross. Akad. Nauk Ser.
Mat. [Russian Acad. Sci. Izv. Math.], Vol. 66, Nº 3, 2002.
[74] Z UDILIN, W., An Apéry-Like Difference Equation for Catalan’s Constant, Elect. J. Comb., Vol.
10, Nº 1, pp. 1–10, 2003.
[75] Z UDILIN, W., Well-poised hypergeometric service for diophantine problems of zeta values, Actes
des 12èmes rencontres arithmétiques de Caen (June 29–30, 2001), J. Théorie Nombres Bor-
deaux, pp. 593–626, 2003.
[76] Z UDILIN, W., Arithmetic of linear forms involving odd zeta values, J. Théorie Nombres Bor-
deaux, pp. 251–291, 2004.
[77] Z UDILIN, W., Binomial sums related to rational approximations to ζ (4), Math. Notes, Vol. 75,
pp. 594–597, 2004.
[78] Z UDILIN, W., Arithmetic hypergeometric series, Russian Math. Surveys, pp. 163–216, 2011.
Sobre el autor:
Nombre: Anier Soria Lorente
Correo electrónico: [email protected]
Institución: Universidad de Granma, Bayamo, Cuba.
1 de abril de 2014
Resumen
El siguiente relato que trata de poner de manifiesto cómo tras la “ignorancia” matemática
puede haber atrevimiento por el que se puede llegar a cualquier resultado por extraño que
éste resulte, todo ello llevado en un tono narrativo ameno y divertido.
Abstract
The following story tries to show that behind the “ignorance” you may find daring and
reach whatever result no matter how strange this could result, everything carried in a fun and
entertaining narrative tone.
Erase una vez hace mucho tiempo en un remoto país, que un joven soldado llegó a un cuartel
llevando tras de sí una reata de 28 caballos. Como es preceptivo, en la entrada le dieron el alto.
El soldado mostró su credencial y el centinela llamó al cabo de guardia, que siguiendo la cadena
de mando lo acompañó al despacho del sargento al que, después de pedir permiso, se presentó
con el saludo de rigor.
- A sus órdenes mi sargento. Traigo conmigo los 28 caballos solicitados por este regimiento.
Están amarrados y en estado de revista en la entrada al cuartel.
- Bien, soldado. Los estábamos esperando pues no tenemos ninguno. Llévalos al patio y
después repártelos cabalmente entre las siete cuadras iguales que hay en él.
- Con gusto cumpliría su orden, mi sargento, pero el problema es no sé cómo hacerlo, pues
soy analfabeto y de números, sólo sé contar con los dedos, y gracias.
El sargento, se rascó la cabeza bajo la gorra, torciendo el gesto.
- Es el caso, soldado, que aunque yo sí entiendo de cuentas, sólo sé sumar, algo que me
enseñó a hacer el capitán, y esa operación requiere otros conocimientos superiores. Ve a su
despacho y que él te resuelva el problema.
45
José Miguel Bel Martínez Cuentos Matemáticos
El soldado se presentó de igual manera al capitán, quién le dijo que él, enseñado a su vez
por el coronel, aparte de sumar y restar, sólo alcanzaba a multiplicar, operación muy útil, pero
insuficiente para la complejidad de aquel problema. Por ello, lo envió al coronel, para quién las
matemáticas no tenían secreto alguno ya que incluso sabía dividir.
Tras solicitar su permiso, entró en su despacho, donde le puso al corriente de las diligencias
hechas con sus subordinados.
El coronel se retrepó en la silla sonriendo con suficiencia.
- No te apures soldado, estás de suerte puesto que además de ser la máxima autoridad del
regimiento, soy el único que domina con soltura la difícil operación de dividir, necesaria para
este menester, así que yo te diré cuantos caballos has de meter en cada cuadra para que queden
equitativamente repartidos.
Dicho esto, cogió papel y lápiz y se puso a la tarea. Escribió, como debe ser, el dividendo a
la izquierda y el divisor a la derecha encerrado en su caja. Entre dientes, razonó:
- Dos entre siete, no cabe, por tanto cojo el ocho, que sí cabe. Ocho entre siete, a uno. Uno
por siete, siete. De siete a ocho, uno y me bajo el dos, lo que nos da veintiuno, que entre siete
son tres. Como el resto es cero, el resultando de la división es exactamente trece.
28 7
21 13
0
- Ese es el número de caballos que has de meter en cada cuadra. Ve, hazlo y después di al
sargento que te den vino, un chusco y un catre donde descansar. Puedes retirarte.
- A las órdenes de usía mi coronel -se despidió el soldado con un marcial taconazo.
Al pasar por la puerta del despacho del capitán éste le llamó interesándose por el resultado
del problema, deseando impresionar con sus habilidades matemáticas al soldado.
- Si el coronel te ha dicho que son 13 los caballos que debes encerrar en cada cuadra, así será.
No obstante, no estará de más que yo lo compruebe para asegurarnos de que esté bien hecha.
Si multiplicamos los trece caballos que has de meter en cada cuadra por el número de estas, o
sea, siete, forzosamente nos han de salir los veintiocho que traes. Veamos pues.
Con no menos pericia que el coronel, fue recitando los pasos de la delicada operación.
- Siete por uno, son siete. Vamos con el otro. Siete por tres, veintiuno, que sumado al siete
de antes, hacen veintiocho: Justo y cabal, soldado. El coronel, como siempre, está en lo cierto.
Cumple pues su orden.
13
7×
7
21 +
28
Tal y como le había ordenado el coronel, el joven se dirigió al despacho del sargento, mas no
le hizo falta entrar, pues antes de llegar a la puerta, ya salió aquel a su encuentro ansioso por
poner en práctica su dominio de la suma e impresionar también al soldado. Una vez que hubo
conocido el resultado obtenido por el jefe y avalado por el oficial, dijo:
- Si el coronel y el capitán han calculado que debes meter trece caballos en cada cuadra,
así habrá de ser, pero por asegurarnos y también porque veas la utilidad de la suma que, al
fin y al cabo, no deja de ser una multiplicación más trabajada, pasa y observa como hago la
comprobación para que vayas sabiendo de cuentas, por si con los años medras en la milicia y
llegas a ser clase de tropa o hasta suboficial, como yo.
Dicho lo cual, sacó lápiz y papel y dispuso, como debe ser en una suma bien hecha, los siete
treces en una columna para después, sin encomendarse a Dios ni al diablo y sin hacer distingos
entre izquierda, derecha, arriba o abajo, sumó de corrido sin dejar ni uno todos los números que
tenía delante.
- Así que tenemos. . . uno más tres, mas uno más tres, mas uno más tres, mas uno más tres,
mas uno más tres, mas uno más tres, mas uno más tres. . . que hacen en total de. . . Mmmm. . .
¡Veintiocho!
13
13
13
13
13
13
13
28
- ¿Ves, como las matemáticas nunca fallan? Mete tranquilo las trece bestias por cuadra y ven
después a que te facilite acomodo para que pases la noche.
El soldado llevó los caballos hasta los establos, abrió las puertas de aquellos reducidos espa-
cios y contándolos con cuidado para no equivocarse, empezó a introducirlos una a uno.
Uuuno, dooos, treees, cuaaaatro, ciiinco, seeeis. . .
Hasta seis llegaron a entrar; y ni uno más, pues los animales, espantados al verse hacina-
dos en aquel cubículo desconocido y oscuro, organizaron tal barahúnda de coces, relinchos y
mordiscos, que las paredes de la cuadra amenazaban con quebrarse. El joven, asustado, los hizo
salir de nuevo, calmándolos después como pudo.
Veamos: tres hombres sabios no pueden estar equivocados y las matemáticas esas, de las
que tan bien hablan todos, tampoco -se dijo-. Piensa, Rufino-que así se llamaba el soldado-, y
cumple bien la orden que te han dado si es que quieres hacer carrera en la milicia.
Trece, trece, trece. . . -repetía angustiado para sí-. Recordó haber visto escrito aquel número
en todas las operaciones: Trece. Y de golpe, una luz iluminó su entendimiento. ¿Qué es un trece
sino un uno y un tres? Esperanzado, metió en la primera cuadra un caballo atravesado al fondo
y tres perpendiculares a él en la parte delantera. Un uno y un tres. O sea, un trece. Trece caballos
cómodamente ubicados. Cerró la puerta y repitió la misma operación en las otras seis cuadras,
comprobando aliviado que no sobraba ni faltaba ninguno.
Mas poco duró la tranquilidad al bueno de Rufino pues el coronel, alarmado por el alboroto
que se había organizado hacía un momento, bajó a ver qué sucedía.
- No se preocupe, mi coronel, que consciente de mi error, lo he enmendado y cada cuadra
está ocupada por los trece caballos que usía indicó. Ahora duermen tranquilos en ellas.
- Eso parece, soldado, pero ya que he bajado, quiero comprobar que lo que dices es cierto.
Abre las cuadras y contemos los animales que pernoctan en ellas.
- Es el caso, mi coronel, que ya están cerradas con llave y además los caballos han hecho un
largo viaje. Es lástima que haya que despertarlos.
- Razón llevas muchacho, mas no hará falta tal cosa puesto que por suerte yo estoy aquí y
con una simple división podremos contar los equinos sin necesidad de abrir ninguna puerta.
Échate al suelo y por el hueco que hay bajo una de ellas, cuenta las patas que veas.
El soldado obedeció al instante y, no sin trabajo, las contó.
- Cabalmente, cuento dieciséis, mi coronel.
- Bien, pues dividamos dieciséis por las cuatro patas que tiene un caballo y el resultado nos
dará el número de ellos que hay dentro.
A falta de papel y lápiz, el coronel se agachó y con el dedo en la tierra del suelo hizo la
consabida cuenta que recitó también en voz alta para que lo viera el ignorante soldado.
16 4
12 13
0
- Uno entre cuatro, no cabe, pasemos pues al seis. Seis entre cuatro, a uno; uno por cuatro es
cuatro; cuatro a seis, dos y me bajo el uno. Doce entre cuatro, tres. Tres por cuatro, doce, al doce,
cero.
Podemos dormir tranquilos, muchacho, y jurar ante Dios que trece, ni uno más ni uno me-
nos, son los caballos que ahora mismo duermen en cada cuadra.
Y ahora, querido lector, ¿no crees que la Matemática, además de ciencia exacta, puede ser
benévola con la ignorancia y también dejar lugar para la fantasía?
Sobre el autor:
Nombre: José Miguel Bel Martínez
Correo Electrónico: [email protected]
Institución: Ingeniero Técnico en Topografía (Colegiado Nº 492), Instituto Geográfico Nacional ,
España.
1 de octubre de 2013
Resumen
La idea de arquitectura paramétrica tiene una vertiente instrumental en el uso de
software orientado a la gestión de la información del proyecto arquitectónico, en particular
cuando se emplean geometrías complejas. Sin embargo, la idea de que la forma, sin
importar su geometría, esté matemáticamente definida mediante un número limitado de
parámetros puede ser de gran alcance en lo que se refiere a las posibilidades de generación
de variaciones por modificación de dichos parámetros, manipulación de la forma en
procesos de optimización y por la relación entre las estructuras matemático-formales
producidas y la idea clásica de unidad formal de la obra de arte.
Abstract
The idea of parametric architecture has an instrumental side in the use of software
oriented to the management of the information of the architectural project, in particular
when complex geometries are involved. But the idea of defining form, regardless its
geometry, by means of a limited number of parameters can be far reaching on account of
the possibilities of generating variations by the modification of these parameters, the
manipulation of form in optimization processes and the relation between the mathematical-
formal structures generated and the classical concept of formal unity in a work of art.
1. Introducción:
Una de las consecuencias del uso del ordenador en la práctica del diseño
arquitectónico ha sido el desarrollo de la llamada “arquitectura paramétrica”, la cual consiste,
49
Óscar del Castillo Sánchez Investigación
a grandes rasgos, en el uso de modelos informáticos del proyecto tanto en la concepción como
en la representación y materialización de edificios y sus componentes.
En este momento, la actividad en el ámbito de la arquitectura paramétrica abarca, por
un lado, la reflexión teórica acerca de las posibilidades de este nuevo paradigma, como puede
encontrarse por ejemplo en [1], y por otro, en el desarrollo de proyectos arquitectónicos que,
haciendo uso del software especializado correspondiente, trata de explotar las posibilidades
formales que, tanto a nivel de diseño como de construcción, estas nuevas herramientas
permiten.
Hay, sin embargo, un aspecto poco explorado de las posibilidades de la idea de
arquitectura paramétrica: la analogía entre la idea clásica de unidad formal –idea de unidad
de naturaleza esencialmente matemática-, y la forma concebida desde un punto de vista
paramétrico y algorítmico.
En este sentido, el nuevo paradigma proporcionaría una poderosa herramienta para
pensar la forma, no solo las superficies de doble curvatura de las últimas arquitecturas sino
también las formas ortogonales que conforman la mayoría de los edificios que se construyen.
Esta auto-limitación en la aplicación del método a formas aparentemente convencionales
pretende demostrar la relevancia de la idea de arquitectura paramétrica para toda actividad
de concepción formal, con independencia de la naturaleza de las formas utilizadas.
Lo que aquí se propone es ensayar y apuntar, mediante algunos ejemplos, líneas
concretas de desarrollo de estas ideas, en una especie de serialismo formal que busca
resultados análogos a los obtenidos en otros ámbitos, en especial el de la música
contemporánea.
2. Retículas espaciales
Observación: en el apéndice al final del presente trabajo se incluye la representación gráfica en
planta y alzado de los proyectos de los cuales se han tomado las redes espaciales utilizadas aquí como
ejemplo, mostrando la correspondencia entre la retícula espacial descrita y la realidad arquitectónica de
la que esta se abstrae.
Cuando en el proyecto arquitectónico se hace uso de una geometría ortogonal y de las
proyecciones planta / alzado / sección para su representación, la forma se define con
frecuencia mediante una configuración de ejes o planos orientados en las direcciones de los
tres ejes coordenados, una retícula o malla tridimensional que pauta el espacio y en referencia
a la cual se disponen los elementos que integran la forma. Esta retícula, que en principio tiene
un carácter meramente instrumental, facilitando tanto el control del proyecto como la
transmisión de la información necesaria para su construcción, se convierte en última instancia
en una herramienta de concepción de la forma [2], por cuanto permite definir a priori
condiciones acerca de las relaciones dimensionales y formales entre distintas partes del
proyecto. La malla espacial es por tanto un instrumento para pensar la forma.
Desde el punto de vista de la arquitectura paramétrica, las dimensiones que definen la
retícula espacial del proyecto pueden ser consideradas como otros tantos parámetros que
definen la forma y por tanto el universo de las posibles retículas espaciales constituye un
espacio de configuraciones.
Este tipo de retícula es muy habitual en edificios de toda índole; se pueden citar como
ejemplos el proyecto de Mies van der Rohe para el Convention Hall de Chicago [3], la iglesia
Bargsvaerd en Copenhague, de Jorn Utzon [4], o el Pabellón de Cristal de la Casa de Campo
de Madrid, de Francisco de Asís Cabrero et. al. [5].
La definición de esta forma requiere al menos cuatro parámetros (nx, ny, nz, K)
correspondiéndose estos al número de celdas en las direcciones x, y, z y a la dimensión de la
celda cúbica unidad, respectivamente. El universo de las retículas cúbicas posibles constituye
por tanto un espacio de configuraciones de 4 dimensiones. Así, en el caso mencionado del
Convention Hall de Mies van der Rohe, la retícula espacial sería:
nx = 24; ny= 24; nz= 4; K= 30’ (914,4 cm).
Este es un edificio de 720 x 720 pies en planta (24 tramos de 30’ cada uno en cada
dirección) y 110 pies en alzado (≈4 tramos de 30’). En la iglesia Bagsvaerd: (nx, ny, nz, K) = (36,
10, ≈7, 220 cm) y en el Pabellón de la Casa de Campo (solo la caja de vidrio): (nx, ny, nz, K) =
(51, 29, ≈8, 250 cm). A los efectos de este estudio se descartan, para simplificar en primera
aproximación, las pequeñas desviaciones de nuestro esquema en la retícula, como es el caso
de dimensiones distintas en tramos extremos. Una sección en tres dimensiones del espacio de
configuraciones para K = constante para los tres casos anteriores sería tal y como se muestra
en la fig. 3.
Figura 3. Sección tridimensional del espacio de configuraciones de redes espaciales cúbicas para un valor de K (lado de
la celda unidad) determinado. Posición en el espacio de configuraciones de cada uno de los ejemplos dados.
Si en lugar de una celda cúbica se parte de una celda ortoédrica con sus tres lados
distintos, cada configuración quedará definida por seis parámetros: (nx, ny, nz, Kx, Ky, Kz),
correspondiendo éstos al número de celdas en las direcciones x, y, z más la dimensión de la
celda unidad en cada una de las direcciones anteriores. Por tanto, el espacio de
configuraciones correspondiente será de dimensión 6, y al igual que en el caso anterior podrán
realizarse secciones tridimensionales de dicho espacio fijando algunos parámetros.
Un ejemplo de este tipo de retícula sería la biblioteca de Rodovre en Copenhague, de
Arne Jacobsen [6], en que: (nx, ny, nz, Kx, Ky, Kz)= (17, 6, 1, 440, 770, 385), con valores de K en
centímetros (fig. 4).
Figura 4. Red espacial correspondiente a la biblioteca de Rodovre; ejemplo de repetición de celdas ortoédricas en las tres
direcciones del espacio con lados de dimensiones diversas. Sección del espacio de configuraciones para (Kx, Ky, Kz)= (770, 440,
385).
clasicismo para c son: 2 (doble cuadrado en alzado), φ ≈ 1,618 (sección áurea), √ , √ , o bien
relaciones fraccionarias como 3:2, 5:3, etc., de modo que la retícula queda definida del modo:
(nx, ny, nz, K, K, c*K) o en el caso más general (nx, ny, nz, K, cy*K, cz*K), aunque no se produce
aquí reducción en el número de parámetros.
Como ejemplo de lo anterior puede tomarse, de entre muchos otros, la biblioteca
universitaria en Salamanca, obra de Carlos Puente [7], en que la estructura principal se
inscribe en una cuadrícula de 615 x 615 cm y una distancia entre ejes de forjados de 357 cm.
Por tanto: (nx, ny, nz, K, K, c*K)= (8, 5, 5, 615, 615, 357), donde c= 615 / 357 = 1,737 ≈ √ (fig. 5).
En este caso, del uso de la proporción √ permite que los rectángulos definidos en alzado por
la estructura porticada sean subdivisibles en tres rectángulos verticales cuya proporción es a
su vez √ , ya que el rectángulo de proporción 1:√ es subdivisible en n rectángulos
semejantes a él (a este respecto, ver el alzado mostrado en el apéndice).
Figura 5. Red espacial correspondiente a la biblioteca de Salamanca; ejemplo de repetición de celdas ortoédricas en las tres
direcciones del espacio con lados de dimensiones ligadas entre sí. Sección del espacio de configuraciones para (Kx, Ky, Kz)= (615,
615, 357), con: (nx, ny, nz)= (8, 5, 5).
Esto en cuanto a la repetición pura de celdas iguales en las tres direcciones del espacio;
sin embargo, también son posibles otras combinaciones como es la alternancia de dos medidas
en cada dirección del espacio, lo que podría denotarse: (nx, ny, nz, Kx1, Kx2, Ky1, Ky2, Kz1, Kz2), en
cuyo caso se necesitan 9 parámetros para definir la retícula. Tal es el caso del proyecto de
Louis Kahn para el centro comunitario de Trenton [8], en que: (nx, ny, nz, Kx1, Kx2, Ky1, Ky2, Kz) =
(20 ½ , 10 ½ , 4, 10’, 20’, 10’, 20’, 8’); los valores fraccionarios de n dan cuenta del hecho de que
la retícula comience y termine con la misma dimensión. En este caso no se produce alternancia
en dirección OZ (fig. 6).
Este es un tipo de retícula muy característico de la arquitectura Beaux-Arts, en que los
espacios de circulación y servicio tienden a ocupar las bandas de dimensión menor de la
retícula. Una derivación de esto es la retícula del tipo “tartán” (en referencia a los tejidos de las
faldas escocesas), basada en una alternancia de dimensiones que por sus propiedades métricas
es a priori adecuada para albergar los usos a incluir en el edificio.
Figura 6. Red espacial correspondiente al centro comunitario en Trenton: (nx, ny, nz, Kx1, Kx2, Ky1, Ky2, Kz)= (20 ½ , 10 ½ ,4, 10’,
20’ 10’, 20’, 8’); espacio de configuraciones para los K dados.
{60, 120, 180, 240, 360, 540, 600, 660, 720, 780, 840, 900, 960...} (unidades en cm). Partiendo de
esta escala se generará automáticamente, a modo de ejemplo, una serie de 16 redes espaciales
que constituyen posibles estructuras geométricas para un mismo proyecto arquitectónico (fig.
7).
Para ello se mantienen constantes los valores (nx, ny, nz, Kz) donde (nx, ny) = (7, 5), y el
valor de (nz, Kz) se supondrá irrelevante. Los parámetros que se hacen variar son (Kx1, Kx2, Ky1,
Ky2), si bien Kx1 y Kx2 por un lado y Ky1 y Ky2 por otro son variables ligadas entre sí de modo
que Kx1 + Kx2= 1200 cm y Ky1 + Ky2= 1200 cm. Por tanto, aunque el espacio de las variaciones
resultantes es en principio de dimensión 4, dado que las medidas en cada dirección están
ligadas dos a dos este espacio puede representarse en un plano (fig. 8).
De este modo es posible generar sistemáticamente las posibles redes espaciales
asociadas a un proyecto determinado y a la escala dimensional que se considera adecuada al
problema en cuestión. Una vez generadas será igualmente posible manipular estas redes
espaciales mediante algoritmos predefinidos.
Figura 7. Redes espaciales generadas mediante trayectorias en el espacio de configuraciones, para (nz, Kz) = constante.
Se representa solo la planta.
Figura 8. Espacio de variaciones del caso anterior, y sección del espacio de configuraciones para (Ky1, Ky2)= constante,
mostrando la trayectoria en dicho espacio descrita en la generación de las variaciones de la figura x. Dicha trayectoria consiste en
una recta debido a la condición impuesta: Kx1 + Kx2= 1200 cm.
Como ejemplo de esto, las redes de la fig. 9 (lado izquierdo) se han generado del
modo siguiente: se dispone una estructura básica consistente en un rectángulo horizontal de
3000 x 1440 cm, dividido en tres bandas horizontales de 240, 480 y 720 cm. A continuación,
para cada banda se genera una serie cuyos términos son elementos de la escala dimensional
de 7 términos E= {120, 240, 360, 600, 720, 840, 960}, tomados aleatoriamente. Esta serie,
consistente en múltiplos de 120 cm al igual que la dimensión de 3000 cm, es tal que siempre es
posible descomponer la dimensión total de 3000 cm en una suma de sumandos de esta serie,
sin resto.
respecto, puede consultarse [10] para una discusión acerca de las limitaciones del LCG y
encontrar listados en lenguaje C con generadores alternativos.
En las fig. 9 se muestran las redes generadas según lo anterior (izquierda) y algunos
ejemplos de aplicación del mismo procedimiento sobre un volumen cúbico (derecha).
Para el análisis de este ejemplo se ha partido del alzado publicado en [12]. A falta de
un plano con todas las cotas, se ha medido sobre el dibujo escalado, por lo que puede haber
un mínimo error en las mediciones. Pequeñas imprecisiones aparte, el alzado parece estar
pautado según un módulo de 30 cm, de modo que todas las pilastras tienen una anchura de 60
cm y los huecos entre ellas responden a las siguientes medidas:
A= 30 cm / B= 60 cm / C= 90 cm / D= 120 cm / E= 150 cm / F= 180 cm / G= 300 cm.
Figura 11. Alzado oriental de la ampliación del ayuntamiento de Murcia, obra de Rafael Moneo. Dibujo del autor.
Figura 12. Ejemplos de secuencias aleatorias de sólido/vacío utilizando la misma escala de dimensiones que en la figura 11.
Existen otros ejemplos de asignación aleatoria de atributos distintos del par sólido /
vacío, como la Colorium Tower en Dusseldorf, obra de Alsop Architects [13], en que los
atributos asignados son colores de una paleta predefinida; sin embargo, la discusión se
limitará aquí al par mencionado.
Figura 13. “Ciudad – Alea”; imagen obtenida mediante generación de una malla tridimensional consistente en una alternancia a
la que se aplica una eliminación aleatoria de planos.
Como ilustración de esto, se propone el siguiente ejemplo: se parte de una red espacial
R1 (nx, ny, nz, Kx, Ky, Kz)= (6, 4, 4, 720, 720, 360), representada en la fig. 15. A continuación se
define el catálogo de atributos a asignar a las celdas con caras contenidas en la envolvente de
la red. En este caso dichos atributos consistirán en una serie de huecos que forman ellos
mismos una secuencia tal que garantiza no solo la unidad formal del catálogo mismo sino
también una relación de autosemejanza con la red definida, a la vez que se satisfacen otros
requisitos puramente arquitectónicos como la legibilidad de las dimensiones del objeto.
Figura 15. Red espacial R1 (nx, ny, nz, Kx, Ky, Kz)= (6, 4, 4, 720, 720, 360).
En esta red las celdas tienen todas las mismas dimensiones en alzado (720 x 360 cm),
lo que facilita la definición del catálogo de atributos. En principio, el número de atributos a
emplear será tal que responda a la simetría del patrón de aplicación a utilizar, de modo que
no se produzcan adyacencias no deseadas de celdas con atributos iguales. Por otro lado,
según el teorema de los cuatro colores (ver por ejemplo [15]), para colorear un mapa
cualquiera sin que aparezcan regiones adyacentes del mismo color, basta utilizar cuatro
colores distintos.
Forma un tanto enrevesada de definir en una sola expresión y con un solo parámetro
(n) la dimensión y la orientación alternante de los huecos, haciéndose depender esta última de
la paridad de n (n par → hueco vertical / n impar → hueco horizontal), de modo que:
n=1 → (lado_horizontal, lado_vertical)= (120 cm, 60 cm).
n=2 → (lado_horizontal, lado_vertical)= (120 cm, 240 cm).
n=3 → (lado_horizontal, lado_vertical)= (480 cm, 240 cm).
Como se puede comprobar, esto se corresponde con las dimensiones de los huecos tal
como se indica en la fig. 16. En esta notación, lado_horizontal corresponde a la medida del
lado horizontal del hueco y lado_vertical con la del vertical.
El último paso consiste en la definición del patrón de asignación de atributos.
Suponiendo cada cara del ortoedro de partida dividido en cuatro regiones, un patrón posible
sería el mostrado en la fig. 17. De entre las varias configuraciones posibles, en esta se evitan
caras adyacentes con el mismo atributo y el patrón resultante posee interesantes propiedades
adicionales en lo que a su simetría se refiere.
La subdivisión del ortoedro de partida (fig. 15) puede realizarse de modo parecido a
como se han generado redes espaciales de forma sistemática, solo que por tratarse de un
espacio de solo tres parámetros, relativamente reducido, podrían obtenerse todas las
combinaciones posibles fácilmente. En la fig. 18 se muestran tres ejemplos de subdivisión de
dicho ortoedro de cara a la asignación de atributos.
Figura 18. Subdivisión del ortoedro de partida en regiones. El parámetro que define las subdivisiones en cada cara del ortoedro se
define a partir del número en base 2 obtenido asignando valor 0 a las celdas que no se corresponden con subdivisión y 1 a las que
sí se corresponden con ésta.
un atributo n+2 en la cara opuesta del ortoedro). Estas simetrías proporcionan a la forma un
grado superior de unidad, y son una buena razón para la aplicación de este tipo de
procedimientos.
Figura 19. Aplicación, sobre la subdivisión (1, 4, 1)10 de la figura 18, de los atributos definidos en la figura 16 asignados según el
patrón descrito.
Con este ejemplo queda por tanto descrito el proceso completo de definición del
proyecto desde la retícula espacial hasta la definición de su estructura de llenos y vacíos
mediante el uso de algoritmos. El alzado de la fig. 19, con la secuencia de huecos obtenida,
prácticamente define la estructura interna del objeto. La forma obtenida es deliberadamente
convencional en orden a demostrar cómo la simetría entendida en sentido clásico se relaciona
con la idea de forma generada algorítmicamente, solo que precisamente por su carácter
algorítmico es fácil generar variantes ya sea con una intención heurística, ya sea con la de
someter las formas obtenidas a procesos de optimización.
3.3 Autosemejanza
La idea matemática de autosemejanza tiene su equivalente en arquitectura en la idea
de coherencia formal como uno de los recursos más extendidos para la unidad formal y la
simetría. La repetición de un único tema formal en los distintos niveles del proyecto es
análoga a la idea de autosemejanza que caracteriza a los objetos fractales, bien entendido que
se trata de una analogía instrumental, casi una metáfora, ya que en el objeto arquitectónico,
con sus requerimientos materiales, funcionales y contextuales, es difícil aplicar rigurosamente
esta idea, y además el concepto de fractal implica el paso al límite para que aparezca la
dimensión fraccionaria, lo que evidentemente no ocurre en arquitectura.
Es evidente que pueden emplearse otros iniciadores, tanto más simples (un plano, por
ejemplo) como más complejos (U’s, Z’s); sin embargo se ha obtenido una gran variedad de
objetos solo con modificar ligeramente algunos de los parámetros que lo definen, por lo que
se limitará aquí el ámbito de búsqueda a este iniciador tipo.
Figura 21. Iniciador en “L” y giros del iniciador en torno a los ejes coordenados.
Figura 22. Representación isométrica de la configuración L (lx, ly, lz, e, gx, gy, gz, vx, vy, vz, n)= (100, 100, 100, 5, -60, -60, -60, 90,
90, 0, 20, 20, 20, 3). El origen de coordenadas se sitúa en el centro de la figura.
(160, 160, 100, 5, -70, -70, (100, 100, 100, 5, -35, -35, (100, 200, 300, 10, -60, -60, (100, 100, 100, 5, 0, 0, 0, 0, (100, 100, 100, 5, -20, -20, -
-70, 90, 90, 0, 0, 0, 40, 3) -35, 90, 90, 0, 0, 0, 0, 3) -60, 0, 90, 90, 0, 0, 0,3) 90, 0, 0, 50, 4) 20, 0, 90, 0, 0, 50, 4)
Figura 23. Ejemplo de generación de formas mediante la modificación de los parámetros que las definen.
4. Conclusión
El presente trabajo ha pretendido ser una exploración inicial, mediante una serie de
ejemplos, en el uso de algoritmos matemáticos para la generación de formas arquitectónicas, y
en la consiguiente aplicación a la forma arquitectónica de una línea experimental ya presente
en otras actividades artísticas como la música, escultura o pintura. Las formas aquí empleadas
han sido deliberadamente convencionales en orden a demostrar tanto el carácter matemático
del concepto tradicional de unidad formal como la aplicabilidad de estas técnicas al común de
la invención arquitectónica que de ello se deriva.
L1 (1870, 1870, 1155, 55, -1065, -880, -860, 90, 90, 0, 0, 0, 0, 3); L2 (1045, 1072, 1100, 55, -220, -220, -220, 90, 90, 0, 0, 0,
0, 3); L3 (1471, 1045, 660, 11, -982, -778, -798, -90, 90, 0, 0, 0, 0, 3).
5. Apéndice
Figura 26.Convention Hall. Arquitecto: Ludwig Mies van der Rohe. Chicago (EE.UU), 1953-54. Dibujo del autor según [3].
Proyecto para el Convention Hall de Chicago. Ejemplo de malla cúbica con K= 30‘ y
(nx, ny, nz, K)= (24, 24, 4, 30’). La planta baja tiene una altura de 20’, por lo que no encaja en la
malla de 30’, aunque sí encajaría en una malla cúbica de K= 10’ que incluiría el nivel
constructivo (por ejemplo, el despiece de fachada).
Figura 27.Iglesia en Bagsvaerd. Arquitecto: Jorn Utzon. Copenhague (Dinamarca), 1973-76. Dibujo del autor según [4].
Figura 28. Pabellón de cristal de la Casa de Campo. Arquitecto: Francisco Cabrero. Madrid (España), 1964. Dibujo del autor
según [5].
Ejemplo de malla cúbica en que K= 250 cm. Esta malla regula sobre todo el despiece de
las fachadas y el intereje entre correas de cubierta, siendo la distancia entre pórticos de 500 cm,
es decir, 2*K.
En los planos puede comprobarse que: (nx, ny, nz, K)= (51, 29, 8, 250).
Figura 29. Biblioteca en Rodovre. Arquitecto: Arne Jacobsen. Copenhague (Dinamarca), 1961-69. Dibujo del autor según [6].
Ejemplo de malla anisotrópica en que Kx ≠ Ky ≠ Kz. En este caso: A= Kx= 440 cm, B= Ky=
770 cm y C= Kz= 385 cm, si se considera que nz=1; sin embargo, de este modo no entrarían en
consideración ni la cubierta sobreelevada del salón de actos ni el sótano (no representado). Si
se toma C= 55 cm, correspondiente a la altura de las hiladas del despiece de piedra de la
fachada, se tendrá que la altura del sótano es de 275 cm= 5 * C, la altura del cuerpo principal
será de 385 cm = 7 * C, y la altura de la cubierta del salón de actos será de 165 cm hasta el
borde superior de las ventanas (3 * C) y de aproximadamente 4,5 * C= 245,5 cm hasta el borde
superior de la cubierta.
La configuración considerada en este trabajo es: (nx, ny, nz, Kx, Ky, Kz)= (17, 6, 1, 440,
770, 385).
Figura 30. Biblioteca universitaria en Salamanca. Arquitecto: Carlos Puente. Salamanca (España), 1997-2000. Dibujo del autor
según [7].
Figura 31. Proyecto de centro comunitario en Trenton. Arquitecto: Louis Kahn. Nueva Jersey (EE.UU), 1954-59. Dibujo del
autor según [8].
Referencias
[1] ORTEGA, Lluis (editor). La digitalización toma el mando. Editorial Gustavo Gili, Barcelona,
2009.
[2] LURÇAT, André. Formes, composition et lois d’harmonie (vol. V), p. 207, Éditions Vincent,
Fréal & Cie., París, 1957.
[3] CARTER, Peter. Mies van der Rohe at work, pp. 101-107, Phaidon, Londres, 1999.
[4] NIETO, Fuensanta y SOBEJANO, Enrique (editores). Jorn Utzon, pp. 87-93, Centro de
Publicaciones - Secretaría General Técnica del Ministerio de Obras Públicas, Transportes
y Medio Ambiente, Madrid, 1995.
[6] SOLAGUREN-BEASCOA, Félix. Arne Jacobsen, pp. 178, Editorial Gustavo Gili, Barcelona,
1991.
[7] REVISTA ARQUITECTURA COAM nº 323 (1º trimestre 2001), pp. 68-88, Revista de
Arquitectura y Urbanismo del Colegio Oficial de Arquitectos de Madrid, Madrid, 2001.
[8] RONNER, Heinz y JHAVERI, Sharad. Louis I. Kahn: Complete work 1935-1974, pp. 82-91,
Birkhäuser, Basilea, 1987.
[10] PRESS, William. J., TEUKOLSKY, Saul A., VETTERLING, William T. y FLANNERY,
Brian P. Numerical Recipes: the Art of Scientific Computing, pp. 340-377, Cambridge
University Press, Nueva York, 2007.
[11] MONEO, Rafael. Apuntes sobre 21 obras, p. 445, Editorial Gustavo Gili, Barcelona, 2010.
[12] REVISTA EL CROQUIS nº 98 “Rafael Moneo: 1995-2000”, pp. 84, El Croquis editorial, El
Escorial, 2000.
[13] POWELL, Kenneth. Will Alsop 1990-2000, pp. 234-235, Laurence King Publishers Ltd.,
Londres, 2002.
[14] ZEVI, Bruno. El lenguaje de la arquitectura moderna, pp. 17-24 (en particular imagen pág.
19), Editorial Poseidón, Barcelona, 1978.
[15] WILSON, Robin. Four colors suffice, Princeton University Press, Princeton, 2003.
Sobre el autor:
Nombre: Óscar del Castillo Sánchez
Correo Electrónico: [email protected]
Profesión: Arquitecto
1 de abril de 2014
Resumen
Abstract
Fractional Calculus provides an opportunity to spread the concepts of derivative and in-
tegral to not integer orders. In this context, we may consider a generalization of Newton’s
second Law in which the second derivative is replaced by other of α order (1, 5 < α ≤ 2). In
this research some classical systems (pendulum, projectile and spring) are discussed under
this new approach.
1. Introducción
Desde el punto de vista matemático, el estudio del movimiento de un cuerpo a partir de las
fuerzas que actúan sobre él se reduce a aplicar la ecuación fundamental de Newton F = ma,
teniendo en cuenta que F es la fuerza resultante o suma vectorial de las fuerzas que intervienen
sobre el cuerpo. Se trata, en suma, de resolver la ecuación diferencial de segundo orden
F = mx ′′ , (1.1)
77
Antón Lombardero Ozores Investigación
Puede ser interesante, como experimento físico-matemático, seguir este proceso desde una
perspectiva más amplia, que utilice como punto de partida una versión generalizada de la se-
gunda Ley de Newton. Esto es lo que se pretende con este trabajo, y para ello haremos uso del
Cálculo Fraccionario.
El Cálculo Fraccionario es una activa rama del Análisis Matemático que nace de una idea
muy básica: de la misma forma que una función se puede derivar o integrar un número entero
de veces, ¿podemos hallar su derivada de orden 1/2? ¿Y su integral π-ésima? Surgen así opera-
dores diferenciales e integrales de orden racional, real e incluso complejo. Con estos cimientos
se puede construir un nuevo edificio del Cálculo, que, sin perder de vista al Cálculo tradicional
(al que no sustituye, sino que extiende), posee sus propios conceptos fraccionarios. Un ejemplo,
entre otros muchos, es el de ecuación diferencial fraccionaria.
Como en cualquier área matemática viva, el interés del Cálculo Fraccionario no se restringe
al plano teórico, siendo este notable. Los campos de aplicación práctica se incrementan año tras
año. En particular, los modelos de derivadas fraccionarias proporcionan mejores representacio-
nes que los modelos tradicionales en aquellos materiales con amortiguación interna.
Volviendo a la segunda Ley de Newton (1.1), autores como Ebaid [4] proponen una genera-
lización fraccionaria de la misma en la forma
F = mD α x (1.2)
en la que D α x se entiende como la derivada de orden α de x (1, 5 < α ≤ 2). A partir de la ecua-
ción alternativa (1.2), seguiremos el procedimiento habitual para el estudio de la dinámica de
algunos sistemas físicos clásicos: péndulo, proyectil y resorte. Este nuevo enfoque fraccionario
de la dinámica conduce a resultados, ecuaciones e interpretaciones diferentes a las tradicionales,
que pueden resultar de interés en el futuro.
2. Historia
El Cálculo Fraccionario (CF) puede ser considerado una disciplina tanto antigua como nueva.
Antigua porque, de la mano del propio Leibnitz, ya da sus primeros pasos a finales del s. XVII,
casi paralelamente al nacimiento del Análisis. Nueva, porque no es hasta la década de los setenta
del pasado siglo cuando ha sido objeto intensivo de conferencias especializadas, tratados y
artículos académicos.
Hasta tiempos recientes, el CF tenía reputación de teoría matemática sin aplicaciones, de un
carácter un tanto esotérico. Pero en las últimas décadas se ha producido una explosión de la
actividad investigadora en este campo, debido a que encuentra aplicación en las más diversas
ramas científicas.
Sin embargo, fueron los padres del Análisis Matemático quienes dieron los primeros pasos.
Quizá estimulados por la propia notación de Leibnitz para la derivada de orden n,
dn
f (x) (2.1)
dx n
son éste y L’Hopital quienes, en una correspondencia del año 1695, se plantean por primera vez
el sustituir n por el valor no entero 1/2.
La denominación adoptada, Cálculo Fraccionario, perdura por razones históricas sin ser la
más adecuada: las derivadas no sólo serán racionales, sino también reales, e incluso complejas,
por lo que una designación más apropiada sería la de cálculo de orden arbitrario.
Tras los primeros tanteos de Leibnitz y L’Hopital, el CF siguió siendo abordado por algu-
nas de las mentes matemáticas más privilegiadas. Euler (s. XVIII), Laplace, Lacroix, Fourier,
Riemann (s. XIX), Hardy y Littlewood (s. XX) se ocuparon del tema, pero la gran cantidad de
problemas abiertos en esas épocas en el Análisis Matemático ordinario quizás fue lo que llevó a
que no profundizaran en este tipo de planteamientos alternativos.
Destacan, a principios del s. XIX, las figuras de Abel y Liouville, que propician el salto cua-
litativo desde la Prehistoria a la Historia del CF. Abel estaba interesado en el problema de la
tautócrona, es decir, en determinar la curva tal que el tiempo requerido por una partícula en
deslizarse por ella hasta su punto más bajo fuera independiente del punto desde el que se sol-
tase la partícula. Abel llegó a que dicha curva g( x ) debía satisfacer la ecuación integral
Z x
−1/2
( x − t) g(t)dt = k (2.2)
0
El término izquierdo de la ecuación es, excepto por la constante multiplicativa 1/Γ( 12 ) la inte-
gral fraccionaria de orden 1/2 de la función g( x ). Con ello Abel establece la primera aplicación
práctica del CF.
Por su parte, Liouville publica una serie de artículos en los que lleva a cabo el primer gran
intento de definir formalmente derivadas e integrales fraccionarias. Dichos trabajos hacen de
Liouville el artífice de las bases fundamentales del CF. Su punto de partida consistió en extender
la fórmula de la derivada de orden entero de una función exponencial a derivadas de orden α
arbitrarias:
D α e ax = aα e ax (2.3)
De los trabajos de Liouville y otros matemáticos posteriores surgen las integrales y deriva-
das fraccionarias de Riemann-Liouville (1870-1884). Estas definiciones son, junto con la inte-
gral fraccionaria de Gründwald-Letnikov (1867-1868), las generalizaciones más popularizadas
del CF, una disciplina en la que hay casi tantas propuestas alternativas de derivada e integral
fraccionaria como estudiosos del tema. Otras fórmulas integro-diferenciales fraccionarias son
la integral fraccionaria de Weyl (1917), la integral fraccionaria de Riesz (1936), o la derivada
fraccionaria de Caputo (1967).
Por fin, tras más de tres siglos en los que el CF se desarrolla como un campo puramente teó-
rico casi exclusivo de los matemáticos, varios autores señalan recientemente que las derivadas e
integrales no enteras son especialmente apropiadas para la modelización matemática de mate-
riales como los viscoelásticos y otros. Se ha demostrado que los modelos de orden fraccionario,
con sus propiedades de no localidad y memoria, son, en ciertos casos, más adecuados que los
clásicos de orden entero en áreas como Teoría de Materiales, Teoría del Transporte, Electromag-
netismo, Teoría de la Señal, Teoría del Caos o Fractales. Desde hace cuatro décadas empiezan
a aparecer los primeros tratados consagrados exclusivamente al tema, y el flujo de artículos,
libros, conferencias internacionales, paquetes de software y grupos de investigación dedicados
al CF se refuerza año a año.
El primer libro dedicado netamente a la materia es publicado en 1974 por Oldham y Spanier
[11]. Posteriormente, son destacables el libro enciclopédico de Samko, Kilbas y Marichev [16],
las obras de Miller y Ross [10], Podlubny [14] y el más reciente de Kilbas, Srivastava y Trujillo
[9].
Entre los retos actuales de los investigadores se encuentra el de dar una mayor consistencia,
claridad y rigor a la teoría matemática sobre operadores fraccionarios. Se hace necesaria una
unificación de los muchos conceptos que ahora hay dispersos, estableciendo algunas definicio-
nes estándar. Queda también abierto el encontrar interpretaciones físicas y/o geométricas a las
derivadas e integrales fraccionarias, de la misma forma que las tienen las de orden entero. En
esta cuestión hay propuestas interesantes, sobre las que destaca el estudio de Podlubny [15].
3. Operadores fraccionarios
Las derivadas e integrales fraccionarias son generalizaciones de las habituales que las inclu-
yen como casos particulares. Sin embargo no hay un consenso generalizado sobre qué son la
derivada y la integral fraccionaria de una función, y nos encontramos con una gran variedad de
propuestas. Definimos en esta sección algunos de los operadores fraccionarios más extendidos,
deteniéndonos especialmente en la integral y derivada de Riemann-Liouville. Estos operadores
facilitan la extensión con continuidad al semiplano {z ∈ C/ℜ (z) > 0} de las nociones clásicas
de derivadas e integrales enteras.
Para la simplificación del cálculo de la integral n-ésima de una función f , es bien conocida
la llamada Fórmula de Cauchy de la Integral Repetida,
Z x Z t1 Z t n −1
( Ian f )( x ) = dt1 ...
dt2 |{z} f (tn )dtn
a a a
n
1
Z x
= ( x − t)n−1 f (t)dt (n ∈ N) (3.1)
( n − 1) ! a
fácilmente demostrable por inducción. Dado que (n − 1)! = Γ (n) ∀n ∈ N, siendo Γ (n) la
función Gamma, observamos que el miembro derecho de (3.1) puede tener significado para
valores no enteros de n. Por lo tanto es natural definir la integral fraccionaria de orden no entero
como sigue.
Definición 1. Sea f ∈ L1 (a, b) (−∞ < a < b < ∞). La integral fraccionaria de Riemann-Liouville
de orden α ∈ C (ℜ (α) > 0) de f se define por
1
Z x
f (t)
( Iaα f ) (x ) := dt (x > a) . (3.2)
Γ (α) a ( x − t )1− α
1
Z b
f (t)
Ibα− f (x ) := dt ( x < b) (3.3)
Γ (α) x ( t − x )1 − α
lı́m ( Iaα f ) ( x ) = f (x )
α →0+
Para la definición de la derivada fraccionaria de orden α, una posibilidad podría ser tomar
la fórmula de la integral (3.2) y sustituir en ella α por −α. Esta solución plantea dos problemas.
El primero, muy evidente, es que para valores negativos de α la integral de (3.2) puede no ser
convergente, además de que para los enteros negativos Γ (α) no está definida. En segundo lu-
gar, debemos tener en cuenta que la derivación y la integración ordinarias h noison operaciones
n n (
recíprocas: es cierto que ( /dx) Ia f ( x ) = f ( x ), pero sin embargo Ia f n) ( x ) 6= f ( x ) ya
d n
que difieren en un polinomio de orden n − 1. La derivada entera es por tanto la inversa por la
izquierda de la integral, y pretendemos que los operadores fraccionarios hereden esta propie-
dad: D α I α = id . Para ello, la derivada fraccionaria de Riemann-Liouville se definirá como la
derivada entera de una cierta integral fraccionaria.
Definición 11. Sea f ∈ L1 (a, b) (−∞ < a < b < ∞). La derivada fraccionaria de Riemann-Liouville
de orden α ∈ C (ℜ (α) ≥ 0) de f se define, si existe1 , por
n
d
( Daα f ) (x ) : = Ian−α f ( x )
dx
n Z x
1 d f (t)
= dt ( x > a) (3.4)
Γ (n − α) dx a ( x − t ) α − n +1
n
d
donde n = [ℜ (α)] + 1 y dx es la n-ésima derivada usual.
1
Z x
d f (t)
Daiθ f (x ) = dt ( θ ∈ R \ {0 } ; x > a ) (3.5)
Γ (1 − iθ ) dx a (x − t)iθ
1 En el Teorema 19 se dan condiciones suficientes de existencia.
Al igual que en el caso entero, la derivada fraccionaria de una función f ( x ) existe bajo unas
condiciones más restrictivas que la correspondiente integral. No basta, como ocurría con Iaα f ,
con que f ( x ) sea de L1 . Antes de enunciar el resultado principal sobre existencia, necesitamos
la noción de continuidad absoluta.
Observación 17. Ser absolutamente continua es una propiedad bastante fuerte: implica continui-
dad uniforme (y por tanto continuidad). Sin embargo es menos restrictiva que ser Lipchiziana:
toda función de Lipschitz es absolutamente continua.
Observación 18. Denotamos por AC n [ a, b] (n ∈ N ) al espacio de funciones f en [ a, b] tales que
f ( n−1) ∈ AC [ a, b]. Es claro que AC1 [ a, b] = AC [ a, b] y que el espacio AC n [ a, b] consiste en todas
las funciones representables por una integral de Lebesgue n-múltiple sumada a un polinomio
de orden n − 1.
Corolario 20. Si f ∈ AC [ a, b], entonces Daiθ f existe en casi todo punto de [ a, b] para cualquier θ ∈
R \ {0} y se representa en la forma (3.5).
El siguiente resultado refleja que, análogamente a lo que ocurre en el caso entero, la derivada
fraccionaria de Riemann-Liouville es el operador inverso izquierdo de la integral fraccionaria
de Riemann-Liouville.
Proposición 23. Sea f ∈ L1 ( a, b) (−∞ < a < b < ∞) y α, β ∈ C (ℜ (α) , ℜ ( β) > 0). Entonces las
igualdades
( Daα Iaα f ) (x ) = f ( x ) (3.6)
β β−α
Daα Ia f (x ) = Ia f (x ) , si α ≤ β (3.7)
β α− β
Daα Ia f (x ) = Da f (x ) , si α ≥ β (3.8)
son ciertas en casi todo punto de [a,b].
con ci ∈ C ∀i = 1, 2, ..., n − 1.
y
10
0
0.25 0
0.5
0.25
0.75
8 1
0.5
1.25
1.5
0.75
1.75
2
6 1
1.25
4 1.5
1.75
2 2
0
0 0.4 0.8 1.2 1.6 2 2.4 2.8
Figura 1. Evolución de las derivadas fraccionarias D0α x2 , con α variando entre 0 (parábola) y 2 (función constante)
Observación 27. Como caso particular de (3.11), se tiene que cuando β = α − 1, α − 2, ..., α −
n (n = [ℜ (α)] + 1) la derivada de orden α-ésimo se anula, por lo que Daα anula las potencias
( x − a)α−1 , (x − a)α−2 , ..., (x − a)α−n .
Observación 28. Para el cálculo de la derivada de una función constante f ( x ) = k hacemos
tender β a 0 por la derecha en (3.11) e introducimos la constante multiplicativa k:
Γ ( β + 1) k
Daα k = k lı́m Daα ( x − a) β = k lı́m ( x − a) β−α =
β →0 β →0 Γ ( β + 1 − α) Γ (1 − α ) ( x − a ) α
n −1
f (k ) ( a )
C
Daα f ( x ) = ( Daα f ) ( x ) − ∑ (x − a)k−α ( x > a) (3.13)
k =0
Γ (k + 1 − α )
Observación 32. Por lo tanto, para órdenes de derivación no enteros, las derivadas de Caputo y
Riemann-Liouville coincidirán cuando se cumpla
f ( a ) = f ′ ( a ) = · · · = f ( n −1 ) ( a ) = 0 (3.14)
Riemann-Liouville
2 Caputo
-1
-2
0 4 8 12 16 20 24
x
Figura 2. Comparativa de derivadas
con ci ∈ C ∀i = 1, 2, ..., n − 1.
La derivada de Caputo puede considerarse por tanto como una regularización de la derivada
de Riemann-Liouville efectuada mediante la sustracción de un polinomio de Taylor, a través
de la cual se obtienen ciertas ventajas importantes, como las citadas condiciones iniciales en-
teras, unos requisitos menos restrictivos para el cumplimiento de la propiedad de semigrupo
C D αC D β f = C D α+ β f , y la no menos significativa característica que enunciamos a continuación.
a a a
Proposición 36. Sea α ∈ C (ℜ (α) > 0) y sea n = [ℜ (α)] + 1. Se verifica que C Daα ( x − a)k =
0 ∀k = 0, 1, · · · , n − 1.
Esta propiedad, demostrada en Kilbas, Srivastava y Trujillo [9], evidencia que el valor de la
derivada fraccionaria de Caputo operada en una constante es nulo.
Al contrario que los operadores de Riemann-Liouville o Caputo, que obtienen sus definicio-
nes de una integral repetida, el enfoque de Grünwald-Letnikov toma como punto de partida la
derivada, en concreto de la siguiente fórmula para la derivada n-ésima, fácilmente deducible:
n
d 1 n n
f ( x ) = lı́m n ∑ (−1)k f ( x − kh) . (3.17)
dx h → 0 h k =0 k
Esta expresión puede ser generalizada para valores no enteros de n, teniendo en cuenta que
para α > 0 podemos interpretar (αk ) como
α Γ ( α + 1)
= . (3.18)
k k!Γ (α + 1 − k)
Entonces, si consideramos una serie infinita se obtiene
α
d 1 ∞ k α
f (x ) ≈ lı́m α ∑ (−1) f ( x − kh) .
dx h → 0 + h k =0 k
Por otra parte,dado un h determinado, no tienen sentido términos de la serie para valores de
k superiores a x − a
h , ya que en ese caso f ( x ) tomaría valores fuera del intervalo de definición
[ a, b]. Se llega entonces a la definición que sigue.
Definición 37. Sea f una función acotada en [ a, b] (−∞ < a < b < ∞). La derivada fraccionaria
de Grünwald-Letnikov de orden α ∈ R (α > 0) de f se define, si existe, por
[ x−h a ]
α) 1 α
fa ( x ) := lı́m+ α ∑ (−1)k f (x − kh) ( x > a) (3.19)
h →0 h k =0
k
Todos los operadores fraccionarios que hemos tratado se definían sobre un intervalo acota-
do de la recta real [ a, b]. Las formas integrales de Riemann-Liouville y diferenciales de Caputo,
con esta característica, serán las que utilizaremos en los cálculos a lo largo de este trabajo. Pe-
ro existen otros tipos de operadores, definidos sobre intervalos no acotados de la recta real, o
incluso sobre el plano complejo. Otros están diseñados para funciones de varias variables.
El operador integral de Weil (también llamado de Liouville) es una extensión de la integral
fraccionaria de Riemann-Liouville al semieje no acotado (−∞, x ]:
1
Z x
α f (t)
( I+ f ) (x ) := dt (−∞ < x < ∞) .
Γ (α) −∞ ( x − t )1 − α
Esta forma, y su correspondiente forma diferencial, son especialmente adecuadas para realizar
integración y derivación fraccionaria de funciones periódicas.
Pero hay una gran variedad de operadores fraccionarios. Otros registrados en la literatura
especializada son el de Riesz, para derivadas e integrales fraccionarias de funciones de varias
variables, Erdélyi-Kober, Hadamard, Bessel, Chen o Dzherbashyan.
ejemplo, es conocido que una EDF de orden α necesita [ℜ (α)] + 1 condiciones iniciales para ser
resuelta de manera única. Sin embargo surgen también notables diferencias, la primera de las
cuales es el hecho de que una misma EDF tendrá diferentes significados según los operadores
fraccionarios implicados sean los de Riemann-Liouville, Caputo u otros.
Problemas análogos a los de Cauchy y Dirichlet para ecuaciones diferenciales surgen tam-
bién en el Cálculo Fraccionario. Aquí hemos de diferenciar explícitamente en función del tipo de
operador fraccionario, ya que hay diferencias importantes en la formulación de las condiciones
iniciales. Veremos los problemas de Cauchy para derivadas de Riemann-Liouville y Caputo.
Definición 44. Si utilizamos derivadas fraccionarias de Riemann-Liouville y queremos encon-
trar la solución y (x ) de la EDF
h i
α
( Daα y) ( x ) = f x, y (x ) , Daα1 y ( x ) , Daα2 y ( x ) , ..., Da m−1 y ( x ) (4.3)
no verifica que α − n tenga parte real positiva. Debe ser interpretada por tanto como la integral
fraccionaria
Ian−α y ( a) = bn .
donde n = [ℜ (α)] + 1, decimos que estamos ante un problema tipo Cauchy con derivadas fraccio-
narias de Caputo.
1
E0,1 (z) =
1−z
1
Z c +i∞
L−1 { F (s)} = f (t) = est F (s) ds
2πi c − i∞
n −1
L { Daα f (t)} = sα F (s) − ∑ s k f ( α − k −1 ) ( 0 ) (derivada Riemann − Liouville) (4.7)
k =0
n o n −1
L C
Daα f (t) = sα F (s) − ∑ sα−k−1 f (k) (0) (derivada Caputo ) (4.8)
k =0
Este resultado proporciona una nueva expresión para la derivada fraccionaria, válido tanto para
Riemann-Liouville como para Caputo siempre que las condiciones iniciales se anulen.
Una última herramienta útil para la búsqueda de las soluciones de una EDF es la convolu-
ción de dos funciones y su relación con la transformada de Laplace.
Definición 57. Dadas f , g : [0, ∞)→ C, f , g ∈ L1 [0, ∞), se define su convolución f ∗ g como la
nueva función Z x
( f ∗ g) ( x ) := f (t) g ( x − t) dt
0
Teorema 58 (de convolución). Dadas f , g : [0, ∞)→ C, f , g ∈ L1 [0, ∞) tales que existan las respec-
tivas transformadas de Laplace F (s) , G (s). Entonces
o, de forma equivalente,
L−1 { F (s) G (s)} = ( f ∗ g) ( x ) . (4.12)
Teorema 61 (para Riemann-Liouville). Dado el problema tipo Cauchy con derivadas fraccionarias
de Riemann-Liouville (4.3)-(4.4), si la función f [ x, y, y1, y2 , ..., ym−1] cumple la condición de Lipschitz
(4.13) para las variables y, y1 , y2 , ..., ym−1, y además se da que f [ x, y, y1, y2, ..., ym−1 ] ∈ L1 ( a, b) para
cualquier (y, y1, ..., ym−1 ) ∈ R m , se verifica que existe una solución única y (x ) al problema en el espacio
Lα ( a, b).
Teorema 62 (para Riemann-Liouville, caso lineal). Sea la EDF lineal (4.2) cumpliéndose que b ( x ) ∈
L1 (a, b). Si ak ( x ) ∈ L∞ ( a, b) o si ak ( x ) están acotadas en [ a, b] (para todo k = 0, 1, ..., m − 1), entonces
el siguiente problema tipo Cauchy con derivadas fraccionarias de Riemann-Liouville
m −1
α
( Daα y) (x ) = a0 (x ) y (x ) + ∑ ak (x ) Da k y ( x ) + b ( x )
k =1
Daα−k y (a) = bk , bk ∈ C (k = 1, 2, ..., n)
En este apartado resolveremos explícitamente algunos problemas tipo Cauchy cuyas solu-
ciones van a ser utilizadas después en el estudio fraccionario del péndulo y del proyectil. Se
trata de sencillas ecuaciones diferenciales fraccionarias de orden α ∈ (1, 2]. Por simplicidad a
la hora de traducir matemáticamente las condiciones iniciales, optamos de aquí en adelante
por la derivada fraccionaria de Caputo, en conjunción con la integral fraccionaria de Riemann-
Liouville. Como método para resolver las EDFs mostraremos dos caminos distintos: uno basado
en (3.16) y en las propiedades elementales de los operadores fraccionarios; otro haciendo uso
de la transformada de Laplace de forma semejante a la habitual para resolver ecuaciones dife-
renciales ordinarias.
Dada una función continua f : [0, ∞) → C y α ∈ (1, 2] vamos a resolver el problema tipo
Cauchy con derivadas fraccionarias de Caputo dado por
C
D0α y (x ) = f ( x ) ( x > 0) (4.14)
Como 1 < α ≤ 2 se tiene que n = [α] + 1 = 2, por lo que necesitamos dos condiciones
iniciales en el límite inferior de integración 0:
y (0) = b0 (4.15)
y′ (0) = b1 (4.16)
Si aplicamos el operador integral I0α sobre ambos miembros de la ecuación (4.14) obtenemos
I0α C D0α y ( x ) = ( I0α f ) ( x )
y (x ) = c0 + c1 x + ( I0α f ) ( x ) (4.17)
y (x ) = b0 + b1 x + ( I0α f ) ( x )
1
Z x
= b0 + b1 x + ( x − t)α−1 f (t) dt (4.18)
Γ (α) 0
de la EDO
y′′ (x ) = f ( x )
con las mismas condiciones iniciales (4.15)-(4.16).
Vamos a tratar ahora una ecuación algo más complicada, que será una generalización de la
conocida ecuación diferencial de oscilación y ′′ ( x ) = λy (x ) + f ( x ). Sustituyendo la segunda
derivada por una derivada de Caputo de orden α ∈ (1, 2] obtenemos una EDF que denomi-
naremos ecuación de oscilación fraccionaria. Para resolver el problema usaremos el método de la
transformada de Laplace. Este sistema tiene la ventaja de que haciendo uso de la transformada
podemos introducir las condiciones iniciales desde el primer momento, llegando directamente
a la solución particular.
Dada una función continua f : [0, ∞) → C, λ ∈ C y α ∈ (1, 2] vamos a resolver el problema
tipo Cauchy con derivadas fraccionarias de Caputo dado por
C
D0α y ( x ) = λy (x ) + f ( x ) ( x > 0) (4.19)
Como 1 < α ≤ 2 se tiene que n = [α] + 1 = 2, por lo que necesitamos dos condiciones
iniciales en el límite inferior de integración 0:
y (0) = b0 (4.20)
y′ (0) = b1 . (4.21)
s α −1 s α −2 F (s)
Y (s) = b0 α
+ b1 α
+ α .
s −λ s −λ s −λ
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 93
Antón Lombardero Ozores Investigación
Se trata ahora de aplicar la transformada inversa de Laplace a los dos miembros de la ecua-
ción, para recuperar y (x ). Utilizando la linealidad de L−1 ,
α −1 α −2
−1 s −1 s −1 F ( s)
y (x ) = b0 L + b1 L +L . (4.22)
sα − λ sα − λ sα − λ
Calculemos por separado las tres transformadas inversas. El resultado (4.10) proporciona el
valor de las dos primeras de forma directa:
α −1
−1 s
L = Eα,1 (λx α ) (4.23)
sα − λ
α −2
−1 s
L α
= xEα,2 (λx α ) (4.24)
s −λ
1
Para el cálculo de L−1 sα − λ F ( s ) utilizamos la propiedad (4.12),
−1 1 −1 1
L F (s) = L ∗ L−1 ( F (s)) (4.25)
sα − λ sα − λ
con lo que reducimos el problema a la obtención de L−1 (1/sα −λ). De nuevo volvemos a recurrir
a (4.10), estableciendo los valores k = 0, β = α y a = λ, con lo que resulta
−1 1
L = x α−1 Eα,α (λx α ) ,
sα − λ
que junto con (4.25) proporciona
−1 1 h i
L α
F (s) = x α−1 Eα,α (λx α ) ∗ f ( x )
s −λ
Z x
(x − t)α−1 Eα,α λ ( x − t)α f (t) dt
= (4.26)
0
Observación 63. Las funciones Eα,1 (λx α ), xEα,2 (λx α ) forman, por tanto, un sistema fundamental
de soluciones de la correspondiente EDF homogénea
C α
D0 y (x ) − λy (x ) = 0
y teniendo en cuenta que las funciones de Mittag-Leffler verifican Eα,1 (0) = Eα,2 (0) = 1,
Eα,α (0) = 1/Γ(α), llegamos a la solución (4.18) esperada:
1
Z x
y (x ) = b0 + b1 x + ( x − t)α−1 f (t) dt
Γ (α) 0
Observación 65. Gorenflo y Mainardi [5] han obtenido una solución alternativa equivalente para
este problema (con λ = −1). En su trabajo, la solución se expresa en términos de derivadas e
integrales de funciones de Mittag-Leffler de un solo parámetro Eα (z):
Z x d
y (x ) = b0 Eα (− x α ) + b1 I01 Eα (− x α ) + Eα (−tα ) f (x − t) dt
0 dt
El péndulo simple es un sistema físico constituido por una partícula de masa m que, sus-
pendida de un punto fijo O por medio de una varilla de longitud L, puede oscilar en un plano
vertical fijo por efecto de la fuerza de la gravedad. La posición de la partícula en el instante t se
especifica mediante el ángulo θ que la varilla forma con la vertical en ese momento.
El estudio del péndulo constituye uno de los problemas clásicos de la dinámica elemental, y
todas sus componentes están perfectamente determinadas desde el punto de vista matemático.
Lo que se pretende con este apartado es realizar un análisis del movimiento pendular enfocado
desde la óptica más extensa del Cálculo Fraccionario, y para ello se va a generalizar la idea
clásica de aceleración (derivada segunda de la posición) a una derivada fraccionaria de Caputo
de un orden comprendido entre 1.5 y 2.
Obviamente el péndulo simple es un sistema idealizado, al que para su análisis vamos a
exigir una serie de hipótesis simplificativas:
El sistema no pierde energía por efecto de la resistencia del aire ni por fricción alguna.
Para determinar la función del movimiento del péndulo fraccionario empezamos buscando la
ecuación diferencial fraccionaria que modeliza dicho movimiento. La ecuación buscada tiene
por incógnita la función θ (t), que determina en radianes el ángulo de la varilla en el instante t.
La partícula se mueve sobre un arco de circunferencia bajo el efecto de dos fuerzas: su propio
peso (mg, donde g = 9,80665 m/s2 es la fuerza gravitatoria ejercida por la Tierra cerca de su
superficie) y la fuerza de tensión T ejercida por la varilla. Si descomponemos el peso en sus
componentes tangencial y normal, observamos que la componente normal de esta fuerza se ve
contrarrestada por la fuerza de tensión de la varilla. Por lo tanto, la única fuerza actuante en lo
que concierne al movimiento del sistema será la componente tangencial del peso, que tendrá
signo negativo por ir siempre en dirección opuesta al movimiento:
F = ma, (5.2)
L
T
θ
mg sen θ
mg cos θ
mg
Figura 1. Fuerzas que intervienen en el péndulo
siendo a la aceleración tangencial del movimiento, y se llega así a la ecuación tradicional del
péndulo. Aquí vamos a sustituir (5.2), en la que la aceleración a (t) es la derivada segunda de la
posición r (t), por una fórmula alternativa más general haciendo uso de derivadas fraccionarias:
F (t) = m C D0α r (t) , (5.3)
C Dα θ
donde 0 (t) será la aceleración angular.
Finalmente, de (5.1), (5.3) y (5.4) obtenemos
C
−mg sen (θ (t)) = mL D0α θ (t)
5.1.2. Solución
La EDF (5.5) no tiene fácil solución, al no tratarse de una ecuación lineal debido al término
no lineal sen (θ (t)). Ahora, recordemos que sen (θ ) ≈ θ cuando θ es pequeño; en particular, θ y
sen (θ ) coinciden en las dos primeras cifras decimales cuando θ < π/12.
Por lo tanto, si nos restringimos a oscilaciones de pequeña amplitud (como las que describe,
por ejemplo, el péndulo de un reloj), parece razonable simplificar nuestro modelo matemático
Esta EDF encaja dentro del tipo de las estudiadas en el Método de la transformada de La-
place. Si establecemos las condiciones iniciales del péndulo, es decir, ángulo de partida θ0 y
velocidad angular inicial v0 ,
θ (0 ) = θ 0 (5.7)
θ ′ (0) = v0 (5.8)
entonces utilizando la fórmula (4.27), y dado que en este caso f = 0, la solución particular, que
determina la posición del péndulo en función del tiempo, es
g g
θ (t) = θ0 Eα,1 − tα + v0 tEα,2 − tα . (5.9)
L L
Por otra parte, en este sistema es habitual considerar nula la velocidad inicial de la partícula.
Esto es, en origen el péndulo se deja oscilar sin ejercer sobre él ningún impulso. En esta situación
la segunda condición inicial es
θ ′ (0) = 0 (5.10)
y la solución queda simplificada a
g g
θ (t) = θ0 Eα,1 − tα = θ0 Eα − tα . (5.11)
L L
En el caso particular α = 2, teniendo en cuenta que
k r 2k r
g ∞ − Lg t2 ∞
(−1)k g g
E2,1 − t2 = ∑ Γ (2k + 1) = ∑ (2k)! t = cos t
L k =0 k =0
L L
2 g
se obtiene la solución tradicional del péndulo linealizado ddt2θ (t) + L θ (t) = 0 con condiciones
iniciales (5.7)-(5.10): r
g
θ (t) = θ0 cos t
L
5.2. Proyectil
Otro problema dinámico clásico es el del proyectil . Se trata de un sistema en el que una
partícula es proyectada formando un determinado ángulo con la superficie de la Tierra. El pro-
yectil, bajo el efecto de la fuerza de la gravedad, describirá una curva de tipo parabólico. Esta
trayectoria está perfectamente determinada matemáticamente incluso para el caso, un poco más
complicado, en el que consideremos otras fuerzas intervinientes en el sistema además de la gra-
vedad, como son las de rozamiento. Áreas como la balística han estudiado en profundidad este
tema. Aquí, siguiendo trabajos como los de Ebaid [4] u Otero [12], vamos a estudiar el mo-
vimiento de un proyectil desde el enfoque fraccionario. Para ello, tal y como hicimos con el
péndulo, utilizaremos como hipótesis inicial una fórmula generalizada de la segunda Ley de
Newton.
Vamos a partir de unas hipótesis que simplifican el problema, de forma que la única fuerza
actuante sea la gravedad:
3 Esto es lo que se hace en el estudio de la dinámica del péndulo clásico, y es razonable para valores pequeños del
ángulo debido a que la solución constante cero es estable. En el caso fraccionario, la estabilidad de la solución cero,
aunque parece obvia, es todavía un problema abierto (ver [3]).
El medio no ofrece oposición al avance del proyectil ni por resistencia del aire ni por
ninguna otra fricción.
No hay curvatura de la superficie terrestre, que es plana y sin rugosidades.
El movimiento se traza dentro de un plano vertical fijo OXY.
La fuerza de la gravedad es uniforme, y no disminuye con la altura del proyectil.
No se tiene en cuenta la fuerza de Coriolis, debida al movimiento de rotación de la Tierra.
v0 sen φ
H
v0
φ x
v0 cos φ
mg
Supongamos que se lanza el proyectil, de masa m, con una velocidad inicial de módulo
v0 y formando un ángulo φ con la horizontal. Escogemos el plano OXY coincidiendo con el
plano de la trayectoria, de forma que el origen O se corresponda con la posición inicial de la
partícula. Podemos descomponer el movimiento en sus componentes horizontal y vertical, que
serán independientes una de la otra. Como la única fuerza considerada, la fuerza gravitatoria,
solo tiene componente vertical, la trayectoria será la composición de un movimiento rectilíneo
uniforme horizontal con un movimiento rectilíneo uniformemente acelerado vertical. Es decir,
si x (t) e y (t) son las respectivas componentes horizontal y vertical del movimiento, lo anterior
junto con la segunda Ley de Newton F = ma nos lleva a las ecuaciones diferenciales clásicas
del proyectil:
d2 x
m 2 (t) = 0 (5.12)
dt
d2 y
(t) = −mgm (5.13)
dt2
que con las condiciones iniciales comentadas,
g 2
y (t) = v0 t sen φ − t . (5.17)
2
Para aproximarnos al problema desde el punto de vista fraccionario sustituimos, tal y como
hicimos al estudiar el péndulo, la fórmula clásica de Newton
F = ma
por la generalización
C
F (t) = m D0α x (t) , (5.18)
con 1,5 < α ≤ 2. Así, aplicando el nuevo enfoque (5.18) a (5.12) y (5.13), y simplificando el factor
común m, obtenemos las EDFs del proyectil fraccionario:
C α
D0 x ( t ) = 0 (5.19)
C
D0α y (t) = − g (5.20)
siendo 1,5 < α ≤ 2.
5.2.2. Soluciones
Las ecuaciones (5.19) y (5.20) son sencillas EDFs lineales del tipo de las estudiadas en el
Método 1 (operadores fraccionarios) de la página 92. Si fijamos las condiciones iniciales (5.14)-
(5.15) y aplicamos la fórmula obtenida para este tipo de ecuaciones (4.18) llegamos a las solu-
ciones
x (t) = v0 t cos φ (5.21)
gtα
y (t) = v0 t sen φ − (5.22)
Γ (α + 1)
Rt
donde hemos utilizado que 0 (t − s)α−1 ds = tα/α y que αΓ (α) = Γ (α + 1). Estas son las ecua-
ciones paramétricas que rigen el movimiento del proyectil en el plano OXY partiendo de (5.18)
como principio físico. Como es natural, al fijar α = 2 se convierten en las soluciones tradiciona-
les (5.16)-(5.17).
Eliminando la variable tiempo y unificando las dos ecuaciones paramétricas se llega a la
ecuación algebraica que relaciona x e y, que, al contrario que en el caso α = 2, donde se trata de
una parábola, es un polinomio fraccionario de grado α:
− gx α
y= + x tan φ (5.23)
Γ (α + 1) (v0 cos φ)α
A continuación se estudian, para los casos fraccionario y clásico, las cuatro magnitudes más
relevantes en el análisis del movimiento de proyectiles: alcance, ángulo óptimo, altura máxima
y tiempo de vuelo. Cada uno de estos valores dependerá de ciertos parámetros, e indicare-
mos estas relaciones con notación de dependencia funcional. Todos los resultados parten de la
hipótesis de un lanzamiento de proyectil (5.23) con condiciones iniciales (5.14)-(5.15). Las de-
mostraciones de los resultados que siguen pueden encontrarse en Ebaid [4].
5.2.3. Alcance
Definición 66. Se denomina alcance ( L) a la distancia horizontal recorrida por el proyectil entre
el instante del lanzamiento y la toma de contacto con la superficie. Es decir, al valor x (t) en el
tiempo t de impacto con la superficie.
7.5 1.5
2
2.5
3
2.5
0
0 2 4 6 8 10 12 14
x
Figura 3. Comparación de trayectorias para diferentes valores de α (φ = π/4, v0 = 10)
100 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Cálculo Fraccionario y dinámica newtoniana Antón Lombardero Ozores
Proposición 71. Fijados α y v0 , el alcance máximo (que viene determinado por el ángulo óptimo) es
r 1
α − 1 Γ (α + 1) v0α α−1
Lmax (α, v0 ) = √
α g α
v20
Lmax (2, v0 ) = .
g
Definición 72. La altura máxima ( H ) es el valor de y (t) en el punto más elevado de la trayectoria
del proyectil.
v20 sen2 φ
H (2, φ, v0 ) = .
2g
Proposición 74. Fijados α y v0 , las alturas máximas fraccionaria y clásica coinciden solamente para un
ángulo de lanzamiento concreto
h i 1
2 α −1
2− α
g 2− α Γ (α ) π
H (α, φ, v0 ) = H (2, φ, v0 ) ⇐⇒ φ = arc sen 0<φ<
v0 2
Definición 75. El tiempo de vuelo ( T ) es el tiempo total en que el proyectil está en movimiento,
desde su lanzamiento hasta que toma tierra.
2v0 sen φ
T (2, φ, v0 ) =
g
Ciertos dispositivos elásticos, como los muelles o las tiras de goma, tienen la particularidad
de recuperar su longitud inicial después de ser estirados. Este comportamiento elástico se debe
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 101
Antón Lombardero Ozores Investigación
a la acción de una fuerza recuperadora que se opone al estiramiento. Supongamos que dispo-
nemos de un muelle que se encuentra suspendido de un extremo, y cuya longitud es l0 . Al
colgar del otro extremo un cuerpo de masa m, el muelle se estira hasta alcanzar la longitud l1 ,
quedando entonces inmóvil. En ese momento el sistema está en lo que se denomina la posición
de equilibrio, y esa posición del cuerpo se toma como origen coordenado de un sistema unidi-
mensional. Denotemos por x (t) la distancia del cuerpo a la posición de equilibrio en el instante
t, considerando x (t) > 0 cuando el muelle está estirado (el cuerpo se mueve hacia abajo) y
x (t) < 0 cuando está comprimido (el cuerpo se mueve hacia arriba). En este caso sí que vamos
a tener en cuenta la fuerza de rozamiento ejercida por el medio (aire) sobre la masa suspendida,
pero fijamos las siguientes hipótesis de partida:
Vamos a analizar qué fuerzas intervienen en el sistema. De acuerdo con la Ley de Hooke,
la fuerza recuperadora FH que un resorte ejerce sobre la masa es directamente proporcional
a la distancia en la que el resorte se ha estirado o comprimido. Puesto que ésta es igual al
desplazamiento de la masa m respecto de la posición de equilibrio, se deduce que
La constante de proporcionalidad k > 0 se llama constante elástica del muelle, y depende de sus
características físicas. Obsérvese que FH y x tienen signos opuestos.
Además, un cuerpo que se mueve en el aire está sometido a fricciones que se traducen en la
aparición de una fuerza que tiende a frenarlo. Suele admitirse que el rozamiento del aire, a ve-
locidades no excesivamente altas, influye sobre los cuerpos oponiendo una fuerza proporcional
y de sentido contrario a la velocidad del cuerpo, es decir
con 1,5 < α ≤ 2. De esta manera obtenemos la nueva EDF que gobierna la dinámica del muelle:
m C D0α x (t) = FH (t) + FR (t) ⇒
m C
D0α x (t) + cx ′ (t) + kx (t) = 0 (5.27)
102 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Cálculo Fraccionario y dinámica newtoniana Antón Lombardero Ozores
a la que podemos añadir unas condiciones iniciales que determinen la posición inicial de la
masa,
x (0) = b0 , (5.28)
x ′ (0) = b1 . (5.29)
5.3.2. Soluciones
∞
(−k/m )n αn ∞ Γ (n + j + 1) (−c/m) j t(α−1)j
x1 ( t ) = ∑ n!
t ∑
Γ (αn + 1 + αj − j) j!
+
n =0 j =0
c ∞
(−k/m)n αn+α−1 ∞ Γ (n + j + 1) (−c/m ) j t(α−1) j
+ ∑ n! t ∑ Γ (αn + α + αj − j) (5.30)
m n =0 j =0
j!
∞
(−k/m )n αn+1 ∞ Γ (n + j + 1) (−c/m ) j t(α−1)j
x2 ( t ) = ∑ n! t ∑ Γ (αn + 2 + αj − j) j!
+
n =0 j =0
c ∞
(−k/m )n αn+α ∞ Γ (n + j + 1) (−c/m) j t(α−1)j
+ ∑ t ∑ (5.31)
m n =0 n! j =0
Γ (αn + 1 + α + αj − j) j!
No hay garantía de que estas dos soluciones formen un sistema fundamental de soluciones
de la ecuación. El encontrarlas se convierte en una cuestión bastante complicada. Por ejemplo,
si α es un racional p/q con p < q, necesitaríamos p soluciones para configurar un sistema funda-
mental. Para un α real no racional, el proceso se complica más.
6. Conclusiones
De la misma forma que las propiedades de no localidad y memoria de los operadores frac-
cionarios han demostrado ser de utilidad en áreas como la mecánica de fluidos, la teoría de
viscoelasticidad, el procesado de imagen y otras muchas, es razonable conjeturar que en sis-
temas físicos de la mecánica clásica pueda, bajo ciertas circunstancias, ocurrir algo parecido.
Un siguiente paso en la investigación consistiría en determinar, a nivel práctico, si existen con-
diciones (de fricción del medio, climatológicas, gravitatorias, etc.) para las que los modelos
fraccionarios expuestos representen la realidad de forma más precisa que los modelos clásicos.
En un plano más teórico, el Cálculo Fraccionario, como rama matemática reciente que es,
proporciona innumerables líneas futuras de investigación. En relación a los aspectos tratados
en este artículo señalamos el problema abierto de determinar la estabilidad de la solución trivial
0 en el péndulo fraccionario, y sobre todo la búsqueda de la solución general para sencillas
ecuaciones diferenciales fraccionarias lineales, como la obtenida en el problema del resorte.
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 103
Antón Lombardero Ozores Investigación
Referencias
[1] A GRAWAL, O. P., A new lagrangian and a new Lagrange Equation of Motion for fractionally dam-
ped systems, Journal of Applied Mechanics, vol. 68, (2001) 339-341.
[2] B ALEANU, D., A SAD, J. H., P ETRAS, I., Fractional-order two-electric pendulum, Romanian Re-
ports in Physics , vol. 64(4), (2012) 907-914.
[3] C HEN, F., N IETO, J. J., Z HOU, Y., Global attractivity for nonlinear fractional differential equations,
Nonlinear Analysis: Real World Applications, 13, (2012) 287-298.
[4] E BAID, A., Analysis of projectile motion in view of fractional calculus, Applied Mathematical
Modelling, 35, (2011) 1231-1239.
[5] G ORENFLO, R., M AINARDI, F., Fractional Calculus: Integral and differential equations of fractional
order, CISM courses and lectures, vol. 378, (1997) 223-276.
[6] H AUBOLD, H. J., M ATHAI, A. M., S AXENA, R. K., Mittag-Leffler functions and their applications,
Journal of Applied Mathematics, vol. 2011, (2009).
[7] H EYMANS, N., P ODLUBNY, I., Physical interpretation of initial conditions for fractional differential
equations with Riemann-Liouville fractional derivatives, Rheologica Acta, vol. 45(5), (2006).
[8] K ILBAS, A. A., Some aspects of differential equations of fractional order, Rev. R. Acad. Cienc.
Exact. Fís. Nat., vol. 98(1), (2004) 27-38.
[9] K ILBAS, A. A., S RIVASTAVA, H. M., T RUJILLO, J., Theory and Applications of Fractional Diffe-
rential Equations, Elsevier, Amsterdam, 2006.
[10] M ILLER, K. S., R OSS, B., An Introduction to the Fractional Calculus and Fractional Differential
Equations, Wiley and Sons, New York, 1993.
[11] O LDHAM, K. B., S PANIER, J., The Fractional Calculus, Academic Press, New York, 1974.
[12] O TERO, O., Ecuaciones diferenciales de orden fraccionario (Trabajo de Fin de Máster), Tutor: Juan
José Nieto Roig, Universidad de Santiago de Compostela, Facultad de Matemáticas (2012).
[13] P IERANTOZZI, T., Estudio de generalizaciones fraccionarias de las ecuaciones estándar de difusión
y de ondas (Tesis doctoral), Tutor: Vázquez, L., Universidad Complutense de Madrid, Facultad
de Matemáticas (2006).
[14] P ODLUBNY, I., Fractional Differential Equations, Academic Press, San Diego, 1999.
[15] P ODLUBNY, I., Geometric and physical interpretation of fractional integration and fractional dif-
ferentiation, Fract. Cal. Appl. Anal, 5(4), (2002) 367- 386.
[16] S AMKO, S. G., K ILBAS, A. A., M ARICHEV, O. I., Fractional Integrals and Derivatives: Theory
and Applications, Gordon and Breach Science Publishers, Suiza, 1993.
[17] S HIMA, H., How far can Tarzan jump?, European Journal of Physics, vol. 33(6) (2012) 1687-
1696.
[18] T ENREIRO M ACHADO, J., K IRYAKOVA, V., M AINARDI, F., Recent history of fractional calculus,
Commun Nonlinear Sci Numer Simulat, vol. 16(3), (2011) 1140-1153.
[19] VÁZQUEZ M ARTÍNEZ, L., Una panorámica del cálculo fraccionario y sus aplicaciones, Rev. R.
Acad. Cienc. Exact. Fís. Nat., vol. 98(1), (2004) 17-25.
104 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Cálculo Fraccionario y dinámica newtoniana Antón Lombardero Ozores
[20] V ELASCO C EBRIÁN, M. P., Modelos diferenciales y funciones especiales en el ámbito del cálculo
fraccionario (Trabajo de Fin de Máster), Tutores: Trujillo, J. J. (Universidad de La Laguna) y
Vázquez, L. (Universidad Complutense de Madrid) (2008).
Sobre el autor:
Nombre: Antón Lombardero Ozores
Correo electrónico: [email protected]
Institución: Departamento de Análise Matemática, Universidade de Santiago de Compostela,
España.
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 105
Investigación
Programación en paralelo GPGPU del método en
diferencias finitas FTCS para la ecuación del calor
FTCS finite difference scheme GPGPU parallel
computing for the heat conduction equation
Vicente Cuellar Moro, Miguel Martín Stickle
y Manuel Pastor Pérez
Revista de Investigación
1 de abril de 2014
Resumen
En el presente artículo se muestran las ventajas de la programación en paralelo
resolviendo numéricamente la ecuación del calor en dos dimensiones a través del método
de diferencias finitas explícito centrado en el espacio FTCS. De las conclusiones de este
trabajo se pone de manifiesto la importancia de la programación en paralelo para tratar
problemas grandes, en los que se requiere un elevado número de cálculos, para los cuales la
programación secuencial resulta impracticable por el elevado tiempo de ejecución. En la
primera sección se describe brevemente los conceptos básicos de programación en paralelo.
Seguidamente se resume el método de diferencias finitas explícito centrado en el espacio
FTCS aplicado a la ecuación parabólica del calor. Seguidamente se describe el problema de
condiciones de contorno y valores iniciales específico al que se va a aplicar el método de
diferencias finitas FTCS, proporcionando pseudocódigos de una implementación secuencial
y dos implementaciones en paralelo. Finalmente tras la discusión de los resultados se
presentan algunas conclusiones.
Palabras Clave: Computación en paralelo, GPU, Diferencias finitas, ecuación del calor.
Abstract
In this paper the advantages of parallel computing are shown by solving the heat
conduction equation in two dimensions with the forward in time central in space (FTCS)
finite difference method. Two different levels of parallelization are consider and compared
with traditional serial procedures. We show in this work the importance of parallel
computing when dealing with large problems that are impractical or impossible to solve
them with a serial computing procedure. In the first section a summary of parallel
computing approach is presented. Subsequently, the forward in time central in space
(FTCS) finite difference method for the heat conduction equation is outline, describing how
the heat flow equation is derived in two dimensions and the particularities of the finite
107
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
difference numerical technique considered. Then, a specific initial boundary value problem
is solved by the FTCS finite difference method and serial and parallel pseudo codes are
provided. Finally after results are discussed some conclusions are presented.
1.1. Introduction
Traditionally, computer software has been written for serial computation. To solve a
problem, an algorithm is constructed and implemented as a serial stream of instructions.
These instructions are executed on a central processing unit (CPU) on one computer. Only one
instruction may execute at a time. After that instruction is finished, the next is executed.
Parallel computing, on the other hand, uses multiple processing elements simultaneously
to solve a problem. This is accomplished by breaking the problem into independent parts so
that each processing element can execute its part of the algorithm simultaneously with the
others. The processing elements can be diverse and include resources such as a single
computer with multiple processors, several networked computers, specialized hardware, or
any combination of the above.
Although the diversity of parallel machines is too big to describe them here, they can be
classified into two big groups. On one side we have the distributed memory systems, on the
other the shared memory systems. The distributed memory systems are those in which every
computing unit has a private memory which only he can access directly. Shared memory
systems are those in which each computing node can access all the memory of the system. It is
important to note that these systems are sometimes combined as the computing unit in a
distributed system is usually a parallel machine on its own, whether a multicore processor or
any other type, which in itself is a system of the shared memory kind.
This distinction between models is critical as it describes how each system handles the
fundamental problem of parallel computing, which is the data locality problem. Except in the
simplest algorithms, often called “embarrassingly parallel problems” parallel machines have
to deal with the access to information between its computing units, as some units need
information calculated by different units.
Shared Memory Systems can have simplified solution compare to the distributed ones, as
the whole system memory is available to every unit it is just a matter of synchronization and
coherence rather than communication. Shared memory systems need to care only about
avoiding race conditions over the data, they have to make sure the calculations of unit are
completed before other units access that information. The main limitation of the Shared
Memory Systems is the poor scaling, to go beyond certain sizes they need to be arranged as
nodes in a Distributed Memory Model.
Distributed Memory Systems have additional needs, they need to establish and manage
the communication of the information across the different nodes to solve the problem. This
means additional complexity and requires a careful design not only of the system itself but
also in the algorithms to avoid the communications from becoming a bottleneck.
108 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
FTCS finite difference scheme Vicente Cuellar Moro,
GPGPU parallel computing for the heat conduction equation Miguel Martín Stickle y Manuel Pastor Pérez
The GPU is in itself a Shared Memory parallel machine. In general a multi-GPU system
would constitute a Distributed Memory System, although it must be indicated that latest
developments in some of the GPU languages as CUDA can make a multi-GPU system appear
to the programmer as a Shared memory one at the cost of some performance loss.
Nowadays a GPU (graphics processing unit) can be used together with a CPU to
accelerate general-purpose scientific and engineering applications. Among the different
existing procedure to perform parallel computing, CPU + GPU is a powerful combination
because CPUs consist of a few cores optimized for serial processing, while GPUs consist of
thousands of smaller, more efficient cores designed for parallel performance. In this section
the principal aspects concern the General Purpose Graphic Processing Unit (GPGPU) for
parallel computing are outline.
Simply known as GPU Computing or GPGPU (General Purpose Graphic Processing Unit)
is a relatively new field in the parallel programming landscape that uses graphic processors
for its calculations.
The usual applications of GPGPU include all sort of physical simulations, image and video
processing, biological and chemical analysis, applications that require many linear algebra
operations, MRI image reconstruction, and all sorts of financial calculations.
The closer model of computers that existed so far was vector computers, although the
differences between vectors computers and GPUs are substantial, the algorithms designed for
vector computers can be easily translated into a GPU version.
The first apparent big difference between GPU computing and “traditional” parallel
computing is the huge amount of execution threads in GPU computing. While conventional
parallel computing rarely goes beyond few hundred threads, GPU computing goes several
orders of magnitude beyond, in some cases launching kernels with millions of threads,
although in reality only few thousands will run concurrently at any given time.
The development of programming languages for GPU computing is very recent. In 2001
Nvidia Corporation develops the first GPU that can easily be programmed. Shortly after that
Thompson et al (2002) write the first academic paper exploring the use of GPU as vector
computers for arithmetic floating point calculations. Ian Buck et al (2004)[5] created Brook, the
first GPGPU language, as an extension of C. Shortly after that the first proprietary language
emerges, CUDA in 2006. OpenCL the main open language is released in 2008 by the Khronos
Group.
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 109
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
CUDA is probably the most advanced and developed language, currently in its 5.0 version
with a pretty solid 5.5 version as release candidate. CUDA development cycle brings a new
major version every year. It is currently accessible as language extensions to C/C++ (Nvidia)
and Fortran (Portland Group), as well as the compiler-directive language OpenACC.
Nevertheless there are bindings and tools to connect CUDA with almost any language.
The latest architecture for GPU computing with CUDA is the Kepler one, it is a 28nm
process, the most powerful card on said architecture is the Tesla K20X, having 2688 stream
processors (CUDA cores) grouped in 14 SMX (Streaming multiprocessors). Even when we are
running our calculations on a modest Geforce GT640 with only 2 SMX (384 sps) the algorithm
provides substantial speedups compared to the serial version.
The CUDA GPUs are able to work with three different levels of parallelism. At the bottom
level we find the execution thread, which will be executed by a Streaming Processor (SP).
Above the single thread we have the execution block, which consists of a group of threads
whose number go from one to tens of thousands (although rarely is less than 32 threads), a
Block will be executed by the SPs grouped in a Streaming Multiprocessor (SMX), in the latest
architecture the SMX holds 192 SP, every SMX is capable of running several execution blocks
concurrently if it has enough hardware resources to do so. Finally there is a third level of
parallelism, although is not frequently used, which is the streams, each stream can schedule
many GPU functions (kernels) even memory transfer operations to and from the CPU, which
depending on the hardware resources may run concurrently, streams are especially useful in
multi-GPU configurations usually in combination with MPI (message passing interface).
Compared to other GPGPU languages CUDA provides both simplicity and depth, it can
provide a great control of the hardware without requiring too many lines of code as opposed
to OpenCL. CUDA enjoys a strong support from Nvidia, in form of many libraries, very
extensive documentation and code examples available for free, despite CUDA proprietary
nature. In the present article we focus only in the first two levels of paralelism, leaving for the
third level a mere indication of how to scale the problem for multi-GPU environments.
2.1. Introduction
There are many physical phenomena that are able to be properly represented by partial
differential equations (PDEs). For instant heat flow through metallic bar, evolution of water
wave motion or the static deformation of an elastic bar.
In general, a PDE in several spatial variables x, y,z , and time t is an equation G involving
an unknown function u = u x, y,z,t , sometimes called the state variable, and some of its
partial derivatives
110 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
FTCS finite difference scheme Vicente Cuellar Moro,
GPGPU parallel computing for the heat conduction equation Miguel Martín Stickle y Manuel Pastor Pérez
where I is a given interval of time, often define as positive time t 0 , and is a region
in the space. A PDE model or an initial boundary value problem is a PDE supplemented with
initial and boundary conditions that specify the state initially and on the boundary.
By a solution of the PDE (1) it is meant a function u = u x, y,z,t defined on the space-time
domain x, y,z , t I that satisfies, upon substitution, the equation (1) identically on the
domain.
In most practical cases coming from engineering practice PDEs are almost always solved
numerically on a computer because most real world problems are too complicate to be solved
analytically. By a numerical method we generally mean that the continuous PDE problem is
replaced by a discrete problem that can be solved on a computer in finitely many steps. The
result is a discrete solution known approximately at only finitely many points.
In this section we focus attention on the forward in time central in space (FTCS) finite
difference method for the heat conduction equation. First we will set up a PDE model for the
heat flow in three and two dimensions. As the heat equation describes a diffusive-like
behaviour the finite difference procedure to solve the diffusion equation is outline
subsequently.
Many PDE models come from a basic balance law or conservation law. A conservation
law, whose foundation goes back to Lavoisier’s principal [2], is just a mathematical
formulation of the basic fact that the rate at which a quantity changes in a given domain must
equal the rate at which the quantity flows across the boundary plus the rate at which the
quantity is created or destroy within the domain.
To set up a PDE model for the heat flow in three dimensions let be a region in space
where heat is flowing, and let u = u x, y,z,t be the temperature at time t at the point x, y,z
in . We assume that the region is homogeneous and is characterized by a constant specific
heat c and a constant density . Now let B be an arbitrary sphere contained in (Figure 1).
We shall apply an energy balance principal to B that requires [1] that rate of change of the
total amount of heat energy in B must equal the rate at which heat flows into B across its
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 111
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
boundary plus the rate at which heat energy is generated by any sources in B. The total
amount of heat in a small volume element dV = dxdydz is c udV , and thus the total amount
of heat energy in B is given by three-dimensional integral
c udV
B
(2)
We assume that the heat sources (or sinks) are given by a point function f = f x, y,z,t ,
where fdV is the rate at which heat is generated in dV ; thus, the rate at which heat is
generated is the whole of B is
fdV
B
(3)
Next, we introduce the heat flux vector x, y,z,t ; its direction corresponds to the
direction of heat flow at position x, y,z at time t . The rate at which heat flows across a
surface element d A oriented by a unit vector n is n dA (Figure 2). Consequently, the net
rate that heat energy flows across the boundary of B, denoted by B , is the surface integral
n dA
B
(4)
So, the conservation law or equivalently the energy balance law [1] is
d
dt B
c udV n dA fdV (5)
B B
Now making use of the divergence theorem, under sufficient differentiability requirement
on the vector field , we can write (5) as
Bringing the time derivative under the integral on the left side and finally rearrange all the
terms under one volume integral we get
c u div f dV 0
B
t
(7)
As this balance law must hold for every sphere B in , the integrand must vanish [1],
giving the partial differential equation
c ut div f (8)
For all t and all x, y,z . The equation (8) is the local form of the conservation law. It
still contains two unknowns, the scalar temperature u and the vector heat flux . A
constitutive relation can be postulate to connect the two. One such relation is the three
dimensional version of Fick’s law, or Fourier’s heat conduction law, which states, consistent
with the laws of thermodynamics, that heat flows down the gradient, or in symbols,
K grad u K ux ,uy ,uz (9)
112 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
FTCS finite difference scheme Vicente Cuellar Moro,
GPGPU parallel computing for the heat conduction equation Miguel Martín Stickle y Manuel Pastor Pérez
where the proportionality constant K is the thermal conductivity. Substituting (9) in (8)
and using the identity div grad u uxx uyy uzz we obtain the heat equation for the
temperature u = u x, y,z,t in three dimensions:
c ut K uxx uyy uzz f (10)
1
ut u f (11)
c
ut u (12)
Of course, these calculations can be made in two dimensions in order to describe heat flow
in the plane. The two-dimensional analogue of (11) has the same expression where now
u = u x, y,t , f = f x, y,t , and u u yy u zz .
Figure 2. Heat flux through a surface element dA oriented by its unit normal vector n
The auxiliary conditions to complete the PDE model for heat conduction includes an
initial temperature condition
and conditions on the boundary. There are two main types of boundary conditions that
often occur in heat conduction problems. We can prescribe the temperature on the boundary
of , that is,
where g is given; this condition is called Dirichlet condition. We could also specify the
heat flux on the boundary
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 113
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
In order to illustrate the forward in time central in space (FTCS) finite difference method
we applied it to the following one dimension diffusion equation problem
ut Duxx , 0 x 1, t 0 (16)
The general idea of the finite difference method is to replace the partial derivatives in the
equation by their difference-quotient approximations and then let the computer solve the
resulting difference equation. The forward in time central in space scheme is a specific choice
for the difference-quotient approximation [2,3].
The first step consist in discretize the region of space-time where we want to obtain a
solution. In this case the region is 0 x 1, 0 t T . We have put a bound on time because in
engineering practice we only solve a problem up until a finite time. Discretizing means
defining a lattice of points, also known as nodes, in this space-time region by
xj jh, tn nk, j 0,1,..., J ; n =0,1,..., N (19)
where the fix numbers h and k are the spatial and temporal step sizes, respectively.
Here, h 1 J while k T N . The integer J is the number of subintervals in 0 x 1 and the
integer N is the number of time steps to be taken (Figure 3). At each node x ,t
j n of the
lattice we seek an approximation, which we call U , to the exact value n
j u x ,t
j n of the
solution. Note that the superscript n refers to time, and the subscript j refers to space. We
regard U nj as two dimension array, or matrix, where n is a row index and j is a column
index.
0
0
Figure 3. Discrete lattice and nodes involved in the FTCS finite difference method
114 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
FTCS finite difference scheme Vicente Cuellar Moro,
GPGPU parallel computing for the heat conduction equation Miguel Martín Stickle y Manuel Pastor Pérez
To obtain equations for the U nj , we replace the partial derivatives in the PDE by their
difference approximation. From calculus, using the Taylor series expansion in t , we might get
[1, 2] the difference quotient
u x j ,tn 1 u x j ,tn
ut x j ,tn k
(20)
u x j 1 ,tn 2 u x j ,tn u x j 1 ,tn
uxx x j ,tn h 2
(21)
The first of these formulas is the usual forward in time difference approximation for a first
derivative, and the second is central in space difference approximation for a second
derivative. Making use of (20) and (21), the PDE (16) at the point x j ,tn is replace by the
difference equation
U nj 1 U nj U nj 1 2U nj U nj 1
D (22)
k h2
U nj 1 U nj
Dk n
h2
U j 1 2U nj U nj 1 (23)
The difference equation (23) gives the approximate solution at the node x j ,tn 1 in terms
of approximations at the three earlier nodes x j 1
,tn , x j ,tn , x j 1 ,tn , see Figure 3. These
four points geometric arrangement form the stencil for the FTCS finite difference method in
one dimension. Now we see how to fill up the lattice approximate values. We know the values
at t 0 from the initial condition (18). That is
U 0j = f x j , j 0,1,...,J (24)
We can observe that if j 0 then U n1 appears in equation (23) and for j J then U Jn1
also enters in equation (23). However, following (19) there are no nodes place at j 1
neither at j J 1 . To solve this drawback [3] we consider the boundary conditions (17),
namely we only consider equations in (23) with subscript j running from 1 to J 1 while
U 0n = 0 , U Jn = 0, n 1, 2,..., N (25)
The difference formula (23) with j 1,...,J 1 can now be applied at all the interior lattice
points, beginning with the values at the t 0 level, to compute the values at the t t1 level,
then using those to compute the values at the t t2 level, and so on. Therefore, using the
difference equation (23), we can march forward in time, continually updating the previous
temperature profile.
A finite difference scheme like (23) is called explicit scheme because it permits the explicit
calculation of an approximation at next time step in terms of values at a previous time step.
Because of the error, called truncation error, that is present in (23) due to replacing derivatives
by differences, the scheme is not always accurate. If the time step k is too large, then highly
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 115
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
inaccurate results will be obtained. It can be shown [3] that for a one dimension problem we
must fulfil the stability condition
kD 1
(26)
h2 2
for the scheme to work; this condition limits the time step. This is why the forward
scheme (23) is known to be conditionally stable. To obtain a method that is unconditionally
stable we should consider an implicit difference method. However, these schemes involved
the solution of a linear system to advance in time. The solution of algebraic linear systems
scales poorly with the problem size, usually as an O(n3) problem, therefore in many cases it is
often avoided when the objective is performance on large problems. From now on we will
focus on the forward in time central in space scheme (23).
ut
1
u u
16 yy zz
(27)
Where is the unit square in xy-plane. By the second section of the present work, this
PDE model might be regarded as a heat conduction problem in two dimensions where the
diffusivity has been set to 1 16 with no heat source. Homogeneous Dirichlet boundary
condition has been considered while initial temperature is given by(29).
After applying separation of variables method [4] to the initial boundary value problem
define by (27)-(29) we get the following analytical solution
The same PDE model is numerically solved by the forward in time central in space (FTCS)
finite difference method up until T 1 s . To make a significant comparison between serial and
different parallel computing implementations, the number J of subintervals for the spatial
discretization, equal for x and y directions, varies in powers of two from 32 to 512 . The
number of time steps is set to N 100000 . Regarding the stability condition for two dimension
k 1
case, being h 1 J and k T N , we get 2 . Therefore stable numerical simulations are
h 4
performed for each value of J [3].
The FTCS finite difference method in two dimensions is associated with a five-point
stencil. Therefore, calculation of every grid node value at time tn 1 requires reading the values
of its four neighbour nodes at time tn , xi , y j 1 ,tn , xi , y j 1 ,tn , xi 1 , y j ,tn , xi 1 , y j ,tn and the
value at time tn of the node itself xi , y j ,tn , as it is shown in Figure 4.
116 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
FTCS finite difference scheme Vicente Cuellar Moro,
GPGPU parallel computing for the heat conduction equation Miguel Martín Stickle y Manuel Pastor Pérez
Bearing in mind this five-point stencil structure, a serial pseudo code solution would
contain code like:
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 117
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
For a naive GPU implementation the pseudo code solution would contain code like:
3: calculate node value U for xi , y j at time step tn 1
4: write result into memory.
Another kernel is launched with one thread per boundary node
5: apply boundary conditions at time step tn 1
As we can see the most aparent change in the pseudo code comes from the fact that
instead of having a nested loop to traverse the grid in the x and y directions, in the parallel
version we simply launch a kernel (GPU subroutine) with as many threads as grid nodes.
Although not all the threads run concurrently, high level of paralelism is achieved. Moreover
the GPU will schedule and execute the threads in a manner that while some threads are
waiting for the values they need from memory, other threads are being executed, this is called
“latency hidding” wich is another reasons GPUs can run so fast calculations. Effective latency
hidding requires launching more threads than Streaming Processors are available at the GPU.
The particular tunning of the parameters like threads per block and total number of threads
must be adapted for every particular algorithm.
We can improve over the naive version when we consider the second level of parallelism
as we can use the on-chip shared memory that every SMX has. In the second version we start
using better the GPU hardware, in particular as we are working at the SMX level we can
employ the on-chip shared memory.
Shared memory presents one of the best opportunities to obtain a large increase in
performance if the algorithm can work within the restriction imposed by its limited size.
Fetching values from shared memory is much faster than reading global memory. With
shared memory we have a latency of around 8 cycles, while uncached global memory has a
latency of up to 600 cycles. That is almost two orders of magnitude.
The main obstacle in the use of shared memory is its limited size. In the best of cases we
have available around 48 Kb of shared memory per SMX which means around 6.000 values in
double precision or 12.000 in single precision.
The benefit of using shared memory is better understood if we consider that in the case of
this stencil operation we can reduce the memory bandwidth in use to a fifth of the naive
implementation.
Any GPU subroutine can be classified from a performance point in two rough groups.
“Memory Bandwidth Bound” or simply Memory-Bound are algorithms whose performance is
limited by the latencies of memory Access, that is the number of floating point operations they
execute cannot compensate for the time each thread has to wait to fetch the values is using, in
other words a faster clock rate in the GPU would not yield the result any faster.
118 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
FTCS finite difference scheme Vicente Cuellar Moro,
GPGPU parallel computing for the heat conduction equation Miguel Martín Stickle y Manuel Pastor Pérez
“Computationally Bound” kernels in the other hand, have enough computations so the clock
rate in the SP is the limiting factor, in this case increasing the clock frequency would result in
the increase of performance. A given implementation can be either memory or
computationally bound, as the computation threshold is usually several times higher that
would be ideally the limit we would like to reach.
Most implementations, specially the naive ones, are “Memory-Bound”. The main reason
for that is how the hardware is built. As the computational power of a chip depends of the
number of transistors in the chip that means is directly proportional to the “effective
computational area” of the chip, while the bandwidth is proportional to the perimeter of the
chip. This may change in the future with the advent of 3D-stack chips, but for the time being is
the main challenge to current chip architectures.
So our second implementation tries to reduce the amount of times we need to fetch values
from global memory, as it can be seen that we are reading five times every value from global
memory, one for the calculation of the node himself and another four for the calculation of its
neighbours.
The new implementation is now divided in two fundamental phases, the first one consist
of loading the node value into shared memory, the second one consists in reading the
necessary values from shared memory to perform the calculation.
There are two very important aspects to consider in this implementation. First given the
limitation of shared memory size we have to divide the grid into “tiles” that fit shared
memory, this also requires some overlap between the tiles so all the threads have access to all
the necessary values. This introduces a small redundancy as the “halo” threads are only used
for memory loads and not for calculation.
The second aspects is a mere practical one, we need to introduce a synchronization barrier
at block level. We need to make sure that all the threads belonging to the same block have
finished loading their node values into shared memory before we make the calculation.
Figure 5. Simplified tile division of the grid for each SMX. Calculation nodes are in blue while “halo” nodes are yellow.
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 119
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
t=0s t = 0.25 s
t = 0.5 s t=1s
120 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
FTCS finite difference scheme Vicente Cuellar Moro,
GPGPU parallel computing for the heat conduction equation Miguel Martín Stickle y Manuel Pastor Pérez
In order to verify that numerical solutions are reproducing adequately the heat
conduction phenomena, the approximate solution achieve with FTCS scheme is compared
with exact solution (30). To this end the error profile for t 0, 0.25, 0.5, 1 s , measure as the
absolute value of the difference of the approximate and the exact solution is displayed in
Figure 7 for J 32 and N 1000
In this last figure it is clear that numerical and analytical solutions are almost coincident.
Obviously at initial time no difference exists. As the time integration proceeds, some error
appears in the numerical solution, being higher at nodes (0.25 , 0.25), (0.25 , 0.75), (0.75 , 0.25),
(0.75 , 0.75), which are the nodes where maximum and minimum temperature is attained.
Once it has been shown numerical solutions reproduce correctly the heat conduction
problem in 2D, running time comparisons for the different implementations are analyzed
subsequently.
We decided to run the parallel version in the last three GPU architectures developed by
Nvidia. These are, Tesla, Fermi and Kepler which is the latest one. Although we decide to
report only the two most diverse cards the Quadro FX3800 and the Geforce GT640, as the
results for the GT420 card didn’t add any interesting insights that were not provided by the
other two (Table 1)
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 121
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
Memory Power
Arquitecture Clock Band. Draw
Processor Name Chipset Process Mhz Gb/sec SMX SP (W)
Intel Xeon Bloomfield W3520 45 nm 2660 130
Quadro
FX3800 Tesla GT200 55 nm 600 51.2 24 192 108
Geforce
GT640 Kepler GK107 28 nm 728 28.5 2 384 50
We also provide the speedup per streaming multiprocessor to account the difference
between each card, as the Tesla one was a high end card at the time of its release (Quadro
FX3800) while the other two are low-end cards (Geforce GT420 and GT640).
In Figure 8 we can observed the running time in milliseconds for different number of
subintervals for the spatial discretization and the different implementations considered. It is
clear that an important running time improvement is achieved when parallel instead of serial
implementation is considered.
5
x 10
9
Intel Xeon - Serial
8 Quadro FX 3800 - Naive parallelization
Quadro FX 3800 - Improved parallelization
7 Geforce GT640 - Naive parallelization
Geforce GT640 - Improved parallelization
Running time [ms]
0
32 64 128 256 512
Number of subintervals for the spatial discretization
A different way to cluster data of Figure 8 is to consider the relative speed up between
serial and parallel computations. This is shown in Table 2. Again, it is clear the important
speed up achieve when parallel instead of serial implementation is considered.
122 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
FTCS finite difference scheme Vicente Cuellar Moro,
GPGPU parallel computing for the heat conduction equation Miguel Martín Stickle y Manuel Pastor Pérez
From Table 2 it can be observed how time scales with the grid size. In principle this is an
O(n2) problem, meaning the execution time (per time step) should scale roughly with the
square of the grid size (spatial steps), so a 64 x 64 grid should take four times as much as the
32 x 32 grid. This is observed in the serial version of the algorithm, with some irregularities
maybe caused by cache sizes.
In the GPU implementations we see a better scaling up to sizes of 512 x 512. This can be
observed in Table 3. In the case of the smaller grids it seems to scale almost free of cost up to
128 x 128. The reason for this advantage has to do with the way GPUs operate. To hide
latencies properly the GPU needs grids with a lot of blocks and threads. In the small problems
we only launch few blocks as the block size is 256 treads organized as 16 x 16 tiles. That means
the 32 x 32 problem only gets four blocks in the naive implementation. As the size of the
problem continues to increase we use all of the scheduling advantage of the GPU and
execution times and problem size scale at the expected rate.
Table 3. Scaling of execution times per time step compared to the previous problem size.
If we analyze the results only considering parallel implementation, we observed that the
improved version performs slightly worse than the naive implementation for the Fermi and
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 123
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
Kepler architectures (Figure 9). The reason is very simple. In this algorithm the use of shared
memory is limited while Kepler and Fermi include caches which in practice is the same
memory used as shared. So the naive implementation is in practice using the same memory
access but launching a bit less blocks and threads. Additionally the naive implementation
suffers no delay from the synchonization barriers.
That is not to say the improve algorithm is worse, whith a different memory access
pattern (a differente type of stencil) or being the problem a 3D instead of a 2D, we would see
much better results. If shared memory had to hold intermediate values like in reductions
algorithms we would see greater performance improvements using the shared memory
4
x 10
4
Quadro FX 3800 - Naive parallelization
3.5 Quadro FX 3800 - Improved parallelization
Geforce GT640 - Naive parallelization
3 Geforce GT640 - Improved parallelization
Running time [ms]
2.5
1.5
0.5
0
32 64 128 256 512
Number of subintervals for the spatial discretization
Another good comparison is the speedup per Streaming Multiprocessor (Table 4). The
number of SMX seems to be more or less constant between architectures, with high end cards
having between 13 to 15, or even twice that if they are dual cards. In this case we can clearly
see that the kepler architecture achieves more than 27x speedup per SMX over the serial
version, so even if we don’t account for the bandwith improvement, which is also a big
difference between low end and high end cards, we could see speedups exceding the 380x
with the high end cards of the Kepler architecture such as the Geforce Titan or the K20 cards.
124 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
FTCS finite difference scheme Vicente Cuellar Moro,
GPGPU parallel computing for the heat conduction equation Miguel Martín Stickle y Manuel Pastor Pérez
Finally if we adjust the speedup with the ratio of power draw between the CPU and the
GPU we obtain a table (Table 5) that brings the attention to the energy efficiency of the new
architectures of GPUs. The performance adjusted for energy consumption is simply fantastic,
these cards can be the building block of low cost, low energy cosumption clusters, there are
already hardware solutions like the Carma and Kayla development kits that point into this
direction. This is from our point of view one of the most interesting results, very good
speedups with very cheap hardware at very low operation costs.
5. Concluding remarks
A PDE model of a heat conduction problem in two dimensions has been numerically solve
by the forward in time central in space (FTCS) finite difference method. Three different
implementations have been developed one under serial computing and the other two under
parallel procedure. The parallel implementations have been run over Quadro FX3800 and
Geforce GT640 GPU architectures.
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 125
Vicente Cuellar Moro, Miguel Martín Stickle y Manuel Pastor Pérez Investigación
Emphasis has been devoted to GPU parallel computing with CUDA language. This
combination of architecture and language is outstanding. It permits a huge amount of
execution threads, compared to conventional parallel computing, while not many lines of
code are required to achieve a great control of the hardware due to CUDA language.
The analysis of the results shows the importance of parallel computing when dealing with
large problems that are impractical or impossible to solve them with a serial computing
procedure as the execution times grow so fast.
Additionally we can see the advantages the new architectures of GPUs. On one side
inclusion of caches has improved the performance of the simpler implementations. On the
other side the increase of performance from one architecture to the next is remarkable, but the
performance adjusted for energy consumption is even more impressive.
6. References
[1] LOGAN JD, Applied Partial Differential Equations, Springer, New York, 2004.
[3] BURDEN RL, FAIRES JD, Numerical Analysis, Brooks/Cole-Thomson Learning, 2001.
126 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Investigación
Estimación de sucesos raros mediante Redes
Bayesianas
Rare events estimation using Bayesian Networks
Francisco Soler Flores y José Ángel Olivas Varela
Revista de Investigación
1 de abril de 2014
Resumen
La teoría y aplicaciones de sucesos de baja probabilidad de ocurrencia o sucesos raros es-
tá en boga en los últimos años debido a su importancia práctica en muy diferentes campos,
como seguros, finanzas, ingeniería o ciencias del medio ambiente. Este trabajo presenta una
metodología para predecir sucesos raros basada en Redes Bayesianas que permite a su vez
estudiar situaciones alternativas.
Abstract
Theory and applications of rare events, very important in recent years due to its practical
importance in very different fields such as insurance, finance, engineering or environmental
science. This paper presents a methodology for predicting rare events based on Bayesian net-
works that in turn allows to study alternative scenarios.
1. Introducción
El tratamiento de sucesos raros, sucesos que ocurren con una probabilidad baja, es un pro-
blema complejo y amplio cuyo tratamiento se enmarca en el ámbito de modelización de la
incertidumbre, teoría de la decisión y donde el estudio del riesgo, definido en términos de teo-
ría de decisiones como las pérdidas promedio o las pérdidas que se pronostican cuando algo
malo sucede, es importante. La ’Ley de Eventos Raros’, demostrada por Poisson, fundamenta
matemáticamente el concepto de suceso raro. Esta ley que lleva su nombre, la ley de los sucesos
raros ([1]) es también llamada por Bortkiewicz ’ley de los pequeños números’ ([2]).
En las últimas décadas se han desarrollado numerosas técnicas de análisis y modelización
de datos en distintas áreas de la estadística ([3],[4] ) y la inteligencia artificial ([5] ) que se han
127
Francisco Soler Flores y José Ángel Olivas Varela Investigación
aplicado al estudio de sucesos raros . La Minería de Datos (MD) ([6] ) es un área moderna in-
terdisciplinar que engloba a aquellas técnicas que operan de forma automática (requieren de la
mínima intervención humana) y, además, son eficientes para trabajar con las grandes cantida-
des de información disponibles en las bases de datos de numerosos problemas prácticos. Estas
técnicas permiten extraer conocimiento útil (asociaciones entre variables, reglas, patrones, etc.)
a partir de la información cruda almacenada, permitiendo así un mejor análisis y comprensión
del problema. En algunos casos, este conocimiento puede ser también post-procesado de forma
automática permitiendo obtener conclusiones, e incluso tomar decisiones de forma casi auto-
mática, en situaciones prácticas concretas (sistemas inteligentes). La aplicación práctica de estas
disciplinas se extiende a numerosos ámbitos comerciales y de investigación en problemas de
predicción, clasificación o diagnosis.
Con el avance y la extensión de los modelos de inteligencia artificial, tales como las redes ba-
yesianas, se plantea, en este trabajo, la utilización de estos últimos para estimar la probabilidad
de ocurrencia de un suceso raro.
El objetivo de este trabajo es presentar una metodología que, a partir de los datos brutos
de los que se dispone habitualmente a la hora de estudiar cualquier problema en el marco
del análisis de datos, sistematiza el estudio de sucesos raros y es capaz de estudiar diferentes
situaciones que modifican la probabilidad de ocurrencia de estos.
El trabajo se divide en tres secciones principales:
Minería de datos y redes bayesianas: Se describen los conceptos básicos en los que se
apoya la metodología.
Metodología propuesta: Se describe en detalle el procedimiento propuesto.
Conclusiones.
Entre las diferentes técnicas disponibles en minería de datos, las redes probabilísticas o re-
des bayesianas permiten modelizar de forma conjunta toda la información relevante para un
128 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Estimación de sucesos raros mediante Redes Bayesianas Francisco Soler Flores y José Ángel Olivas Varela
n
∏ p(xi | pa(xi )) (1)
i =1
Normalmente las redes bayesianas consideran variables discretas o nominales, por lo que si
no lo son, hay que discretizarlas antes de construir el modelo. Aunque existen modelos de redes
bayesianas con variables continuas, éstos están limitados a variables gaussianas y relaciones
lineales. Los métodos de discretización se dividen en dos tipos principales: no supervisados y
supervisados ([19]).
El concepto de causalidad ([20]) en una red bayesiana se traduce en un caso particular de es-
tas denominado red causal ([21]). Las redes bayesianas pueden tener una interpretación causal
y aunque se utilizan a menudo para representar relaciones causales, el modelo no tiene por que
representarlas de esta forma, Naive-Bayes (apartado 2.2) es un ejemplo de esto, sus relaciones
no son causales.
Las redes probabilísticas automatizan el proceso de modelización probabilística ([22]) utili-
zando la expresividad de los grafos. Los modelos resultantes combinan resultados de la teoría
de grafos (para representar las relaciones de dependencia e independencia del conjunto de va-
riables) y de la probabilidad (para cuantificar estas relaciones). Esta unión permite realizar de
forma eficiente tanto el aprendizaje automático del modelo, a través del cálculo de parámetros
([23]) que para el caso de variables binarias se modela con una distribución Beta y para varia-
bles multivaluadas mediante su extensión, que es la distrubución Dirichlet (tabla 1),como la
inferencia a partir de la evidencia disponible. La base de conocimiento de estos sistemas es una
estimación de la función de probabilidad conjunta de todas las variables del modelo, mientras
que el módulo de razonamiento es donde se hace el cálculo de probabilidades condicionadas. El
estudio de esta técnica proporciona una buena perspectiva global del problema del aprendizaje
estadístico y la minería de datos.
estimador expresión
N
Máxima verosimilitud. Multinomial θk∗ = k
N
Nk + ak
Estimación bayesiana. Dirichlet θk∗ =
N + ∑ri=1 ai
Una de las formas más sencillas que se pueden idear al considerar la estructura de una
red Bayesiana con objetivos clasificatorios, es el denominado Naive-Bayes o Bayes ingenuo ([24]).
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 129
Francisco Soler Flores y José Ángel Olivas Varela Investigación
X1 X2 X3 X4 Xn
El modelo Naive Bayes es muy utilizados debido a que presentan, entre otras, ciertas venta-
jas ([6]):
Dado que suponemos que las variables predictoras son condicionalmente independeientes
dada la variable C, se obtiene que (ecuación 3)
n
P(C = θi | X1 = x1j , . . . , Xn = xnj ) ∝ P(C = θi ) · ∏ P( Xr = xrj |C = θi ) (3)
r =1
130 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Estimación de sucesos raros mediante Redes Bayesianas Francisco Soler Flores y José Ángel Olivas Varela
conocida para el estudio de sucesos raros. De esta forma y conocidos los valores de las varia-
bles que se utilizan como predictoras se pueden estudiar las diferentes situaciones posibles y
observar cuándo es más probable la ocurrencia de un suceso raro.
3. Metodología propuesta
A continuación se detallará la metodología presentada en este trabajo. El objetivo de la mis-
ma es sistematizar el estudio de sucesos raros, proporcionar una herramienta para su estudio
partiendo de los datos brutos, a priori, de cualquier problema.
Así las fases de este metodología son las siguientes:
Normalmente las redes bayesianas consideran variables discretas o nominales, por lo que si
no lo son, hay que discretizarlas antes de construir el modelo. Aunque existen modelos de redes
bayesianas con variables continuas, éstos están limitados a variables gaussianas y relaciones
lineales. Los métodos de discretización se dividen en dos tipos principales: no supervisados y
supervisados ([19]).
En este apartado de preparación de los datos, éstos son formateados de forma que las he-
rramientas informáticas que se uilicen puedan manipularlos. A su vez, este apartado consiste
también en la selección de variables y discretización de los datos en el caso de variables con-
tinuas. En este paso se identifica también la variable de la cual se quieren estimar los sucesos
raros.
En este caso se propone un modelo concreto, el Naive-Bayes (apartado 2.2), adecuado para
el proceso de clasificación y para el que la variable de la cual se quieren estimar los sucesos
raros será la variable denominada padre(([18])).
Se va a optar para el desarrollo de la metodología por el modelo Naive Bayes. El modelo
Naive Bayes es muy utilizados debido a que presentan, entre otras, ciertas ventajas:
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 131
Francisco Soler Flores y José Ángel Olivas Varela Investigación
Las inducciones de son extremadamente rápidas, requiriendo sólo un paso para hacerlo.
Una vez determinada la estructura de la red se calculan las tablas de probabilidades que
permiten realizar la descomposición de la distribución de probabilidad y posteriormente reali-
zar inferencia basándose en las evidencias obtenidas y la distribución de probabilidad asociada
a la red.
Esta información que nos proporciona la red construida se completa aplicando el modelo
Naive-Poisson descrito a continuación (apartado 3.3).
Suposición Poisson
λx
P ( X = x ) = e−λ (4)
x!
Para el cálculo del parámetro que determina la distribución e Poisson, la media, tomaremos
su estimador máximo verosímil que viene dado por la ecuación 5, donde xi son los valores
discretos de la accidentalidad y p( xi ) los valores de probabilidad que proporciona la red cons-
truida mediante el algoritmo Naive-Bayes , siendo la variable discreta a clasificar C, la variable
de la que queremos estudiar los valores que consideramos sucesos raros y el resto, las variables
utilizadas para tal fin el resto de las variables.
n
λ= ∑ xi · p( xi ) (5)
i =1
Para cada uno de los valores de los estratos a1 , a2 , . . . , an que proporciona la discretización
de la variable a estudiar los sucesos raros, se tomarán los valores reales de la frecuencia del
suceso raro a estimar excepto para el an para el que se asigna un valor dado por la media de los
valores x i superiores al mayor estrato an−1 (ecuación 6 ).
132 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Estimación de sucesos raros mediante Redes Bayesianas Francisco Soler Flores y José Ángel Olivas Varela
0.25 0.25
caso 1 caso 1
caso 2 caso 2
caso 3
caso 4
caso 5
0.2 0.2
0.15 0.15
Probabilidad
Probabilidad
0.1 0.1
0.05 0.05
0 0
0 2 4 6 8 10 12 14 16 0 2 4 6 8 10 12 14 16
De esta forma descrita obtenemos las distribución de poisson asociada a cualquiera de las
situaciones (j = 1 . . . , m) que estudiemos (ecuación 7) con cualquier conjunto de valores de las
diferentes variables seleccionadas (en la figura 2 se observan dos ejemplos distintos), con lo cual
es posible determinar en que situación es más alta la probabilidad de ocurrencia de los sucesos
de baja probabilidad una vez detectado cuales son los valores de la variable estudiada que dada
su distribución son de probabilidadmás baja (figura 3).
λ xj
Pj ( X = x ) = e−λ j (7)
x!
0.25
caso 1
caso 2
caso 3
caso 4
caso 5
0.2
0.15
Probabilidad
0.1
0.05
sucesos raros
0
0 2 4 6 8 10 12 14 16
Figura 3. Sucesos raros y sus distribuciones según los valores de las variables predictoras
4. Conclusiones
En este trabajo se propone una metodología para que, de forma sencilla, a partir cualquier
conjunto de datos y una vez seleccionada la variable objeto de estudio, se analicen los sucesos
de baja probabilidad o sucesos raros, además de comparar diferentes situaciones para poder
controlar así la probabilidad con la que se producen. La metodología presentada permite de
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 133
Francisco Soler Flores y José Ángel Olivas Varela Investigación
manera sistemática el estudio de la ocurrencia de los sucesos objeto del estudio y plantear alter-
nativas para controlar su ocurrencia.
Las líneas que abre este trabajo son entre otras la del desarrollo de software que permita su
utilización de manera intuitiva, el estudio de las medidas de bondad de ajuste que completen
su aplicación y el contraste con otros modelos de inteligencia artificial.
Referencias
[1] Simon Denis Poisson. Recherches sur la probabilité des jugements en matière criminelle et en matière civile,
précédées des règles générales du calcul des probabilités. Bachelier, 1837.
[2] Ladislaus Bortkiewicz. Das gesetz der klienen zahlen (the law of small numbers.). Leipzig, Germany:
Teubner, 1898.
[3] Micahael Tomz, Gary King, and Langche Zeng. Relogit: Rare events logistic regression. Journal of
statistical software, 8(i02), 2003.
[4] F. Soler-Flores, J. M. Pardillo Mayora, and R. Jurado Piña. Tratamiento de outliers en los modelos de
predicción de accidentes de tráfico. VIII CIT, 02/07/2008-04/07/2008, La Coruña, España., 2008.
[5] Chris Seiffert, Taghi M. Khoshgoftaar, Jason Van Hulse, and Amri Napolitano. Mining data with
rare events: a case study. In Tools with Artificial Intelligence, 2007. ICTAI 2007. 19th IEEE International
Conference on, volume 2, pages 132–139. IEEE, 2007.
[6] D. E. Holmes, J. Tweedale, and L. C. Jain. Data mining techniques in clustering, association and
classification. Data Mining: Foundations and Intelligent Paradigms, pages 1–6, 2012.
[7] D. J. Hand. Data mining: statistics and more? The American Statistician, 52(2):112–118, 1998.
[8] E. Castillo, J. M. Gutiérrez, and A. S. Hadi. Expert systems and probabilistic network models. Springer
Verlag, 1997.
[9] J. Bromley, NA Jackson, OJ Clymer, AM Giacomello, and F. V. Jensen. The use of hugin to develop
bayesian networks as an aid to integrated water resource planning. Environmental Modelling and
Software, 20(2):231–242, 2005.
[10] K. J. Cios, W. Pedrycz, and R. W. Swiniarski. Data mining: A knowledge discovery approach. Springer
Verlag, 2007.
[11] Usama Fayyad, Gregory Piatetsky-Shapiro, and Padhraic Smyth. From data mining to knowledge
discovery in databases. AI magazine, 17(3):37, 1996.
[12] S. Acid and L. de Campos. Approximations of causal networks by polytrees: an empirical study.
Advances in Intelligent ComputingâĂŤIPMU’94, pages 149–158, 1995. ID: 149.
[13] Gary M. Weiss. Mining with rarity: a unifying framework. ACM SIGKDD Explorations Newsletter,
6(1):7–19, 2004.
[14] Seong-Pyo Cheon, Sungshin Kim, So-Young Lee, and Chong-Bum Lee. Bayesian networks based rare
event prediction with sensor data. Knowledge-Based Systems, 22(5):336–343, 2009.
[15] A. Ebrahimi and T. Daemi. Considering the rare events in construction of the bayesian network
associated with power systems. In Probabilistic Methods Applied to Power Systems (PMAPS), 2010 IEEE
11th International Conference on, pages 659–663. IEEE, 2010.
[16] J. H. Kim and J. Pearl. A computational model for causal and diagnostic reasoning in inference
systems. In Proceedings of the 8th International Joint Conference on Artificial Intelligence, pages 190–193.
Citeseer, 1983.
[17] J. Pearl. Probabilistic reasoning in intelligent systems: networks of plausible inference. Morgan Kaufmann,
1988.
[18] E. Castillo, J. M. Gutierrez, and A. S. Hadi. Expert systems and probabilistic network models. Springer
Verlag, 1997.
[19] James Dougherty, Ron Kohavi, and Mehran Sahami. Supervised and unsupervised discretization of
continuous features. In ICML, pages 194–202, 1995.
134 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Estimación de sucesos raros mediante Redes Bayesianas Francisco Soler Flores y José Ángel Olivas Varela
[20] Cristina Puente Agueda. Causality in sciencie. Pensamiento Matemático, (1):12, 2011.
[21] Judea Pearl. Causality: models, reasoning and inference, volume 29. Cambridge Univ Press, 2000.
[22] Laura Uusitalo. Advantages and challenges of bayesian networks in environmental modelling. Eco-
logical Modelling, 203(3):312–318, 2007.
[23] Luis Enrique Sucar. Redes bayesianas.
[24] R. O. Duda, P. E. Hart, and D. G. Stork. Pattern classification. 2000. NY Wiley. ID: 129.
[25] P. Langley, W. Iba, and K. Thompson. An analysis of bayesian classifiers. In Proceedings of the National
Conference on Artificial Intelligence, pages 223–223. JOHN WILEY and SONS LTD, 1992.
[26] José Antonio Gámez and José Miguel Puerta. Sistemas expertos probabilísticos, volume 20. Univ de
Castilla La Mancha, 1998.
[27] Sang-Bum Kim, Hee-Cheol Seo, and Hae-Chang Rim. Poisson naive bayes for text classification with
feature weighting. In Proceedings of the sixth international workshop on Information retrieval with Asian
languages-Volume 11, pages 33–40. Association for Computational Linguistics, 2003.
[28] F. Soler-Flores. Naive-poisson, a mathematical model for road accidents frequency estimation. Con-
ference of Informatics and Management Sciences, pages 384–391, 2013.
[29] F. Soler-Flores. Expert system for road accidents frequency estimation based in naive-poisson. Global
Virtual Conference, page 646, 2013.
Sobre el autor:
Nombre: Francisco Soler Flores
Correo electrónico: [email protected]
Institución: Universidad Politécnica de Madrid, España.
Nombre: José Ángel Olivas Varela
Correo electrónico: [email protected]
Institución: Departamento de Tecnologías y Sistemas de Información, Escuela Superior de Infor-
mática, Universidad de Castilla La Mancha, España.
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 135
Juegos Matemáticos
Números triangulares cuadrados (la cuadratura del
triángulo)
Triangular squares (the quadrature of the triangle)
Dionisio Pérez Esteban
Revista de Investigación
1 de abril de 2014
Resumen
Entre los números triangulares, se encuentran camuflados algunos cuadrados. En este ar-
tículo se les desenmascara, y en el proceso nos topamos con algunos personajes interesantes:
fracciones continuas, ecuación de Pell, aproximación de soluciones, etc.
Abstract
Among triangular numbers, some are squares in disguise. We uncover them, and in so
doing we come across some interesting people such as continuous fractions, Pell’s equation
and aproximation of solutions.
Un problema geométrico clásico, que trajo de cabeza a los antiguos griegos (esos “colegas
de otras universidades”, como les llamó Hardy), es el de la cuadratura del círculo, que no con-
siste en dibujar un círculo cuadrado, sino en construir un cuadrado cuya área sea la misma que
la de un círculo dado, y hacerlo usando sólo la regla y el compás. La imposibilidad de cua-
drar el círculo no se logró demostrar hasta el siglo XIX, cuando se probó que el número π es
trascendente.
La cuadratura del triángulo, por el contrario, es muy sencilla de realizar, y puede verse como
un agradable entretenimiento. Aquí vamos a abordar un problema aritmético en lugar de uno
geométrico: qué números triangulares resultan ser también números cuadrados.
137
Dionisio Pérez Esteban Juegos Matemáticos
Recordemos que los números triangulares son los que se obtienen sumando varios enteros
consecutivos a partir de 1. La sucesión de los números triangulares comienza así:
x2 − 2y2 = −1
En el segundo caso, se tiene algo parecido, con los papeles de n y n + 1 cambiados: n + 1 es el
cuadrado de x en tanto que n es el doble del cuadrado de y (y en consecuencia
x2 − 2y2 = 1
Recogiendo los dos caso en uno solo, escribimos
x2 − 2y2 = ±1
No es difícil hallar algunas soluciones “pequeñas”de la nueva ecuación por tanteo, proban-
do con números bajos a ver si el doble de su cuadrado se diferencia de otro cuadrado en una
unidad (por exceso o por defecto). Así, encontramos las soluciones
x = 1, y = 1; x = 3, y = 2; x = 7, y = 5; x = 17, y = 12; . . .
138 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Número triangulares cuadrados (la cuadratura del triángulo) Dionisio Pérez Esteban
cada una de las cuales nos facilitan un número triangular cuadrado: por ejemplo, la solución
x = 7, y = 5 de la ecuación
x2 − 2y2 = −1
n ( n +1 )
corresponde a n = x 2 = 49, y el correspondiente número triangular es 2 = 49 · 25, que es
el cuadrado de 35; y la solución x = 17, y = 12 de la ecuación
x2 − 2y2 = 1
288
corresponde a n = 2y2 = 288 y nos proporciona el número triangular (y cuadrado) 289 · 2 =
41616, que es el cuadrado de 204 = 17 · 12.
No es difícil caer en la cuenta de que conociendo una solución x, y de la ecuación, se obtiene
otra haciendo x ′ = x + 2y, y′ = x + y. Así, podemos generar una sucesión de soluciones
(1, 1), (3, 2), (7, 5), (17, 12), (41, 29), (99, 70), . . .
que a su vez nos proporciona una sucesión de números triangulares cuadrados;
Ya podemos contestar una de las preguntas que nos planteábamos: ¿cuántos números trian-
gulares hay que sean también cuadrados? La respuesta es una cantidad infinita (y además sa-
bemos generarlos y determinar el lugar que ocupan en la lista).
√
3. Una dosis de irracionalidad: aparece 2. Y se acabó
La ecuación
x2 − 2y2 = ±1
que hemos resuelto “a ojo”, es una ecuación diofántica de un tipo bien estudiado, tanto que
√ como ecuación de Pell, y se resuelve fácilmente mediante el desarrollo
tiene nombre: se la conoce
en fracción continua de 2. Los convergentes de ese desarrollo proporcionan las soluciones de
la ecuación. Como ese desarrollo es
[1; 2, 2, 2, 2, 2, . . .]
x2 − 2y2 = ±1
proporcionan medidas para los lados de unos triángulos isósceles que son “casi”rectángulos:
así, el triángulo de lados 12, 12 y 17 es ligeramente obtusángulo, y el de lados 29, 29 y 41 es
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 139
Dionisio Pérez Esteban Juegos Matemáticos
levemente acutángulo; hace falta mucha agudeza para detectar la diferencia . Estimando el resto
en la aproximación dada por el primer polinomio de Taylor, podemos concluir que la diferencia
(en radianes) entre el ángulo que forman los dos “cuasi-catetos ”(de longitud l) y un recto es
casi exactamente igual a l12 , lo que supone aproximadamente 4 minutos de arco para el último
triángulo (no digamos en triángulos de lados 408, 408 y 577, o 5741, 5741 y 8119). Desde el punto
de vista del dibujo técnico estaríamos ante un triángulo rectángulo a todos los efectos.
Partiendo de la pregunta acerca de los números que son a la vez cuadrados y triangulares,
hemos realizado una excursión que nos ha llevado desde algunos problemas clásicos de los
antiguos geómetras griegos a los polinomios de Taylor, pasando por ecuaciones diofánticas y
fracciones continuas, en una pálida muestra de la variedad e interconexión de la Matemática.
Sobre el autor:
Nombre: Dionisio Pérez Esteban
Correo electrónico: [email protected]
Institución: ETSI de Caminos, Canales y Puertos. Universidad Politécnica de Madrid, España.
140 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Juegos Matemáticos
Una Sucesión Matemática Curiosa
A Curious Mathematical Sequence
1 de abril de 2014
Resumen
Existe una Sucesión de números naturales, que inician en cero y va aumentando cada vez
el entero correspondiente al orden del elemento de la misma, que tiene un comportamiento
a primera vista curioso, mas al realizar un estudio de la misma se entiende a cabalidad este
comportamiento. La sucesión genera además ciertos resultados interesantes, como el poder
determinar resultados para la potencia de cualquier número entero, en función de sus ele-
mentos.
Abstract
There exists a sequence of natural numbers, starting at zero and increases every time the
integer corresponding to the order of the same element which have curious properties that
are fully understood after a study of the same . The sequence also generates some interesting
identities such as one that determine the powers of any integer as a function of the elements
of the sequence.
1. Introducción
Las sucesiones matemáticas se construyen partiendo de elementos base y generando los
sucesivos con alguna regla y con operaciones básicas que generalmente responden al orden de
cada elemento, por tanto siempre existirá una relación entre los elementos, más en la mayoría
de estas sus elementos no presentan algún resultado que genere interés alguno.
Se ha descubierto una sucesión de números naturales, que siendo de fácil construcción, ge-
nera elementos que sí brindan resultados interesantes y que exigieron un estudio para que se
141
Marco Vinicio Vásquez Bernal Juegos Matemáticos
entienda la naturaleza de este hecho, a sabiendas de que la naturaleza de las matemáticas siem-
pre brindará explicaciones lógicas y científicas a esos resultados que aparentemente surgen del
azar.
Esta sucesión presenta resultados por demás interesantes, relacionando sus elementos con
potencias de todos los números naturales, en este trabajo se explica esta sucesión y además la
relación entre sus elementos y las potencias.
2. Presentación de la Sucesión
que permite hacer los cálculos para obtener los elementos de la sucesión:
a0 = 0
a1 = a0 + 1 = 0 + 1 = 1
a2 = a1 + 2 = 1 + 2 = 3
a3 = a2 + 3 = 3 + 3 = 6
a4 = a3 + 4 = 6 + 4 = 10
a = a + 5 = 10 + 5 = 15
5 4
a6 = a5 + 6 = 15 + 6 = 21
a7 = a6 + 7 = 21 + 7 = 28
a8 = a7 + 8 = 28 + 8 = 36
a9 = a8 + 9 = 36 + 9 = 45
a 10 = a9 + 10 = 45 + 10 = 55
...
Resultado 1
n = a n − a n −1
142 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Una Sucesión Matemática Curiosa Marco Vinicio Vásquez Bernal
Resultado 2
15 + 21 = 36 = 62
21 + 28 = 49 = 72
28 + 36 = 64 = 82
36 + 45 = 81 = 92
45 + 55 = 100 = 102 . . .
n2 = a n + a n −1
Demostración
Utilizaremos inducción sobre n.
Supondríamos que se cumple ( p − 1)2 = a p−1 + a p−2 , que sería nuestra hipótesis de induc-
ción.
Recordemos además que la definición de la Sucesión establece:
a p = a p−1 + p, ∀ p ≥ 1, (1)
a p = a p−2 + p − 1 + p = a p−2 + 2p − 1, ∀p ≥ 2
de donde:
a p−2 = a p − 2p + 1
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 143
Marco Vinicio Vásquez Bernal Juegos Matemáticos
( p − 1)2 = a p−1 + a p − 2p + 1
n2 = a n + a n −1
n2 = an + an−1 ⇒ n2 = an + an − n ⇒ n2 + n = 2an
que se reduce a:
n2 + n n( n + 1)
an = =
2 2
Como consecuencia lógica se tendrá que:
(n − 1) · ((n − 1) + 1) n( n + 1)
an = =
2 2
resultado que ya se obtuvo y que muestra la razón en la que se sustenta la curiosidad de los
elementos de la sucesión.
4. Generación de Potencias
Partiendo de los hechos ya demostrados:
n2 = a n + a n −1 y n = a n − a n −1
n3 = a2n − a2n−1
144 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Una Sucesión Matemática Curiosa Marco Vinicio Vásquez Bernal
Potencia Resultado
n1 a n − a n −1
n2 a n + a n −1
n3 a2n − a2n−1 = ( a n − a n −1 ) · ( a n + a n −1 )
n4 ( a n − a n −1 ) · ( a n − a n −1 ) = ( a n − a n −1 ) 2 · ( a n + a n −1 )
2 2
Se debe notar que los exponentes de los resultados son menores, significativamente que los
de la potencia, así para expresar n14 , se utilizan las potencias de 5 y 4, evidenciándose que los
factores del resultado son menores o iguales a la tercera parte de la potencia.
Con lo que se observa que se tiene un resultado general, para todo s, numero natural:
s s
3s
n = ( a n + a n −1 ) · ( a n − a n −1 )
n 3s + 1 = ( a n + a n −1 ) · ( a n − a n −1 ) s +1
s
3s+2
n = ( a n + a n −1 ) s +1 · ( a n − a n −1 ) s
Demostración
La demostración la realizaremos por inducción.
Veamos para s = 0.
0 0 0
n = ( a n + a n −1 ) · ( a n − a n −1 ) = 1
n1 = ( a n + a n −1 ) 0 · ( a n − a n −1 ) 1 = ( a n − a n −1 ) = n
2
n = ( a n + a n −1 ) 1 · ( a n − a n −1 ) 0 = n2
Ahora supongamos que lo expuesto se cumple para s = k, por lo que tendremos las siguien-
tes igualdades:
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 145
Marco Vinicio Vásquez Bernal Juegos Matemáticos
Con los resultados (5), (6) y (7), podemos concluir que lo propuesto se cumple con s = k + 1,
por lo que podemos concluir con la demostración.
Y como es conocido:
n
n( n + 1)
∑i= 2
;
i =0
146 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Una Sucesión Matemática Curiosa Marco Vinicio Vásquez Bernal
entonces:
n( n + 1)
an = ;
2
a n + a n −1 = n2 , y a n − a n −1 = n
Utilizando la fórmula de la suma para una sucesión de números naturales consecutivos que
inician en 1 se tendría que:
!2
n n n n
(i − 1) · i
2
∑i = ∑ i + 2i · = ∑ i 2 + i · (i − 1) · i = ∑ (i 2 + i 3 − i 2 ) ⇒
i =1 i =1
2 i =1 i =1
!2
n n
⇒ ∑i = ∑ i3
i =1 i =1
Consecuentemente:
n
a2i = ∑ i3
i =1
n3 = a2n − a2n−1
Referencias
[1] O LMAN, Bernard, y H ILL, David R., Álgebra Lineal, México, Pearson Educación, 2006, ISBN
970-26-0696-9.
[2] J OHNSONBAUGH, Richard., Matemáticas Discretas, México, Pearson Educación, 2005, ISBN
970-26-0637-3.
[3] B RASSARD, G. y B RATLEY, P., Fundamentos de Algoritmia, Madrid, Prentice Hall, 1997, ISBN
84-89660-00-X.
[4] K ENNETH, H. Rosen., Discrete mathematics and its applications, McGraw Hill, 2003. ISBN 0-
07-123374-1.
[5] K ENNETH, H. Rosen y M ICHAELS, John G., Handbook of discrete and combinatorial mathematics,
CRC, 1999. ISBN 0-8493-0149-1.
[6] V OROBIOV, N. N., Números de Fibonacci, Editorial Mir, Moscú, 1974, Colección Lecciones Po-
pulares de Matemáticas. Traducción al español de Carlos Vega, catedrático de Matemáticas
Superiores y candidato a doctor en ciencias físico-matemáticas.
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 147
Marco Vinicio Vásquez Bernal Juegos Matemáticos
Sobre el autor:
Nombre: Marco Vinicio Vásquez Bernal
Correo electrónico: [email protected]
Institución: Instituto Tecnológico Andrés F. Córdova, Quito, Ecuador.
148 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Críticas
Informe sobre el libro: “El enigma de Fermat. Tres
siglos de desafío a la matemática”, Albert Violant
A report of the book: “El enigma de Fermat. Tres
siglos de desafío a la matemática”, Albert Violant
Javier Rodrigo Hitos
Revista de Investigación
1 de abril de 2014
Resumen
En este artículo se presenta un informe sobre el libro “El enigma de Fermat. Tres siglos
de desafío a la matemática”, perteneciente a la colección “El mundo es matemático”, de la
que ya se comentó otro libro en un artículo anterior del autor.
Abstract
This paper presents a report of the book “El enigma de Fermat. Tres siglos de desafío a
la matemática”, that belongs to the collection “El mundo es matemático”. Another book
from this collection was commented by the author in a previous paper.
1. Ficha técnica
Título: El enigma de Fermat. Tres siglos de desafío a la
matemática.
Autor: Albert Violant.
De la colección “El mundo es matemático”
Editado por: RBA
149
Javier Rodrigo Hitos Críticas
2. Crítica
El libro del que se hace el informe pertenece a una colección de divulgación de las
matemáticas de la que ya se comentó otro libro (Los números primos. Un largo camino al
infinito. Ver [2]). Destaca principalmente por su esfuerzo para hacer asequibles para el público
general temas de mucha dificultad matemática, sin perder el rigor necesario en un texto
científico.
Es especialmente notable en el texto que nos ocupa el tratamiento que se hace de las
teorías de Kummer, explicando de forma amena y brillante conceptos revolucionarios en su
día como el de los números ideales ó el de los primos regulares, nada elementales y difíciles
de asimilar. Se pone muy bien de manifiesto la importancia de estas teorías como avances en
el intento de demostrar el teorema de Fermat.
Como aspectos mejorables, destacar el desarrollo desigual de los diferentes capítulos del
libro: la introducción “prehistórica” al teorema de Fermat desarrollada en los primeros
capítulos, en la que se trata el teorema de Pitágoras, es quizás algo larga, remontándose a
tiempos demasiado antiguos y, aunque es de gran interés cultural, desvía la atención del tema
principal del libro.
A este respecto, se puede decir que se dedican unas 128 páginas a los preliminares y sólo
15 páginas a la prueba del teorema. Por la gran dificultad de la misma, quizás el capítulo que
la trata es uno de los más débiles del libro, deslizándose algún error que se comenta
posteriormente en este informe.
A pesar de esto, se puede decir que es un libro altamente recomendable incluso para
aquellos que no tengan noticias previas del teorema de Fermat, ni una formación matemática
elevada (ver [3] y [4] como referencias de libros sobre el teorema de Fermat más avanzados).
- Cicloide: lugar geométrico que describe un punto de la circunferencia al girar sobre una
recta
- Braquistócrona: curva que da la forma más rápida posible de ir de un punto “A” a otro
“B” sometidos únicamente a la acción de la gravedad
- Números perfectos: números que son suma de sus divisores propios
150 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Informe sobre el libro: “El enigma de Fermat” Albert Violant Javier Rodrigo Hitos
- Números poligonales: aquellos que pueden ser representados con puntos dispuestos en
forma de polígono regular, empezando por el 1
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 151
Javier Rodrigo Hitos Críticas
- 60 es el número más pequeño que se puede expresar de 6 formas distintas como suma
de 2 primos
- Si n es un entero positivo tal que sus divisores menores que n son consecutivos,
entonces n es primo ó el doble de un primo ó uno de los números 1, 8, 12, 24, 60.
- Los números con sólo factores primos 2, 3 ó 5 cumplen que sus inversos siempre tienen
un número finito de cifras decimales en el sistema sexagesimal
- La longitud de un arco de cicloide es 8 veces el radio de la circunferencia que gira.
- Si se deja caer un objeto sobre un arco de cicloide invertida, no importa desde qué
altura, siempre tarda lo mismo en llegar a la base
- Principio de Fermat: la luz sigue trayectorias de tiempo mínimo (relacionado con la ley
de refracción de Snell)
- El área de un triángulo rectángulo de lados enteros no puede ser nunca un número
cuadrado
- Un número par es perfecto si y sólo si es de la forma 2n1 2n 1 , con 2n 1 primo
(primo de Mersenne, por tanto. Demostrado por Euler)
- No hay números perfectos impares menores que 10 300
- Un polígono regular es construible con regla y compás si la factorización en primos de
su número de lados contiene doses ó números de Fermat distintos.
- 243112609 1 es el mayor primo de Mersenne que se conoce (hasta el 12 de Junio de 2009) 1
- Una ecuación de primer grado con coeficientes enteros tiene soluciones enteras si y sólo
si el máximo común divisor de los coeficientes divide al término independiente.
- Todo número natural se puede expresar como suma de cuatro cuadrados (demostrado
por Lagrange)
- Todo número natural se puede expresar como suma de como mucho tres números
triangulares (demostrado por Gauss)
- Todo número natural se puede expresar como suma de como mucho n números n-
gonales (generalización de las anteriores. Demostrado por Cauchy)
- Todo número primo de la forma 4 k 1 puede expresarse como suma de dos cuadrados
(demostrado por Euler)
- Teorema fundamental de la aritmética: Todo número se descompone de manera única
como producto de factores primos
- Si el producto de dos números primos entre sí es un cuadrado, cada uno de ellos es un
cuadrado
57885161
1Actualmente el mayor número primo de Mersenne es el 2 1 , descubierto el 25 de enero de 2013 por Curtis
Cooper con el proyecto GIMPS -Great Internet Mersenne Prime Search - ("Gran búsqueda de números primos de
Mersenne por Internet") un proyecto colaborativo de voluntarios que utilizan los programas gratuitos Prime95 y
MPrime con el fin de buscar números primos de Mersenne.
152 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Informe sobre el libro: “El enigma de Fermat” Albert Violant Javier Rodrigo Hitos
- Dividir un número cuadrado en suma de dos cuadrados (problema 8 del libro II.
Inspiró el teorema de Fermat)
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 153
Javier Rodrigo Hitos Críticas
1 8 1 6
x2 x2
2 5 2 5
1 8 1 12
- En la página 94, pone x3 por x3
2 5 2 5
1 8 1 9
x4 x4
2 5 2 5
- En la página 119, pone las las curvas elípticas por las curvas elípticas
- En la página 123, pone desanimar hasta el más por desanimar hasta al más
- En la página 125, pone dada una una solución por dada una solución
- En el final de la misma página, pone el último teorema no tiene soluciones. Sería más
correcto decir: la ecuación de Fermat no tiene soluciones no triviales
- En la página 126, pone estadounindense por estadounidense y asisistido por asistido
- En la página 127, se dice que la conjetura conlleva una relación necesaria entre la
conjetura de Taniyama-Shimura y el último teorema: si la primera era cierta, la otra no podía
serlo. En realidad es: si la primera era cierta, el otro lo sería
- En la página 137, pie de foto: y el interviene por y en él interviene
- En la página 138, un documento de complejidad por un documento de la complejidad
- En la página 143, falta la interrogación en ¿quién de nosotros…?
154 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Informe sobre el libro: “El enigma de Fermat” Albert Violant Javier Rodrigo Hitos
Referencias
[1] VIOLANT, Albert. El enigma de Fermat. Tres siglos de desafío a la matemática, RBA. Colección
El mundo es matemático, Madrid, 2010.
[2] GRACIÁN, Enrique. Los números primos. Un largo camino al infinito, RBA. Colección El
mundo es matemático, Madrid, 2010.
[3] EDWARDS, Harold M. Fermat’s last theorem, Springer Verlag, New York, 1977.
[4] RIBENBOIM, Paulo. 13 lectures on Fermat’s last theorem, Springer Verlag, New York, 1979.
Sobre el autor:
Nombre: Javier Rodrigo Hitos
Correo Electrónico: [email protected]
Institución: Universidad Pontificia Comillas, Madrid, España.
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 155
Entrevista
Pablo Amster: La Música de la Matemática
Pablo Amster: The Music of Mathematics
Rosa María Herrera
Revista de Investigación
1 de abril de 2014
Resumen
“Soy matemático. En una época era músico, o al menos eso creía. Después empecé a decir que era
buen matemático entre los músicos y buen músico entre los matemáticos. Pero ahora me asumo como
matemático, aunque siga haciendo algo de música. Quizás la diferencia es que ahora que no me siento
músico la disfruto más”. Este es el autorretrato que arranqué a Pablo al empezar este trabajo,
conociéndole un poco esto es muy escaso, pero es una herramienta válida para poner en con-
tacto al lector y al matemático.
Abstract
“I’m a mathematician. Before I was a musician, or so I thought. Then, I began to say I was a good
mathematician among musicians and a good musician among mathematicians. But now I assume my-
self as a mathematician, although I continue making some music. Perhaps the difference is that now
I do not feel musician I enjoy it much more”. Here the self-portrait I have extracted from Pablo
when we started this conversation; for those who know him a bit this is not much, but it is a
valuable tool to introduce this mathematician to reader.
1. Introducción
Pablo Amster nació en Buenos Aires, pero al hablar con él enseguida se descubre un ciuda-
dano del mundo. Me impresionó su talento de escritor de matemática como cultura antes de
que nos hiciéramos amigos. Su biografía es de una riqueza extraordinaria: matemático de pri-
mer nivel perteneciente a los círculos internacionales en el ámbito del análisis no lineal (ejerce
como investigador en el CONICET y es profesor de la universidad de Buenos Aires); la música
le dio el mundo, o quizá sería mejor decir que tengo la impresión de que aprendió la vida en la
música a través del mundo.
Como el hilo conductor de la serie de entrevistas que la directora del grupo me propuso, me
gustaría que fuese la belleza matemática, las relaciones y las estructuras. La primera pregunta
que hice a nuestro matemático fue la siguiente:
157
Rosa María Herrera Entrevista
- Pablo, ¿eres un músico que hace matemática o un matemático que hace música?, ¿qué Pablo es el
que es más tú el músico o el matemático?
Su respuesta concisa y simpática me ha servido para escribir el resumen de apertura. Con-
fieso que no es la primera vez que se lo oigo decir, recuerdo un noviembre en Madrid, en que
hablando a un público severo, sin soporte, sin guión, sin pizarra ni otras apoyaturas, mantuvo
una charla en la que con naturalidad iba intercambiando los papeles protagonistas, del mate-
mático al músico poeta para finalizar acordándose de su guitarra, le faltaba entre los dedos.
Me complazco en las fronteras y como dice un amigo mío citando, más o menos, a cierto
manchego: lo interesante de las puertas es que tienen bisagras.
- ¿Por qué escribes historias poético literarias de música y matemáticas?
La literatura es una de mis grandes pasiones. Todavía no me animo a escribir tantas “histo-
rias” y me quedo más bien en la divulgación, pero siento que ya iré encontrando el espacio para
algo más. Escribo sobre matemática porque es lo mío y me gusta el desafío de la transmisión.
Me interesan especialmente los cruces de discursos. Y la música surge también bastante natu-
ralmente en lo que escribo, al menos en los últimos años, sobre la matemática tiene la ventaja
de que a todo el mundo le gusta, aunque cuando uno se mete a hablar de teoría musical se ve
que no son cosas muy distintas.
- Me emociona, creo que como a ti, la combinación de ideas de mundos diferentes porque sirve para
construir, me gustaría saber, Pablo, ¿por qué te interesan los cruces de discursos?
Lo del “cruce” es más o menos lo que encontrás2 en mis textos: discursos en apariencia
disímiles que confluyen . . . en mi caso, el hilo conductor es sin duda la matemática, pero me
gusta hacer dialogar distintas voces, como me gusta la intertextualidad. En general pienso los
escritos en varias dimensiones. Alguna vez se me ocurrió la idea de escribir como lector, lo que
provocó algunos efectos interesantes: entre otros, el de un texto infinito.
1 http://mate.dm.uba.ar/~pamster/
2 En español argentino en el original.
158 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Pablo Amster: La Música de la Matemática Rosa María Herrera
Con esta respuesta lo primero que siento es deseos de leer más textos tuyos, es fácilmente
comprensible para quien ha tenido la fortuna como yo de escuchar lo que dices y leer lo que
escribes, aunque discrepemos en ocasiones, nos entendemos bien.
3. Hablemos de Matemáticas
Llegados a este punto creo que podemos abordar otro de los motivos que nos ha traído hasta
aquí, en nuestra conversación, el mundo matemático de Pablo.
Figura 2. Portadas de los libros de divulgación escritos por Pablo Amster “La Matemática como una de las bellas artes” y
“Matemática maestro”3
- La hermosa frase “pensar los discursos en varias dimensiones” leída sin detenerse parece un recurso
poético prestado de las matemáticas y me confunde realmente, pero quizá sea una sencilla expresión del
matemático interior . . . En cualquier caso me parece de una riqueza extraordinaria, y un signo del hombre
actual más avanzado, ya que hemos tomado el camino de los discursos cruzados, del infinito y de las varias
dimensiones, háblame un poco del análisis no lineal, de las ecuaciones diferenciales no lineales (que son
las que nos interesan a los físicos) y de tu matemática.
Creo que mi actividad matemática también tiene algo del “cruce de discursos” pues siempre
me gustó leer de distintos temas e intentar combinarlos. En general me siento un poco ajeno en
la mayoría de los encuentros sobre “ecuaciones diferenciales”: si bien se supone que es mi área
de investigación, es muy amplia y contiene algunas cosas que no me interesan mayormente.
En tal sentido, mi visión choca bastante con la de mucha gente que se dedica a esto. En todo
caso, mi aporte puramente matemático es más bien modesto, nunca me sentí capaz de grandes
cosas. Sí creo que logro transmitir bastante bien las ideas a los estudiantes. Siempre me gustó
particularmente destacar la experiencia del error; el error como instancia de aprendizaje. Me
parece que convivimos demasiados años con un sistema que reprime el error, y eso termina
produciendo una gran inhibición y pone límites a la creatividad.
3 http://articulo.mercadolibre.com.ar/MLA-471166749-la-matematica-como-una-de-las-bellas-artes-pablo-amster-
_JM; http://articulo.mercadolibre.com.ar/MLA-476669568-pablo-amster-matematica-maestro-ed-siglo-xxi-_JM
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 159
Rosa María Herrera Entrevista
El lector no tiene por qué saber que esta conversación entre Pablo y yo se desarrolla entre dos
continentes, varios países, viajes, cambios de horario . . . y ese hecho común en nuestros días era
impensable hace pocas décadas, nos contamos de este modo los pequeños y no tan pequeños
avatares de nuestra vida cotidiana, así es que nuestra conversación amistosa se intercala en
la entrevista. Se trata de un aspecto no menor que nos permite el mundo híbrido, el mundo
combinado de la proximidad (de internet) y la lejanía. Pablo está preparando un nuevo libro
dedicado a estudiantes, pero que también sirve de provecho a estudiosos no especialistas y me
cuenta estas cosas desde un encuentro de matemáticos.
Pablo: [El congreso] Más allá de las charlas, me sirvió para ver a la editora y confirmar que
pronto me enviarán el diseño de tapa y comenzará el proceso de edición/corrección. Es un tex-
to claramente diferente a los otros que has leído, aunque hay algunos pasajes en los que pude
deslizar algo de mi estilo. En todo caso, es un libro académico4 , no de divulgación, aunque
tiene toda la intención de ser más o menos sencillo, para que pueda ser entendido por estudian-
tes avanzados. Pero cuando avanza la topología, especialmente con la construcción del grado
topológico, no se puede mentir demasiado.
- ¿No crees, Pablo, que es un poco escasa la buena divulgación y de nivel en matemática, con las
honrosas excepciones que tenemos todos en mente y cuya labor es casi heroica, porque más pronto que
tarde (como ocurre en la construcción del grado topológico) “no se puede mentir demasiado”?
En realidad no es tan escasa, al menos en los últimos años ha crecido muchísimo . . . pero lo
del nivel es, por cierto, discutible. No todo lo que se divulga tiene la misma calidad; aun así,
valoro el esfuerzo de tanta gente que empieza a encontrarle el gusto a transmitir ideas mate-
máticas a personas alejadas, pero que se van acercando tímidamente. Es verdad que siempre
se llega a un punto en el que es difícil avanzar sin que el discurso se vuelva más complejo, pe-
ro entiendo que existe divulgación a distintos niveles. El ejemplo que mencionas, el del grado
topológico, es difícil explicarlo a gente que no tiene ninguna idea de matemática (a grandes
rasgos, puedo decir que se trata de una especie de conteo de la cantidad de ceros de una fun-
ción... aunque en ciertos casos también es preciso decir qué significa “función” y entonces hay
que ir todavía más atrás). Pero puedo introducir bastante bien (y sin mentir demasiado, como
decíamos) la idea intuitiva de grado a estudiantes que apenas conocen unas pocas nociones de
análisis matemático. Si además saben algo de variable compleja, entonces puedo dar una idea
todavía más precisa. La parte que requiere más trabajo es la construcción rigurosa, para la que
hacen falta herramientas más sofisticadas. Sea como fuere, en algún sentido la esencia de la
matemática es mentir, así que por ese lado estoy tranquilo.
- Me gusta mucho tu idea sobre la experiencia del error como instancia de aprendizaje, estoy de acuerdo
contigo, el reprimir el error es una manera de limitar. Me gustaría que desarrollaras un poco más esta
propuesta, ¿en qué consiste una buena transmisión de pensamiento matemático a los estudiantes?
Hay teoremas que en la pizarra se demuestran en pocos minutos y sin embargo llevaron
vidas enteras de trabajo, conjeturas, construcciones conceptuales . . . y especialmente errores de
esos que hacen volver todo al comienzo. Ante los alumnos trato de transmitir parte de esa ex-
periencia que uno vive como investigador, de seguir una línea que al final se revela equivocada
pero a su vez permite comprender muchas cuestiones. Pensar un problema tiene sus momen-
tos de entusiasmo, pero también de desaliento. En general me gusta mostrar los teoremas que
se trabajan en la clase desde mi punto de vista: ¿cómo haría yo para convencerme de que el
4 “Topological Methods in the Study of Boundary Value Problems” Springer 2013, ya publicado en el momento de edición
de esta publicación.
160 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Pablo Amster: La Música de la Matemática Rosa María Herrera
enunciado es cierto? Y en general lo primero que haría seguramente estaría equivocado, pero
luego uno rearma las ideas y finalmente llega al resultado. Y más aun, un punto que me parece
interesante y que cuesta mucho trabajo a los alumnos es el de entender cómo se llega a cier-
tos teoremas. En general, los estudiantes de matemática pueden entender las demostraciones
pero suelen decir que nunca se les ocurriría algo así. En muchos casos hay un hilo natural de
pensamiento que hace que uno llegue a ciertos resultados. Es muy lindo conjeturar enuncia-
dos aunque después resulten falsos; la idea de establecer una conjetura es que tiene que haber
razones para ello.
- Por cierto, aunque para ti sea una obviedad, para los lectores que no están familiarizados con ningún
desarrollo de la topología, sería factible que contaras en pocas palabras por qué es especialmente difícil
“dulcificar” la construcción del grado topológico.
Algo de esto ya te contesté antes. Como te dije, me siento capaz de que estudiantes con una
base mínima de análisis entiendan la idea intuitiva del grado: primero se trata de contar los
ceros de una función “buena” y luego ver cómo eso puede generalizarse. La parte delicada es
la construcción bien rigurosa, pero uno se puede saltar los detalles técnicos y trabajar con las
propiedades del grado sin mayor cargo de conciencia . . . de hecho, creo que mucha gente lo
hace: emplear la teoría porque tiene múltiples aplicaciones muy poderosas, sin molestarse en
entenderla del todo. No digo que esté bien pero, al fin y al cabo, es lo mismo que hacemos con
montones de cosas que empleamos a diario.
- Intuyo que eres un buen profesor, porque creo que combinas el sosiego de quien piensa y escribe con
el dominio escénico del músico intérprete, ¿de dónde emana ese equilibrio que parece convivir en ti con
tanta naturalidad?
Respecto de mi “dominio escénico”, no es más que intuición. En general trato de ponerme
en el lugar del que escucha, algo parecido a lo que te dije respecto de la escritura. Creo que hay
una lógica de la escena, aunque en mi caso es más bien una improvisación. Digamos que me
gusta sorprenderme con lo que yo mismo digo.
A nadie se le escapa el ajetreo que conlleva una vida activa en nuestros días, en la de Pablo,
el continuo trasiego entre países y otros “ires y venires” no le impide realizar además tareas
más personales, como el cuidado de sus pequeños hijos, la atención a la familia, el ocuparse de
los amigos.
No es infrecuente que las personas bien dispuestas asuman cada vez mayores cantidades de
retos. . . ; algo me cuenta Pablo de una traducción al inglés de sus trabajos divulgativos, ¡buena
noticia!, al menos para los lectores de habla inglesa . . . que los hispano hablantes tienen la
oportunidad al alcance de la mano con facilidad.
Tampoco quisiera dejar de señalar aquí su participación en la serie de conciertos Noches
de música y ciencia que se desarrolla en la capital argentina, la “complicidad” del pianista y
el científico conforma, seguro, una atmósfera de gran belleza. Sería grato que el evento estu-
viese disponible para otras ciudades también, y que de esta colaboración se beneficiasen sus
ciudadanos. Algo habrá que pensar al respecto, lanzo el reto . . .
Cabría hablar de las tesis que dirige, de las clases que imparte, de otros asuntos ineludibles
menos académicos propios del quehacer universitario, y es que el lector debe tener presente
que Pablo es un profesor; pero, quién sabe, si este asunto no surgió en la conversación es que
no hace falta tratarlo.
Me gustaría decir que de Pablo he aprendido un poco a mirar los hechos de la matemática
de otra manera. Hay quienes enseñan matemáticas, algunos otros muestran el camino para
Volumen IV, Número 1, Abr’14, ISSN 2174-0410 Revista “Pensamiento Matemático” | 161
Rosa María Herrera Entrevista
aprender a mirarlas, con otros una se da cuenta de que la matemática casi no es lo que ella
pensaba, y así . . .
Para finalizar este bosquejo quería agradecer a Pablo su amabilidad, no se nos escapa a
ninguno la complejidad del mundo global y también del mundo individual y local, que la hacen
más valiosa.
Nos hemos divertido redactando estas notas, el cambio continuo de ambientes, el intercam-
bio de opiniones en medio de la vorágine de nuestras actividades ha sido un aliciente más, el
detalle simpático y dulce. Y para cerrar volvemos a Buenos Aires y a Madrid, que es el punto
de partida. Hasta pronto.
Referencias
[1] A MSTER , P., La matemática como una de las bellas artes, Siglo XXI editores, Buenos Aires, 2008.
[2] A MSTER , P., ¡Matemática maestro! Un concierto para números y orquesta, Siglo XXI editores,
Buenos Aires, 2010.
[3] Plataforma ConCIENCIA Musical, http://concienciamusical-plataforma.blogspot.com.es/
Sobre la autora:
Nombre: Rosa María Herrera
Correo electrónico: [email protected]
Institución: APYCE.
162 | Revista “Pensamiento Matemático” Volumen IV, Número 1, Abr’14, ISSN 2174-0410
Este material está registrado bajo licencia Creative Commons 3.0 Reconocimiento - No Comer-
cial - Compartir Igual, por lo que tienes que tener en consideración que:
Atribución Debes reconocer y citar la obra de la forma especificada por el autor o el licen-
ciante.
Licenciar Igual Si alteras o transformas esta obra, o generas una obra derivada, sólo puedes
distribuir la obra generada bajo una licencia idéntica a ésta.
Al reutilizar o distribuir la obra, tienes que dejar bien claro los términos de la licencia de esta
obra.
Alguna de estas condiciones puede no aplicarse si se obtiene el permiso del titular de los dere-
chos de autor.