Primitivas Elementales

Revista digital
Matemtica, Educacin e Internet

(http://tecdigital.tec.ac.cr/revistamatematica/).
Vol 15, No 2. Marzo Agosto 2015.
Artculo de seccin
ISSN 1659 -0643
Se puede saber si una funcin tiene

primitiva elemental?
Walter Mora F.
[email protected]
Escuela de Matemtica
Instituto Tecnolgico de Costa Rica
Resumen. Este artculo presenta una introduccin elemental a la teora de Lioville (en versin antigua
y moderna) y un teorema de Chebyshev, sobre la integracin en trminos finitos. Se presenta la parte
histrica, la parte computacional y la parte algebraica.
Palabras clave: Primitivas, funciones elementales, integracin en trminos finitos
Abstract. This article presents an elementary introduction to Lioville theory (in old version and
modern version) and also a Chebyshev theorem, on integration in finite terms. It is presented the
historical part, the computational part and the algebraic part.
KeyWords: Primitive, elementary functions, integration in finite terms
1.1
Introduccin
En los cursos de clculo elemental, el problema de clculo de integrales consiste

en: "dada una funcin
Z
elemental f , calcular otra funcin elemental (si existe) g tal que g0 = f (o
f = g )". La solucin par-
cial de este problema se presenta como un proceso heurstico que emplea una "caja de trucos" (tablas
de integrales, sustituciones, integracin por partes y descomposicin en fracciones parciales). Solo en
el caso de funciones racionales P/Q aparece un algoritmo finito (fracciones parciales). Los mtodos
empleados no son mtodos de decisin, solo de clculo: Si no hay xito en el clculo de la primitiva de
Se puede saber si una funcin tiene primitiva elemental? . Walter Mora F.
Derechos Reservados 2015 Revista digital Matemtica, Educacin e Internet (http://tecdigital.tec.ac.cr/revistamatematica/)
Revista digital Matemtica, Educacin e Internet (http://tecdigital.tec.ac.cr/revistamatematica/). Vol 15, No 2. Marzo Agosto 2015.
una de una funcin elemental f , no se sabe si esta primitiva se puede expresar o no, en trminos de
funciones elementales.
Si uno observa varios ejemplos de clculo de primitivas, no parece que hubiera un patrn regular entre
el integrando y la antiderivada.
Ejemplo 1.1
Por ejemplo, usando Wolfram Mathematica obtenemos,
In[]:=
Integrate[1/(x^2 + 1), x]
Integrate[1/Sqrt[1 - x^2], x]
(* = ArcTan[ x ] *)
(* = ArcSin[ x ] *)
h
i
Integrate[Log[x+Sqrt[x^2-1]], x] (* = 1 + x2 + xLog x + 1 + x2 *)
Sin embargo, si introducimos logaritmos y el nmero i se puede empezar a ver un patrn que, como
vamos a ver, es persistente.
Ejemplo 1.2
Usando Wolfram Mathematica obtenemos,
In[]:=
1
1
Integrate[1/(x^2 + 1), x] (* = i Log( x + i ) Log( x i )
2
2
Integrate[1/Sqrt[1 - x^2], x] (* = i Log( 1 x2 + i x )
*)
*)
h
i
Integrate[Log[x+Sqrt[x^2-1]], x] (* = 1 + x2 + xLog x + 1 + x2 *)
J. Liouville fue el primero que prob un teorema que permite tener una base terica para decidir si
una funcin elemental (las funciones usuales del clculo y otras ms) tiene primitiva elemental basado
en la forma que debera tener la funcin y su primitiva.
En general, si F es una funcin algebraica con primitiva elemental, entonces bajo algunas condiciones
se tiene,
Z
F ( x, y1 , y2 , ..., yn ) dx = V0 +
c j Log(Vj )
j =1
Aqu, F ( x, y1 , y2 , ..., yn ) es una expresin algebtraica F ( x, log x, ex , ...).

Por medio de este principio, Liouville pudo probar que varias familias de funciones no tienen primitiva
elemental, en especial ciertas integrales elpticas, que era uno de los temas que dominaba el universo
matemtico de su tiempo.
Este teorema tiene su generalizacin en trminos puramente algebraicos, sigue dominado la forma en
la que deben aparecer las primitivas elementales en el caso general. Sorpresivamente, el algoritmo para
la integracin de funciones racionales P/Q que veremos ms adelante, es similar al que se usa en el
caso general.
Organizacin. Este artculo se divide en dos partes. En la primeria parte vamos a ver los teoremas de
integracin en trminos finitos que fueron establecidos en el siglo XIX.

a) El algoritmo de J. Bernoulli (1703) que es el mtodo de integracin por descomposicin en fracciones parciales que se ensea todava en clculo. Este algoritmo tiene valor terico solamente.
b) El teorema de Chebyshov (1853) que nos permite decidir de manera fcil si funciones de la forma
x m ( a + bx n ) p con a, b, n, m, p Q, tienen primitiva elemental o no.
c) El teorema (fuerte) de Liouville y dos casos especiales, para aplicarlo a ciertas familias de funciones (la teora la desarrollo Lioville entre 1833 y 1841).
En la segunda inicia con el algoritmo para integrar funciones racionales (mtodo de Hermite y el
algoritmo de Rothstein/Trager) pues esta es una gua para los otros algoritmos. Para considerar una
introduccin algoritmo general, se incluye una introduccin mnima al lgebra diferencial y el teorema
de Liouville desde esta perpectiva. Se hace la prueba de un caso especial del teorema de Liouville (solo
un caso, porque la maquinaria algebraica es muy extenuante y a veces muy avanzada). La segunda
parte estar dedicada a una introduccin mnima al algoritmo de Risch; este un algoritmo general para
calcular primitivas o decidir si la primitiva no es elemental. Solo consideramos el algoritmo para casos
de funciones con extensiones logartmicas.
En general, solo necesitamos hacer clculos con lgebra de polinomios: Divisin de polinomios con
resto, mcd de dos polinomios y el algoritmo extendido de Euclides. Como la perpectiva es computacional, todos los ejemplos de prueba para los algoritmos se harn con Wolfram Mathematica.
1.2
Primitivas en general
Para efectos computacionales, una funcin y = f ( x ) tiene primitiva si existe una funcin g( x ) tal que
g0 ( x ) = f ( x ). Como es usual, se escribe
Z
f ( x ) dx = g( x ) + constante
Interludio: Primitivas en anlisis matemtico

Toda funcin tiene una primitiva?. Aunque estamos acostumbrados a una clase de
funciones que si tiene primitiva, en realidad las primitivas son escasas (en el mismo
sentido que "casi" no hay racionales en R ). Histricamente la integracin fue definida
como el proceso inverso de la derivada, es decir, la integral de f es entendida como una
funcin F tal que F 0 = f . Desde el punto de vista del anlisis matemtico y algunas
aplicaciones en fsica e ingeniera, esto es muy limitado porque hay funciones que no son
una derivada (de otra funcin). Las funciones continuas sin son derivadas (de alguien),
pero las funciones con "discontinuidades de salto" no. Alrededor de 1850 A. Cauchy y B.
Riemman construyen un concepto de integral independiente del concepto de derivada,
usando la nocin de "rea bajo la curva". Finalmente un funcin acotada en [ a, b] es
integrable en el sentido Riemann si su conjunto de discontinuidades tiene medida cero (un
0
conjunto numerable de discontinuidades). Hay funciones
R 0 f cuya derivada f tienen una
cantidad no numerable de discontinuidades, as que f = f no aplica, al menos para la
integral de Riemann (en otra nocin de integral, llamada "generalizada de Riemann", esta
identidad si aplica [1, Seccin 8.1]).
Una vez que la derivada y la integral han sido definidas de manera independiente, el
"Teorema fundamental del clculo" revela la relacin inversa entre ambos conceptos. En
realidad nos interesa "el segundo" teorema fundamental
del clculo: Si f : [ a, b] R es
Z
x
integrable (en el sentido Riemann), entonces F ( x ) =

continua en x = c [ a, b], entonces
F 0 (c)
= f ( c ).
f ( x ) dx es continua y si f es
Como las funciones continuas son integrables, tenemos: Si f : [ a, b] R es continua,

entonces
f tiene una primitiva (es derivada de alguien). Si x [ a, b] , la funcin
Z
x
F(x) =
f ( x ) dx es derivable y satisface F 0 ( x ) = f ( x ) en [ a, b].
No toda funcin es derivada (de alguien). El teorema del valor intermedio para derivadas
(Teorema de Darboux) dice que una funcin con una discontinuidad de salto no puede
ser derivada de otra funcin. Por ejemplo, la funcin f ( x ) = signo(x) toma los valores
1, 0, 1 si x < 0, x = 0 o x < 0 , respectivamente; como el conjunto de discontinuidades
es finito, es Riemann integrable, pero no es derivada de otra funcin. Esta es una de las
bondades de que la integral no dependa de la derivada.
Como una curiosidad adicional, las funciones continuas tienen primitiva... pero la mayora
de funciones continuas no tienen derivada en ningn punto!. Raro?. Esto es una consecuencia del "Teorema de Categora de Baire". Desde el punto de vista de este teorema,
los conjuntos "flacos" son los conjuntos "nunca densos" (su complemento es denso).
Los conjuntos que son unin contable de este tipo de conjuntos se llama "magro" o de
"primera categora". Los conjuntos de segunda categora son los que no son de primera
categora, son los conjuntos "gordos", como R . Resulta que las funciones continuas que
son derivables al menos en un punto es un conjunto de "primera categora", es decir,
las funciones diferenciables son excesivamente atpicas en general ([1, Seccin 8.2]), algo
como lo que pasa con los familiares nmeros racionales y los irracionales en R . Como
consecuencia, las primitivas son funciones bastante raras.
Desde el punto de vista computacional, un comando como "Integrate[f,x]" devuelve una funcin
(si tiene xito) g tal que g0 = f , y en general no se especifica el dominio. Como es sabido, las primitivas de f difieren solo en una constante y pueden tener aspecto diferente, as que pueden aparecer
diferencias de apariencia.
Ejemplo 1.3
Vamos a usar Mathematica para calcular algunas primitivas. Este software asume que algunas
funciones son de variable compleja por lo que a veces hay que especificar que la variable x es
real.
a) Primitiva de f ( x ) =
1
(en el software Log[ x ] corresponde a ln x ).
x
In[]:= Integrate[1/x, x]
In[]:= D[Log[x],x]
Out[]=
Out[]= 1/x
Log[x]
b) Tambin podemos calcular g( x ) =
Z x
1
1
dt
In[]:=
F=Integrate[1/t, {t,-1,x}, Assumptions -> {t Reals, x Reals}]
Out[]= ConditionalExpression[Log[-x], -1 < x < 0]
In[]:= D[F,x]
Out[]= ConditionalExpression[1/x,-1<x<0]
1.3
Primitivas elementales
Para efectos computacionales, una funcin f tiene primitiva elemental si existe una funcin g elemental
tal que g0 = f . Como es usual, se escribe
Z
f ( x ) dx = g( x ) + constante
Las funciones "elementales" son las funciones del clculo usual: Funciones racionales, funciones algebraicas (explcitas e implcitas), exponenciales y logartmicas (estas son las "trascendentales elementales") y las fuciones definidas por uan combinacin finita de estas tres clases anteriores.
q
2
Por ejemplo, tan(ex ) ln(1 + x ) + x x + ln(arcsen x ) es una funcin elemental. La expresin x x
se interpreta como ex ln x .
Las funciones racionales se definen como el resultado de un nmero finito de adiciones y multiplicaciones en la variable x y pueden ser expresadas como un cociente de polinomios en x (con exponentes
1 + ex
1 + e x
enteros y coeficientes en R o C). Por ejemplo,
es una funcin racional de ex y
es
2x
1 2e
1 + 3e2x
x
1+e
una funcin racional de ex si se reescribe como x
.
e + 3e3x
Las funciones algebraicas son funciones y que satisfacen una ecuacin del tipo
yn + f n1 yn1 + .... + f 0 = 0
p
donde cada f i es una funcin racional. Por ejemplo, y = x + x + x es algebraica porque satisface
la ecuacin y4 4xy2 4xy x = 0.
Esta "definicin" de funcin elemental es suficiente para los propsitos de este artculo.
Si admitimos funciones de variable compleja, las funciones elementales se reducen a funciones racionales,
exponenciales y logartmicas y las combinaciones algebraicas mencionadas. Sea z = x + i y. La exponencial compleja se define como
ez = ex (cos y + i sen y) para todo z C

Si x = 0 tenemos, ei y = cos y + i sen y, y de aqu deducimos sen y =
ei y ei y
2ii
y cos y =
ei y + ei y
La extensin compleja de estas funciones es,

sen z =
ei z ei z
2ii
y cos z =
ei z + ei z
2
, para todo z C.
Usando log z = ln |z| + i arg(z), z 6= 0; las funciones trigonomtricas inversas se obtienen como

1
log i z + (1 z2 )1/2
i

1
b) arccos(z) = log z + (z2 1)1/2
i
a) arcsen(z) =
c) arctan(z) =
1.4
1
z +i
i
z +i
log
= log
.
2ii
z i
2
z i
El algoritmo de Johan Bernoulli (1703)
El problema de calcular primitivas de funciones racionales es viejo, tanto como

el problema de calcular la derivada de una funcin. Segn Ostrogradsky ([18,
1845]) Newton prefera desarrollar estas funciones en serie e integrar trmino a
trmino (para evitar funciones trascendentes). Leibniz calculaba primitivas de funciones racionales obteniendo una factorizacin irreducible del denominador sobre
los reales, y luego integrando la descomposicin en fracciones parciales. Esto no
inclua el manejo de todos los casos en los cuales los denominadores son cuadrticos.
En el siglo XVIII Johan Bernoulli (1703) perfeccion el mtodo de descomposicin en fracciones parciales y
complet el mtodo de Leibniz ([20, pg 353]). La prueba definitiva del algoritmo solo se pudo completar hasta que la prueba rigurosa del teorema fundamental del clculo fue establecida por Gauss y
Weiertrass. Este parece ser el algoritmo ms viejo que se ha registrado sobre este tema. Es el mtodo
que se ensea hoy en los textos de clculo. El problema principal del mtodo es, por supuesto, la factorizacin completa del denominador. Ostrogradsky ([18]) present en 1845 un nuevo algoritmo que
calcula la parte racional de la integral sin factorizar. Sin embargo este mtodo solo se conoca en Rusia.
En 1872 Hermite ([19]) public un algoritmo diferente para calcular la parte racional de la integral
tambin sin factorizar (ver seccin 1.12).
El algoritmo de Bernoulli (1703). Este mtodo no es computacionalmente eficiente debido al costo de
factorizar en R[ x ] pero tiene importancia terica.

