Sociedad Colombiana de Matemáticas Lecturas Matemáticas

XV Congreso Nacional de Matemáticas Volumen Especial (2006), páginas 285–298

2005 Apuntes

Ecuaciones diofánticas y engranajes cilı́ndricos

Carlos Alberto Casas Figueroa

Universidad Javeriana, Cali, Colombia

Abstract. The problem of designing a cylindrical gear in order to

transmit fast wheel work starting from an existent ratio between the
number of teeth of two sprockets, can be approached as integer positive
solutions of Diophantine equations of the form ax + by = c.
Key words and phrases. Diophantine equations, Gears, Continued frac-
tions, Convergents.
2000 AMS Mathematics Subject Classification.Primary: 54H25. Sec-
ondary: 47H10.
Resumen. El problema de diseñar un engranaje cilı́ndrico en el cual
se transmite un movimiento rápido de rodaje o desfogue del engranaje
a partir de una relación existente entre el número de dientes de dos
piñones, se puede abordar como las soluciones enteras positivas de una
ecuación diofántica de la forma ax + by = c .

1. Introducción
Las ruedas dentadas sirven para la transmisión
de potencia o movimientos entre dos o más ejes.
La absorción de los choques en un engranaje
cilı́ndrico es igual en cada uno de los dientes,
si la relación de transformación es de la forma
a/b, donde a y b son números primos relativos
entre sı́. La obtención de esta relación de trans-
formación donde a y b son el número de dientes
del engranaje cilı́ndrico, se puede apoyar en las
ecuaciones diofánticas y el concepto fundamental
de las reducidas. El proceso abordado con las
ecuaciones diofánticas y las reducidas de una
286 Carlos Alberto Casas Figueroa

fracción es un método alternativo para el cálculo de las relaciones de transfor-

mación que conlleva a disminuir el error en las relaciones de transformación.
Inicialmente se desarrolla el punto de vista teórico del problema como so-
porte, para finalmente aplicarlo a los engranajes cilı́ndricos.

2. Ecuaciones diofánticas

Una ecuación lineal con dos incógnitas no es suficiente para determinar el

valor de éstas, si deseamos encontrar una solución única. Pues dos cantidades x
e y no están unı́vocamente determinadas más que en el caso de haber entre ellas
dos ecuaciones independientes. Pero si alguna de las incógnitas x e y, se les
exige que sean números enteros positivos, esta condición puede reemplazar la
segunda ecuación. En este caso la ecuación se llama diofántica, en recuerdo del
matemático Diofanto, quien vivió en Alejandrı́a hacia el año 275 de nuestra
era y creó el análisis indeterminado o análisis diofántico, que es un método para
buscar soluciones enteras de ecuaciones y de sistemas de ecuaciones algebráicas
con coeficientes enteros.
Reciben el nombre de ecuaciones diofánticas de primer grado con dos incóg-
nitas, las ecuaciones de la forma ax + by = c, donde las incógnitas x, y son
números enteros positivos.
Soluciones posibles: Aún en el caso de ser posible hallar valores de x e y
que, satisfaciendo la condición de ser enteros y positivos, verifiquen estas ecua-
ciones, existirá siempre, una notable diferencia entre las ecuaciones diofánticas
de primer grado y las ecuaciones ordinarias de primer grado con dos incógnitas.
En efecto ocurre que mientras:

Una ecuación lineal con dos incógnitas Una ecuación diofánticas de primer
grado
Tiene infinitas soluciones en los números 1. No tiene solución
reales 2.Tiene una o un número finito de solu-
ciones
3. Tiene infinidad de soluciones

Condición para la existencia de soluciones:

Una ecuación Diofántica ax + by = c, en la cual los coeficientes de las
incógnitas tienen un factor común, que no divide exactamente al segundo miem-
bro, no tiene solución.
Si se divide la ecuación por el factor común a esos coeficientes enteros, en el
primer miembro, se obtienen coeficientes enteros para las incógnitas, mientras
que en el segundo miembro, resulta una fracción. Como un par de valores
enteros x e y hacen que el primer miembro sea entero, serı́a preciso, para que
existiera solución, que un número entero fuera igual a una fracción. Esto no es
posible.
 Ecuaciones Diofánticas y engranajes cilı́ndricos 287

Ejemplo: Sea la ecuación 20x + 15y = 7

