Introducción A Las Ecuaciones en Diferencias
Introducción A Las Ecuaciones en Diferencias
Introducción A Las Ecuaciones en Diferencias
+ = A si B
A
A
U A n U
n
n
1 , ) 0 ( ) ( = + = A si nB U n U
Si 1 = A , esta concuerda con la ecuacin (5), si U(0) se escribe como C.
xB C x U + = ) ( .
13
Si 1 A , se escribe en la forma
A
B
A
B
U A n U
x
+
)
`
=
1 1
) 0 ( ) ( .
Si hacemos C
A
B
U =
1
) 0 ( , entonces
1 ,
1
) (
+ = A
A
B
CA n U
n
Encontrando as la solucin general de la ecuacin (7) con una constante
arbitraria C.
Definicin 1.11: Si 0 = C entonces ) (n U se reduce a la funcin constante
A
B
1
. Este se llama valor de equilibrio de ) (n U y se denota por ) (
*
n U ;
A
B
n U
=
1
) (
*
.
Por consiguiente la solucin general se escribe entonces de la forma
) ( ) (
*
n U CA n U
n
+ = .
Ejemplo 1: 5 2 ) ( 3 ) 1 ( = + n n U n U .
En este caso, se sugiere que
2 1
*
) ( K n K n U + = , como el segundo miembro
de la ecuacin. Entonces
{ } { } ) 2 ( 2 3 ) 1 ( ) ( 3 ) 1 (
2 1 1 2 1 2 1
* *
K K n K K n K K n K n U n U + = + + + = +
comparando el coeficiente de n y el termino constante de esta expresin con
los correspondientes de 5 2 n , se tiene
5 2 , 2 2
2 1 1
= + = K K K .
De donde se obtiene 2 1
2 1
= = K y K . Por lo tanto 2 ) (
*
+ = n n U es
solucin particular. Una solucin general es
2 3 ) ( + = n C n U
n
.
Ejemplo 2:
n
n U n U 3 ) ( 3 ) 1 ( = + .
Puesto que 3 = = a , se supone que
n
n K U 3
1
*
= . Entonces
n n n
K n K n K n U n U 3 3 3 3 ) 1 ( ) ( 3 ) 1 (
1
1
1
1
1
* *
= + = +
+
.
Por lo tanto, 1 3
1
= K , esto es,
3
1
1
= K y
n
n n U 3
3
1
) (
*
= . Una solucin
general es
n
n C n U 3
3
1
) ( |
.
|
\
|
+ = .
14
1.2.2 Ecuaciones en diferencias y razones comunes no constantes
Teorema 1.4:
Sea ,... 2 , 1 , 0 0 ) ( ) ( ) 1 ( = = + n n U n A n U (8)
para los ) (n A tales que no se hace 0 para ningn valor de n.
Demostracin: Si se busca la solucin sucesivamente por el mtodo de
iteracin, partiendo de U(0), se tiene
),.... 2 ( ) 1 ( ) 0 ( ) 0 ( ) 2 ( ) 2 ( ) 3 (
), 1 ( ) 0 ( ) 0 ( ) 1 ( ) 1 ( ) 2 (
), 0 ( ) 0 ( ) 1 (
A A A U A U U
A A U A U U
A U U
= =
= =
=
Normalmente se adopta el smbolo para representar un producto de
trminos sucesivos de una sucesin { } ) (n A , y se define
=
+ =
i
m k
i A i A m A m A k A ) ( ) 1 ( )... 1 ( ) ( ) ( ,
para i m .
Entonces
, ) ( ) 0 ( ) 3 (
, ) ( ) 0 ( ) 2 (
, ) ( ) 0 ( ) 1 (
2
0
1
0
0
0
=
=
=
=
=
=
k
k
k
k A U U
k A U U
k A U U
A partir de estas expresiones se infiere una solucin general de la ecuacin
(8),
=
= = =
1
0
,... 2 , 1 ), 0 ( ), ( ) (
n
k
n U C k A C n U .
Esta se generaliza en
=
+ + = = =
1
,... 2 , 1 ), ( ), ( ) (
n
k
n U C k A C n U
para una sucesin que empiece desde = n .
15
1.3 Mtodos de solucin por Diferencias y Sumas
1.3.1 Diferencias
Definicin 1.12: Dada una funcin U y una constante h talque h n +
pertenezca al dominio de dicha funcin, denominaremos primera diferencia
de U a aquella funcin cuyo valor en n viene dado por
) ( ) ( ) ( n U h n U n U + = ,
designando as a esta funcin por U , y su valor particular en n lo
simbolizaremos por ) (n U .
El smbolo es llamado operador diferencia y colocado antes de la funcin
indica que esta se ha trasformado en la funcin diferencia. El numero h
recibe el nombre de intervalo de diferencia.
En nuestro caso 1 = h entonces la ecuacin seria:
) ( ) 1 ( ) ( n U n U n U + = .
1) La primera diferencia correspondiente a la funcin n n U = ) ( es:
n n n U + = ) 1 ( ) ( sea 1 = n .
2) La diferencia correspondiente a un valor constante C esta definida
por:
0 ) ( = = C C C U , entonces 0 = C .
3) La diferencia correspondiente a una funcin
n
a esta dada por:
n n n n
a a a a a U ) 1 ( ) (
1
= =
+
,
por lo que se tiene que
n n
a a a ) 1 ( = .
Ejemplo 1:
Si 3 2 ) ( + = n n U
, 2 ) 3 2 ( } 3 ) 1 ( 2 { ) ( = + + + = n n n U
Ya que . 3 ) 1 ( 2 ) 1 ( + + = + n n U
Adems K , 9 ) 3 ( , 7 ) 2 ( , 5 ) 1 ( , 3 ) 0 ( = = = = U U U U
Entonces K , 2 7 9 ) 2 ( , 2 5 7 ) 1 ( , 2 3 5 ) 0 ( = = = = = = U U U
La grafica de 3 2 + n es una recta. La diferencia 2 se llama pendiente de la
recta, que es precisamente el coeficiente de n .
1.3.2 Diferencias de polinomios
16
Definicin 1.13: El calculo en diferencias de un polinomio se reduce a
calcular la diferencia de
k
n . Sin embargo, la diferencia no resulta expresada
con una formula muy sencilla.
Ejemplo 2: . 1 3 3 ) 1 (
2 3 3 3
+ + = + = n n n n n
. 1 4 6 4 ) 1 (
2 3 4 4 4
+ + + = + = n n n n n n
Una funcin que es fundamental para calcular la diferencia y suma de
polinomio es la funcin factorial
) (k
n , definida como
) 1 (
)! (
!
) 2 )( 1 (
) 1 (
1
) (
) 3 (
) 2 (
) 1 (
) 0 (
k n
n
n
n n n n
n n n
n n
n
k
=
=
=
=
=
M
Ejemplo 3:
. 2 3 ) 2 )( 1 (
) 1 (
2 3 ) 3 (
2 2
) 1 (
n n n n n n n
n n n n n
n n
+ = =
= =
=
Para poder calcular la diferencia de la ecuacin (1), entonces tenemos para
1 k :
.
))! 1 ( (
!
)! 1 (
!
)! 1 (
! ) 1 ( )! 1 (
)! (
!
)! 1 (
)! 1 (
) 1 (
) 1 (
) ( ) ( ) (
=
=
+
=
+
+ +
=
+
+
=
+ =
k
k k k
kn
k n
n
k
k n
n
k
k n
n k n n
k n
n
k n
n
n n n
Por lo tanto.
) 2 ( ) , 2 , 1 ( ,
) 1 (
K = =
k kn n
k
.
17
Ejemplo 4: . 3 3
2 ) 1 3 ( 3
n n n = =
. 4 4
3 ) 1 4 ( 4
n n n = =
La funcin
|
|
.
|
\
|
k
n
definida por ,...). 2 , 1 , 0 (
!
) (
= =
|
|
.
|
\
|
k
k
n
k
n
k
Se llama coeficiente binomial. De la ecuacin (2)
.
)! 1 ( ! !
) 1 ( ) 1 ( ) (
= =
k
n
k
kn
k
n
k k k
Por lo tanto .
1
|
|
.
|
\
|
=
|
|
.
|
\
|
k
n
k
n
1.3.3 Sumas
Definicin 1.14: Sean ) ( ) ( n u y n U dos funciones talque ) ( ) ( n u n U =
entonces ) (n u se llama diferencia de ) (n U . Recprocamente, ) (n U se
llama una suma indefinida de ) (n u .
As la solucin general de ) ( ) ( n u n U = , se denota por ) (
1
n u
.
Entonces si ) ( ) ( n u n U = tenemos que
C n U n u + =
) ( ) (
1
(3)
donde C es una constante arbitraria
2
.
Ejemplo 5:
3
4
4
n
n
=
4
) 4 (
) 3 ( 1
.
Las operaciones para calcular la diferencia y obtener una suma de ciertas
funcin, se representa escribiendo los smbolos y
1
, respectivamente,
antes de dicha funcin. En este sentido y
1
se llaman operadores.
Si se suma la diferencia ) (n U , se obtiene )) ( (
1
n U
. Esta expresin se
escribe en la forma ) (
1
n U
y
1
se llama producto de los operadores
y
1
. El producto
1
se define de manera semejante. El producto de dos
operadores significa que se realizan dos operaciones sucesivas.
2
Nota: si se denota una solucin particular de la ecuacin en diferencias con
la propia ) (n U , la solucin general de la ecuacin en diferencias es
C n U + ) ( . Pero si se permite que la solucin n tome valores no enteros,
entonces C debe ser una funcin arbitraria con periodo 1.
18
Aunque tanto
1
como
1
son el producto de y
1
, no son
operadores iguales, es decir, no representan la misma operacin. Si se
sustituye la ecuacin (3) en el primer miembro de la segunda ecuacin para
eliminar ) (n u , entonces resulta
C n U n U + =
) ( ) (
1
. (4)
Puesto que, por definicin, ) (
1
n u
es una solucin de ) ( ) ( n u n U = ,
entonces
) ( ) (
1
n u n u =
(5)
Como ) (n u es arbitraria se puede remplazar con ) (n U . Entonces de las
ecuaciones (4) y (5),
C n U n U =
) ( ) (
1 1
.
Por consiguiente, en general ) ( ) (
1 1
n U n U
.
1.3.4 Algunas propiedades de los operadores
Teorema 1.5 : Ley conmutativa respecto a las constantes. Si C es una
constante, la diferencia de la funcin ) (n CU es igual al producto de C por
la diferencia de ) (n U , sea:
| | ) ( ) ( n U C n CU =
Demostracin:
| |
| |
) (
) ( ) 1 (
) ( ) 1 ( ) (
n U C
n U n U C
n CU n CU n CU
=
+ =
+ =
Ejemplo 6:
) 1 2 ( 5 5 ) 5 (
, 3 3 ) 3 (
2 2
+ = =
= =
n n n
n n
Teorema 1.6: Ley distributiva. La diferencia de la suma de dos funciones es
igual a la suma de sus diferencias, es decir si
2 1
U y U son dos funciones, se
tendra:
| | ) ( ) ( ) ( ) (
2 1 2 1
n U n U n U n U + = +
19
Demostracin:
| | | | | |
| | | |
) ( ) (
) ( ) 1 ( ) ( ) 1 (
) ( ) ( ) 1 ( ) 1 ( ) ( ) (
2 1
2 2 1 1
2 1 2 1 2 1
n U n U
n U n U n U n U
n U n U n U n U n U n U
+ =
+ + + =
+ + + + = +
Colorario: Dadas k funciones
k
U U U U K ; , ,
3 2 1
y k constantes
arbitrarias,
k
C C C C K ; , ,
3 2 1
se tendr
| | ) ( ) ( ) ( ) ( ) ( ) (
1 21 2 1 1 1 21 2 1 1
n U C n U C n U C n U C n U C n U C
k k k k
+ + + = + + + K K
para todo valor entero positivo k .
1.3.5 Aplicacin de la sumatoria.
Definicin 1.15: Suma de series. Si en la ecuacin
=
+ + = = =
1
,.... 2 , 1 ), ( ), ( ) (
n
k
n U C k R C n U
=
=
1
) ( ) ( ) (
n
k
k u U n U
.
Se puede considerar que esta formula dice que la suma de la serie del
segundo miembro esta representada con las sumas de ) (n u del primer
miembro, debido a que si C n U n u + =
) ( ) (
1
, entonces
) ( ) ( ) ( ) (
1 1
U n U u n u + =
, lo cual implica que
=
=
1
1 1
) ( ) ( ) (
n
k
u n u k u
.
Hagamos n n = 1 y m = , donde m y n son enteros no negativos.
Entonces la ecuacin anterior se convierte en
=
+ =
n
m k
m n m u n u k u ) ( ) 1 ( ) (
1 1
.
El segundo miembro de la ecuacin se denota por | |
1
1
) (
+
n
m
n u y se llama
suma definida de ) (n u de m a ) 1 ( + n .
Soluciones generales de ecuaciones en diferencias de primer orden lineales.
El problema de resolver una ecuacin en diferencias de primer orden lineal
se puede reducir a un problema de sumatoria. Por simplicidad, se considera
una ecuacin en diferencias con coeficientes constantes,
0 ) ( ) ( ) 1 ( = + a n R n aU n U . (6)
20
Una solucin particular de la ecuacin homognea 0 ) ( ) 1 ( = + n aU n U ,
que resulta de la ecuacin (6) haciendo que ) (n R sea cero, es
n
a . Sea ahora
) ( ) ( n V a n U
n
=
y obtngase una ecuacin en diferencias para ) (n V . El primer miembro de
la ecuacin (6) resulta ser
{ } ) ( ) 1 ( ) ( ) 1 ( ) ( ) 1 (
1 1
n V n V a n V a a n V a n aU n U
n n n
+ = + = +
+ +
dividiendo ambos miembros de la ecuacin (6) por
1 + n
a , se tiene
) ( ) (
1
n R a n V
n
= .
Por consiguiente, { } C n R a n V
n
+ =
) ( ) (
1 1
.
Si se multiplica esta expresin por
n
a , se tiene
)
`
+ =
n
n n
a
n R
a Ca n U
) (
) (
1 1
, (7)
que es la solucin general de la ecuacin (6). El primer termino es
precisamente la solucin general de 0 ) ( ) 1 ( = + n aU n U .
En el caso que el coeficiente no es constante, basta reemplazar
1 n
a ,
n
a y el
primer termino con
=
1
0
) (
n
k
k A ,
=
n
k
k A
0
) ( y la ecuacin
=
= = =
1
0
,... 2 , 1 ), 0 ( ), ( ) (
n
k
n U C k A C n U ,
respectivamente.
Ejemplo: Resolver 1 ) ( ) 1 ( = + a b n aU n U .
Sustityase b por ) (n R en la ecuacin (7). Entonces
a
b
Ca
a
a
b a Ca
a
b a Ca
a
b
a Ca n U
n
n
n n
n
n n
n
n n
+ = |
.
|
\
|
+ =
|
.
|
\
|
+ = |
.
|
\
|
+ =
+
1
1
1
1
1
1
) (
1
1 1 1 1
.
1.3.6 Diferencias de segundo orden y orden N.
Definicin 1.16: A partir de la funcin ) (n U , se define la funcin ) (n U ,
y si a su vez se puede definir su diferencia | | ) (n U . Esta se llama
diferencia de segundo orden de ) (n U y se escribe ) (
2
n U .
2
representa el
cuadrado del operador diferencia . De manera semejante se puede
21
definir ) (n U
N
por medio de | | ) (
1
n U
N
para cualquier numero natural N.
por ejemplo, si
2
) ( n n U = ,
{ }
. 0 , 0 2
, 2 ) 1 2 ( , 1 2
2 5 2 4 2 3
2 2 1 2 2
= = = = =
= + = + = + =
K n n n
n n n n n n
De la definicin anterior es obvio que
N
es un operador lineal. Esto es,
{ } ) ( ) ( ) ( ) ( n V b n U a n bV n aU
N N N
+ = + .
Ahora bien, dado que ) ( ) 1 ( ) ( n U n U n U + = ,
{ } { }
), ( ) 1 ( 3 ) 2 ( 3 ) 3 (
) ( ) 1 ( 2 ) 2 ( ) 1 ( ) 2 ( 2 ) 3 (
) ( ) 1 (
2
n U n U n U n U
n U n U n U n U n U n U
n U n U
+ + + + =
+ + + + + + + =
+ =
y asi sucesivamente.
Una tabla donde se enlistan los valores de K ), ( ), ( ), (
2
n U n U n U , se llama
tabla de diferencias. Se acostumbra colocar el valor ) (a U entre los
renglones U(a) y U(a+1).
n 0 1 2 3 4 5
) (n U
) (n U
) (
2
n U
) (
3
n U
0 1 4 9 16 25
1 3 5 7 9
2 2 2 2
0 0 0
Tabla 1.1 Diferencias de
2
x .
Ya sabemos como expresar un polinomio ) (n U de grado N en trminos de
funciones factoriales. Ahora se demuestra que los coeficientes se
representan utilizando ) 0 ( ),..., 0 ( ), 0 ( U U U
N
. Sea
N
k
n A n A n A n A A n U + + + + + = ... ) (
3
3
2
2
1
1 0
y calculamos sucesivamente las diferencias
. ! ) (
) 2 )( 1 ( ... 1 2 3 ) (
) 1 ( ... 2 3 1 2 ) (
... 3 2 ) (
) 0 (
) 3 (
3
3
) 2 ( ) 1 (
3 2
2
) 1 ( ) 2 (
3
) 1 (
2 1
n A N n U
n A N N N A n U
n A N N n A A n U
n NA n A n A A n U
N
N
N
N
N
N
N
N
=
+ + =
+ + + =
+ + + + =
M
22
se sustituye 0 = n en las formas anteriores, puesto que 0 0 =
N
, se tiene
N
N
A N U A U A U A U A U ! ) 0 ( ,..., ! 3 ) 0 ( , ! 2 ) 0 ( , ) 0 ( , ) 0 (
3
3
2
2
1 0
= = = = = .
De aqu se determinan de inmediato
N
A A A ,..., ,
1 0
.
Si ) (n U es un polinomio en la variable n, entonces
) 0 (
!
... ) 0 (
! 3
) 0 (
! 2
) 0 ( ) 0 ( ) (
) (
3
) 3 (
2
) 2 (
) 1 (
U
N
n
U
n
U
n
U n U n U
N
N
+ + + + + = .(8)
Y esta ecuacin se llama formula de interpolacin de Newton. Si se emplean
coeficientes binomiales, se escribe en la forma
) 0 ( ... ) 0 (
3
) 0 (
2
) 0 (
1
) 0 ( ) (
3 2
U
N
n
U
n
U
n
U
n
U n U
N
|
|
.
|
\
|
+ +
|
|
.
|
\
|
+
|
|
.
|
\
|
+
|
|
.
|
\
|
+ = . (9)
Si se usa una tabla de diferencias, la ecuacin (8) sirve para representar un
polinomio en trminos de funciones factoriales.
Ejemplo: representar 2 5
3
+ n n en trminos de funciones factoriales.
N 0 1 2 3
) (n U
) (n U
) (
2
n U
) (
3
x U
2 -2 0 14
-4 2 14
6 12
6
Tabla 1.2 Tabla divisin sinttica.
) 3 ( ) 2 ( ) 1 (
) 3 ( ) 2 (
) 1 (
3 4 2
! 3
6
! 2
6 4 2 ) (
n n n
n n
n n U
+ + =
+ + =
Dada una funcin, utilizando la formula de Newton, se puede construir un
polinomio que adquiera los mismos valores que tiene la funcin cuando
k n ,..., 2 , 1 , 0 = . Se puede considerar al polinomio como una aproximado de
la funcin en un punto n que no sea entero. En este sentido, las ecuaciones
(8) o (9) se les llama algunas veces formula de interpolacin de Newton.
23
2. ECUACIONES EN DIFERENCIAS LINELAES DE ORDEN
SUPERIOR
2.1 Ecuaciones en diferencias lineales
Definicin 2.1: La ecuacin
) ( ) ( ) ( ) 1 ( ) ( ..., ) 1 ( ) ( ) ( ) (
1 1 0
n R n U n A n U n A N n U n A N n U n A
N N
= + + + + + +
se llama ecuacin lineal en diferencias de una funcin incgnita ) (n U .
Definicin 2.2: Si tanto ) (
0
n A como ) (n A
N
son funciones que nunca se
hacen cero, se dice que la ecuacin es de N-simo orden o una ecuacin de
orden N. ) (n R se llama trmino no homogneo, y si es la constante es 0, la
ecuacin se llama homognea. Si todas las funciones
) ( ) ( ), ( ),..., ( ), (
1 1 0
n R n A n A n A n A
N N
=
=
.
Por lo tanto, dados y arbitrariamente, si se determinan
1
C y
2
C de las
ecuaciones anteriores, se cumple que
= + = = + =
2 2 1 1 2 1
) 1 ( , ) 0 ( e C e C U C C U ,
3
Denominadas races mltiples.
27
esto es, se satisface ciertamente la condicin inicial. De esta manera, se ha
demostrado que la ecuacin (7) satisface 2. En el caso en el que la raz es
mltiple, se puede seguir el mismo procedimiento (Verificar).
2.1.2 Soluciones particulares de ecuaciones en diferencias lineales de
segundo orden con coeficientes constantes
Definicin 2.6: Una solucin general de la ecuacin no homognea
0 ), ( ) ( ) 1 ( ) 2 ( = + + + + b n R n b n aU n U , (9)
es, como se establece en (ii) la suma de una solucin general de la
correspondiente ecuacin homognea (una solucin homognea) y una
solucin particular de la ecuacin (9). La ecuacin caracterstica y las races
caractersticas de (4) son tambin llamadas as con respecto a la ecuacin
(9).
Para encontrar una solucin particular de (9), el mtodo de los coeficientes
indeterminados tambin es muy apropiado y til. Por lo que se refiere a la
forma supuesta de una solucin particular, tambin es valida una regla
semejante a la que se cumple para as ecuaciones de primer orden. Sean
1
e y
2
e las races caractersticas de la ecuacin (9) y ) (n P
N
y ) (n Q
N
polinomio
en n de grado N. en este caso, si ) ( ) ( n P n R
N
n
= , se puede encontrar una
solucin particular ) (
*
n U de (9) haciendo
= = =
= = =
=
. ), ( ) (
, ), ( ) (
, ), ( ) (
2 1
2 *
1 2 2 1
*
2 1
*
e e si n Q n n U
e e o e e si n nQ n U
e y e si n Q n U
N
n
N
n
N
n
Si las constantes son nmeros complejos, la clasificacin anterior incluye
el caso en que ) (n R es el producto de un polinomio en n por una suma de
funciones trigonometritas. Si
) ( ) cos ( ) ( n P r n n sen n R
N
n
+ = ,
se puede hacer
) ( ) cos ( ) (
2 1
*
n Q r n k n sen k n U
N
n
+ =
y as sucesivamente, segn la clasificacin dada anteriormente. Pero incluso
si ) (n R es de la forma anterior, a veces es conveniente considerarla como la
parte real de una funcin de valor complejo y usar el teorema de De Moivre.
2.1.3 Ecuaciones en diferencias lineales de tercer orden o superior
28
Definicin 2.7: Los mtodos de solucin ya mencionados tambin se
pueden aplicar a ecuaciones en diferencias lineales de orden tres o mayor.
Puesto que la ecuacin caracterstica de este caso es una ecuacin algebraica
de grado N, el numero de races es N. Segn sea e una raz simple p una de
multiplicidad k, defnase una funcin mediante
n
Ce
n k
k
e n C n C n C C ) ... (
1 2
3 2 1
+ + + + , (10)
cada una de las cuales corresponde a un trmino de la ecuacin (7) y (8),
respectivamente. Si se hace la suma de dichas funciones, se tiene una
funcin que contiene N constantes arbitrarias. Aunque se omite la
demostracin, se puede probar que esa funcin es una solucin general de la
ecuacin homognea, como la del caso de segundo orden.
Si una raz compleja es de multiplicidad k, su conjugada compleja tambin
es raz de multiplicidad k. Por lo tanto, en este caso, basta remplazar los
coeficientes
n
e C
1
de
1 i
n en la ecuacin (10) con la expresin de la formula
) ( ' ) cos( ) ( n sen r C n Cr n U
n n
+ = .
Pero en muchos casos, a diferencia del caso de segundo orden, resulta muy
laborioso resolver las ecuaciones caractersticas.
Para encontrar una solucin particular, si el trmino no homogneo de la
ecuacin es ) (n P
m
n
y si coincide con una raz e de multiplicidad k, se
puede hacer
) ( ) (
*
n Q n n U
m
k n
= .
2.2 Sistemas lineales de ecuaciones en diferencias
Definicin 2.8: Las ecuaciones en diferencias tratadas hasta ahora consisten
en una sola ecuacin en una funcin incgnita. Un conjunto de N ecuaciones
en diferencias de N funciones incgnitas se llama sistema de ecuaciones en
diferencias simultneas. Para resolver el sistema se trata de eliminar
) 1 ( N funciones de las n ecuaciones para formar una ecuacin en
diferencias en una sola funcin.
Ejemplo:
= +
= +
. 0 ) ( 4 ) 1 (
, 0 ) ( ) 1 (
n U n V
n V n U
(1)
La primera ecuacin se escribe como 0 ) 1 ( ) 2 ( = + + n V n U . De esta y la
segunda ecuacin, se tiene
0 ) ( 4 ) 2 ( = + n U n U .
Si luego se remplaza, por ejemplo, la segunda ecuacin del sistema con la
ecuacin en diferencias de ) (n U anterior, se obtiene el sistema
29
= +
= +
. 0 ) ( 4 ) 2 (
, 0 ) ( ) 1 (
n U n U
n V n U
(2)
Siguiendo el mismo procedimiento, de este nuevo sistema se obtiene la
segunda ecuacin del sistema original. Por lo tanto, se puede elegir
cualquiera de los dos sistemas.
La segunda ecuacin del nuevo sistema es la ecuacin en diferencias lineal
de segundo orden cuyas races caractersticas son 2 . Por consiguiente, una
solucin general de la ecuacin es
n n
C C n U ) 2 ( 2 ) (
2 1
+ = .
En seguida se obtiene
1
2
1
1
) 2 ( 2 ) 1 ( ) (
+ +
+ = + =
n n
C C n U n V
de la primera ecuacin. Dados los valores iniciales ) 0 ( U y ) 0 ( V , se tiene
2 1
) 0 ( C C U + = , y
2 1
2 2 ) 0 ( C C V = . De estas ecuaciones, se determinan
1
C y
2
C . Adems, dado que ) 1 ( U se determina por ) 0 ( V , la condicin de
unicidad para ) (n U se cumple. Por consiguiente ) (n V resulta determinar de
manera nica de la primera ecuacin del sistema. De esta manera dicha
pareja de funcione es una solucin general del sistema.
2.2.1 Sistemas de ecuaciones lineales homogneos
Primeramente se considera una solucin general del sistema homogneo
) ( ) ( ) 1 (
) ( ) ( ) 1 (
n dV n cU n V
n bV n aU n U
+ = +
+ = +
(3).
Supngase que
0 , ' ) ( , ) ( = = e e C n V Ce n U
n n
. (4)
Donde al menos una de las constantes C y ' C no es 0, sustituyendo estas
expresiones en (3) y dividiendo por
n
e , se tiene
0 ' ) (
0 ' ) (
= +
=
c d e cC
bC C a e
(5)
Eliminando C y ' C de este sistema, se obtiene
0 ) )( ( ) )( ( = c b d e a e . (6)
Recprocamente, si se satisface (6) entonces (5) tiene una solucin tal que
no ambas C y ' C son cero. De este hecho se ve que si e satisface (6),
entonces la pareja ) (n U , ) (n V es una solucin del sistema (3). La ecuacin
(6) se llama ecuacin caracterstica y sus races se llaman races
caractersticas. Si se escribe la ecuacin (6) en forma de determinante
0 =
d e c
b a e
, (7)
30
la correspondencia entre (6) y el sistema (3) se ve ms fcilmente. Adems,
obsrvese que si el termino constante ) ( bc ad de (6) es 0, existe solamente
una solucin no nula.
Una vez encontradas las races caractersticas, la razn de los coeficientes
C y ' C en (4) se determina de (5) utilizando la primera o la segunda
ecuacin. Se supone que b y c no son ambas cero. Si las races
caractersticas
1
e y
2
e no son iguales, usando los coeficientes determinados
en esta forma, se puede hacer
n n n n
e C e C V e C e C n U
2
'
2 1
'
1 2 2 1 1
() , ) ( + = + = (8)
Si las races caractersticas
1
e y
2
e son iguales, en este caso la ecuacin
tiene races mltiples, se hace
n n
e n C C n V e n C C n U
1
'
2
'
1 1 2 1
) ( ) ( , ) ( ) ( + = + = . (9)
Sin embargo, se debe hacer notar que C y ' C no son necesariamente
proporcionales. Para determinar los coeficientes se debe sustituir (9) en el
sistema (3).
Ejemplo:
). ( ) ( ) 1 (
), ( ) ( ) 1 (
n V n U n V
n V n U n U
+ = +
+ = +
Las races de la ecuacin caracterstica , 0 1 ) 1 (
2
= + e son i e + =1
1
y
i e =1
2
, de (5),
1
'
1
iC C = y
2
'
2
iC C = . Por lo tanto
. ) 1 ( ) 1 ( ) (
, ) 1 ( ) 1 ( ) (
2 1
2 1
n n
n n
i iC i iC n V
i C i C n U
+ =
+ + =
para expresar estas funciones en forma real se escribe
)) 45 ( ) 45 (cos( 2 1
o o
isen i = ,
Entonces
{ }
{ }. ) 45 ( ) 45 cos( ' ) 2 ( ) (
) 45 ( ' ) 45 cos( ) 2 ( ) (
n Csen n C n V
n sen C n C n U
o o n
o o n
=
+ =
Este resultado se puede obtener tambin como sigue. Puesto que
) (n U satisface la ecuacin de segundo orden cuya ecuacin caracterstica es
igual a la dada anteriormente,
31
{ } ) 45 ( ' ) 45 cos( ) 2 ( ) ( n sen C n C n U
o o n
+ = ,
de la primera ecuacin del sistema,
{ }
{ }, ) 45 ( ' ) 45 cos( ) 2 (
) 45 45 ( ' ) 45 45 cos( ) 2 (
) ( ) 1 ( ) (
1
n sen C n C
n sen C n C
n U n U n V
o o n
o o n
+
+ + + =
+ =
+
2.2.2 Sistemas no homogneos
Definicin 2.9: Si un sistema lineal es no homogneo, entonces una
solucin general del sistema es la suma de una solucin particular del
mismo y una solucin general del correspondiente sistema homogneo. Para
encontrar una solucin particular, se puede emplear el mtodo de los
coeficientes indeterminados.
Ejemplo:
. 0 ) ( 4 ) 1 (
, ) ( ) 1 (
= +
= +
n U n V
n n V n U
Dado que las races caractersticas son distintas de 1. se puede hacer
. ) (
, ) (
'
2
'
1
2 1
n k k n V
n k k n U
+ =
+ =
sustituyendo estas expresiones en el sistema se tiene
. 0 4 4 ) 1 (
, ) 1 (
2 1
'
2
'
1
'
2
'
1 2 1
= + +
= + +
n k k n k k
n n k k n k k
comparando los coeficientes de las potencias de n , resulta
. 0 4 . 1
, 0 4 , 0
2
'
2
'
2 2
1
'
2
'
1
'
1 2 1
= =
= + = +
k k k k
k k k k k k
cuya solucin es
3
4
,
9
8
,
3
1
,
9
5
'
2
'
1 2 1
= = = = k k k k .
Por lo tanto una solucion general es
.
3
4
9
8
) 2 ( 2 ) (
,
3
1
9
5
) 2 ( 2 ) (
1
2
1
1
2 1
n C C n V
n C C n U
n n
n n
+ =
+ =
+
32
3. ECUACIONES EN DIFERENCIAS Y ECUACIONES
DIFERENCIALES
3.1 Analogas entre el clculo de diferencias y el clculo diferencial
Al definir la derivada de una funcin como limite de cociente de diferencias
se deducen interesantes analogas entre el clculo de diferencias finitas y el
clculo diferencial.
Definicin 3.1: Dada una funcin y, se denomina derivada de dicha funcin
a una nueva funcin Dy, cuyo valor en x es
h
x y
h
x y h x y
x Dy
x x
) ( ) ( ) (
) (
lim lim
0 0
=
+
=
D es el operador de diferenciacin que aplicado a una funcin da lugar a la
derivada de dicha funcin. El valor de
| |
h
x y ) (
es la pendiente de la recta
que une los puntos de la representacin grafica de y correspondientes a unas
abcisas x y h x + ; el valor de ) (x Dy es la pendiente de la tangente
geomtrica en x . Por ejemplo, dado que
2 2
2 h xh x + =
Entonces tenemos
x h x
h
h xh
Dx
h h
2 ) 2 (
2
lim lim
0
2
0
2
= + =
|
|
.
|
\
| +
=
las sucesivas diferenciaciones de una funcin se indican con potencias
sucesivas del operador D; as, la segunda derivada de y se representa por
y D
2
y proviene de ) (Dy D , la tercera derivada se indica por y D
3
, etc.
Consideremos ahora la operacin inversa a la diferenciacin. Si una funcin
Y es tal que y DY = , se dice que Y es la funcin primitiva de y . As el
operador inverso de diferenciacin es
1
D y escribiremos y D Y
1
= o
\
|
9. si C y y y Y + = =
1
, siendo C una
funcin peridica.
1.
h
x y
x Dy
h
) (
) (
lim
0
2. K , 2 , 1 ) (
1
= =
n y D D y D
n n
3. CDy Cy D = ) (
4.
2 2 1 1 2 2 1 1
) ( Dy C Dy C y C y C D + = +
5. Si y es un polinomio de grado n,
y D
n
es constante y las derivadas de
orden superior son nulas.
6.
1
=
n n
nx Dx
7. vDu uDv v u D + = ) (
8.
2
v
uDv vDu
v
u
D
= |
.
|
\
|
9. si C Y y y DY + = =
, siendo C
una funcin peridica.
Tabla 3.1 Resultados del clculo diferencial y analogas
Correspondientes al clculo de diferencias.
3.2 La ecuacin en diferencias como aproximacin a la ecuacin
diferencial.
Dada la estrecha relacin existente entre el operador en diferencias finitas
y el operador derivada D no es extrao encontrar muchas conexiones entre
las ecuaciones en diferencias y las ecuaciones diferenciales. Consideraremos
la ecuacin diferencial lineal de primer orden con coeficientes constantes.
34
Sea y una funcin definida para valores reales de x en un intervalo
b x a , tal que satisfaga una ecuacin diferencial. Entonces,
b x a B x Ay
dx
x dy
x Dy + = = ) (
) (
) ( (1)
siendo Ay B constantes arbitrarias y adems 0 A . Supongamos que el
valor de y para a x = sea dado, es decir, que debemos hallar una funcin
y que satisfaga a la ecuacin diferencial con condicin inicial
0
) ( y a y =
donde
0
y es una constante dada.
Para lograr una aproximacin a esta ecuacin diferencial a travs de una
ecuacin en diferencias ser preciso un conjunto discreto de valores de x
para los que pueda definirse la ecuacin en diferencias. Por lo tanto
escogemos un numero entero positivo cualquiera n y dividimos el intervalo
| | b a, en n partes iguales, cada una de longitud
n
a b
h
=
siendo los puntos de subdivisin
b nh a x h a x h a x a x
n
= + = + = + = = K , 2 , ,
2 1 0
Simplifiquemos la notacin haciendo
) ( ) (
0
kh x y x y y
k k
+ = =
adems
h
y
x Dy
k
h
k
=
0
lim ) (
parece bastante razonable sustituir la ecuacin diferencial por la ecuacin en
diferencias
1 , , 2 , 1 , 0 = + =
n k B Ay
h
y
k
k
K (2)
o bien
1 , 2 , 1 , 0 ) 1 (
1
= + + =
+
n k Bh y Ah y
k k
K
La condicin inicial indica que el valor
0
y es dado.
35
La ecuacin en diferencias es una ecuacin lineal de primer orden con
coeficientes constantes que puede resolverse fcilmente aplicando la
funcin de valor de equilibrio entonces tenemos
A
B
Ah
Bh
y =
+
=
) 1 ( 1
*
y por tanto
1 , 2 , 1 , 0 ) 1 (
0
= |
.
|
\
|
+ = n k
A
B
A
B
y Ah y
k
k
K
Si partiendo de la ecuacin en diferencias (2) hacemos que
x kh x x y h
k
+ =
0
0
es evidente que dicha ecuacin se convierte en la ecuacin diferencial (1).
Ser por tanto de inters averiguar si la solucin de la ecuacin en
diferencias tiende a una solucin de la ecuacin diferencial al imponer los
mismos lmites. Demostracin,
( )
( )
( )
+ |
.
|
\
|
+ =
+ |
.
|
\
|
+ =
|
.
|
\
|
+ + =
=
A
B
Ah
A
B
y
A
B
Ah
A
B
y
A
B
A
B
y Ah
y x y
h
x x
k
k
k
0
1 lim
1 lim
1 lim
lim ) (
0
0
0
escribiendo e Ah = obtenemos
( )
) (
1
0
1 lim
x x A
e e
+
definimos a e como el siguiente limite
e
e
e e
1
0
) 1 ( lim + =
al pasar el limite, haciendo que 0 h y siendo A constante, el producto
e Ah = tendera tambin a 0. Por tanto,
( )
) (
) (
1
0
0
1 lim
x x A
x x A
e e e
+
36
entonces
A
B
e
A
B
y y x Y
x x A
k
|
.
|
\
|
+ = =
) (
0
0
lim ) (
Puede demostrarse que esta funcin y es, efectivamente, una solucin de la
ecuacin diferencial (1) calculando su derivada y probando que la igualdad
(1) se satisface para todo valor de x comprendido en el intervalo b x a
(Verificar).
En particular para 0 = B tenemos la ecuacin diferencial
) (
) (
) ( x Ay
dx
x dy
x Dy = = (3)
es decir, que el tipo instantneo de cambio de la funcin y respecto a x es
proporcional al valor de la funcin y . Entonces
1 , 2 , 1 , 0 ) 1 (
1
= + =
+
n k y Ah y
k k
K
de solucin
k
k
Ah y y ) 1 (
0
+ =
bajo las condiciones restrictivas ya descritas la solucin correspondiente a la
ecuacin diferencial (3) viene dada por la funcin exponencial
) (
0
0
) (
x x A
e y x y
=
siendo
0
y el valor dado para y cuando
0
x x = . Si y toma solamente
valores positivos se dice que y sigue un crecimiento exponencial si 0 > A o
que sigue un decrecimiento exponencial cuando A<0 cuando x crece a
partir del valor inicial
0
x x = . El crecimiento exponencial viene a ser, por
tanto, el anlogo en forma continua del crecimiento en progresin
geomtrica, que es discreto.
3.3 Un modelo probabilstico sobre el aprendizaje
Supongamos que deseamos estudiar como varia el comportamiento de un
sujeto sometido a unas condiciones dadas de experimentacin. Imaginemos
para ello una serie de sucesos que empiezan con la percepcin de un
estimulo, al que sigue la realizacin de una respuesta y acaban con la
aparicin de un hecho circunstancial externo.
37
El comportamiento de un individuo se mide mediante la probabilidad p de
que de una respuesta determinada, durante un intervalo de tiempo
especificado, una vez que la secuencia de experimentos se ha iniciado. El
lneas generales, p , indica el nivel de disposicin del sujeto para dar una
respuesta; su valor crece y decrece, segn que los factores circunstanciales
favorezcan o se opongan a la repeticin de aquella.
Si imaginamos un experimento en el que se somete repetidamente a un
sujeto a esta serie de sucesos, podemos dividir dicho experimento en varias
etapas, siendo cada una de las cuales una prueba en la que el sujeto recorre
la serie completa. La actuacin del sujeto puede considerarse como una
funcin del nmero de la prueba n lo que nos permite definir
n
p como la
probabilidad de obtener una respuesta dada la prueba nmero n . El nmero
0
p se considerara como el valor inicial que indica la disposicin del sujeto
a emitir la respuesta esperada en el momento en que se inicia por primera
vez el experimento. A partir de
n
p se define la funcin de p para el
dominio de valores de K , 2 , 1 , 0 = n . Por ser p una funcin de probabilidad
tendr que ser
K , 2 , 1 , 0 1 0 = n p
n
siendo los valores extremos 1 y 0 los de certeza absoluta con respecto a que
se de o no la respuesta.
Suponemos, en primer lugar, que
1 + n
p depende solamente de
n
p y no de
los valores anteriores a
n
p . En otras palabras, la actuacin del sujeto en la
prueba 1 + n , que se considera dependiente de lo acontecido en la prueba n ,
lo que no es sino decir que el modelo que acabamos de exponer es
markoviano o que goza de la propiedad de Markov. Siguiendo a Bush y
Mosteller, supondremos para simplificar que la dependencia entre
1 + n
p y
n
p es lineal, es decir, que al representar grficamente
1 + n
p a partir de
n
p
resultara una lnea recta. La forma explicita de esta relacin lineal ser la
ecuacin.
K , 2 , 1 , 0
1
= + =
+
n mp a p
n n
(4)
siendo m el coeficiente angular o pendiente de la recta esto es la cantidad
en que varia
1 + n
p por unidad de incremento
n
p y la ordenada en el origen o
sea el valor de
1 + n
p cuando 0 =
n
p .
Esta relacin lineal podemos escribirla en forma denominada de prdidas y
ganancias. Definimos un parmetro b as
b a m =1
38
con ello, la relacin (4) puede escribirse de la forma
K , 2 , 1 , 0 ) 1 (
1
= + =
+
n bp p a p p
n n n n
. (5)
Si la probabilidad de que el sujeto actu en forma determinada en la prueba
numero n es
n
p , con
n
p 1 se indicara e incremento mximo posible que
puede tener lugar partiendo de
n
p y
n
p se designara el decremento
mximo posible que puede sufrir
n
p al efectuarse la prueba 1 + n . Ello es
consecuencia de que el valor de
n
p ha de quedar comprendido ente 0 y 1.
La ecuacin (5) puede interpretarse diciendo que la variacin en la
probabilidad de obtener respuesta
n n n
p p p =
+1
es proporcional al
incremento o ganancia mxima posible y al decremento o perdida mxima
posible. Las constantes de proporcionalidad son, respectivamente, a y b ,
midindose, por tanto, con el parmetro a aquellas circunstancias que
incitan a que se de la respuesta (premio) y con b aquellas otras que
desanimas al sujeto para que no de la respuesta (castigo).
Los parmetros a y b presentan unas limitaciones originadas por el hecho
de que
n
p y
1 + n
p quedan comprometidos entre 0 y 1. Cuando 0 =
n
p , se
tendr que a p
n
=
+1
, por lo que 1 0 a mientras que si 1 =
n
p , ser
b p
n
=
+
1
1
, por lo que ) 1 ( b tendr que estar comprometido entre 0 y 1,
lo que indica que 1 0 b , estas son las condiciones necesarias que han de
satisfacer a y b para que
1 + n
p quede comprometido entre 0 y 1.
Cuando 0 = a , la ecuacin describe una situacin en la que no se da
compensacin o premio alguno cuando se produce la respuesta, mientras
que cuando 0 = b significa que no se somete al sujeto a ningn castigo si da
la respuesta. Si b a = , la ecuacin indica que la compensacin y el castigo
son iguales.
Bush y Mosteller el cambio progresivo que se da en la probabilidad de
obtener una respuesta en un experimento como el del corredor de Graham
Gagne, o la caja de Skinner, en los que se hace aparecer las mismas
circunstancias cada vez que se obtiene la respuesta esperada
consideremos:
Ejemplo:
Si 4 . 0 = a y 1 . 0 = b , la ecuacin (5) se convierte en
n n n n
p p p p 1 . 0 ) 1 ( 4 . 0
1
+ =
+
o sea
39
K , 2 , 1 , 0 4 . 0 5 . 0
1
= + =
+
n p p
n n
(6)
Suponiendo que el valor inicial 2 . 0
0
= p , podemos calcular sucesivamente
5 . 0 4 . 0 ) 2 . 0 ( 5 . 0
1
= + = p , K , 65 . 0 4 . 0 ) 5 . 0 ( 5 . 0
2
= + = p . Para resolver la
ecuacin (6) y obtener con ello
n
p para un valor cualquiera n , utilizaremos
el valor de equilibrio con 5 . 0 = A y 4 . 0 = B , con lo que
8 . 0
5 . 0 1
4 . 0
1
*
=
=
A
B
p
por lo que obtenemos con valor inicial 2 . 0
0
= p ,
8 . 0 ) 8 . 0 2 . 0 ( ) 5 . 0 ( + =
n
n
p
o bien
K , 2 , 1 , 0 8 . 0 ) 5 . 0 ( 6 . 0 = + = n p
n
n
.
Esta solucin expone la forma en que la probabilidad de obtener la respuesta
varia a medida que se suceden las pruebas. Dado que 1 0 < < a y
*
0
p p < , la
sucesin { }
n
p ser montona creciente y de limite 8 . 0
*
= p . Por tanto, en
pruebas repetidas en las que la proporcin entre compensaciones y castigos
se establezca en 4 a 1, las probabilidades
n
p y
n
p 1 , de que se de y no se
de respuesta, tienden a los valores limites
*
p y
*
1 p , que guardan entre si
la misma proporcin de 4 a 1.
Volviendo al caso general, observamos que la ecuacin (5) puede escribirse
de la forma explicita siguiente
K , 2 , 1 , 0 ) 1 (
1
= + =
+
n a p b a p
n n
(7)
remplazando A b a = 1 y B a = , tenemos
b a
a
A
B
p
+
=
=
1
*
siempre que a y b no sean ambos nulos a la vez. Se obtiene, por tanto, la
solucin
K , 2 , 1 , 0 0
0 ) 1 (
0
0
= = + =
+
+
+ |
.
|
\
|
+
=
n con b a si p p
b a si
b a
a
b a
a
p b a p
n
n
n
(8)
Teniendo en cuenta las limitaciones 1 0 a y 1 0 b se observa que la
constante b a A =1 queda comprendida entre -1 y 1, valores extremos
que adopta A solamente en el caso en que a y b sean ambos iguales a 1 o
a 0. En el caso en que 1 = = b a , ser 1 = A , con lo que { }
n
p ser una
sucesin oscilante finita entre los valores
0
p y
0
1 p . Pero en los dems
40
casos la sucesin { }
n
p ser convergente, de limite
0
p cuando 0 = = b a , o
de limite
*
p para los dems valores a y b . Si es 1 0 < + < b a , resultara
1 0 < < A y la sucesin { }
n
p ser montona decreciente de limite
*
p si
*
0
p p > , montona creciente de limite
*
p si
*
0
p p < y constante igual a
*
p si
*
0
p p = . Si 2 1 < + < b a , tendremos que 0 1 < < A , y { }
n
p ser una
sucesin de limite
*
p . El caso particular en que 0 = + b a da lugar a una
sucesin constante en la que cada elemento es
*
p .
Consideremos dos casos mas, primero 0 = a , se supone que no se ofrece
compensacin alguna cuando se produce respuesta. La ecuacin en
diferencias (5) se convierte en
n n
p b p ) 1 (
1
=
+
(9)
de solucin
K , 2 , 1 , 0 ) 1 (
1
= =
+
n p b p
n
n
n
que describe el decrecimiento continuo que sufre la probabilidad de obtener
respuesta a partir de la probabilidad de obtener a partir de la probabilidad
inicial
0
p . Representado grficamente p como funcin de n, se obtiene la
curva de extincin experimental.
En el caso dos b a = , se dan por igual compensacin y castigo al sujeto cada
vez que da la respuesta. Dejando aparte los casos extremos en que 0 = = b a
y 1 = = b a , se tiene que cuando b a = la cantidad 0 ) 1 (
n
b a cuando
n y aplicando la relacin (8) se obtiene que
*
p p
n
, que es igual a
0.5. Esto es, a medida que aumenta el nmero de pruebas la respuesta tiende
a aparecer en la mitad de las mismas. El equilibrio entre compensacin y
castigo da lugar, a largo plazo, a un comportamiento indeciso del sujeto,
que, que por igual, da respuesta o se abstiene de darla.
41
4. Conclusiones
El estudio de las ecuaciones en diferencias contribuye a que el estudiante
aborde con menos dificultad los sistemas discretos. Adems estos sistemas
tienen la ventaja de ser modelos mas ajustados a la realidad. Aunque a nivel
Matemtico no se le ha dado gran importancia a el estudio de los sistemas
discretos en otras profesiones como Ingeniera y Economa entre otras se
hace necesario su estudio.
Los sistemas continuos son una alternativa para la solucin de los problemas
prcticos que no tienen respuesta en sistemas discretos, la investigacin que
se puede hacer directamente en ecuaciones en diferencias por ejemplo puede
permitir encontrar nuevas soluciones comnmente representadas en modelos
continuos.
42
Bibliografa:
Takahashi, Takehito.
Ecuaciones en diferencias con aplicaciones.
Mxico: grupo editorial iberoamericana.
Goldberg, Samuel.
Ecuaciones en diferencias finitas. Barcelona:
Marcombo S.A. 1964.
Dennis G, Zill.
Ecuaciones diferenciales con aplicaciones.
Grupo editorial iberoamericana.
Kenneth, Bogart.
Matemtica discretas.
Varona Malumbres, Juan Luis.
Mtodos clsicos de resolucin de ecuaciones en diferencias ordinarias.
Universidad de la Rioja
http://personales.ya.com/casanchi/mat/varona01.htm.
Maha, Ramn
Formulacin y solucin de ecuaciones lineales en diferencias con
coeficientes constantes en el contexto del anlisis.
Junio 1998.
http://www.uam.es/otroscentros/klein/doctras/doctra9804.pdf.
Neuman, Carlos Enrique.
Ecuaciones en Diferencias.
Septiembre 14, 2000.
http://www.intec.ceride.gov.ar/CN/EcDife5.pdf