Supongamos que f = P/Q con P, Q R[ x ]. Despus de dividir los polinomios (divisin con resto),
podemos suponer que f = C + A/Q con C, A, Q R[ x ], mcd( A, Q) = 1 y grado ( A) <grado ( Q). Se
factoriza Q como
n
i =1
j =1
Q = ( x i ) ni ( x 2 + b j x + c j ) m j
con las funciones cuadrticas irreducibles. La descomposicin parcial de f es
ni
j
m
Bjk x + Cjk
Aik
+
( x )k ( x2 + b x + c )k
i
j
j
j =1 k =1
i =1 k =1
f =S+
El algoritmo para calcular la integral es ([17, pg,36]),

a) Calcular
C (no ofrece problema, pues C es un polinomio).
b) Calcular cada integral
Aik ( x i )1k
Aik
(1 k )
=
( x i )k
Ai1 log( x i )
si
k > 1,
si
k=1
c) Calcular las integrales con los denominadores cuadrticos. Recordemos que b2j 4c j < 0. Sea
j = 4c j b2j > 0. Hay dos casos.
Si k = 1,
2C
b
B
Bj1
2x
+
b
Bj1 x + Cj1
j
j1
j j1
q
=
log( x2 + b j x + c j ) +
arctan q
2
( x2 + b j x + c j )
j
j
(1.1)
Si k > 1, se usa la siguiente frmula recursiva, hasta que k 1 = 1, y se integra como el caso
anterior,
Bj1 x + Cj1
( x2 + b j x + c j )k
Z
(2Cjk b j Bjk ) x + b j Cjk 2Cj Bjk
(2k 3)(2Cjk b j Bjk )
+
2
k
1
( k 1) j ( x + b j x + c j )
( k 1 ) j ( x 2 + b j x + c j ) k 1
(1.2)
Descomposicin en fracciones parciales?. La descomposicin en fracciones parciales se usa en inte-
gracin de funciones racionales y en el clculo de la inversa de la transformada de Laplace. Pero "una

descomposicin en fracciones parciales" depende de la factorizacin que se tenga. Hay varios mtodos de
descomposicin en fracciones parciales basados por ejemplo en la descomposicin "libre de cuadrados
del denominador" (ver [17, Seccin 2.8 ],[25]).
Casos simples. Para polinomios que factorizan de manera simple, hay varias tcnicas rpidas y muy
eficientes para obtener los coeficientes de la descomposicin en el estilo del algoritmo de Bernoulli (ver
[21]). Por ejemplo,
1
1
a)
=
[ D ( x ) + a ][ D ( x ) + b ] b a
b) Sea
Q( x ) = x2 + bx + c
1
1
D(x) + a
D(x) + b

, a 6= b
1
1
con Q( a) 6= 0. Entonces
=
2
Q( a)
( x a)( x + bx + c)
1
x+a+b
2
xa
x + bx + c
c) Si P( x )/Q( x ) es una fraccin propia con la cuadrtica Q( x ) = ax2 + bx + c irreducible, entonces

P( x )
A2 Q0 ( x ) + B2
A1 Q0 ( x ) + B1
An Q0 ( x ) + Bn
=
+
+
...
+
Q( x )
Qn ( x )
( ax2 + bx + c)n
Q2 ( x )
d) Todas las races simples: Cuando Q( x ) tiene n races simples 1 , 2 , , n y grado P < grado Q,
se puede usar los residuos en la expansin de Laurent de la funcin racional P/Q en los polos i :
n
P(k )
P( x )
ck
=
con ck = 0
Q ( x ) k =1 x k
Q (k )
Este mtodo se conoce como "mtodo de los residuos" (por su origen el clculo de residuos en
funciones de variable compleja) y se puede generalizar a races de multiplicidad m > 1 (ver [17,
pg 54]). Si mcd( Q, Q0 ) = 1 entonces todas las races de Q son simples1 .
Ejemplo 1.4
Calcular
1
.
3 + x + x 2 + x 3
1
3 + x + x 2 + x 3
1
( x 1)( x2 + 2x + 3)

Z
1
1
x+3
=
2
, usando la frmula b)
6 x 1 ( x + 2x + 3)

1
1
4
2x + 2
=
log( x 1) log( x2 + 2x + 3) arctan
6
2
8
8
(se us la frmula 1.1, con 1 = 8)
Ejemplo 1.5
Calcular
x+1
.
15 13x + 3x2 + x3
(* Seguimos el algoritmo usando Mathematica *)

Clear[P, Q, Qp]
P[x_] = x + 1;
Q = -15 - 13 x + 3 x^2 + x^3;
(* Derivada de Q *)
Q tiene una raz de multiplicidad m > 1 entonces Q = ( x )m R( x ) . Por lo tanto Q0 = m( x )m1 R( x ) + ( x )m R( x )
con lo cual mcd( Q, Q0 ) 6= 1.
1 Si
10
Qp[x_] = D[Q, x];

(*Q solo tiene races simples *)
PolynomialGCD[Q, Qp];
(* = 1 *)
(* Clculo de las 3 races con Solve *)
{a1, a2, a3} = Flatten[x /. {Solve[Q == 0, x]}]; (*Races
(* Integrar *)

Z
P[a1]
P[a2]
P[a3]
1
1
1
+
+
dx
Qp[a1] x - a1
Qp[a2] x - a2
Qp[a3] x - a3
Z
1
1
1
-1
+ 0 +
dx *)
(* =
8 x+5
8 x-3
= {5, 1, 3}*)
Out[]= 1/8 Log[-3 + x] - 1/8 Log[5 + x]
Caso de polinomios irreducibles en Q[ x ]. Hay muchos polinomios de apariencia sencilla que son ir-
reducibles en Q[ x ] o cuya factorizacin en Q[ x ] no se adapta al algoritmo de Bernoulli. Algunos

ejemplos los podemos encontrar en [22]. Por ejemplo,
a) x5 + x + 1 = ( x2 + x + 1)( x3 x2 + 1)
b) x3 x2 + 1
c) x7 + x + 1
d) x6 5x4 8x3 2x2 + 2x + 1
No todo est perdido. Como veremos ms adelante, integrar funciones racionales no requiere siempre
obtener la factorizacin del denominador en el estilo del algoritmo de Bernoulli. Hay reducciones que
no ocupan factorizar y simplifican, muchas veces, el clculo de las primitivas.
Ejemplo 1.6
Consideremos la integral racional
emos
In[]:=
x7
1
. Al hacer el clculo con Mathematica obten+x+1
1
dx
x7 +x+1
Out[]= RootSum[1 + #1 + #17 &, Log[x - #1]/(1 + 7 #16 ) &]

Este Q = x7 + x + 1 no factoriza en Q[ x ] y tampoco las races se pueden expresar en trminos
de radicales en R. Pero sus races son simples pues mcd( Q, Q0 ) = 1; as que podemos dar una
respuesta simblica.
11
Si las races de Q son 1 , 2 , , 7 , entonces, aplicando el mtodo de los residuos,

Z
1
dx
x7 + x + 1
Q0 (k ) x k dx
k =1
7
log( x k )
+C
7k 6 + 1
k =1
Esta es la interpretacin de la salida "extraa" que nos da Mathematica. Esta frmula es til
para aproximar integrales definidas una vez calculadas la races aproximadas ([5]).
En la prctica computacional, no se factoriza el denominador totalmente para obtener la descomposicin en factores parciales, en vez de eso se usa una factorizacin "libre de cuadrados". En todo caso
no siempre es posible factorizar: Como se sabe, para polinomios de grado mayor o igual a cinco no
siempre es posible encontrar las races del polinomio en trminos de radicales.
Aunque Q( x ) = x6 5x4 8x3 2x2 + 2x + 1 no factoriza en Q[ x ], se puede hacer reducciones que

solo ocupan operaciones con polinomios. En el captulo sobre "integracin de funciones racionales"
Z
6x5 + 6x4 8x3 18x2 + 8x + 8
(seccin 1.11) vamos a ver como se calcula la integral
. Por ahora,
x6 5x4 8x3 2x2 + 2x + 1
veamos en el siguiente ejemplo, como va a quedar el resultado.
Ejemplo 1.7
Vamos a calcular
Rothstein/Trager:
6x5 + 6x4 8x3 18x2 + 8x + 8

con el mtodo de Hermite y el mtodo de
x6 5x4 8x3 2x2 + 2x + 1
El clculo solo requiere operaciones con polinomios y calcular las races del polinomio
R(z) = (2 2z + z2 )3 (ver seccin 1.11).
R 6x5 + 6x4 8x3 18x2 + 8x + 8

x6 5x4 8x3 2x2 + 2x + 1
1+
h
i

3 Log 1 + x + 3x + 3x2 x3
h
i

3 Log 1 x + 3x + 3x2 + x3
12
1.5
Principio de Laplace, 1812.

Laplace, en el primer libro de la "Thorie analytique des probabilits" [7, pg 5],
estableci en 1812 por primera vez una teora consistente de funciones "elementales" y la integracin en trminos finitos, en forma de esbozo general. Por ejemplo,
Z
dx
p
Laplace indica que con su razonamiento, que la integral
no
1 + x2 + x4
tienen primitiva elemental, sin embargo Liouvlle seala que "la demostracin no
se ha publicado y nosotros ignoramos de qu naturaleza esta demostracin podra
ser" ([23, pg 145]).
En todo caso, Liouville reconoce que este esbozo sobre el tema de integracin en trminos finitos es
notable porque lo cita casi completo en su "Premier Memoire sur la dtermination des intgrales dont
la valeur est algbrique" ([23, pg 144]).
La cita de Lioville que aqu presentamos, en traduccin libre, sigue la versin en ingls de Luetzen ([20,
pgs 357-58]): Habiendo introducido las funciones algebraicas, exponencial y logartmicas, Laplace contina:
Estas cantidades son esencialmente distintas: la exponencial a x , por ejemplo, nunca puede ser idntica
a una funcin algebraica de x. De hecho, cada funcin algebraica es reducible a una serie descendente
0
de la forma kx m + k0 x nn + ...; pero es fcil demostrar que cuando a es mayor que la unidad y x es
x
infinita, a es infinitamente mayor que kx n , por muy grande que k y n puedan ser. Del mismo modo,
es fcil ver que cuando x es infinito, x es infinitamente mayor que k(log x )n .
Por lo tanto, las funciones exponenciales, algebraicas y logartmicas no pueden reducirse de una a la
otra; las cantidades algebraicas estn en el medio entre las exponenciales y los logaritmos, de hecho,
cuando la variable es infinita, podemos considerar los exponentes como infinitos en las exponenciales,
finito en las algebraicas e infinitamente pequeas en las cantidades logartmicas.
Adems, uno puede en principio establecer que una funcin radical de una variable no puede ser idntica a una funcin racional de la misma variable o con otra funcin radical. Por lo tanto, (1 + x3 )1/4
es esencialmente diferente de (1 + x3 )1/3 y de (1 + x )1/2 . Estos principios, que se basan en la naturaleza de las funciones, pueden ser de gran utilidad en la investigacin analtica, ya que nos indican
la forma que necesariamente deben tener este tipo de funciones que andamos buscando y mostrar
que esta forma no existe en un gran nmero de casos. Sin embargo, uno debe entonces estar absolutamente seguro de no omitir ninguna de las formas posibles. Por lo tanto, puesto que la diferenciacin
permite a la exponencial y las cantidades radicales subsistir y slo hace que las cantidades logartmicas desaparezcan cuando estn multiplicadas por las constantes, se puede concluir que la integral de
una funcin diferencial no puede incluir ninguna otra exponencial y ningn radical que los que ya
figuran en el presente funcin. De este modo, me he dado cuenta que no se puede obtener la integral
R
dx2
como una funcin finita explcita o implcita. Del mismo modo, he demostrado que
4
1+x + x
las ecuaciones diferenciales en derivadas parciales lineales de segundo orden y de tres variables, no
pueden ser integradas de forma finita, y esto me ha llevado a un mtodo general de integrar stas
en esta forma cuando sea posible. En los otros casos, no se puede obtener una integral finita sino por
medio de integrales definidas. [Laplace 1812, pp 4-5]
13
Ms adelante va a ir quedando claro como la "forma" del integrando y de la primitiva es crucial en la

decisin de si una funcin tiene o no primitiva elemental.
Acerca de esto, en esta primera memoria de Liouville sobre el tema ([23]), l inicia probando que la
Z
dx
integral
no tiene primitiva elemental. Luetzen ([20, pgs 370-372]) indica que el camino
1 + x4
que Liouville sigui (de acuerdo a su manuscritos) fue el de expresar esta integral en "la forma"
dx
1 + x4
M + P 1 + x4
=
Q
con M, P y Q polinomios en Z[ x ]. La justificacin elude a los principios de Laplace: La primitiva solo

podra tener el mismo tipo de radicales que aparece en el integrando (no cbicas, ni de otro orden, ni
exponenciales ni logaritmos, etc.). Luego prueba que es imposible encontrar M, P y Q que satisfagan
esta identidad.
1.6
Teorema de P. Chebyshov, 1853.

Nuestro inters est ahora en un teorema de P. Chebyshov porque, aunque es posterior (y complementaria) de la teora establecida por Liouville, es un teorema de
fcil aplicacin. P. Chebyshov estaba interesado en casos especiales de funciones
algebraicas y su trabajo sobre la integracin de este tipo de funciones est conectado de manera muy cercana con el trabajo de Abel, Liouville y Ostrogradsky. En el
artculo "Sur lintegration des diffrentielles irrationelles" (1853, [3]) P. Chebyshov
resuelve el problema de obtener la parte logartmica de la integral
Z
P( x ) dx
p
,
Q( x ) m R( x )
donde P, Q y R son polinomios y m es un entero positivo. Pero el artculo es ms bien conocido

por la solucin completa del problema de la integracin del binomio diferencial
Z
x m ( a + bx n ) p dx.
Ya Goldbach y Euler haban mostrado que este tipo de integral es expresable en trminos de funciones
elementales en los casos en que p es un entero, (m + 1)/n es un entero o (m + 1)/n + p es un entero.
P. Chebyshov demostr, con mtodos analticos, que estos son los nicos casos en los que esta integral
se puede expresar en trminos de funciones elementales.

14
Figura 1.1: P.L. Tchebychef. "Sur lintegration des diffrentielles irrationelles" (en Oeuvres, 1899), [3]
Para la presentacin de este teorema definitivamente lo mejor es seguir la exposicin del libro "Clculo
diferencial e Integral" de N. Piskunov [10, pgs. 375-378].
Teorema 1.1 (Chebyshov, 1853).

Si m, n y p son nmeros racionales, entonces la integral
Z
x m ( a + bx n ) p dx
se puede expresar en trminos de funciones elementales si y slo si al menos uno de los nmeros
m+1
m+1
o
+ p , es entero.
p,
n
n
En el caso de que la integral se puede expresar en trminos de funciones elementales, la manera de

obtener la primitiva es va sustituciones [10, pgs. 375-378]. Estas sustituciones convierten la integral
en una funcin racional que, por el algoritmo de J. Bernoulli, siempre tiene primitiva elemental. La
otra parte de la prueba del teorema de Chebyshov no es elemental, se puede ver en [8, pgs. 37-39] y
en espaol en [9, pgs. 45-48].
Aplicacin. El teorema se puede usar para decidir si una funcin de la forma x m ( a + bx n )
diferencial") tiene primitiva elemental y, tambin para el clculo.

1 1 1
u n du y
n
Z
Z
1
m+1
x m ( a + bx n ) p dx =
uq ( a + bu) p du con q =
1
n
n
Para calcular, primero se hace la sustitucin x = u1/n , entonces dx =
("binomio
a) Caso I: Si p es entero y q es racional. Si q =

b) Caso II: Si
c) Caso III. Si
Si q =
h
, se hace la sustitucin u = td .
d
h
m+1
es entero y p es racional. Si p = , se hace la sustitucin a + bu = td .
n
d
m+1
+ p es entero y p es racional. Primero se transforma la integral,
n