Factor común para 20 y 15 es 5 y no lo es para 7 luego la ecuación serı́a:
4x + 3y = 7/5
Por lo tanto: Para que una ecuación Diofántica ax+by = c tenga soluciones,
es preciso que los coeficientes a, b de las incógnitas sean números primos rela-
tivos. Es decir factores comunes para (a, b) = 1.
Ejemplo: Sea la ecuación 12x − 7y = 4.
Cumple la condición factores comunes (12, 7) = 1. Algunas posibles solu-
ciones son:
x 12 19
y 20 32
12(12) − 7(20) = 4 12(19) − 7(32) = 4
Verifiquemos:
144 − 140 = 4 228 − 224 = 4
¿Pero cómo se obtuvieron estas soluciones?
Veamos como hallar una solución cualquiera de la ecuación 127x − 52y = −1
127
Sea la relación donde a = 127 y b = 52.
52
Ejecutemos la división entre el numerador y el denominador de la fracción,
esto origina un cociente incompleto y un residuo, en seguida se efectuan di-
visiones a partir del divisor original ası́: primera división, el divisor original
dividido por el primer residuo y cada nueva división, el nuevo residuo dividirá
al anterior residuo, hasta que se logre un residuo cero.
Sean las divisiones, donde:
ai : Los cocientes incompletos
ri : Los residuos
a1 a2 a3 a4 a5
a b r1 r2 r3 r4
r1 r2 r3 r4 r5
La expresión obtenida da como resultado:

127 23 1 1 1
=2+ =2+ =2+ =2+
52 52 52/23 2 + 6/23 1
2+
23/6
1 1 1
=2+ =2+ =2+
1 1 1
2+ 2+ 2+
3 + 5/6 1 1
3+ 3+
6/5 1 + 1/5
288 Carlos Alberto Casas Figueroa

La expresión obtenida se llama fracción continua limitada o finita.

Suprimiendo el último término de esta fracción, que es 1/5, transformamos
la fracción continua que acabamos de obtener en una fracción ordinaria y la
restamos de la fracción inicial.
1 1 1 1 4 22
2+ =2+ =2+ =2+ =2+ =
1 1 8+1 9 9 9
2+ 2+
1 4 4 4
3+
1
127 22 1143 − 1144 1
− = =
52 9 52 × 9 52 × 9
Reduciendo la expresión hallada a un denominador común, se obtiene:
127 × 9 − 52 × 22 = −1 ⇒ 127(9) − 52(22) + 1 = 0
127X − 52Y + 1 = 0
Comparando esta igualdad obtenida con la ecuación. Concluimos que: x = 9
e y = 22 son solución de esta ecuación.
La validez de esta solución se apoya en el concepto de las reducidas de una
fracción. Llámase reducida de una fracción al proceso de suprimir el último
término en una fracción continua, esto significa que se obtiene de una fracción
continua, tantas reducidas como cocientes incompletos se haya podido efectuar.
Observemos en nuestro ejemplo el número de reducidas que se pueden plantear.
ai
2 2 3 1 5 Recuerde:
127 52 23 6 5 1 ai : conciente incompleto
23 6 5 1 0 ri : residuo
ri
Luego:

ai 2 2 3 1 5

La expresión obtenida da como resultado:

127 1
=2+
52 1
2+
1
3+
1
1+
5
Tenemos cinco cocientes incompletos luego podemos plantear: n = 5 reduci-
das ası́:
Sea Ri una reducida cualquiera
 Ecuaciones Diofánticas y engranajes cilı́ndricos 289

1 127 pn p5 1 22 pn−1 p4
R5 = 2 + = = = R4 = 2 + = = =
1 52 qn q5 1 9 qn−1 q4
2+ 2+
1 1
3+ 3+
1 1
1+
5
1 17 pn−2 p3 1 5 pn−3 p2
R3 = 2 + = = = R2 = 2 + = = =
1 7 q n−2 q3 2 2 qn−3 q2
2+
3
pn−4 p1
R1 = 2 = =
qn−4 q1

El cálculo de cada una de ellas se vuelve un proceso dispendioso.

¿Cómo se pueden obtener?

ai 2 2 3 1 5

R1 = a1 = 2 ⇒ p1 = 2; q1 = 1

2 × 2 + 1 (a1 ) × (a2 ) + 1 5
R2 = = = ⇒ p2 = 5; q2 = 2
2 a2 2