Z
Z
1
a + bu p
m+1
1
uq ( a + bu) p du =
uq+ p
du con q =
1
n
n
u
n
h
a + bu
, se hace la sustitucin
= td .
d
u
Ejemplo 1.8
1
dx
=
x 1/3 (1 + x )1/2 dx. En este caso m = ,
3
3
x x+1
m+1
2
m+1
2 1
= y
+ p = . Como ninguno de estos nmeros es entero, la integral no se
n
3
n
3 2
Z
Consideremos la integral
puede expresar en trminos de funciones elementales.
Ejemplo 1.9
dx
x 2 )3
x 3/2 (4 x )3/2 dx. Como
m+1
1
+p=
n
2
(4x
3
= 2, la integral s se puede expresar en trminos de funciones elementales. Usando las
2
indicaciones de sustitucin mencionadas, se obtiene
Z
x 1/3 (1 + x )1/2 dx
1
x2
p
+K
4 x (4 x )
Ejemplo 1.10
a)
sen x dx no es elemental: Si u = sen x obtenemos
15
16
sen x dx =
u1/2 (1 u2 )1/2 du
que, de acuerdo al teorema de Chebyshov, no se puede expresar en trminos de funciones

elementales.
Z
b)
tan x dx es elemental. Si u = tan x , u2 = tan x, la integral queda
Z
tan x dx =
2u2 (u4 + 1)1 du
que si pasa el test de Chebyshov pues p = 1.

1
(1 x n )1/k dx. Como m = 0 y p = , entonces la integral es elemental
k
1
1
1
es entero o 1/k es entero o + es entero, es decir, solo si k = 1, o n = 1
solo si
n
n k
1
1
o k = n = 2 o n = k (en otro caso + < 1 ).
n
k
k
Si k > 0 , esto solo nos deja 1 x n , k 1 x, k 1 1/x, 1 x2 y 1 x k con primitiva
elemental.
c) Consideremos
Longitud de arco. La longitud (de arco) de una curva suave de ecuacin y = f ( x ) desde x = a hasta
x = b, se calcula con la frmula

s=
Z bq
a
1 + ( f 0 ( x ))2 dx.
Como se ve, son pocas las integrales de este tipo que se pueden calcular usando el teorema fundamental del clculo por la escasez de primitivas elementales.
Por ejemplo, si f ( x ) =
1.7
1
entonces, s =
x
Z p
1 + x 4 dx no es elemental.
Primitivas elementales y funciones inversas
Usando integracin por partes, la integral de una funcin f se puede expresar en trminos de la integral de su inversa f 1 . Adems de ser un recurso de clculo, tambin se puede usar como mtodo de
decisin en el caso de que se conozca que f o f 1 tiene o no tiene primitiva elemental.
Si f es invertible en [ a, b] y si f y f 1 tienen derivada continua, entonces aplicando integracin por
partes con u = f ( x ) y dv = dx, se obtiene
f ( x ) dx
= x f (x)
= x f (x)
x f 0 ( x ) dx
f 1 ( f ( x )) f 0 ( x ) dx
Aplicando la sustitucin t = f ( x ) en la integral de la derecha, obtenemos la primitiva G (t) =

Sustituyendo queda
Z
17
f 1 (t) dt .
= x f ( x ) G ( f ( x ))
f ( x ) dx
A veces, abusando del lenguaje, se escribe
f ( x ) dx = x f ( x )
Z f (x)
f 1 (y) dy
Ejemplo 1.11
Z
a)
ln x dx es elemental pues
Z ln x
b)
As,
ln x dx
= x ln x
ey dy es elemental:
ln x dx
ln x
= x ln x
Z ln x
ey dy
ey dy
ln x dx no es elemental porque (como veremos ms adelante)
e x dx no es
elemental.
c)
1
x
dx =
ln x
ln x
As,
Z 1/ ln x
e1/y dy
e1/y dx no es elemental porque (como veremos ms adelante)
elemental.
1
dx no es
ln x
18
1.8
El principio de Liouville
El estudio sistemtico del problema de decidir cuando una integral tiene primitiva
elemental inicio en el siglo XIX. Una dcada despus de Laplace, J. Liouville
establece y prueba un teorema ms general y preciso: Hablando informalmente,
Liouville establece que si la primitiva de una funcin elemental f es una funcin
elemental, entonces esta primitiva se puede expresar usando solo las funciones
que aparecen en el integrando y una combinacin lineal de logaritmos de tales
funciones. El teorema de Liouville nos da un criterio de integracin de funciones
elementales en trminos de la solucin de una ecuacin diferencial de primer orden
con una funcin racional.
Ejemplo 1.12
Z
a)
Z
1 ix
2e
+ 12 ei x dx = 21 i ei x + 21 ei x + K
p
1 x2 + i x ) + K
1 x2
Z
p
p
p
c)
log( x + x2 1) dx = x log( x + x2 + 1) x2 1 + K
b)
dx = i log(
Para enunciar el teorema "fuerte" de Liouville necesitamos ubicarnos con las funciones algebraicas y
trascendentes.
Funciones algebraicas. Recordemos que una funcin algebraica (ver por ejemplo [13]) y = y( x ) es una
funcin que satisface una ecuacin de la forma

Rn ( x )yn + Rn1 ( x )yn1 + + R0 ( x ) = 0,
con Rn ( x ) 6 0.
donde cada Ri ( x ) una funcin racional. Multiplicando a ambos lados por el mximo comn denominador de los denominadores nos queda
P( x, y) = an ( x )yn + an1 ( x )yn1 + + a0 ( x ) = 0,
con an ( x ) 6 0.
donde cada ai ( x ) es un polinomio con coeficientes en algn conjunto. Desde el punto de vista computacional interesa que ai ( x ) Q[ x ].
Las funciones racionales son algebraicas: Si y =
19
P( x )
entonces y es algebraica porque satisface la
Q( x )
ecuacin
Q( x )y P( x ) = 0.
Una expresin que puede ser formada a partir de x y ciertas constantes aplicando un nmero finito de
veces las operaciones de suma, resta, producto, divisin y extraccin de races, se llama una funcin
algebraica explcita. Por ejemplo,
s
r
b
a + x4
1
3
y=a 2 +
+
a + bx
+
x
1 + x3
a + b + dx3
puede ser obtenida a partir de x, a, b, d aplicando las operaciones indicadas.

Se sabe que las funciones algebraicas explcitas satisfacen la definicin de funciones algebraicas, pero
como las ecuaciones polinomiales de grado > 4 no pueden ser, en general, resueltas por radicales, no
todas las funciones algebraicas puede ser explcitas.
La suma, resta, producto, divisin y "extraccin de races" de funciones algebraicas vuelve a ser una
funcin algebraica. Si y es algebraica y es invertible, su inversa es algebraica. Adems las funciones
algebraicas son continuas en su dominio y tienen un nmero finito de ceros.
Funciones trascendentes. Las funciones que no son algebraicas se llaman trascendentes. Por ejemplo,
sen x, cos x, tan x, ln x y sus inversas son trascendentes. Veamos un par de pruebas.
sen x es trascendente: Si fuera algebraica, tendra un nmero finito de ceros, pero sen(k ) = 0 para
todo k Z. Tambin se puede probar que sen x es trascendente usando el hecho de que P( x, sen x ) es
una funcin continua (ver apndice ??).
y = ln x es trascendente (por tanto tambin su inversa y = e x ). La prueba de este hecho usa una idea
que vamos a encontrar varias veces. La idea de la prueba es esta: Si ln x satisface la ecuacin polinomial, ln x debe ser racional y eso no puede pasar. Iniciamos probando que ln x no puede ser una
funcin racional.
Si ln x =
1
QP0 PQ0
P( x )
con P y Q coprimos, entonces =
, es decir,
Q( x )
x
Q2
Q2 = x ( QP0 PQ0 ).
Por tanto x es un factor de Q. Sea Q( x ) = x m R( x ), con R(0) 6= 0 y m 1. Sustituyendo y dividiendo

por x m obtenemos,
x m R2
= x [ x m RP0 P(mx m1 R + x m R0 )]
= xRP0 mPR xPR0
Pasando a restar y sacando x a factor se tiene,
20
mP( x ) R( x )
= x [ RP0 PR0 x m1 R2 ]
Entonces x divide a PR, como no divide a R, entonces divide a P, lo que no puede ser pues P y Q
son coprimos.
Ahora supongamos que y = ln x es algebraica y que
a N ( x )(ln x ) N + a N 1 ( x )(ln x ) N 1 + + a0 ( x ) = 0,
con a N ( x ) 6 0,
donde la ecuacin se escogi con el menor grado N posible. Entonces
(ln x ) N +
a N 1 ( x )
a (x)
(ln x ) N 1 + + 0
= 0,
a N (x)
a N (x)
con N 2
Derivando obtenemos,
N
(ln x ) N 1 +
x
a N 1 ( x )
a N (x)
0
(ln x ) N 1 + = 0, con N 2
que es un polinomio de grado menor que N , por lo que lo nico que queda es que los coeficientes
(polinomios) se deben anular para satisfacer la ecuacin. En particular,
N
+
x
Pero integrando obtenemos que ln x =
no puede ser como ya probamos.
a N 1 ( x )
a N (x)
0
= 0.
a N 1 ( x )
, es decir, ln x sera una funcin racional. Pero esto
Na N ( x )
El teorema fuerte de Liouville es enunciado aqu en trminos de una funcin f ( x, y1 , ..., yn ) que se expresa en trminos de funciones elementales yi ( x ) pero que requiere que las derivadas yi0 sean algebraicas.
Teorema 1.2 (Lioville, 1835).

dyi
es funcin algedx
f ( x, y1 , y2 , ..., yn ) dx es elemental si y slo si
Sea f ( x, y1 , y2 , ..., yn ) una funcin algebraica donde cada yi = yi ( x ) y cada

braica de x, y1 , ..., yn . Entonces
Z
f ( x, y1 , y2 , ..., yn ) dx = U0 +
Cj ln(Uj )
j =1
donde U0 es funcin algebraica de x y cada Cj es constante y cada Uj1 es funcin algebraica

de x, y1 , ..., yn .
21
dy
Si f ( x, y1 , y2 , ..., yn ) es una funcin racional y cada i es funcin racional de las funciones
dx
x, y1 , ..., yn , entonces U0 es una funcin racional de x y cada Uj , j 1, deben ser funcin
racional de x, y1 , ..., yn .
La notacin es un poco fuerte, pero es en realidad sencilla.
Ejemplo 1.13
Sea y1 = ex , y2 = sen x y y3 = cos x .
Sea f ( x, y1 , y2 , y3 ) = 2x +
y3
2 + 2y1
+
. En este caso, f es racional
y2
y1 + x
dy1
dy2
dy3
= y1 ,
= y3 , y
= y2 , se cumple que f ( x, ex , sen x, cos x ) es algebraica y
dx
dx
dx
dyi
son funciones algeraicas de x, y1 , y2 , y3 .
que las derivadas
dx
Como
El teorema de Liouville dice que

Z
2x +
f ( x, y1 , y2 , y3 ) dx es elemental si y slo si
y3
2 + 2y1
dx = U0 +
+
y2
y1 + x
Cj ln(Uj )
j =1
con U0 = U0 ( x ) algebraica y las Cj son constantes y las Uj1 son funciones algebraicas de
x, y1 , , y2 y3 .
En este caso la integral efectivamente es elemental y
Z
2x +
cos x
2e x + 2
+
dx = x2 + 2 ln ( x + e x ) + ln(sen( x )) + K
sen x
x + ex
y como se ve, U0 = x2 , U1 = x + y1 = x + e x y U2 = y2 = sen( x ).
El teorema de Liouville tiene valor terico y prctico y ha sido la base de el trabajo desarrollado en
la rama de integracin "en trminos finitos". Hay dos casos especiales del teorema que nos dan herramientas para decidir si ciertas funciones importantes tienen o no primitiva elemental.
22
1.9
Caso especial del teorema de Lioville, 1835.
x
ZConsideremos la funcin elemental f ( x )e
con f Q( x ). De acuerdo al teorema de Liouville,
x2
f ( x )e dx es elemental si y slo si
f ( x )ex dx = U0 +
Cj ln(Uj )
(1.3)
j =1
2
con U0 ( x ) una funcin racional y cada Uj1 = Uj ( x, ex ).

Derivando a ambos lados en (1.3) obtenemos
f ( x )ex = U00 +
Uj0
Cj Uj
j =1
Si U0 = P/Q se obtiene, con un poco de esfuerzo, una ecuacin diferencial que se debe satisfacer:
f ( x ) = U00 + 2xU0 . Esta ecuacin diferencial es la clave para demostrar, eso s de manera algebraica,
2
que f ( x )ex tiene o no tiene primitiva elemental. Tratar de resolver la ecuacin diferencial (va factor
integrante) no ayuda en nada.
En general,
Teorema 1.3 (Caso especial Liouville, 1835)

Si f y g son funciones racionales y g no es constante, entonces
Z
f ( x ) eg( x) dx es elemental si y slo si existe R( x ) racional, tal que f ( x ) = R0 ( x ) + R( x ) g0 ( x )
Hay una parte fcil en la prueba del teorema. Una direccin del teorema dice que si existe una funcin
racional R tal que f ( x ) = f ( x ) = R0 ( x ) + R( x ) g0 ( x ), entonces
Z
f ( x ) eg( x) dx = R( x )eg( x) .
En efecto, derivando obtenemos

d
R ( x )e g ( x )
dx
23
R0 ( x ) e g( x ) + e g( x ) R ( x )
= eg(x) ( R0 ( x ) + g0 ( x ) R( x ))
f ( x ) e g( x )
La otra direccin de la prueba es la que no es elemental.
Ejemplo 1.14
Z
ex dx no tiene primitiva elemental (y por tanto
En efecto, de acuerdo al teorema 1.3,

que
1
Sea R( x ) =
ln x dx tampoco, segn la seccin 1.7).
1 ex dx es elemental si y slo si existe R( x ) racional tal
R0 ( x ) + R( x ) 2x
(1.4)
P( x )
con mcd( P, Q) = 1. Entonces tenemos
Q( x )
1
P0 Q Q0 P
P
P0 Q Q0 P + 2xPQ
+
2x
=
Q
Q2
Q2
de donde
Q2
= Q0 P + P0 Q + 2xPQ
Como Q|( P0 Q + 2xPQ) entonces Q| Q0 P y como mcd( Q, P) = 1 se concluye que Q| Q0 . Pero la

derivada Q0 tiene menor grado que Q, por lo tanto la nica posibilidad es que Q sea constante.
Digamos Q q
Entonces 1.4 se puede reescribir como
P0 + 2xP
Pero esto es imposible pues P0 + 2xP tiene grado 1 mientras que la constante tiene grado
cero.
24
Ejemplo 1.15
Z
x2n eax dx, con n N y a 6= 0, no tiene primitiva elemental. La prueba sigue un argumento
similar al del ejemplo anterior (en [14] puede ver todos los detalles).
Ejemplo 1.16
Z
ex
dx no tiene primitiva elemental.
x
En efecto, de acuerdo al teorema 1.3,

que
1
x
Sea R( x ) =
1 x
e dx es elemental si y slo si existe R( x ) racional tal
x
R0 ( x ) + R( x )
(1.5)
P( x )
con mcd( P, Q) = 1. Entonces tenemos
Q( x )
1
x
P0 Q Q0 P
P
P0 Q Q0 P + PQ
+
=
2
Q
Q
Q2
de donde
Q( Q xP0 xP)
= xQ0 P
(1.6)
Si Q tiene una raz 6= 0 de multiplicidad k 1 entonces Q = ( x )k S( x ) con S() 6= 0.

Entonces, como P no comparte races con Q, la ecuacin 1.6 se puede reescribir como
( x )k+m N ( x ) = ( x )k1 D ( x ), con N () 6= 0, D () 6= 0 y m 0.
Como las multiplicidades son distintas, Q no tiene races no nulas. Por tanto, si Q tiene una raz,
debera ser = 0. Sea Q = Cx k con C constante y k 1. En este caso la ecuacin 1.6 se puede
reescribir como
( x )k+m N ( x ) = ( x )k1 D ( x ), con N () 6= 0, D () 6= 0 y m 0.
se puede reescribir como

x k+m N ( x )
x k D ( x ), con N (0) 6= 0, D (0) 6= 0 y m 1.
25
De nuevo las multiplicidades no coinciden, por lo que Q es constante, Q q. Sustituyendo en

1.6 nos queda
q = xP0 + xP
que es imposible pues el lado derecho tiene grado 1.
Se puede probar que otras integrales no tienen primitiva elemental usando identidades, sustituciones
o el mtodo de integracin por partes.
Ejemplo 1.17
1
dt no tiene primitiva elemental pues
ln t
x = ln t.
Z
a)
b)
ex
dx =
x
1
dt usando la sustitucin
ln t
ex ln x dx no tiene primitiva elemental pues, integrando por partes con u = ex y
dv =
1
dx , se obtiene
x
Z
1.10
ex
dx = ex ln x
x
ex ln x dx
Caso especial Liouville-Hardy, 1905.
El siguiente caso especial fue obtenido por G. Hardy ([6]) en 1905 y se aplica a funciones de la forma
f ( x ) ln x.
Teorema 1.4 (Liouville-Hardy, 1905)

Si f ( x ) es una funcin racional, entonces
f ( x ) ln x dx es elemental si y slo si existe una
funcin racional R( x ) y una constante C tal que

f (x) =
C
+ R0 ( x )
x
26
La prueba de este teorema tiene una parte fcil: Si f ( x ) =

mental. En efecto,
Z
f ( x ) ln x dx
=
=
=
Z
Z
C
+ R0 ( x )
x
C ln x
dx +
x
C
+ R0 ( x ) entonces
x
f ( x ) ln x dx es ele-
ln x dx
ln x R0 ( x ) dx
C ln2 ( x )
+ ln x R( x )
2
1
R( x ) dx, integrando por partes con u = ln x y dv = R0 ( x )dx.
x
R( x )
R( x )
dx es elemental porque
racional. Por lo tanto
x
x
f ( x ) ln x dx es elemental.
La otra parte de la prueba no es elemental.
Ejemplo 1.18
En el teorema 1.4, para probar que una integral es elemental, se debe encontrar dos cosas: una
constante C y una funcin racional R tal queda
R0 ( x ) = f ( x )
C
x
Integrando tenemos
R=
f ( x ) dx C ln x
por lo tanto, si C = 0 y f ( x ) es un polinomio, el teorema se cumple.

Z
P( x ) ln x dx es elemental si P es un polinomio.
Esto es algo fcil de verificar (y calcular) usando integracin por partes.
Ejemplo 1.19
Sea Q es un polinomio con todos sus ceros simples y no nulos, digamos Q( x ) = ( x 1 ) ( x n ),
Z
ln x
entonces
dx no es elemental.
Q( x )
27
n
1
ln( x i )
dx = C ln x +
, que no es una
Q( x )
Q 0 ( i )
i =1
funcin racional para ningn valor de C (recuerde que dijimos que cada i 6= 0).
Esto es as pues sera R = C ln x +
A la luz de este resultado,

Z
ln x
dx, etc.
x2 + 1
ln x
dx es elemental pero
x
ln x
dx no es elemental, ni tampoco
x+1
Ejemplo 1.20
Z
arcsec2 ( x ) dx no es elemental pues, aplicando integracin por partes,

Z
1.11
arcsec2 ( x ) dx = x arcsen2 ( x ) + 4
ln t
dt
t2 + 1
Algoritmo para Integrar una Funcin Racional
Como ya sabemos por el algoritmo de Bernoulli, toda funcin racional P/Q tiene primitiva elemental.
En el algoritmo que sigue, el clculo de esta primitiva no se hace buscando una factorizacin completa
para Q, que de por s puede ser imposible, sino que se divide la integral en una parte racional y una
parte logartmica. La parte racional se obtiene con el mtodo de Hermite y no necesita factorizaciones.
La parte logartmica de la integral se hace con el algoritmo de R/T y usa una factorizacin libre de
cuadrados y el clculo de races de un polinomio especial.
Este algoritmo aparece de manera muy similar en el caso general.
Conocimientos previos: Factorizacin libre de cuadrados. Desde el punto de vista computacional, para
integrar una funcin racional, se trata de factorizar el denominador pero de una manera menos demandante que la factorizacin total. Esta factorizacin se llama libre de cuadrados. Esta factorizacin
empaqueta los factores de Q que tiene la misma multiplicidad, en un solo factor y solo requiere
divisin de polinomios y clculo del mximo comn divisor.
Un polinomios es libre de cuadrados si no tiene factores repetidos. Si un polinomio tiene factores
repetidos, entonces estos factores vuelven a aparecer en la derivada: Si P = A2 B entonces P0 =
28
2AA0 B + A2 B0 = A(2A0 B + AB0 ). Por tanto, P es libre de cuadrados si mcd( P, P0 ) = 1.
Definicin 1.1
a) Un polinomios es libre de cuadrados si no tiene factores repetidos
b) Una factorizacin libre de cuadrados de P es una factorizacin P = ik=1 ( Pi )i (potencias
crecientes) donde cada Pi es libre de cuadrados y mcd( Pi , Pj ) = 1 si i 6= j. Algunos Pi
pueden ser 1 .
Teorema 1.5
Un polinomio P Q[ x ] es libre de cuadrados si y slo si mcd( P, P0 ) = 1
Prueba. (=) Si P tiene algn factor repetido, entonces P = Qn D con n > 1. Por tanto P0 = nQn1 D +
D 0 Q. Como Q es un factor comn entre P y P0 , esto contradice la hiptesis de que mcd( P, P0 ) = 1.
(=) Si P es libre de cuadrados (sin factores repetidos) y si P tiene la factorizacin irreducible
P = Q1 Q2 Qn con todos los Qi distintos,
entonces
P0 = Q10 Q2 Qn + Q1 Q20 Q3 Qn + ... + Q1 Q2 Qn1 Q0n
Supongamos que mcd( P, P0 ) > 1 y sea D Q[ x ] no trivial con D | P y D | P0 , entonces como D Q[ x ],
D debe ser un mltiplo de algn Qi . Por tanto podemos asumir que algn Qi divide a ambos P y P0 .
Entonces Qi divide todos los sumandos de P0 donde l aparece y como divide a P0 , deber tambin
dividir al sumando Q1 Qi0 Qn , por tanto debe dividir a Qi0 , pero esto es imposible pues Qi0 tiene
menor grado que Qi .
Note que hay una diferencia entre ser libre de cuadrados y la factorizacin libre de cuadrados.
Algoritmo para obtener la factorizacin libre de cuadrados. Para ver la manera en que se obtienen
los factores Pi libres de cuadrados en la factorizacin libre de cuadrados de P, hagamos un clculo pequeo: Supongamos que la factorizacin libre de cuadrados de P es P = P1 P22 P33 (hasta ahora
desconocida) entonces queremos calcular estos tres Pi . Primero calculamos P0 y factorizamos,
P0 = P2 P32 [ P2 P3 P10 + 2P1 P3 P20 + 3P1 P2 P30 ] .

P
Entonces A1 =mcd( P, P0 ) = P2 P32 y B1 =
= P1 P2 P3 . Por tanto, si C1 = mcd ( A1 , B1 ) = P2 P3 ,
mcd ( P, P0 )
obtenemos
B
P1 = 1
C1
De manera similar obtenemos P2 y P3 .

En general: Primero calculamos A1 = mcd ( P, P0 ) y B1 =
P
mcd ( P, P0 ).
Luego iteramos:
Desde i = 1 hasta completar el grado de P
Ci = mcd ( Ai , Bi )
A i +1 =
Ai
Ci
Bi+1 = Ci
Pi =
Bi
Ci
Ejemplo 1.21
Obtener la factorizacin libre de cuadrados de P( x ) = x8 2x6 + 2x2 1.
(*Usamos Wolfram Mathematica*)
P = x^8 - 2 x^6 + 2 x^2 - 1;
DP = D[P, x];
A1 = PolynomialGCD[P, DP];
B1 = PolynomialQuotient[P, A1, x];
C1 = PolynomialGCD[A1, B1] ;
P1 = PolynomialQuotient[B1, C1, x]
A2 = PolynomialQuotient[A1, C1, x];
B2 = C1;
C2 = PolynomialGCD[A2, B2];
P2 = PolynomialQuotient[B2, C2 , x]
A3 = PolynomialQuotient[A2, C2, x];
B3 = C2;
C3 = PolynomialGCD[A3, B3];
P3 = PolynomialQuotient[B3, C3 , x]
Entonces P( x ) = ( x2 + 1) 12 ( x2 1)3 .
(* P1 = x2 + 1 *)
(* P2 = 1 *)
(* P3 = x2 1 *)
29
30
1.12
Mtodo de Hermite para obtener la parte racional, 1872.
Una variante del algoritmo de Bernoulli consiste en usar una factorizacin completa del denominador
en C[ x ]. Para esta factorizacin, la descomposicin en fracciones parciales no lleva a la integral
f = v + ci log(ui )
i =1
donde v, u1 , ..., un C( x ) y cada ci C. En esta expresin, v se llama la parte racional de la integral y

la suma de logaritmos se llama la parte trascendental. El mtodo de Hermite calcula la parte racional
de la integral respecto a la factorizacin libre de cuadrados del denominador nicamente. La parte
logartmica requiere otro algoritmo. Por ejemplo,
x7 24x4 4x2 + 8x 8
1
3x
6x
dx = +
+
+
2
8
6
4
2
x (2 + x 2 )
2 + x2
x + 6x + 12x + 8x
|
{z
}
parte racional
Figura 1.2: Mtodo de Hermite, 1872
Pues bien, el mtodo de Hermite hace una reduccin del tipo

Z
P
C
=
+
Q
D
A
B
dx
x
31
con P, Q, A, B, C, D Q[ x ], grado( A) <grad( B) y B mnico y libre de cuadrados.

C/D es la parte racional de la integral mientras que
A/B es la parte logartmica (si esta parte no es
nula, requiere extensiones logartmicas).

Algortimo. El mtodo de Hermite procede como sigue: Primero se normaliza P/Q de tal manera que
mcd( P, Q) = 1. Esto se logra dividiendo P y Q por el mximo comn divisor de ambos. Tambin se
requiere que Q sea mnico. Luego aplicamos divisin euclidiana a P y Q y obtenemos L, R Q[ x ]
tal que P = Q L + R con R = 0 o grado( R ) <grado( Q ) . Obtenemos,
P
=
Q
Luego, para
L+
R
Q
R
calculamos la factorizacin libre de cuadrados de Q,
Q
k
Q = Qii
i =1
con cada Qi mnico, libre de cuadrados, mcd( Qi , Q j ) = 1 si i 6= j, y grado( Qi ) > 0.

Ahora calculamos la descomposicin en fracciones parciales de R/Q respecto a al factorizacin libre
de cuadrados:
k i R
R
ij
= j
Q i =1 j =1 Q
i
con Rij Q[ x ], grado( Rij ) <grado( Qi ) si grado( Qi ) > 0 y Rij = 0 si Qi = 1.

Z
k i
R
=
Q i =1 j =1
Rij
se puede reducir usando integracin por partes y el algoritmo extenj

Qi
dido de Euclides de tal manera que cada sumando tenga un denominador libre de cuadrados en vez
de una potencia de un factor libre de cuadrados.
La integral
Consideremos un sumando particular
Rij
en el que j > 11. Como Qi es libre de cuadrados, mcd( Qi , Qi0 ) =
j
Qi
1. Usando el algoritmo extendido de Euclides determinamos polinomio t, s Q[ x ] tal que t Qi + s

Qi0 = 1. Y en particular se pueden obtener polinomios S, T Q[ x ] tal que
S Qi + T Qi0 = Rij
j
donde grado( S ) <grado( Qi ) 1 y grado( T ) <grado( Qi ) . Dividiendo por Qi