1 3 3×2×2+3+2
R3 = 2 + =2+ =
2×3+1 2×3+1 2×3+1
3
3(2 × 2 + 1) + 2 a3 × (a1 × a2 + 1) + a1
= =
3×2+1 a3 × a2 + 1
Observemos detenidamente P2 = a1 × a2 +1; P 1 = a1 ; q 2 = a2 ; q 1 = 1. Luego:
a3 × p2 + p1 3 × 5 + 2 17
R3 = = = ⇒ p3 = 17; q3 = 7 .
a3 × q2 + q1 3×2+1 7
Según lo anterior, entonces R4 y R5 serán:
a4 × p3 + p2 1 × 17 + 5 22
R4 = = = ⇒ p4 = 22; q4 = 9 .
a4 × 3 + 2
q q 1×7+2 9

a5 × p4 + p3 5 × 22 + 17 127
R5 = = = ⇒ p5 = 127; q5 = 52 .
a5 × q4 + q3 5×9+7 52
De lo anterior se puede plantear lo siguiente:
290 Carlos Alberto Casas Figueroa

Teorema 2.1. Utilizando el método de inducción matemática se demuestra

que las relaciones de este mismo tipo:
Pk = ak Pk−1 + Pk−2 ; q = ak qk−1 + qk−2 (2.1)
son válidas para k ≥ 3.

Demostración: Supongamos que las igualdades (??) se cumplen para cualquier

k ≥ 3. De la definición de las fracciones reducidas se desprende directamente,
1
que sustituyendo en la expresión Rk las magnitudes ak por: ak + ak+1 . La
expresión Rk se convierte en Rk+1 . Por inducción tenemos que:
pk ak pk−1 + pk−2
Rk = q = q
k ak k−1 + qk−2
1
Sustituyendo aquı́ ak por ak + ak+1 resulta


1 pk−1
ak + pk−1 + pk−2 ak pk−1 + + pk−2
ak+1 ak+1
Rk+1 =
 =
1 qk−1
ak + qk−1 + qk−2 ak qk−1 + + qk−2
ak+1 ak+1
pk−1 pk−1
ak pk−1 + pk−2 + pk +
ak+1 ak+1
Rk+1 = =
qk−1 qk−1
ak qk−1 + qk−2 + qk +
ak+1 ak+1
(ak+1 ) pk + pk−1
ak+1 (ak+1 ) pk + pk−1 pk+1
Rk+1 = = =
(ak+1 ) q k + qk−1 (ak+1 ) q k + qk−1 qk+1
ak+1
Ası́, pues, de la validez de las igualdades (??) para cualquier k ≥ 3, se deduce
su validez para k + 1.


3. Propiedad de la diferencia entre dos reducidas consecutivas

La diferencia de dos reducidas es:

(k)
(−1)
Rk − Rk−1 = q q (k > 1) .
k k−1

En efecto,
pk pk−1 pk qk−1 − qk pk−1
Rk − Rk−1 = q − q =
qk qk−1
k k−1
 Ecuaciones Diofánticas y engranajes cilı́ndricos 291

Apoyándonos en las igualdades (??) reemplacemos k por k − 1

Pk q−1 − qk Pk−1 = (ak Pk−1 + Pk−2 )qk−1 − (ak qk−1 + qk−2 )Pk−1
= ak Pk−1 qk−1 + qk−1 Pk−2 − ak Pk−1 qk−1 − Pk−1 qk−2
= −(Pk−1 qk−2 − qk−1 Pk−2 ) ,
reduciendo términos. Repitiendo la misma transformación en las expresiones
obtenidas, resulta una sucesión de igualdades:
(Pk qk−1 − qk Pk−1 ) = (−1)(Pk−1 qk−2 − qk−1 Pk−2 )
= (−1)2 (Pk−2 qk−3 − qk−2 Pk−3 )
= (−1)3 (Pk−3 qk−4 − qk−3 Pk−4 ) . . .
= (−1)k−2 (P2 q1 − q2 P1 ) ;
como P2 = a1 a2 + 1, q2 = a2 , P1 = a1 , q1 = 1,
(Pk qk−1 − qk Pk−1 ) = (−1)k−2 ((a1 a2 + 1) × 1 − a2 (a1 )) = (−1)k−2 .
De esto se deduce que:
k−2 k
pk qk−1 − qk pk−1 (−1) (−1)
Rk − Rk−1 = qk qk−1
=
qk qk−1
= 2 .
(−1) q k qk−1
k
(−1)
Rk − Rk−1 = q q . (3.1)
k k−1