S Qi
j
T Qi0
j
Qi
Qi
Rij
j
Qi
Ahora aplicamos integracin por partes a
T Qi0
j
Qi
Rij
j
Qi
Z
T Qi0
j
Qi
j 1
Qi
j 1
( j 1) Q i
T Qi0
j
Qi
con u = T y v =
1
j 1
( j 1) Q i
T0
j 1
( j 1) Q i
. Nos queda
32
Entonces tenemos una disminucin en el grado del denominador:
Rij
j
Qi
T/( j 1)
j 1
S + T 0 / ( j 1)
j 1
Qi
Qi
En resumen,
a) Con este proceso, si j > 1, cada integral

la integral original y la integral
Rij
j
Qi
contribuye con
S + T 0 / ( j 1)
termina.
b) En el caso de que j 1 = 1, la integral
S + T 0 / ( j 1)
integral original.
j 1
Qi
a la parte racional de
podra ser cero y todo el proceso de reduccin
j 1
Qi
T/( j 1)
j 1
Qi
contribuye a la parte logartmica de la
c) Si j 1 > 1 se puede aplicar nuevamente el proceso de reduccin a
S + T 0 / ( j 1)
hasta que
j 1
Qi
los denominadores restantes queden libres de cuadrados y esto termina con la parte racional de
la integral.
Ejemplo 1.22
Para fines ilustrativos vamos a aplicar el mtodo a la integral
Mathematica.
In[]:= P = 4 x^3 + 16 x^5 + 2 x^7;

Q = 2 x^2 + 8 x^4 + 8 x^6;
d = PolynomialGCD[P, Q];
(* Normalizar*)
P = PolynomialQuotient[P, d, x];
Q = PolynomialQuotient[Q, d, x];
PolynomialGCD[P, Q];
(* Dividir *)
(* = 2x2
2x7 + 16x5 + 4x3

usando
8x6 + 8x4 + 2x2
*)
(* = 2x + 8x3 + x5 *)
(* = 1 + 4x2 + 4x4 *)
(* = 1
*)
{L = PolynomialQuotient[P, Q, x], R = PolynomialRemainder[P, Q, x

]};
x
(* P=L+R/Q con L=
y R= 7x
4 +7
4
*)
(* Factorizacin libre de cuadrados de Q*)
FactorSquareFree[Q];
(* Q= (1 + 2x2 )2 *)
(* Expansin en fracciones parciales libre de cuadrados *)
7x
7x
ApartSquareFree[R/Q]; (*
*)
2(1 + 2x2 ) 4(1 + 2x2 )2
R11 = 7 x/2;
R12 = -7 x/4 ;
Q1 = 1 + 2 x^2;
(* Solo tenemos el caso j = 2 *)
(* Calcular S y T con algoritmo extendido de Euclides para Q1 y
Q1 *)
Q1p = D[Q1, x];
{d, {s, t}} = PolynomialExtendedGCD[Q1, Q1p, x];
7x2
{S = s*R12, T = t*R12};
(* T =
*)
8
(* Resultado: Parte racional *)
PRacional=Integrate[L,x] - T/Q1
(* Parte logartmica *)
Integrate[R11/Q1] + Integrate[(S+Tp)/Q1];
Out[]=
7x2
x2
8
8 (2x2 + 1)
Ejemplo 1.23
Aplicar el mtodo de Hermite a la integral
matica.
x6
6x5 + 6x4 8x3 18

usando Mathe 5x4 8x3 2x2 + 2x + 1
In[]:= (*Aplicar Hermite a *)

Clear[P, Q]
P = 6 x^5 + 6 x^4 - 8 x^3 - 18; x^2 + 8 x + 8;
Q = x^6 - 5 x^4 - 8 x^3 - 2 x^2 + 2 x + 1;
(* Dividir *)
L = PolynomialQuotient[P, Q, x];
(* = 1 *)
(* = 0 *)
33
34
R = PolynomialRemainder[P, Q,x];
= 18 8x3 + 6x4 + 6x5 *)
(* P = L + R/Q *)
(*Factorizacin "libre de cuadrados" de Q*)
(*Expansin en fracciones parciales*)
ApartSquareFree[R/Q]; (*=
(*
2(9 4x3 + 3x4 + 3x5 )

*)
1 + 2x 2x2 8x3 5x4 + x6
(* j=1. La integral solo tiene parte logartmica *)
1.13
Algoritmo de Rothstein/Trager para obtener la parte logartmica.
Despus de aplicar el mtodo de Hermite, podemos suponer que en la integral solo quedan fracciones
propias Rij /Qi con Qi mnico y libre de cuadrados. Las integrales de estas fracciones contribuyen a
la parte logartmica de la integral y como tal, deben tener la forma
Z
R
=
Q
ci log(vi )
(1.7)
i =1
donde los vi se pueden asumir libres de cuadrados y primos relativos y los ci son todos distintos.
Esto se puede suponer porque podemos usar las propiedades del logaritmo para separar o agrupar,
c1 log(v1 v2 ) + c2 log(v1 v3 ) = (c1 + c2 ) log v1 + c1 log v2 + c2 log v3 .
Rothstein ([27, 1977]) y Trager descubrieron de manera independiente que en la ecuacin (1.7), cada
vi = mcd( R ci Q0 , Q) y los ci son los nmeros tales que mcd( R ci Q0 , Q) 6= 0, esto es equivalente
a que el resultante res x ( R ci Q0 , Q) = 0. Este resultante es un polinomio en c con grado menor o
igual que el de Q y con races posiblemente repetidas, por tanto se puede factorizar de manera rpida
usando factorizacin libre de cuadrados.
Teorema 1.6
(Rothstein/Trager, 1977). Sean R, Q Q[ x ] con mcd( R, Q) = 1 tal que Q es mnico y libre de
cuadrados. Sea S el conjunto de ceros distintos del polinomio res x ( R cQ0 , Q). Entonces
Z
R
=
Q
c log( mcd ( R cQ0 , Q))
c S

35
y esta es la expresin que usa la mnima cantidad de extensiones algebraicas de Q.
Antes de ver un ejemplo, vamos a ver todos los detalles con los que se llega a este resultado.
Clculo del Resultante: El mtodo de Rothstein/Trager usa el Resultante para calcular la parte logart-
mica de la integral. El resultante apareci en el siglo IXX ligado a la solucin de problemas de sistemas
de ecuaciones no lineales.
Matriz de Sylvester. Sean A( x ), B( x ) F [ x ] con A( x ) = im=0 ai xi y B( x ) = in=0 bi xi . La matriz de
Sylvester de A y B es la matriz (m + n) (m + n)
am
M=
bn
a m 1
am
bn 1
bn
a m 1
bn 1
a1
am
b1
bn
a0
a1
b0
b1
a0
a0
b0
b0
donde las entradas que no se muestran son cero.
Ejemplo 1.24
Si A( x ) = a3 x3 + a2 x2 + a1 x + a0
y
B( x ) = b2 x2 + b1 x + b0 ,
entonces
M=
a3
0
b2
0
0
a2
a3
b1
b2
0
a1
a2
b0
b1
b2
a0
a1
0
b0
b1
0
a0
0
0
b0
Definicin 1.2
El determinante res ( A, B) =det M se llama el resultante de A y B. Si se quiere incluir la variable
x, se escribe res x ( A, B).
36
La primera propiedad importante es que el resultante se puede escribir en trminos de la races de A

y B.
Teorema 1.7
teores Sean A, B D [ x ] donde D es un dominio integral. Si A( x ) = am im=1 ( x i ) y B( x ) =
bm in=1 ( x i ), entonces
n
res ( A, B)
= (1)mn bnm A( i )
i =1
= a n B ( i )
i =1
= a n b n ( i i )
i =1
En particular, si A y B tienen un factor ( x ) (o una raz) en comn, entonces

res ( A, B) = 0 si y slo si grado( mcd ( A, B)) > 0.
Como Q es un campo, Q[ x ] es un dominio integral con divisin ecuclidiana. El algoritmo SubResultante se usa para calcular el resultante ([17, pg 24]).
Frmula de Rothstein/Trager. Primero necesitamos un lema acerca de la forma de R y Q en trminos de
los componentes de la primitiva.
Lema 1.0
n
R
= ci log(vi ) con mcd( R, Q) = 1, Q mnico y libres de cuadrados y adems
Q
i =1
los vi primos relativos, mnicos y libres de cuadrados. Entonces
!
lema1 Sea
Q = v1 v2 v n
R=
i =1
ci vi0 v j
j 6 =i
Prueba. En vez de hacer una prueba formal sobre todos los ndices, vamos a ver la idea de la prueba con
el caso n = 2. La generelizacin despus es inmediata.
37
Sea
R
= c1 log v1 + c2 log v2
Q
(1.8)
donde los v1 , v2 son libres de cuadrados y primos relativos y c1 , c2 son nmeros distintos.
Entonces, derivando a ambos lados de la ecuacin (1.8) se obtiene
R
Q
= c1
=
v10
v0
+ c2 2
v1
v2
c1 v10 v2 + c2 v20 v1
.
v1 v2
de donde,
Rv1 v2 = Q(c1 v10 v2 + c2 v20 v1 )
(1.9)
Recordemos que mcd(v1 , v2 ) = 1 y cada vi no tiene factores repetidos,es decir, mcd(v1 , v10 ) = 1.
Ahora, tanto v1 como v2 dividen a Q(c1 v10 v2 + c2 v20 v1 ), pero como no dividen la suma c1 v10 v2 + c2 v20 v1 ,
deben dividir a Q, por hiptesis entonces v1 v2 | Q. Ahora, como Q| Rv1 v2 , por hiptesis entonces
Q|v1 v2 . Como Q y los vi son mnicos se concluye Q = v1 v2 . Finalmente, cancelando Q en la ecuacin
(1.9) se obtiene la expresin para R. Es decir,
Q = v1 v2
R = c1 v10 v2 + c2 v20 v1
La prueba en el caso general sigue este mismo razonamiento.
Ahora vamos a demostrar, igual con n = 2, que vi = mcd ( R ci Q0 , Q).

Primero hay que recordar que mcd ( A B, B) = mcd ( A, B). Este hecho se usa repetidamente para
calcular v1 y v2 . Usando el lema ?? tenemos
38
mcd ( R ct Q0 , vi )
mcd (c1 v10 v2 + c2 v20 v1 ct (v10 v2 + v1 v20 ), vi )
mcd ((c1 ct )v10 v2 + (c2 ct )v20 v1 , vi )
0
0
mcd ((c1 ct )v1 v2 + (c2 ct )v2 v1 , v1 )
si
i=1
mcd ((c2 ct )v20 v1 + (c1 ct )v10 v2 , v2 )
si
i=2
0
mcd ((c1 ct )v1 v2 , v1 )
si
i=1
mcd ((c2 ct )v20 v1 , v2 )
si
i=2
mcd (0, v1 ) = v1
si
si
mcd (0, v2 ) = v2
t=1
o
t=2
si
t 6= 1
si
t 6= 2
Finalmente, si se conocen los c1 y c2 , entonces

mcd ( R ci Q0 , Q) = mcd ( R ci Q0 , v1 v2 )
mcd ( R ci Q0 , v1 ) mcd ( R ci Q0 , v2 )
v1 1 = v1
si
i=1
1 v2 = v2
si
i=2
Para el clculo de los ci solo hay que observar que estos nmeros c son precisamente los nmeros
tales que
mcd ( R cQ0 , Q) 6= 1
pero, por el teorema ??, estos las races del polinomio res x ( R cQ0 , Q) = 0.
De nuevo, la demostracin del caso general sigue el mismo razonamiento.
Ejemplo 1.25
En el ejemplo 1.22 calculamos la parte racional de la integral
Z
2x7 + 16x5 + 4x3

8x6 + 8x4 + 2x2
con el mtodo de Hermite. Ahora vamos a completar el clculo con la parte logartmica.
Tenamos en el ejemplo 1.22
2x7 + 16x5 + 4x3

x2
7x2
7x/2
+
=
+
6
4
2
8
8 (2x2 + 1)
1 + 2x2
8x + 8x + 2x
|
{z
Z
0
1 + 2x2
}
Parte logartmica
Solo debemos calcular
7x/2
con el algoritmo de Rothsteinh/Trager.
1 + 2x2
In[]:= R=7x/2;
Q1=1+2x^2
Q1p = D[Q1, x];
Resultant[R11 - c Q1p , Q1, x];
(* = 1/4(49 112c + 64c2 )*)
{c1, c2} = Flatten[c /. {Solve[Resultant[R11 - c Q1p , Q1, x] ==
0, c]}];
(* = {7/8, 7/8} *)
(* Races repetidas, c1=c2. Solo calculamos v1 *)
v1 = PolynomialGCD[R11 - c1 Q1p, Q1, x]; (* v1 = 1 + 2x2 *)
(* Parte Logartmica *)
PLogaritmica = c1 Log[v1];
(*
2
= 7/8 log(1 + 2x ) *)
(* Finalmente I = Parte Racional + Parte Logartmica*)
PRacional + PLogaritmica
Out[]=
x2
7x2
+ 7/8 log(1 + 2x2 )

8
8 (2x2 + 1)
Ejemplo 1.26
Ya habamos visto, en el ejemplo 1.23, que la integral
Z
7x1 3 + 10x8 + 4x7 7x6 4x3 4x2 + 3x + 3

dx
x1 4 2x8 2x7 2x4 4x3 x2 + 2x + 1
solo tiene la parte logartmica. Bien, la vamos a calcular ahora.
In[]:= Clear[P, Q1, Q1p, Q, a1, a2];

P = 6 x^5 + 6 x^4 - 8 x^3 - 18 x^2 + 8 x + 8;
Q = x^6 - 5 x^4 - 8 x^3 - 2 x^2 + 2 x + 1;
(* = Q *)
Q1 = Q;
Q1p = D[Q1, x];
39
40
R = Factor[Resultant[P - c Q1p , Q1, x]];(* = 1453248(2 2c + c2 )3 *)

R = -2 - 2 c + c^2;
{c1, c2} = Flatten[c /.{Solve[R == 0, c]}];(* c1 = 1 3, c2 = 1 + 3*)
(* Se debe usar la extensin 3 para obtner el MCD en Q( 3)[ x ] *)

v1 = PolynomialGCD[P - c1 Q1p, Q1, Extension -> {Sqrt[3]}];
(*v1 = 1 x + 3x + 3x2 + x3 *)
v2 = PolynomialGCD[P - c2 Q1p, Q1 , Extension -> {Sqrt[3]}];
(* v2 = 1 + x + 3x + 3x2 x3 *)
(* Parte logartmica *)
c1*Log[v1] + c2*Log[v2]
Out[]=
1+
1.14
h
i
h
i

3 Log 1 + x + 3x + 3x2 x3 + 1 3 Log 1 x + 3x + 3x2 + x3
lgebra diferencial y el algoritmo de Risch
Un hecho notable en integracin simblica es que toda la teora se puede poner en trminos algebraicos, sin usar procesos de lmite. A mediados del siglo XX, J. Ritt invent el lgebra Diferencial
que es un marco apropiado para tratar el problema de integrar una funcin en trminos finitos. En
1968, M. Rosenlicht public la primera versin puramente algebraica del teorema de Liouvile usando
lgebra Diferencial y R. Risch, en 1970, especific un algoritmo para el teorema. El algoritmo reduce
el problema de integracin a un problema de decisin, pero no ha sido implementado totalmente en
ningn sistema por lo complejo que se puede volver en algunos casos.
Definicin 1.3
Un campo diferencial es un campo F de caracterstica 0 en el cual se ha definido una aplicacin
D : F F tal que para todo f , g F se cumple
D ( f + g) = D ( f ) + D ( g)
(1.10)
D ( f g) = f D ( g) + g D ( f )
(1.11)
La aplicacin D se llama operador diferencial
41
Teorema 1.8
Sea D un operador diferencial en un campo diferencial F , entonces se cumple
a) D (0) = D (1) = 0
b) D ( f ) = D ( f ) para toda f F

g D ( f ) f D ( g)
f
=
para toda f , g F con g 6= 0
c) D
g
g2
d) D ( f n ) = n f n1 D ( f ) para toda n Z, f F con f 6= 0
Prueba. La demostracin usa la definicin 1.3 y es enteramente algebraica.

a) D (0) = D (0 + 0) = D (0) + D (0) por (1.10). Por tanto D (0) = 0.
D (1) = D (1 1) = 1 D (1) + 1 D (1) por (1.11). Por tanto D (1) = 0.
b) D (0) = D ( f + ( f )) = D ( f ) + D ( f ) = 0, entonces D ( f ) = D ( f ) .

1
D ( g)
= 2 si g 6= 0.
c) Como 0 = D (1) = D ( g 1/g) = g D (1/g) + 1/g D ( g), entonces D
g
g

f
f D ( g) D ( f )
g D ( f ) f D ( g)
Ahora, D
= f D (1/g) + 1/g D ( f ) =
+
=
g
g
g2
g2
d) Se prueba por induccin sobre n.
Interludio: Extensiones de campo

Para lo que sigue necesitamos recordar un poco de extensiones de campo. Sea F un
campo. Si F es un subcampo de E entonces E es una extensin de campo de F.
E se puede ver como un espacio vectorial sobre F y como tal, puede tener dimensin
finita o infinta. Si la dimensin es finita, E se dice extensin finita de F.
Por ejemplo, el campo Q( 2) = { a + b 2, : a, b Q} es una extensin de Q.

Sea E, es algebraico sobre F si es raz de un polinomio P 6= 0, P F[ x ]. Si no
es algebraico, es trascendental.
42

f ()
Sea tascendente sobre F. Entonces F() =
: f , g F[ x ], g() 6= 0 . F() es el
g()
ms pequeo subcampo de E que contiene a F y a . A F() se le llama una extensin
simple de F .
Por ejemplo, sea = log x. Entonces Q( x, ) es una extensin del conjunto de funciones
racionales Q( x ), y
x + x2 log2 x
x + x2 2
=
Q( x, )
x
x
x + 2 log2 x
x + 2 2
b)
=
Q( 2)( x, )
x + log x
x+
a)
Si 1 , 2 , ..., k son elementos de E, entonces

F F(1 ) F(1 , 2 ) = F(1 )(2 ) ... F(1 , ..., k ) = F(1 , ..., k1 )(k )
Sea F un subcampo de E y sea E algebraico sobre F. El (nico) polinomio mnico

e irreducible P F[ x ] tal que P() = 0, se llama polinomio mnimo de sobre F. Si
grado( P( x )) = n 1 entonces
F () = F[] = {c0 + c1 + + cn1 n1 : ci F.}
Extensin diferencial de un campo diferencial. Si F, E son campos diferenciales con operadores difer-
enciales DF y DG respectivamente, entonces decimos que E es una extensin diferencial de F si E

es una extensin (de campo) de F y si
DF ( f ) = DG ( f ) para todo f F.
Homomorfismos diferenciales. Si F, E son campos diferenciales con operadores diferenciales DF y DG
respectivamente, la aplicacin : F E es un homomorfismos diferencial si es un homomorfismo

de campo y
( DF ( f )) = DG (( f )) para todo f F
Si E es una extensin de F, existe un homomorfismo inyectivo , de E en F tal que ( DF ( f )) =
DG (( f )) para todo f F. Esto dice que existe un subconjunto en E, isomorfo a F y las derivadas
coinciden en este subconjunto.
Campo de constantes. Sea F un campo diferencial con operador diferencial D. El campo de constantes
de F es el subcampo K definido como
K = { c F : D ( c ) = 0}
43
Extensiones. Una funcin racional f Q( x ) es integrable en Q( x ) si existe g Q( x ) tal que D ( g) = f .
Ya sabemos que 1/x Q( x ) pero no existe g Q( x ) tal que D ( g) = 1/x. Por lo tanto, para expresar
la antiderivada de una funcin racional podra ser necesario extender el campo Q( x ) con nuevas funciones. Como ya sospechamos, solo se necesitan logaritmos.
Ejemplo 1.27
Consideremos el campo de funciones racionales Q( x ). Si D ( f ) = f 0 (es decir, la derivada
ordinaria) entonces Q( x ) es un campo diferencial. En este caso, el campo de constantes es
K = Q.
E = Q( x, log x ) es una extensin de F = Q( x ). La derivada DE coincide con la derivada DF en
Q( x ).
En un contexto general, los logaritmos se definen algebraicamente usando el operador D con una
manera que captura la esencia de los logaritmos:
Definicin 1.4
Sea F un campo diferencial y sea E una extensin diferencial de F. Un logaritmo sobre F es un
elemento E con la propiedad de que existe u F tal que
D ( ) =
D (u)
y se escribe = log(u)
u
El problema de integracin en trminos finitos. Dado un campo diferencial F, el proceso de integracin
de f F es determinar si existe o no, una extensin de campo E en la cual existe g E tal que
D ( g) = f .
Consideremos el campo diferencial F = Q( x ). La funcin
de extensiones
1
es integrable pues existe una torre
x2 2
F = Q( x ) Q( 2)( x ) Q( 2)( x, log( x 2) Q( 2)( x, log( x 2, log( x + 2)) = E
44
en la cual existe g E, g =
1
1
log( x 2) log( x + 2), tal que D ( g) = f .
4
4
Como se sabe, para cada f Q( x ) siempre es posible encontrar una extensin diferencial
E = Q(1 , ..., k )( x, 1 , 2 , ..., n )
en la cual exista un elemento g E tal que D ( g) = f . En este caso, como se sabe, 1 , ..., k son nmeros
algebraicos y cada i es un logaritmo.
Ejemplo 1.28
Z
1
1
dx =
Q( x )
x+1
( x + 1)2
1
dx = log( x ) Q( x, log( x ))
x
a)
b)
Z
c)
x2
Z
d)
1
1
1
dx = log( x 2) log( x + 2) Q( 2)( x, log( x 2), log( x + 2))
4
4
2
1 x3
dx
x3 + x
= x + log( x ) 12 log( x i )( x + i )
x +i
2i log
Q(i )( x, log( x ), log( x i ), log( x + i ))
x i
En esta seccin vamos a trabajar en K( x ) = Q( x ) o en K( x ) = Q(1 , ..., n )(y1 , y2 , ..., yn ) donde los
i son nmeros algebraicos
y las yi son otras variables. Adems, por convencin, D ( g) se expresar
Z
como g0 y D 1 ( f ) =
f y adems no se incluir la constante de integracin.
Una herramienta que necesitamos es integracin por partes,
Teorema 1.9
Si F es un campo diferencial con derivada D, entonces para cualquier u, v F,
Z
u D (v) = u v
La prueba consiste en aplicar D en lado derecho.
v D (u)
1.15
45
Integracin de funciones racionales: La parte racional.
Podemos aplicar el mtodo de Hermite sobre el campo diferencial K( x ) si K es de caracterstica 0 y

si D ( x ) = 1.
Recordemos que, siguiendo la presentacin del mtodo en la seccin 1.12, este mtodo usa operaciones
de polinomios para hace la reduccin
Z
P
C
=
+
Q
D
i =1 j =1
Rij
j
Qi
donde Q = ik=1 Qii es la factorizacin libre de cuadrados de Q. Si j > 1 entonces existe S, T K[ x ]

tal que
S Qi + T Qi0 = Rij
y, usando integracin por partes,
Z
a)
T/( j 1)
j 1
Qi
b) Si j 1 = 1,
Z
Rij
T/( j 1)
S + T 0 / ( j 1)
+
=
j
1
j 1
Qi
Qi
Qi
es la contribucin de cada
Z
S + T 0 / ( j 1)
j 1
Qi
Rij
a la parte racional de la integral.
Qi
contribuye a la parte logartmica.
c) Si j 1 > 1, se repite el proceso de reduccin.

Al final del proceso, obtenemos la parte racional de la integral sumando todas las contribuciones
y adems obtenemos, sin calcular, la parte logartmica de la integral.
1.16
Integracin de funciones racionales: La parte logartmica.
Despus de aplicar el mtodo de Hermite, podemos suponer que en la integral solo quedan fracciones
Rij /Qi con Qi mnico y libre de cuadrados.
El mtodo de Rothstein/Trager sigue las mismas lneas presentadas en la seccin 1.13. Se enuncia en
el teorema que sigue.

46
Teorema 1.10
Sea K ( x ) un campo diferencial con campo de constantes K . Sean P, Q K [ x ] con
mcd( P, Q) = 1, Q mnico, libre de cuadrados y grado( P ) <grado( Q. ) Supongamos que
Z
P
=
Q
ci log(vi )
i =1
donde las ci K son n constantes distintas no nulas y los vi K [ x ] son mnicos, libres
de cuadrados y primos relativos dos a dos y de grado positivo. Entonces los ci son las races
distintas del polinomio
R(c) = res x ( P( x ) c Q0 ( x ), Q( x )) K [ x ]
y, una vez calculadas las races ci , los vi son los polinomios
vi = mcd( P( x ) ci Q0 ( x ), Q( x ))
Si K 6= Q, hay que tener el cuidado de calcular cada mcd( P( x ) ci Q0 ( x ), Q( x )) en K [ x ] y no en

Q( x ).
Un mtodo ms eficiente es el algoritmo mejorado de Lazard/Rioboo/Trager, pero para el propsito
de este artculo, el mtodo que hemos visto es suficiente.
1.17
Algoritmo de integracin de Risch
Como ya vimos, integrar una funcin racional P( x )/Q( x ) requiere solo algoritmos para manipulacin
de polinomios: Divisin de polinomios con resto, clculo del mcd, factorizacin de polinomios y solucin de ecuaciones.
En el caso general de una funcin elemental f , que podra involucrar funciones tales como log, exp, sen, cos,
etc., las herramientas requeridas son algoritmos para manipular polinomios en varias variables. En este
caso, las expresiones no racionales que aparecen en f se tratan como smbolos independientes de tal
manera que f se manipula como una funcin racional en varias variables: Los pasos para integrar f
son entonces muy parecidos a los pasos que usamos para integrar funciones racionales.
En nuestra nueva notacin vamos a definir varias cosas familiares.

p
n
f,
47
Definicin 1.5
Sea F un campo diferencial y E una extensin diferencial de F.
a) E es un logaritmo sobre F si existe u F tal que 0 =

b) E es una exponencial sobre F si existe u F tal que
u0
. Se escribe = log(u).
u
0
= u0 . Se escribe = exp(u).
c) E es algebraico sobre F si existe P F[z] tal queda P( ) = 0

d) Si E no es algebraico sobre F, entonces es trascendente sobre F.
e) La extensin diferencial E se llama extensin trascendental elemental de F si E =
F(1 , ..., n ) donde cada i es trascendental y, o es logartmica o es exponencial sobre el
campo Fi1 = F(1 , ..., i1 ).
f) La extensin diferencial E se llama extensin elemental de F si E = F(1 , ..., n ) donde
cada i es logartmica, exponencial o algebraica sobre el campo Fi1 = F (1 , ..., i1 ).
Definicin 1.6
Sea K( x ) un campo diferencial de funciones racionales sobre un campo de constantes K C.
Si F es una extensin trascendental elemental de K( x ) entonces decimos que F es un campo de
funciones elementales trascendentales. Similarmente, si F es una extensin elemental de K( x )
entonces F se llama un campo de funciones elementales.
Ejemplo 1.29
Sea 1 = log( x ). 1 es logartmica sobre Q( x ).
Sea 2 = exp(1 /2). Como
Q( x, 1 ).
2 21 10
20
1
=
= 10 Q( x, 1 ), entonces 2 es exponencial sobre
2
2
2
2 es algebraica sobre Q( x ) pues satisface la ecuacin 22 x = 0
Teoremas de estructura. es un monomio sobre el campo diferencial F si F( ) y F tienen el mismo

campo de constantes y adems es trascendental sobre F y es, o exponencial o logartmica sobre F.
48
Si F( x, ) es una extensin trascendental de F( x ), las operaciones racionales sobre F( x, ) son equivalentes a ejecutar operaciones aritmticas en el campo de funciones racionales F( ) . El problema es
saber si una extensin es trascendental (esto puede ser difcil, como en el caso de campos de constantes).
Los teoremas de estructura (ver [5], [17]) establecen requerimientos explcitos para que una nueva extensin sea independiente de las extensiones previas. Estos teoremas establecen los clculos que se
deben hacer para decidir si existe alguna relacin algebraica entre los objetos presentes en el campo y
el nuevo objeto que se quiere adjuntar. Si no hay relacin, el nuevo miembro es un monomio.
Por ejemplo, si 1 = exp( x ) y 2 = exp(2x ), entonces 2 = 2 , por lo que Q( x, 1 , 2 ) = Q( x, 1 ).
En lo que sigue, vamos a suponer que se han aplicado los teoremas de estructura y que la representacin de f como f K(1 , ..., n )(1 , ..., m ) solo involucra extensiones trascendentes.
Derivadas de funciones elementales. Ahora ya podemos ser ms precisos en varias cosas. Un campo
de funciones elementales es cualquier extensin del campo de funciones racionales K( x ), finitamente

generada, de tal manera que cada extensin sea o logartmica, exponencial o algebraica.
Desde el punto de vista computacional, es conveniente, como decamos antes, que K se de la forma
Q(1 , 2 , ..., n ) con cada i algebraico sobre Q. El nmero de extensiones
i es el requerido para el
2
problema actual. Por ejemplo, en el caso de 1/( x 2) solo se requiere 1 = 2.
El problema de integracin
Z en trminos finitos es: Si f Q(1 , ..., n )(1 , ..., m ), determinar extensiones
adicionales tal que g =
f K(1 , ..., n+h )(1 , ..., m+k ) y determinar explcitamente g, o probar que
no existe tal funcin elemental.
El teorema de Liouville prueba la observacin de Laplace (seccin 1.5) acerca de que la primitiva de
una funcin no puede incluir ninguna otra exponencial ni ningn otro radical de los que ya figuran en
la funcin. Para probar esto necesitamos establecer primero la accin del operador diferencial D sobre
polinomios en , donde es trascendental o logartmica, trascendental o exponencial u algebraica
sobre un campo difeencial F.
Teorema 1.11 (Derivadas)