Cuando el desarrollo de a/b en fracción continua posee n términos, la

n−ésima fracción reducida Rn coincide con a/b. Por lo tanto Rk − Rk−1 =
k
(−1)
si k = n, entonces;
qk qk−1
a
= Rn ,
b
luego
n
a (−1)
− Rk−1 = .
b qn qn−1
Apliquemos lo anterior a la ecuación ax − by = c, con factores comunes para
(a, b) = 1. Luego a/b es continua y a/b = Rn ; entonces
n
a pn−1 (−1)
− =
b qn−1 qn qn−1
haciendo común denominador y simplificando
n
aqn−1 − bpn−1 (−1)
=
bqn−1 qn qn−1
como a/b y pn−1 /q n−1 son reducidas consecutivas, entonces b = qn .
292 Carlos Alberto Casas Figueroa

Por lo tanto
aqn−1 − bPn−1= (−1)n (3.2)
n−1
multiplicando la expresión (??) por (-1) C tenemos, donde C es una cons-
tante entera
a[(−1)n−1 Cqn−1 ] − b[(−1)n−1 CPn−1 ] = (−1)n (−1)n−1 C
= (−1)2n−1 C,
como 2n − 1 es impar
a[(−1)n−1 Cqn−1 ] − b[(−1)n−1 CPn−1 ] = −C
a[(−1)n−1 Cqn−1 ] − b[(−1)n−1 CPn−1 ] + C = 0
a(x) − b(y) + C = 0.
De donde X = (−1)n−1 Cqn−1 ; Y = (−1)n−1 CPn−1 , conforman una solución
de esta ecuación.
Para otras soluciones se utiliza lo siguiente:
Teorema 3.1. Sean a, b números primos entre sı́ y [xo , yo ] es una solución
arbitrarı́a de la ecuación lineal aX − bY + C = 0. Entonces:
X = xo + bt; Y = yo + at,
t = 0, ±1, ±2, ±3, . . . , ±n, serán otras soluciones.

Demostración: Sea [x, y] solución arbitraria de la ecuación, entonces de las

igualdades ax − by + c = 0 y axo − byo + c = 0 tenemos que:
a(x − x0 )
ax − axo − by + byo = 0, y − y0 = , (3.3)
b
como y − yo es un número entero; a y b números primos entre sı́, la expresión
x−xo debe dividirse por b sin resto, es decir, x−xo adquiere la forma x−xo = bt
en la que t es número entero. Por lo tanto, x = xo + bt; reemplazando este
valor en (??) tenemos:
abt
y − y0 = ⇒ y = y0 + at.
b

4. Aplicaciones

Ruedas dentadas y transmisiones por medio de ellas mismas. Las

ruedas dentadas sirven para la transmisión directa de potencia o movimientos
entre dos o más ejes. La transmisión del esfuerzo tiene lugar por medio del
engrane de los flancos de los dientes.
La velocidad tangencial común es
π π π π
v = v1 = v2 = do1 n1 = do2 n2 = ro1 n1 = ro2 n2 [m/s]
60 60 30 30
 Ecuaciones Diofánticas y engranajes cilı́ndricos 293

Si do y ro se miden en metros. De aquı́ se deduce:

n1 do2 ro2
= = ,
n2 do1 ro1
con v1 = ro1 w1 = ro2 w2 .
La relación de transformación vale:
w1 do2 ro2 n1 z2
i= = = = = .
w2 do1 ro1 n2 z1
Para la elección de la relación de transformación hay que tener en cuenta las
siguientes consideraciones: Si se requiere un rápido rodaje o desfogueo del
engranaje, la relación de transformación deberá ser un número entero, porque
ası́ siempre engranan los mismos pares de dientes.
Para i = 1/1 son los mismos en cada vuelta; para i = 2/1 cada diente de
la rueda menor engrana alternativamente con solo dos dientes de la rueda an-
tagonista, pero si presentan oscilaciones periódicas en el esfuerzo a transmitir
(prensas de rodillera, prensas de estampar, . . . ) entonces las transformaciones
por números enteros han de evitarse a toda costa, de no hacerlo, siempre so-
portarı́an las puntas de carga los mismos pares de dientes y se desgastarı́an
prematuramente.
En tales casos, a ser posible, cada diente de la rueda menor debe engranar
sucesivamente con todos los de la mayor, de manera que todos los dientes se
repartan por igual la absorción de los choques. Esto se consigue con relaciones
de transformación que no se pueden simplificar, por ejemplo 19 : 45.