Sea F un campo diferencial y sea F( ) una extensin diferencial de F con el mismo subcampo
de constantes.
a) Si es trascendental o logartmica sobre F, entonces para cada P( ) F[ ] de grado
positivo, P0 ( ) F[ ] y el grado de esta derivada disminuye en 1 si su coeficiente
principal es una constante, sino, se mantiene igual.
49
b) Si es trascendental o exponencial sobre F, entonces para cada P( ) F[ ] de grado

positivo, P0 ( ) F[ ] y el grado de esta derivada se mantiene igual. Adems: P( ) divide
P0 ( ) si y slo si P( ) = H n para alguna H F y n entero.
c) Si es algebraico sobre F con polinomio mnimo P(z) =
N +1
pi zi F[z] mnico, entonces
i =0
0 =
d( )
F( )
e( )
con d, e F[z] y
N
d(z) =
pi0 zi ,
e(z) =
i =0
( i + 1 ) p i +1 z i
i =0
Ahora podemos ver el teorema de Liouville. En este teorema se basan los algoritmos de integracin en
trminos finitos.
Teorema 1.12 (Teorema de Liouville).

Sea F un campo diferencial con campo de constantes K. Supongamos que para f F, existe
una
extensin elemental G de F, con el mismo campo de constantes K, y g G tal que
Z
f = g. Entonces existen v0 , v1 , ..., vn F y constantes c0 , c1 , ..., cn K tal que
f=
v00
v0
+ ci i ,
vi
i =1
m
es decir,
f = v0 + ci log(vi )
i =1
El teorema es fcil de interpretar, a al luz del teorema 1.11: Si la primitiva de f F requiere una
extensin logartmica trascendental , entonces f = (c + d)0 F[ ]. La primitiva de f F no requiere
una exponencial diferente de las que ya hay en f ni requiere nuevas extensiones algebraicas.
Las primitivas pueden variar en apariencia, pero una de las formas debe encajar con el teorema de
Liouville. Por ejemplo, la siguiente primitiva en principio parece que requiere extensiones nuevas
diferentes a las que aparecen en la integral,
Z
2x3 2x2 1
exp( x2 )
( x 1)2
exp( x2 + log( x )/2)

exp( x2 + log( x )/2)
+
2( x 1)
2( x + 1)
50
Las extensiones exp( x2 + log( x )/2) y x no estn en f . Pero, simplificando exp( x2 + log( x )/2) =
exp( x2 ) x y sumando, se obtiene la forma esperada,

Z
2x3 2x2 1
exp( x2 )
( x 1)2
x
exp( x2 )
x1
Prueba del teorema de Liouville. La prueba es por induccin y se divide en tres casos, dependiendo de
si se requieren extensiones trascendentales logartmicas, exponenciales o extensiones algebraicas (ver

[5], [17], [24]).
La hiptesis de la prueba es que existe g0 = f con g G y que existen 1 , 2 , ..., N , donde cada i es
logartmica, exponencial o algebraica sobre Fi1 , tal que G = F (1 , 2 , ..., N ). Adems cada Fi1 tiene
el mismo campo de constante K.
La prueba es sobre el nmero de extensiones N y consiste en demostrar que f tiene la forma establecida por el teorema, ya sea para las extensiones logartmicas o para establecer que no se necesitan
extensiones exponenciales adicionales ni algebraicas de las ya presentes en la integral.
El teorema se cumple para N = 0 : En efecto, en este caso no hay necesidad de extensiones y g F
con lo cual g0 = f y m = 0 y v0 = g.
La hiptesis de induccin es que el teorema se cumple para cualquier nmero de extensiones < N.
Para el caso de N extensiones, se puede ver el campo F(1 , 2 , ..., N ) como F(1 )(2 , ..., N ). De ahora
en adelante ponemos en vez de 1 .
Podemos aplicar el teorema para
f F( )
g F( )(2 , ..., N ).
por hiptesis de induccin, existen v1 ( ), ..., vm ( ) F( ) y constantes c1 , c2 , ...cn K tal que

m
f = v00 ( ) + ci
i =1
vi0 ( )
vi ( )
(1.12)
Como la prueba completa puede ser extenuante, aqu solo vamos a hacer un caso de la prueba para el
caso de extensiones trascendentales y logartmicas.
Prueba del teorema para el caso de extensiones trascendentales y logartmicas . Sea es extensin
trascendental y logartmica sobre F. Si fuera necesario, podemos aplicar la regla log(vi v j ) = log vi +
log v j , y asumir que, para i = 1, ..., m, cada vi ( ) F o que vi ( ) es mnico e irreducible en F[ ],
con grado positivo y adems todos estos vi son diferentes y las ci no son nulas. Bien, como es
trascendental, v0 ( ) F( ). Sea
v0 ( ) =
A( )
con A, B F[ ], mcd ( A, B) = 1 y B mnico.
B( )
51
Ahora factorizamos B( ) = id=1 Bi ( )ri con los Bi ( ) distintos, mnicos e irreducibles en F( ) y los
ri enteros positivos. Desarrollamos v0 ( ) en fracciones parciales,
d
ri
v0 ( ) = A0 ( ) +
Aij
j
i =1 j=1 Bi ( )
con A0 , Aij , Bi F[ ] y grado Aij < grado Bi
Entonces, derivando respecto a x, la ecuacin 1.12 se convierte en

f = A0 ( ) 0 +
Aij ( )0
ri
i =1 j =1
Bi ( ) j
j Aij ( ) Bi ( )0
Bi ( ) j+1
v ( )0
ci vii ( )
(1.13)
i =1
Recordemos que f no depende de (no est presente en f ). Como es logartmica sobre F, entonces,
segn el teorema 1.11, a), para cualquier polinomio P( ), mnico e irreducible en F[ ] de grado positivo, se tiene P0 ( ) F[ ] y el grado de P0 es inferior al grado de P y, por tanto, P( ) no divide a
P0 ( ). Bien, si P es uno de los Bi en 1.13 entonces el lado derecho de esta ecuacin contiene un trmino
con denominador P( )ri +1. Pero este trmino no puede cancelar con nadie, por lo que debe aparecer en
f pero esto contradice el hecho de que no est en f . Por tanto la doble suma en 1.13 no debe aparecer.
De la misma manera P( ) no puede ser uno de los vi ( ) porque no puede cancelar con la derivada
vi ( )0 y hara depender f de . As, los vi no dependen de .
Entonces f tiene la forma,
f = A0 ( ) 0 +
ci
i =1
vi0
vi
(1.14)
y debe ser que A0 ( )0 sea independiente de . Por el teorema 1.11, a), solo podra pasar que A0 ( ) =
c + d F[ ] con lo que f tiene la forma
f = d0 + c
m
v0
u0
+ ci i con d, u, vi F y c, ci K,
u
vi
i =1
como se deseaba.
Los otros casos usan el teorema 1.11 tal y como se puede ver en [17], [5] o [24], por ejemplo.
1.18
Algoritmo de Risch: Caso de extensiones logartmicas
Ahora ya podemos pasar a describir el algoritmo de Risch. Como el tema es tan extenso, solo lo vamos
a hacer en el caso de extensiones logartmicas.
El algoritmo de Risch es un proceso de decisin para integrar funciones elementales que estn a un
campo de funciones trascendentales elementales, es decir,
f K( x, 1 , 2 , ..., n )
52
con K un subcampo de C y cada i es trascendente y, i o es logartmica o es exponencial sobre Ki1 .

El algoritmo determina si
f existe como una funcin elemental, en otro caso, en el camino indica
por qu no existe esta primitiva elemental.

Como indicamos antes (teoremas de estructura), suponemos que K( x, 1 , 2 , ..., n ) es puramente trascendental. Como cada extensin es un smbolo trascendente, el integrando f puede ser manipulado
como una funcin racional de esos smbolos. Por tanto el algoritmo para funciones trascendentes
parece una reminiscencia del algoritmo para integrar funciones racionales. De hecho usamos el mtodo
de Hermite para la parte racional y el algoritmo de Rothstein/Trager para la parte logartmica.
Dada f K( x, 1 , 2 , ..., n ), f puede ser vista como una funcin racional en la extensin = n , es
decir,
f ( ) =
P( )
Fn1 ( ) con Fn1 = K( x, 1 , 2 , ..., n1 ).
Q( )
Adems suponemos que f est normalizada de tal manera que P, Q Fn1 ( ) son primos relativos y
Q( ) mnico. Como antes, la derivacin 0 es respecto a x.
El algoritmo es recursivo. El caso inicial es integracin de funciones racionales, en F0 = K( x ), el cual
ya est resuelto.
Algoritmo de Risch para extensiones logartmicas. Sea logartmica con 0 = u0 /u con u Fn1 .
Aplicando divisin euclidiana a P( ), Q( ) Fn1 [ ] obtenemos polinomios C ( ), R( ) Fn1 [ ] tal

que
P( ) = Q( ) C ( ) + R( ) con R( ) = 0 o grado R( ) < grado Q( ).
Entonces
f =
C ( )
| {z }
parte polinomial
R( )
.
Q( )
| {z }
Z
parte racional
Extensiones logartmicas: Integracin de la parte racional. Procedemos con el mtodo de Hermite: Cal-
culamos la factorizacin libre de cuadrados de Q( ) Fn1 [ ],

k
Q( ) =
Qi ( )i
i =1
con cada Qi ( ) mnico y libre de cuadrados, mcd( Qi ( ), Q j ( )) = 1 si i 6= j y grado Qk > 0.

d
Q ( )) = 1. Si
d i
0
0
x = 1 para x F, entonces esto implica que mcd ( Qi ( ), Qi ( ) ) = 1 ([5, Teorema 12.6]).
R( )
Usamos operaciones de polinomios para calcular la expansin en fracciones parciales de
en la
Q( )
forma
Qi ( ) es libre de cuadrados en el dominio de polinomios Fn1 [ ] si mcd ( Qi ( ),
R( )
=
Q( )
i =1 j =1
53
Rij ( )
j
Qi ( )
donde Rij ( ) Fn1 [ ] y grado Rij ( ) < grado Qi ( ) si grado Qi > 0 y Rij ( ) = 0 si Qi ( ) = 1.
Entonces tenemos,
Z
R( )
=
Q( )
i Z
i =1 j =1
Rij ( )
j
Qi ( )
Ahora procedemos con las integrandos Rij /Qi no nulos.

Caso j > 1 . Como mcd( Qi ( ), Qi ( )0 ) = 1, existe s, t Fn1 [ ] tal que
S( ) Qi ( ) + T ( ) Qi0 ( ) = Rij ( )
(1.15)
con grado S < grado Qi0 y gradoT < gradoQi .

Ahora dividiendo por Qi ( ) j en la ecuacin 1.15 y aplicando integracin por partes, se obtiene
Z
Z
Rij ( )
T ( ) / ( j 1)
S ( ) + T ( ) 0 / ( j 1)
=
+
j
1
j 1
Qi ( )
Qi ( )
Qi ( )
Si j 1 > 1, se repite el proceso de reduccin. Al final del proceso todos los denominadores quedan
libres de cuadrados y, podemos escribir la parte racional de la integral como,
Z
C ( )
R( )
=
+
Q( )
D ( )
A( )
con A, B, C, D Fn1 [ ], grado A < grado B y B mnico y libre de cuadrados
B( )
A( )
con el mtodo de Rothstein/Trager con la variante adecuada para el caso
B( )
de extensiones logartmicas.
Ahora calculamos
Teorema 1.13
Sea F un campo de funciones elementales con campo de constantes K. Sea trascendental y
logartmica en F y supongamos que la extensin trascendental elemental F( ) tiene el mismo
campo de constantes K. Sea A( )/B( ) F( ) con A( ), B( ) F[ ], mcd ( A( ), B( )) =
1, grado A < grado B y B mnico y libre de cuadrados .
A( )
es elemental si y slo si todas las races del polinomio R(z) = res ( A( ) z
B( )
0
B( ) , B( )) F[z], son constantes
Z
a)
54
m
A( )
A( )
v ( )0
es elemental, entonces
= ci i
donde las ci son las races distintas
B( )
B( )
vi ( )
i =1
de R(z) y donde las vi ( ) se definen como
b) Si
vi ( ) = mcd ( A( ) ci B( )0 , B( )) F(c1 , c2 , ..., cm )[ ]
Extensin logartmica: Integracin de la parte polinomial. A diferencia de la parte polinomial en la integral de una funcin racional simple, el caso polinomial en el caso de extensiones logartmicas es ms
laborioso.
Vamos a dar la descripcin del algoritmo siguiendo a [5, pgs 540-543].

Supongamos que la parte polinomial es
C ( ) = p s s + p s 1 s 1 + + p 0 Fn 1 [ ] y p i Fn 1
Por el principio de Liouville y otras consideraciones, C ( ) debe tener la forma
p s s + p s 1 s 1 + + p 0 = v 0 ( ) 0 +
vi0
v
i =1 i
m
con v0 ( ) = qs+1 s+1 + + q0 con qs+1 K y los otros qi Fn1 . Es decir,
p s s + p s 1 s 1 + + p 0 = q s +1 s +1 + + q 0 +
Z m 0
vi
vi
i =1
Ahora, derivando obtenemos el sistema de ecuaciones
0 = q0s+1
ps = (s + 1) qs+1 0 + q0s
..
p1 = 2q2 0 + q10
p = q 0 + Q0 con Q = q +
0
0
1
ci log(vi )
0
0
i =1
55
En este sistema, los pi Fn1 son los conocidos y las incognitas son qs+1 K, qi Fn1 .
Algoritmo. Aplicando integracin en el sistema, obtenemos primero que qs+1 es una constante de inte-
gracin, sea qs+1 = bs+1 (todava desconocida!)

Paso s + 1) Integrando en la segunda ecuacin, obtenemos
ps = (s + 1)bs+1 + qs
Ahora el proceso de integracin se invoca de manera recursiva para integrar
ps . Para resolver
esta
ecuacin para bs+1 K y qs Fn1 , se requiere que las siguientes condiciones se cumplan
Z
ps :
a*)
ps es elemental
b*) A los sumo hay una extensin logartmica de Fn1 apareciendo en la integral y si hay una
debe ser = log(u).
Si
Z una de estas dos condiciones falla, entonces la ecuacin no tiene solucin y se concluye que
ps no es elemental.
Si las condiciones a) y b) se cumplen, entonces si

bs + 1 =
cs
s+1
ps = cs + ds con cs K y ds Fn1 , se obtiene

q s = d s + bs
donde bs K es una constante de integracin no conocida todava.

Paso s) Sustituyendo qs en la tercera ecuacin, reagrupando e integrando se obtiene,
p s 1 s d s 0 = s bs + q s 1
Si las condiciones a*) y b*) se cumplen (sino,
ps no es elemental y termina el algoritmo),
entonces si
Z
se obtiene
ps1 s ds 0 = = cs1 + ds1 con cs1 K y ds1 Fn1 ,
56
bs =
c s 1
s
q s 1 = d s 1 + bs 1
donde bs1 K es una constante de integracin no conocida todava.
..
.
Paso 2) El proceso descrito se continua hasta llegar hasta la penltima ecuacin donde la solucin debe
ser de la forma
b2 =
c1
2
q1 = d1 + b1
donde b1 K es una constante de integracin no conocida todava.
Paso 1) Sustituyendo q1 en la ltima ecuacin, reagrupando e integrando se obtiene

Z
p0 d1
u0
= b1 + Q0
u
Si la integral de la izquierda es elemental (sino,

Z
p0 d1
ps no es elemental y termina el algoritmo) y
u0
= d0
u
entonces b1 es el coeficiente en d0 de = log(u) y

Q0 = d0 b1 log(u).
Con esto hemos calculado (si en el camino no result ser
ps no elemental) la parte polinomial de la
integral
Z
ps = bs+1 s+1 + qs s + ... + q1 + Q0
Ejemplo 1.30
Apliquemos el algoritmo a el clculo de la integral
1 + 2x log3 x + x log4 x
dx. Solo necesitax log2 x
mos la extensin = log x.
In[]:= Clear[t, P, Q, R, C];

(* t = = log x *)
P = -1 + 2 x t^3 + x t^4;
Q = x t^2;
P = -1/x + 2 t^3 + t^4;
Q = t^2;
C = PolynomialQuotient[P, Q, t];
R = PolynomialRemainder[P, Q,
t] ;
Z
(* ---------Parte Racional:
(* = 1 *)
(* = 2t + t2 *)
(*= 1/x*)
R/Q -------------*)
(* --Hermite -- *)
(*Factorizacin "libre de cuadrados" de Q*)
(*= t2 , = Q1 = t
y j = 2 *)
(*Expansin en fracciones parciales*)
ApartSquareFree[R/Q, t] ;
(* = 1/t2 , j = 2*)
R11 = -1/x;
Q1 = t;
Q1p = D[Q1, t];
(* Como j=2 *)
{d, {s, t}} = PolynomialExtendedGCD[Q1, Q1p, t]; (*s=0, t=1*)
(* 0 Q1 + t Qp = 1 *)
(* 0 Q1 + t 1/x = 1/x = T = 1 *)
S = 0;
T = -1;
(* Hermite entrega *)
Hermite = -T/Q1;
Z
(* -- Rothstein/Trager:
(S + D [ T, x ])/Q1 *)
(S + D[T, x])/Q1;
RothsteinTrager = 0;
(* ----------------
(*= 0/Q1*)
Parte polinomial:
C = t^2 + 2 t;
s = 2; p2 = 1; p1 = 2; p0 = 0;
C ------------------*)
57
58
(*Paso
s=3*)
Z
(*
p2 = x = c2 t + d2
c2 = 0 y d2 = x*)
b3 = 0; d2 = x;
q2 = d2 + b2;
(*Paso
s=2*)
Z
(*
p1 2 d2 t 0 =
2 2 = 0 = c1 t + d1 *)
c1 = 0; d1 = 0; b1 = 0;
q1 = d1 + b1;
(*Paso
s=1*)
Z
(*
p0 d1 t 0 =
0 0 = 0 = d0 *)
d0 = 0;
q0 = d0 + b1*Log[x];
(* Resultado *)
Hermite + RothsteinTrager + b3*t^3 + q2*t^2 + q1*t + q0
Out[]=
1
+ 0 + x 2 =
1
log x
| {z }
parte racional
x log2 x
| {z }
parte polinomial
Ejemplo 1.31
El algoritmo decide si una primitiva no es elemental.
1
P( )
1
. El integrando es
= Q( x, ) con = log x.
log x
Q( )
Solo debemos aplicar el mtodo de Rothstein/Trager:

c
c
R(c) = res (1 , ) = 1 Q( x )[c]
x
x
Como R(c) no tiene races constantes, la integral no es elemental.
59
Como indicabamos, los casos del algoritmo de Risch con las otras extensiones se pueden ver en [16],
[17] o [5].
Bibliografa
[1] S. Abbot. "Understanding Analysis". Springer-Verlag. 2000.
[2] M. Spivak "Calculus" W.A. Benjamin Inc. 1967. http://archive.org/details/Calculus_643.
Consultada Julio 6, 2013.
[3] "Oeuvres de P.L. Tchebychef (1899)". http://archive.org/details/117744684_001. Consultada Julio 6, 2013.
[4] Robert H. Risch. "The solution of the problem of integration in finite terms", BulL Amer. Math.
Soc., vol. 76, pp. 605-608.
[5] Keith O. Geddes, Stephen R. Czapor, George Labahn. Algorithms for Computer Algebra. Kluwer Academic Publisher. 1992.
[6] G. H. Hardy. "The Integration of Functions of a Single Variable." 2nd ed., Cambridge University Tracts
in Mathematics and Mathematical Physics, no. 2, Cambridge, England, 1916. Consultada Julio 6,
2013.http://www.gutenberg.org/ebooks/38993.
[7] P. Laplace. "Oeuvres compltes de Laplace". Tome Septime. Gauthier-Villars. 1820. http://
archive.org/stream/theorieanaldepro00laplrich#page/n193/mode/2up.Consultada
Julio 6, 2013.
[8] J. F. Ritt. "Integration in Finite Terms: Liouvilles Theory of Elementary Methods". Columbia University
Press, New York, 1948.
[9] C. Ivorra. "Funciones sin primitiva elemental".
Primitivas.pdf. Consultada Julio 6, 2013.
http://www.uv.es/ivorra/Libros/
[10] N. Piskunov. "Differential and Integral Calculus". Editorial Mir. Mosc. 3ra Ed. 1977.
https://ia801507.us.archive.org/33/items/DifferentialAndIntegralCalculus_
109/N.Piskunov-DifferentialAndIntegralCalculus1969mir.djvu. Consultada Julio 6,
2013.
[11] J. Moses. "Symbolic Integration". ACM Communications, Vol 14. Number 8. 1971. http://wwwinst.eecs.berkeley.edu/~cs282/sp02/readings/moses-int.pdf
[12] "Sign function". The Wolfram Functions Site. http://functions.wolfram.com/
ComplexComponents/Sign/introductions/ComplexComplements/ShowAll.html.
Consultada Julio 6, 2013.
[13] J. Pierpont. "Lectures on the theory of functions of a real variable. Vol 1, Ginn and Company, 1905.
http://archive.org/details/lecturesonthethe032088mbp. Consultada Julio 6, 2013.
[14] E. A. Marchisotto, G. Zakeri. "An invitation to integration in finite terms". College Mathematics
Journal, Vol 25, No. 4. (Sep, 1994).
[15] B. Conrad. "Impossibility theorems for elementary integration." Clay Institute Research
Academy Colloquium Series. http://www.claymath.org/programs/outreach/academy/
LectureNotes05/Conrad.pdf
[16] R. Risch. "The Solution of the Problem of Integration in Finite Terms." Bull. Amer. Math. Soc., 1-76,
605-608, 1970.
[17] M. Bronstein. "Symbolic Integration I:Transcendental Functions". 2nd edition. Springer-Verlag, 2005.
60
[18] M.W. Ostrogradsky. "De integration des fractions rationelles". Bulletin de la Classe PhysicoMathematiques de l"acadmie Impriale des Sciences de St. Ptersbourg, IV: 145-167,286-300, 1845.
http://books.google.es/books/about/Bulletin_physico_math%C3%A9matique.
html?hl=es&id=DegAAAAAYAAJ&utm_source=gb-gplus-share. Consultada Agosto 3, 2013.
[19] E. Hermite. "Sur Integration des fractions rationelles". Nouvelles Annales de Mathematiques (2eme
serie), 11:145-148, 1872. http://books.google.es/books/about/Nouvelles_annales_
de_math%C3%A9matiques.html?hl=es&id=GDITAQAAMAAJ&utm_source=gb-gplusshare. Consultada Agosto 3, 2013.
[20] J. Luetzen. "Joseph Liouville 1809-1882". Volume 15 of Studies in the history of mathematics and
physical sciences. Springer-Verlag, New York, 1990.
[21] T. Wang. "Techniques on Partial Fractions". Proceedings of the AMATYC 33rd Annual Conference, Minneapolis, Minnesota 2007. American Mathematical Association of Two
Year Colleges. http://c.ymcdn.com/sites/www.amatyc.org/resource/resmgr/2007_
Conference_Proceedings/wang07.pdf. Consultada Agosto 3, 2013.
[22] E. Selmer. "The Irreducibility of Certain Trinomials". http://www.mscand.dk/article.php?
id=1472. Consultada Agosto 3, 2013.
[23] J. Liouville. "Premier Memoire sur la dtermination des intgrales dont la valeur est algbrique". Journ. Ec. Polyt. 14 (22. cahier) (1833), 124-148; Mem. Savans Etrangers Acad.
Sci. Paris 5 (1838), 76-102.http://books.google.es/books?id=o8cAAAAAYAAJ&printsec=
frontcover&hl=es#v=onepage&q&f=false.Consultada Agosto 3, 2013.
[24] B. Terelius. "Symbolic Integration." Master"s Thesis.KHT. Sweden. 2009. http://www.nada.
kth.se/utbildning/grukth/exjobb/rapportlistor/2009/rapporter09/terelius_
bjorn_09095.pdf. Consultada Agosto 3, 2013.
[25] E. Horowitz. "Algorithms for partial fraction decomposition and rational function integration."
In Proceedings of SYMSAM"71, pages 441-457. ACM Press, 1971. http://research.cs.wisc.
edu/techreports/1970/TR91.pdf. Consultada Agosto 3, 2013.
[26] M. Rothstein. "A new algorithm for the integration of exponential and logarithmic functions". Proceedings of the 1977 MACSYMA Users Conference, pages 263-274. NASA Pub. CP-2012, 1977.
[27] B. Buchberger et all. Computer Algebra. Simbolic and Algebraic Computation. 2nd Edition. SpingerVerlag. 1982.


Primitivas Elementales

Cargado por

Copyright:

Formatos disponibles

Primitivas Elementales

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Primitivas Elementales

Cargado por

Copyright:

Formatos disponibles

Revista digital

Matemtica, Educacin e Internet

Se puede saber si una funcin tiene

En los cursos de clculo elemental, el problema de clculo de integrales consiste

elemental f , calcular otra funcin elemental (si existe) g tal que g0 = f (o

f = g )". La solucin par-

Integrate[1/Sqrt[1 - x^2], x] (* = i Log( 1 x2 + i x )

Aqu, F ( x, y1 , y2 , ..., yn ) es una expresin algebtraica F ( x, log x, ex , ...).

integracin en trminos finitos que fueron establecidos en el siglo XIX.

Interludio: Primitivas en anlisis matemtico

integrable (en el sentido Riemann), entonces F ( x ) =

Como las funciones continuas son integrables, tenemos: Si f : [ a, b] R es continua,

f ( x ) dx es derivable y satisface F 0 ( x ) = f ( x ) en [ a, b].

b) Tambin podemos calcular g( x ) =

ez = ex (cos y + i sen y) para todo z C

La extensin compleja de estas funciones es,

El algoritmo de Johan Bernoulli (1703)

El problema de calcular primitivas de funciones racionales es viejo, tanto como

factorizar en R[ x ] pero tiene importancia terica.

El algoritmo para calcular la integral es ([17, pg,36]),

C (no ofrece problema, pues C es un polinomio).

b) Calcular cada integral

Descomposicin en fracciones parciales?. La descomposicin en fracciones parciales se usa en inte-

gracin de funciones racionales y en el clculo de la inversa de la transformada de Laplace. Pero "una

c) Si P( x )/Q( x ) es una fraccin propia con la cuadrtica Q( x ) = ax2 + bx + c irreducible, entonces

(* Seguimos el algoritmo usando Mathematica *)

Qp[x_] = D[Q, x];

Out[]= 1/8 Log[-3 + x] - 1/8 Log[5 + x]

reducibles en Q[ x ] o cuya factorizacin en Q[ x ] no se adapta al algoritmo de Bernoulli. Algunos

Out[]= RootSum[1 + #1 + #17 &, Log[x - #1]/(1 + 7 #16 ) &]

Si las races de Q son 1 , 2 , , 7 , entonces, aplicando el mtodo de los residuos,

Aunque Q( x ) = x6 5x4 8x3 2x2 + 2x + 1 no factoriza en Q[ x ], se puede hacer reducciones que

6x5 + 6x4 8x3 18x2 + 8x + 8

R 6x5 + 6x4 8x3 18x2 + 8x + 8

Principio de Laplace, 1812.

Ms adelante va a ir quedando claro como la "forma" del integrando y de la primitiva es crucial en la

con M, P y Q polinomios en Z[ x ]. La justificacin elude a los principios de Laplace: La primitiva solo

Teorema de P. Chebyshov, 1853.

donde P, Q y R son polinomios y m es un entero positivo. Pero el artculo es ms bien conocido

Se puede saber si una funcin tiene primitiva elemental? . Walter Mora F.

Teorema 1.1 (Chebyshov, 1853).

En el caso de que la integral se puede expresar en trminos de funciones elementales, la manera de

diferencial") tiene primitiva elemental y, tambin para el clculo.

Para calcular, primero se hace la sustitucin x = u1/n , entonces dx =

a) Caso I: Si p es entero y q es racional. Si q =

puede expresar en trminos de funciones elementales.

x 3/2 (4 x )3/2 dx. Como

sen x dx no es elemental: Si u = sen x obtenemos

que, de acuerdo al teorema de Chebyshov, no se puede expresar en trminos de funciones

2u2 (u4 + 1)1 du

que si pasa el test de Chebyshov pues p = 1.

x = b, se calcula con la frmula

Primitivas elementales y funciones inversas

Aplicando la sustitucin t = f ( x ) en la integral de la derecha, obtenemos la primitiva G (t) =

A veces, abusando del lenguaje, se escribe

ln x dx no es elemental porque (como veremos ms adelante)

e1/y dx no es elemental porque (como veremos ms adelante)

funcin que satisface una ecuacin de la forma