Las ecuaciones diofánticas nos permite calcular engranajes cilı́ndricos

dada la relación de velocidad o relación de transmisión de ellos.

Ejemplo N◦ 1. Se desean construir dos ruedas dentadas para un engranaje

que den una relación de transmisión: i = 131/85, esto podrı́a realizarse, si una
de las ruedas tuviese 131 dientes y la otra 85; no obstante, no conviene que las
ruedas tengan menos de 10 dientes ni más de 30. Calcular el número de dientes
que debe tener cada rueda para que esto se cumpla.
294 Carlos Alberto Casas Figueroa

Por la propiedad de la diferencia de reducidas consecutivas y el Teorema ??,

tenemos:
n
pn pn−1 (−1)
− = i = 131/85 = 1.541
qn qn−1 qn qn−1
Pn = 131, qn = 85, Pn−1 = y, qn−1 = x. Por lo tanto
n
131 y (−1) n
− = ⇒ 131x − 85y = (−1) . (4.1)
85 x 85x
Como (131, 85) = 1, (??) es una ecuación diofántica. Calculemos las reducidas:
ai
1 1 1 5 1 1 3
131 85 46 39 7 4 3 1
46 39 7 4 3 1 0
ri

1 p1 2 p2
R1 = =
q1 R2 = =
q2
1 1
a3 p2 + p1 3 p3 a4 p3 + p2 17 p4
R3 = q = = R4 = q = =
a3 2 + q1 2 q3 a4 3 + q2 11 q4
a5 p4 + p3 20 p5 a6 p5 + p4 37 p6
R5 = q q
= =
q5 R6 = q = =
a5 4 + 3 13 a5 5 + q4 24 q6
p
a7 6 + 5 p 131 p7
R7 = q = =
a7 6 + q5 85 q7
De donde Pn−1 = 37 y qn−1 = 24. Como la ecuación es de la forma ax − by = c
entonces tenemos:
131(X) − 85(Y ) + (−1)6 = 0.
La solución es:
X = (−1)n−1 CQn−1 = (−1)6 (1) ∗ 24 = 24
Y = (−1)n−1 CPn−1 = (−1)6 (1) ∗ 37 = 37
Pero el ejercicio nos pide que las ruedas no deben tener menos de 10 dientes
ni más de 30. Luego i = 24/37, entonces planteamos una nueva ecuación
 Ecuaciones Diofánticas y engranajes cilı́ndricos 295

diofántica
n
24 y (−1) n
− = ⇒ 24x − 37y = (−1) . (4.2)
37 x 37x
pn−1 11
Calculemos las reducidas: Calculadas sus reducidas, se halla que = ,
qn−1 17

0 1 1 1 5 2
24 37 24 13 11 2 1
24 13 11 2 1 0

de donde, las ruedas deben tener 11 y 17 dientes respectivamente. Luego

i = 17/11 = 1.545.
Error 131/85 − 17/11 = −0.0042 por exceso.
Error relativo: 0.0042/(131/85) = 0.0027.
Ejemplo N◦ 2. Para calcular el menor número de vueltas n0 que se debe dar
al husillo desembragado del carro de un torno, para pasar de una a otra entrada
al roscar un tornillo de filete múltiple de cuatro entradas, de paso 3/2, siendo
la constante del torno 5/8, se llega a la ecuación:
4(5/8)n0 − 4(3/2)n1 = 3/2 ⇔ 5n0 − 12n1 = 3;
entonces 5(n0 )−12(n1 )−3 = 0, que es una ecuación diofántica, con (5, 12) = 1.
Calculemos las reducidas:

0 2 2 2
5 12 5 2 1
5 2 1 0

0 p1 1 p2
R1 = =
q1 R2 = =
q2
1 2
a3 p2 + p1 2 p3 a4 p3 + p2 5 p4
R3 = q q
= =
q3 R4 = q q
= =
q4
a3 2 + 1 5 a4 3 + 2 12

Como no = (−1)3 Cqn−1 = (−1)3 (−3)5 = 15 y n1 = (−1)3 CPn−1 =

(−1)3 (−3)2 = 6 y nos piden el menor valor de no , entonces por el Teorema
??,
no = 15 + 12t =⇒ no = 15 + 12(−1) = 3.
n1 = 6 + 5t =⇒ n1 = 6 + 5(−1) = 1.

Ejemplo No 3. Se necesitan dos ruedas: Una de 237 dientes y otra de 295.

No disponemos más que de ruedas de 10, 15, 20, 25, 30, 35, y 40 dientes.
296 Carlos Alberto Casas Figueroa

Determinar los que debe tener la otra para que con ellas se pueda efectuar el
trabajo.
237
Si i = , entonces los posibles juegos de ruedas son:
295
(10, 15); (10, 20); (10, 25); (10, 30);
(10, 35); (10, 40)
(15, 20); (15, 25); (15, 30); (15, 35); (15, 40)
(20, 25); (20, 30); (20, 35); (20, 40)
(25, 30); (25, 35); (25, 40)
(30, 35); (30, 40)
(35, 40)
n
237 y (−1) n
Tenemos: − = ⇒ 237x − 2957y = (−1) 237 Y (-1)n .
295 x 295x

0 1 4 11 1 1 2
237 295 237 58 5 3 2 1
237 58 5 3 2 1 0

Podemos observar que aquı́ se generarı́an siete (7) reducidas. Luego la

ecuación es 237X − 295Y = (−1)7 y entonces tenemos:
237X − 295Y + 1 = 0 ; (4.3)
237 94
la relación es i = = 0.8033, R6 = , donde Pn−1 = 94 y Qn−1 = 117
295 117
que están muy lejos de lo que se tiene.
Planteamos otra ecuación:
94X − 117Y = (−1)n (4.4)
se calculan sus reducidas

0 1 4 11 2
94 117 94 23 2 1
94 23 2 1 0
 Ecuaciones Diofánticas y engranajes cilı́ndricos 297

45
R4 = , que está lejos de nuestra condición 94X − 117Y = (−1)5 .
56
De nuevo planteamos una nueva ecuación:
45X − 56Y = (−1)n (4.5)
se calculan sus reducidas

0 1 4 11
45 56 45 11 1
45 11 1 1 0

4
Si R3 = estos valores están por debajo de lo deseado, pero 4 y 5 sı́ son
5
factores a la vez de 20 y 25, respectivamente, que es uno de los juegos de ruedas
que se puede combinar. Por tanto Pn−1 = 4 y Qn−1 = 5.
Sustituyamos en la ecuación (??)

45X − 56Y = (−1)4 ⇔ 45X − 56Y = 1

multiplicada por 5
45(5Qn−1 ) − 56(5Pn−1 ) = 5
Por lo tanto
X = 5Qn−1 = 25
Y = 4Pn−1 = 20.
20
entonces i = = 0.80.
25
237 20
Error − = 0.003 por defecto.
295 25
0.003
Error relativo = = 0.004
237/295

Bibliografı́a
[1] J. Shigley. Diseño en Ingenierı́a Mecánica. Editorial Mc. Graw Hill: México, 1995.
[2] R. Norton. Diseño de Maquinaria. Editorial Mc. Graw Hill: Madrid, 2000.
[3] K. H. Rosen. Matemática discreta y sus aplicaciones. Editorial Mc. Graw Hill:
Madrid, 2004.
[4] A. O. Guelfond. Resolución de ecuaciones en números enteros. Editorial Mir:
Moscú, 1979.
[5] R. P. Grimaldi. Matemática discreta y combinatoria. Addison–Wesley Iberoameri-
cana: México, 1997.
[6] K. H. Rosen. Elementary Number Theory and its applications, 2nd. Ed. Addison-
Wesley: Boston, 1998.
[7] D. M. Burton. Elementary Number Theory. (3a Ed.). Wm. C. Brown: Boston, 1994.
298 Carlos Alberto Casas Figueroa

(Recibido en mayo de 2006. Aceptado para publicación en diciembre de 2006)

Facultad de Ingenierı́a
Universidad Javeriana
Cali, Colombia
e-mail: [email protected]