Ecuaciones diferenciales I

4.1 In trod u cció n

Una ecuación diferencial es la relación existente entre una función desconocida

y sus derivadas de cualquier orden. El término “diferencial” se sustenta en la
particularidad propia de incluir diferenciales en la mencionada relación. A modo
de ejem plo, supongamos que se desea encontrar la función de demanda
marshalliana x = / (P) del bien x, donde P es su precio, si se sabe que su
elasticidad precio puntual es r}fl = -1 ( V P > 0). Por definición,

de esta forma, dx___dP^

x P

La expresión anterior constituye una ecuación diferencial en la medida que in­

cluye de manera explícita los diferenciales de x y P. Al integrar ambos lados se
obtiene

ln(x) = -ln ( P ) + C -> x=^

Ecuaciones diferenciales I 111

Se dice que una ecuación diferencial es ordinaria si en su planteamiento existe

una única variable independiente. Por ejemplo,

— ^ = 2x ~~ = sm (y )
dx

son ecuaciones diferenciales ordinarias ya que la función incógnita, es decir

y = y(x), depende de una variable independiente única. Por otro lado, una ecua­
ción diferencial es parcial (o en derivadas parciales) si involucra a más de una
variable independiente. Por ejemplo, el teorema de Euler (para una función
homogénea de primer orden) postulado como

x M + d l =S
dx dy
es una ecuación diferencial parcial para S = S(x, y). En este texto, se desarrollan
aspectos teóricos y prácticos exclusivamente de ecuaciones diferenciales or­
dinarias (EDO, en adelante).

Ejemplo 4.1 Determinar una curva en el plano xy que contenga al punto (0,3) y
cuya pendiente en el punto (x, y) sea 2 x / y 2.

Dado que la pendiente de la recta tangente de la función en (x, y) es la primera

derivada de y con respecto a x, se tiene que

dy = 2x
dx y 2

Despejando ambos términos e integrando se llega a

j y 2dy = J*2x dx —» j-y3 = x 2 + C

Por dato del problema, y(0) = 3, por lo que C = 9. Así, la curva solicitada viene
dada por

y = (3x2 + 27)1/3 *

En general, puede atribuirse cualquier significado económico a las variables,

independiente y dependiente, en una EDO. En el caso de la elasticidad visto
líneas arriba, las variables son precio y cantidad, respectivamente. Sin embargo,
las EDO han encontrado un amplio espectro de aplicaciones económicas, consi­
112 Apuntes de estudio

derando al tiempo como variable independiente. En este sentido, resultan ser

herramientas fundamentales en el análisis dinámico ya que describen el proceso
de cambio (evolutivo o dinámico) de la función incógnita. Por ello, la solución de
este tipo de ecuaciones implica una trayectoria (esto es, una senda temporal o
cronológica), y = y(t), denotada, indistintamente, como yr

Ejemplo 4.2 Considere un flujo de inversión, l t = 100c-0 05'

En un mundo sin depreciación, donde la inversión / equivale al cambio en el

stock de capital en t, se tiene

K, = -2000e~°°5' + C

donde C es una constante arbitraria. La expresión resultante de es una solu­

ción general de su EDO ya que C es indeterminada. Con el fin de encontrar una
trayectoria exacta de Kf es necesario conocer el valor del stock de capital en
algún m omento determ inado. Por ejem plo, si en t = 0, K es 1000, entonces
C = 3000. Estos aspectos serán desarrollados posteriormente. ♦

Ejemplo 4.3 (Modelo de Harrod - Domar) Se sabe que un cambio en los niveles
de inversión afecta tanto a la demanda agregada como a la capacidad productiva
de una economía. El modelo de Harrod - Domar determina la trayectoria en la que
puede crecer una economía, manteniendo la plena utilización de su capacidad
productiva. Si la propensión marginal a ahorrar, s, y la ratio capital producto, k,
son constantes, procederemos a determinar la función de inversión de esta
economía.

En primer lugar, el cambio en la demanda agregada es igual al cambio de la

inversión por un multiplicador l/s, de modo que

dY = l d l
dt s dt
donde Y denota la demanda agregada e /, la inversión. Por su parte, el cambio de
Q, la capacidad productiva del país, es igual a las variaciones del stock de capital
por la inversa de la ratio k. Asimismo, la inversión es igual a los cambios en el
stock de capital,

dQ _ 1 dK dQ =l¡
dt k dt dt k
Ecuaciones diferenciales I 113

En equilibrio, dQ = dY, se tiene

ld ¿ =l ¡ áL - ± d t = 0
s dt k i k

expresión que al ser integrada arroja como resultado

/, C é " kv

donde C es una constante arbitraria (el nivel de inversión cuando t = 0). La

trayectoria obtenida ilustra la predicción esencial del modelo de Harrod - Do­
mar: la inversión aumenta a una tasa constante determinada por s/k, la razón de
ahorros divido por la relación capital producto1. ♦

Una ecuación diferencial implica el uso de derivadas. El orden de una ecuación

diferencial viene dado por el orden de la mayor derivada explícita en esta rela­
ción. Por ejemplo,

dy
(a) (c) — 5-+ cos(r)- e'
dt- dt

de estos casos, (a) es una EDO de primer orden, (b) es una EDO de tercer orden
y (c) es de segundo orden.

Por otro lado, el grado de una ecuación diferencial es la potencia a la cual está
elevada el diferencial de mayor orden. En los casos mostrados, (a) es una EDO
cuadrática (de grado dos) mientras que (b) y (c) son ecuaciones diferenciales
lineales (de grado uno).

Nota. Por conveniencia, en las secciones y capítulos siguientes denotaremos a

la derivada de v con respecto al tiempo como v'. En general,

jt_ d y , „ _ d 2y m_ d i y ínl d" y

dt ' d t1 d t3 ' dt"

1 Para un desarrollo más extenso, así como un contraste em pírico sobre el modelo de
Harrod - Domar, consúltese Jorge Fernándcz-Baca y Janice Seinfeld, Capital humano,
instituciones y crecimiento, CIUP, Lima, 1995, pp. 36 - 42.
114 Apuntes de estudio

En adelante, se discutirán y desarrollarán ciertos criterios analíticos en la resolu­

ción de EDO lineales, es decir, EDO de primer grado. Adicionalmente, se analiza­
rán enfoques alternativos en la resolución de EDO no lineales (de grado supe­
rior a uno). ♦

La discusión anterior puede ser sintetizada en la siguiente definición:

Definición 4.1 Una ecuación diferencial ordinaria (EDO) de orden n satisface

F ( t , y ,y ',y ', y " ,. .. ,y {n)) = 0

donde F denota alguna forma funcional. La ecuación diferencial es ordinaria

en tanto su solución es una función que depende únicamente del tiempo (en
un intervalo dado),
y = X 0 = y,

Alternativamente, una EDO normal de orden n puede expresarse como

y in) = f ( í , y , y ’, y ' , y m, . . . , y {n-'))

Se dice que una EDO es autónoma si no depende explícitamente del tiempo.
Formalmente,
y(n>=

Recordar. La solución de una EDO es una función que la satisface. Así, y = cp es

solución de la EDO si y solo si

F(t, tp, (p', tp", (p,",...,tp<ll)) = 0

4.2 E cu a cio n es d iferen cia les o rd in aria s lin ea les d e prim er orden

Según la Definición 4.1, una EDO de primer orden es toda relación

y' = f ( t , y ) ...(i)

En particular, una EDO lineal de primer orden puede expresarse como

y ' + a ( t) y = b(t) ...(2)

Ecuaciones diferenciales 1 15

La solución de (2) resulta ser una función y t que satisface la EDO. La manera de
enfrentar este problema se basa en la observación de que, si definimos a u = u(t)
= u¡ según

u
luego,
, _ du
= e• — í a (t)dt = u,a(t)
dt d tJ

Así mismo,
d(uy) ,
------- = u y + u y = uy + uay ...(3)
dt

Si (2) es multiplicado po r« (, utilizando (3) se tiene que

uy' + uay = ub

d(uy) _
...(4)

Al integrar ambos lados de (4), se obtiene

v, =u l [jubdt + c ] = e ^ " ln“, \ ^ e ^ l,,‘“b (t)d ty C \^ ...(5)

donde C es una constante de integración arbitraria.

Este procedimiento es resumido en los tres pasos mostrados en el Cuadro No. 4.1

Cuadro No. 4.1

Resolución de una EDO lineal de primer orden

\a{ndt
Paso 1. Calcular,// = e comúnmente denominado factor de integración

Paso 2. Multiplicar ambos miembros de (2) por u y expresar el resultado como

d(uy) / dt = ub

Paso 3 Integrar ambos lados de la ecuación obtenida en el Paso 2 e incluir una

constante de integración
116 Apuntes de estudio

Ejemplo 4.4 Resolver y '- 4 t y = t

Dado que a(t) = - 4 í , el factor de integración viene dado por ut = e ^ i'><" = e 2' .

Al multiplicar la ecuación dada por u, se obtiene d(e~2‘2y )l dt =te~2'2 e inte­

grando ambos miembros, se llega a

y, =Ce2'2 - j

donde C es la constante de integración. ♦

La ecuación presentada en (2) representa el caso más general de una EDO lineal
de primer orden. Cuando a(t) es una función, como el Ejemplo 4.3, se dice que la
EDO tiene coeficiente variable; por su parte si b{t) es una función, la EDO tiene
término variable o móvil. Cuando ambos, a{t) y b{t), son constantes, la EDO
presenta términos constantes.

La mayor parte de las aplicaciones que las EDO de primer orden han encontrado
en economía son, en forma, más sencillas que las presentadas hasta el momento,
aunque muchos modelos y relaciones de comportamiento parten de plantea­
mientos similares a (2). El Cuadro No. 4.2 resume los casos más comunes.

Cuadro No. 4.2

Planteamientos y soluciones de EDO lineales de primer orden

Caso Planteamiento Factor de integración Solución

1. Caso general y ’ + a (t)y = b(t) uí = e l ^ Ver (5)
Casos homogéneos

2. h(t) = 0 y ' + a (t)y = 0 u , = e ^ ,U‘ y t = Ce

a(t) = a€ R
3. y ' + ay = 0 y, = C e “'
*(0 = 0
Casos no homogéneos

4. a(l) = a e R y ' + ay = b(t) u , = e ^ ' = e “- y ,= e - “ [\e°'b (t)dt + c \

a(» = « e R
5. y ’ + ay = b u, = J '" " = e “‘ y, = C e + b /a
b(t) = h e R
6. a(0 = 0 y = b(t) > y, = \ b ( t) d t
Ecuaciones diferenciales I 117

Evidentemente, los casos 2 - 6 del Cuadro No. 4.2 pueden concebirse como simpli­
ficaciones del caso general. Una observación sobre la información contenida en el
cuadro, se centra en la clasificación homogéneo y no homogéneo. Como puede
apreciarse, una EDO es homogénea si b = 0 y no homogénea de otro modo. La
razón de esta clasificación radica en que la solución de una ecuación diferencial
puede ser expresada como la suma de dos soluciones, la solución homogénea o
complementaria y la solución no homogénea o particular. Por ejemplo, en el caso
5 del cuadro, la solución complementaria viene dada por Ce~“' mientras que la
particular corresponde a bla. Si bien es cierto que esta observación no es necesa­
ria en la resolución de EDO de primer orden (ya que el procedimiento del Cuadro
No. 4.1 es suficiente), el principio es útil para el desarrollo de EDO de orden
superior, como se discutirá detalladamente en la sección 4.3.

Hasta el momento se han determinado soluciones que no son únicas debido a la

presencia de una constante arbitraria. En los ejemplos desarrollados, de hecho,
pueden existir infinitas soluciones toda vez que C pueda adoptar cualquier valor
real. Este tipo de solución se llama solución general de una EDO. Uno de los
problemas centrales en la teoría de ecuaciones diferenciales es responder bajo
qué circunstancias existe o no una solución a la que un punto (tQ, y0) e R2
pertenezca. Al determinar dicha solución, la función resultante es única o solu­
ción específica.

El problema de hallar una solución para (1) que pase por el punto (?0, yQ) es
conocido como el problema del valor inicial donde y(ro) = y0 es la condición
inicial. Dada la condición inicial, la constante de integración puede ser determi­
nada, enfatizando el hallazgo de una solución única.

Ejemplo 4.5 Resolver la ecuación y ' - y l t = t dado y; = 8.

Al aplicarlos pasos del Cuadro No. 4.1, se tiene que y = t\n {t)+ C. Si utilizamos
la condición y x = 8, se tiene que C = 8, y así se llega a una función única que
satisface la EDO. ♦

4.2.1 Aplicación: crecimiento exponencial

Muchas variables cambian con el tiempo proporcionalmente con el valor que la

misma toma en el presente. Ejemplos de ello son el número de habitantes de un
país y algunos tipos de inversión. Las ecuaciones diferenciales han probado ser
una herramienta útil en la modelación de este tipo de comportamiento. De hecho,
la formulación matemática de este fenómeno es:
 Apuntes de estudio

y' =a y ...(*)

Se dice que una variable tiene crecimiento (decrecimiento) exponencial en el

tiempo si su tasa de crecimiento (o decrecimiento) es proporcional al valor que la
variable toma en el presente, lo que es expresado en (*).

El número a es denominado constante de proporcionalidad o coeficiente

malthusiano2. Si a > 0, y > 0, por lo que y se incrementa con el tiempo; por el
contrario, si a < 0, y decrece en todo momento. Por este motivo, cuando a es
positiva se dice que es una tasa de crecimiento, mientras que cuando es negati­
va se trata de una tasa de decrecimiento. Alternativamente, el crecimiento o
decrecimiento exponencial puede ser entendido en términos de una tasa relati­
va. De (*) se desprende que
y '/ y = a
de modo que se dan cambios exponenciales en una variable y, si su tasa de
cambio relativa es constante.

Con el fin de determinar la trayectoria de y, (*) puede ser reescrita como

y '- a y = 0

Es decir, como una EDO homogénea similar al caso 3 del Cuadro No. 4.1. Así, su
solución viene dada por
v, = y0e“'

donde y0 es el valor de v en t = 0.

4.2.1.1 Crecimiento poblacional

Suponga que el número de habitantes de un país o región sigue un modelo de
crecimiento exponencial. En un período inicial tQ= 0 la población de este país
ascendió a AQ; posteriormente, en t ¡ la población fue de A! (es claro que t ¡ > t(j y
A ¡ > Afl). Entonces, cabe preguntarse cuál será la población de este país en r, y
en qué período la población será igual a A,.

2 En honor al economista “clásico'’ Thom as Malthus, quien sostenía en 1798 que, dado un
crecim iento poblacional geométrico junto con un crecimiento aritmético en los niveles de
producción (de alimentos), los salarios se verían dramáticamente reducidos a un nivel míni­
mo de subsistencia. Una nota interesante al respecto puede ser hallada en Thom as M.
Humphey, “Algebraic Production Functions and Their Uses Before Cobb-Douglas”, Economie
Quarterly, vol. 83, 1, Federal Reserve of Richmond, 1997.
Ecuaciones diferenciales I 19

De acuerdo con los datos mencionados, el número de habitantes del país se

comporta según
y, = V "

con la información del período r se desprende que

= A0ea'

por lo que a, la tasa de crecimiento, está dada por3

a =[ln(A1) - l n ( / \ 0)]/í,

Así en t2 la población ascenderá a y h = A0e |ll,(/'1»-'«.Aj>i/2/<,

Por su parte, se deduce que en t = [ln(A2) - ln(A0) ] / a la población será igual

aAr

Resolvamos el ejercicio anterior con cifras. Según datos de la Organización de

Naciones Unidas, la población mundial en 1990 fue de 5.3 mil millones de habi­
tantes con una tasa de crecimiento de 2% al año. Si consideramos a 1990 como
t = 0, tendremos que la población mundial obedece a

.y, = 5.3e 002'

miles de habitantes. Así, este planteamiento predice que en el año 2000 (í = 10)
la población fue de y|(J = 5.3e002x|°, que es aproximadamente 6.5 miles de millones
de habitantes. Por su parte, la población mundial será de 12.5 miles de millones
de habitantes en t = ln(12.5/5.3) = 43, es decir, en el año 2033.

Por otro lado, si y tiene crecimiento exponencial, el tiempo necesario para que su
valor actual se duplique se denomina tiempo de duplicación. Si y tiene un mode­
lo de decrecimiento exponencial, el tiempo necesario para que su valor se reduz­
ca a la mitad es denominado vida media dey. Formalmente, si T denota el número
de períodos necesarios para que y se duplique, se tiene

2 v0 = y*e*r
por lo que
T _ ln(2)
a

3 Note que el caso general de esta expresión es a = [In( A t ) - ln( Att)] /(/, - r„) .
120 Apuntes de estudio

Es importante notar que T no depende de las condiciones iniciales del problema

(tn e y 0). El lector puede demostrar que si a< 0, la vida media de una variable con
decrecimiento exponencial es T= -ln(2)/aperíodos .

Con los datos anteriores, puede concluirse que en T= ln(2)/0.02 = 35. Esto es,
manteniéndose un ritmo de crecimiento de 2% anual, la población mundial será
10.6 miles de millones en el año 2025 y 21.2 miles de millones hacia el año 2060.

Figura No. 4.1:

Población mundial

Según el modelo de crecimiento exponencial,

tomando como punto de partida una pobla­
ción de 5.3 miles de millones de personas en
1990 (punto A), en el año 2000 el planeta ha
sido habitado p o r 6.5 miles de millones de
individuos (punto B), mientras que en el año
2033 serán cerca de 12 mil millones (punto
C). Por otro lado, el punto D corresponde a
2024, periodo consistente con el tiempo de
2040 duplicación de la población de 1990.

Por último, como se observa en la Figura No. 4.1 al comparar los resultados del
modelo de crecimiento poblacional desarrollado con las proyecciones (estima­
cio n es) de la O N U 45, podem os apreciar que el prim ero provee valores
sistemáticamente mayores que las segundas. De hecho, es difícil imaginar un
conjunto poblacional que mantenga, por períodos prolongados, un crecimiento
exponencial (note que en cortos períodos, el modelo exponencial presenta un
perfil aceptable, en particular hasta el año 2000). Posteriormente, se expondrán
modelos alternativos que subsanen estas deficiencias.

4.2.1.2 Interés compuesto

En el primer capítulo se introdujo la idea básica detrás del funcionamiento del
interés compuesto. Considere cierta cantidad de dinero invertido a una tasa
anual de interés r. Si las ganancias obtenidas por concepto de intereses de esta
inversión se capitalizan cada año, se dice que el interés es compuesto anual­

4 Una aplicación económica inmediata es el cálculo de la tasa de depreciación de algún

activo. En general, si se calcula que la vida media de este activo es igual a T, su tasa de
deprecación es d = ln(2) / T.
5 Las cifras han sido tomadas de http://www.popin.Org/popl998/default.htm #contents.
Ecuaciones diferenciales 1 121

mente; si se realizan dos veces al año, semestralmente; si son cuatro veces al año,
trimestralmente y así sucesivamente. Un caso extremo se centra en la capitaliza­
ción o composición del interés cada día, cada hora, cada segundo o incluso cada
fracción de segundo. El límite de esta secuencia es que el interés sea compuesto
de manera instantánea, fenómeno conocido como capitalización continua.

Si h denota una fracción de tiempo, el interés redituado por una inversión At en

el período comprendido entre t y t +h es

Ai+h - A' = rhA,

Así, si h —> 0, se tiene que

hm ----- - = rA
h^o h
que equivale a
A ' = rA,

De ello se deduce que el capital invertido será, en el período t

A, = A0e rl
Es decir, en la medida en que transcurre el tiempo, el dinero invertido inicialmen­
te, Aq, tiende a crecer en un factor erl, conocido como tasa de interés efectiva.

4.2.2 Convergencia y estabilidad

La solución de una EDO es estable o converge si tiende a algún valor real y

definido conforme la variable dependiente, típicamente t, aumenta de manera
indefinida. El valor real al que tiende la trayectoria es denominado estado esta­
cionario o, dado el carácter temporal que hemos adoptado en el análisis, equili­
brio de largo plazo.

Teóricamente, el estado estacionario resulta ser una situación en la que la varia­

ble dependiente de la solución de la EDO permanece estática. Es decir, el estado
estacionario de y en (2) se da cuando y = 0. Si se aplica esta condición a la
ecuación (2), se tiene que el valor de largo plazo de y viene dado por
y, = b (t)/a (t)

Luego, una primera condición para que v sea convergente es

lim b(t)/ a(t) = y (S e D c R
122 Apuntes de estudio

En otras palabras, si el estado estacionario de y, yss, no es definido (o es infinito),

y será una función divergente (se cumple que y & 0 para todo valor posible de
t). En caso de existir estado estacionario, definido sobre un conjunto Q. de los
números reales, la posible diferencia existente entre el valor inicial de y y su valor
de largo plazo, es decir y0- y¡s, es denominada desequilibrio inicial.

Si bien la condición mencionada en el párrafo anterior es necesaria, no es sufi­

ciente. Supóngase una ecuación similar al caso 3 del Cuadro No. 4.2. Claramente,
el estado estacionario es definido en ysj = 0. La solución a este problema es yt =
Ce-“'. Si a < 0, el límite de y es infinito; por el contrario, si a > 0, su límite es
definido (e igual a cero). Aunque existe un estado estacionario, queda abierta la
posibilidad de que la trayectoria de y nunca lo alcance, o incluso se aleje de él.
Así, se desprende una segunda condición, dada por

lim l/u (t) = 0

donde u es el factor de integración de la EDO, definido en el Cuadro No. 4.1

Teorema 4.1 Se dice que la solución de una EDO lineal de primer orden
converge a su estado estacionario (o, alternativamente, que la EDO es esta­
ble) si
lim y t = y ss € Q, R

Para que la convergencia se dé, debe cumplirse que

lim b (t)l a(t) = y ss e Q c R

lim l/u (t) = 0

Por otro lado, si

lim y t = ±

y es divergente. En algunos casos se dice que y tiene una trayectoria

explosiva.
Ecuaciones diferenciales I 123

Ejem plo 4.6 Con el fin de esclarecer los conceptos vinculados con la conver­
gencia de una trayectoria, considérese los ejemplos presentados en el Cuadro
No. 4.3 cuyas soluciones son graficadas en la Figura No. 4.2

Cuadro No. 4 3
Algunos ejemplos de EDO de primer orden

E c u a ció n Condición lim ^ U/ y„ ¿Converge?

/ - y = s e n (í) y G r/2 ) = G °o e~‘ No existe No

“
a)
Solución: y t = ~ ( e r“05jr •- s e n ( / ) - c o s ( 0 )

y ' + 4 ty = 6 t y ( 0 ) = 15 1.5 e^ 0 1.5 Si

b)
Solución: y, = 1 3 .5 c '2' + 1.5

y ' - 2 y = 6e~' y ( 0) = 0 0 e ' 2f 0 No

c)
Solución: y t = 2 e 2' —2e

y ' + y = c o s (2 í) y (0 ) = 0 .2 w e l- U l] e‘ 0 No definido No
d)
Solución: y , = 0 .2 ( c o s ( 2 í) + 2 s e n (2 í))

Según lo definido, el caso (a) constituye una EDO inestable por dos motivos: la
trayectoria de su solución no tiene un estado estacionario y el límite de la inversa
del factor de integración asociado con esta EDO no es definido. En el Panel (a) de
la Figura No. 4.2 se aprecia claramente que en este caso y tiende al infinito confor­
me transcurre el tiempo, esbozando así una trayectoria explosiva. Por el contrario,
el caso (b) es estable y, consecuentemente, y converge hacia 1.5, su valor de largo
plazo, representado por líneas punteadas en la Figura No. 4.3. Dado que y en t = 0
es 15, el desequilibrio inicial de esta trayectoria es 1 5 -1 .5 = 13.5.

Figura No. 4.2:

Ejemplos de trayectorias
124 Apuntes de estudio

Los casos (c) y (d) ilustran la importancia de que las dos condiciones discutidas
en el Teorema 4.1 deben cumplirse para asegurar la convergencia de y. En primer
lugar, el problema del caso (c) presenta un estado estacionario (igual a cero); no
obstante ello, la trayectoria de y se aleja sistemáticamente del mismo, compor­
tándose de manera explosiva (no se cumple la condición vinculada al límite de 1/
u). Por su parte, en el caso (d) el límite dey no es definido ya que se trata de una
función permanentemente fluctuante. Así, se tiene una función divergente, mas
no explosiva.

4.2.2.1 Aplicación: estabilidad en economía

4.2.2.1.a Equilibrio macroecónomico

C o n sid e re m o s un m odelo s e n cillo de d e te rm in a c ió n de e q u ilib rio
macroeconómico. Para tales efectos expresaremos las variables como desviacio­
nes de su equilibrio de largo plazo. Así, si X es una variable macroeconómica, su
desviación con respeto a su estado estacionario es
X , = X t - x ss

El modelo está compuesto por las siguientes ecuaciones:

Cj = cYt 7, = bYt Y' = a (C , + /, - Yt ) 0 < a,b,c < 1

donde se asume que tanto el consumo C como la inversión / son funciones

lineales del ingreso Y. Asimismo, la tercera ecuación nos informa que el ingreso
cambia a un índice proporcional a la demanda en exceso. Sustituyendo las dos
primeras ecuaciones en la tercera, se tiene:

Y' = a(c + b - 1)K

que es una ecuación diferencial similar a la estudiada en el modelo de crecimien­

to exponencial. Así,

Yt = Cea(c+b~l"

En t = 0, Y = Y0 - Yss por lo que

K=[E0-YJe“^1
’-'"+Yís
De este modo, se aprecia que conforme 1 —» °°, Y —» Ks solo si c + b < 1. Es decir,
se alcanza el equilibrio macroeconómico si la suma de la propensión marginal al
Ecuaciones diferenciales 125

consumo y la propensión marginal a invertir, c y b, respectivamente, es menor a

uno. Asimismo, si se tiene una situación inicial Yg < Ys¡ (Yg > Ys¡), Y crecerá
(decrecerá) exponencialmente bajo la condición de estabilidad, con el fin de
alcanzar su valor de equilibrio. Por el contrario, si c + b > 1, Y decrecerá (crecerá)
de modo que se trate de una trayectoria explosiva. Por último si c + b = 1, Y = Yg.

4.2.2. l.b Demanda por dinero e inflación

Bajo el enfoque de la teoría cuantitativa del dinero, la demanda por dinero es
motivada exclusivamente por fines transaccionales y puede ser expresada como
M f = kPtQ

donde k es una constante positiva, P es nivel de precios de la economía en el

instante t y Q es la producción real que se asume fija. La inflación, es decir, la
tasa de cambio del nivel de precios, es proporcional al exceso de oferta moneta­
ria (constante), AP, sobre la demanda de dinero,
P' = § ( M S - M f ) # >0

La dinámica que siguen los precios puede ser expresada como una ecuación
diferencial; reemplazando la ecuación de demanda por dinero en la ecuación de
inflación, se tiene
P '+ m q p = m 5

con lo que se deduce que

P, =\Pa - P ]e -*<* + P P = M ' / k(¿

Puesto que Ü kQ > 0 , la convergencia del nivel de precios está asegurada.

Asimismo, puede verificarse que la inflación converge a su valor de estado
estacionario, P' = 0.

Supongamos ahora que la demanda por dinero se ve reducida por las expectati­
vas inflacionarias, Pe', de modo que

M f = kPtQ - gP‘ '

Por simplicidad, asumiremos que la inflación esperada es función positiva de la

inflación observada de modo que P 1' = hP'. Bajo estas especificaciones, puede
verificarse que el nivel de precios obedece a

P = [PQ_ p j e -«KW-W»> +P"

126 Apuntes de estudio

En este caso, P —> Pss si ühg < 1. Este resultado es interesante ya que incluso
cuando h > 1, lo que quiere decir que los agentes esperan que la inflación se
acelere, la economía es estable toda vez que t? y/o g tomen valores suficiente­
mente reducidos.

Figura No. 43 :
Nivel de precios e inflación bajo condiciones de convergencia

4.2.2. l.c Equilibrio microeconómico

Consideremos un mercado donde las funciones de demanda y oferta, Q d y Q \
vienen definidas como

Q< = a , - P llP¡
Q\ = a s + PsPt + e~e‘ eos(t)

La función de demanda resulta convencional. Sin embargo, la función de oferta

presenta un término móvil poco usual cuya presencia es justificada en la medida
en que el productor se vea afectado por shocks periódicos. Ejemplos de ello son
los denominados productos estacionales como uniformes escolares o adornos
navideños. Asimismo, el término móvil puede representar un ajuste dinámico en
el proceso de aprendizaje del productor; en el límite, este ajuste es nulo (dado
que 0 es positivo) y, por ende, el productor fija su nivel de producción sin
distorsiones.

Los precios en este mercado se ajustan proporcionalmente con el exceso neto

de demanda,

p ' = m ‘! - Q / ) ¿>0
Ecuaciones diferenciales I 127

Asumiendo k = 1 y reemplazando las funciones de oferta y demanda en la

ecuación de ajuste de precios se llega a la EDO

P' + (P,, + P ,)P = ( a rf - a s ) - e~9' c o s (f) ...(* )

Según lo establecido en el Teorema 4.1, se tiene

( a d - a s ) - e ~ a cos(t) (orrf - or,)

lim — í-----í-------------------- e R y lim e ~ ^ =0
( P d + P J ÜÍ7+K)
por lo que puede afirmarse que el estado estacionario de los precios es definido
y que la interacción de oferta y demandan aseguran la convergencia hacia el
mencionado precio de equilibrio.

La EDO (*) es en forma similar al caso 4 del Cuadro No. 4.2. La solución es
hallada mediante integración por partes del producto del término móvil por el
factor de integración,

(Pd+ Ps- ú ) c o s ( f)- sen(r)~j

c o í(t)]d t
(P,,+Ps - 0)2+1 J
que es multiplicado (junto con una constante de integración) por la inversa del
factor de integración. De esta forma, la trayectoria de los precios viene dada por

a ,ra s ,-et (Pd + - 0 )c o s (/)-s e n (O

Pr = C e~<íi"+p->' +
Pd + p, L ([C+p7- o)~+ i _

donde C es una constante determinada a partir de alguna condición inicial sobre

el mercado. En la Figura No. 4.4 se presenta la trayectoria de los precios y de las
cantidades demandadas y ofertadas, según la situación inicial del mercado. Si
enf = 0el exceso de demanda es positivo (es decir, Q‘'> Q S) los precios tenderán
al alza de modo que la brecha sea cerrada. Por el contrario, los precios son
decrecientes si en t = 0, Q d < Q \ Asimismo, note que la presencia del término
móvil retarda la consecución del equilibrio ya que tanto P como Q'1y Q 1oscilan
alrededor de sus valores de equilibrio hasta que éstos son alcanzados.
128 Apuntes de estudio

El modelo desarrollado tiene la particular característica de presentar un factor de osci­

lación en las decisiones de oferta que puede interpretarse como un componente estacional
o ajustes por aprendizaje del productor. Los Paneles (a) y (b) presentan la trayectoria del
precio y las cantidades a partir de un exceso de oferta inicial. Claramente, las cantidades
ofertadas son fluctuantes (lo que genera el mismo efecto sobre el resto de variables),
hecho que puede ser apreciado en el Panel (c), donde la situación inicial corresponde a
los primeros puntos de la derecha de la Figura. Los Paneles (d), (e) y (f) presentan la
dinámica del modelo ante un exceso de demanda inicial.

4.2.3 Resolución de un caso especial

En esta sección desarrollaremos brevemente una metodología de resolución de

EDO lineales de la forma especificada en el caso 4 del Cuadro No. 4.2, es decir,
nos centraremos en resolver

y ' + ay = b(t) ...(6)

Esta EDO tiene la particularidad de contar con coeficiente constante, a, y un

término móvil, b(t). La lógica de esta sección es la de introducir al lector en el
método de resolución más común de EDO de orden superior. Por tal motivo, los
aspectos formales detrás de los pasos que a continuación se discuten pueden
ser encontrados en las secciones 4.3.1 y 4.3.2.
Ecuaciones diferenciales I 129

Como se mencionó anteriormente, la solución de una EDO puede ser expresada

como la suma de dos soluciones
y, = y c + y p ...o
donde es denominada solución o integral complementaria (u homogénea) e
y es conocida como integral particular (o no homogénea). La idea de denomi­
nar “integral” a las soluciones está vinculada con el intrínseco carácter diferen­
cial de un planteamiento como (6).

La solución complementaria, y c, es aquella función que satisface la versión ho­

mogénea (conocida como forma reducida) de (6)
y ' + ay = 0 ...(8)

Del análisis de las secciones anteriores (ver Cuadro No. 4.2) se sabe que las
funciones del tipo Cen satisfacen (8). Entonces, cabe preguntarse qué función
en particular satisface (8). Para ello, planteamos una solución yc = Cen. Así,
y ' = rCe". Reemplazando este intento de solución en (8), se tiene

H er,(r + a ) = 0
Luego, para que yc sea solución de (8), r = -a . En otras palabras, _v = Ce““', lo que
es totalmente coincidente con la solución respectiva dada en el Cuadro No. 4.3.

Por su parte, la solución particular, y , es aquella función que satisface (6), la

ecuación no homogénea o forma completa de la EDO. Una manera alternativa al
método del Cuadro No. 4.1 de encontrarla es el método de los coeficientes inde­
terminados, que será detallado en la sección 4.3.2. Los lincamientos básicos de
este método sugieren plantear una función 'F(í), de la familia de funciones de
b(t), que cumpla con (6),
y¥ ' + a' ¥=b( t )

dado que Y es similar a b(t) la expresión anterior constituye una igualdad de

polinomios. Por ello, para que 'F sea solución de (8), sus coeficientes deben ser
determinados (de ahí el nombre del método) de acuerdo con la EDO.

Ejemplo 4.7 Resolver y ' + y = 5 t2

La solución complementaria viene dada por yr = Ce~‘. Con el fin de hallar la

solución particular, se plantea un polinomio de segundo orden, de similares
características a b(t) = 511,
130 Apuntes de estudio

y p = c 2t 2 + c,r + c 0 y'p = 2c2t + c¡

Al reemplazar este intento en la EDO,

c2t 2 + ( lc 2 + cl)t + (cl + c 0) = 5 t2

se obtiene una igualdad de polinomios que debe cumplirse con el fin de que y
satisfaga la EDO. Resolviendo la mencionada igualdad,

c2 = 5
2 c2 + c, = 0 —» c, = - 2 c2 = - 1 0

c, + c0 = 0 —> c0 = -c , = 10

Por consiguiente, la solución particular es y = 5 í 2 - lOí +10 y el problema es

resuelto por
y, =yc + yp = C e x+ 5t2- íor + 10

El lector puede comprobar que esta solución satisface la EDO planteada. ♦

Ejemplo 4.8 Resolver y ' + y = 4e4'

La solución complementaria se deduce inmediatamente, yc- Ce'. Por su parte,

como solución particular se plantea
y p = ce4' y ’p = 4ce4'
Al reemplazar este intento en la EDO,
5ce4' = 4e4'
de modo que c = 415 y la solución de la EDO es

y, =He~‘ + j e 4' ♦

Con los ejemplos mostrados un aspecto queda claro: la solución particular siem­
pre es definida y única. No obstante, la solución complementaria requiere de
una condición inicial para ser determinada, dada la presencia de la constante de
integración C. Aspectos generales sobre el método de coeficientes indetermina­
dos y los posibles problemas en su uso serán abordados posteriormente.

Al adaptar las condiciones de estabilidad a ( 8 ), puede notarse que y es conver­

gente si
Ecuaciones diferenciales 1 131

a >0
lim b(t) = ¿te £2 c R

De modo que el estado estacionario en este caso puntual es ys¡ = bla. Note que
este valor es plenamente coincidente con la solución particular cuando b{t) es
constante.

Ejemplo 4.9 (La curva de aprendizaje): La ecuación

Q' = k ( B - Q t ) k ,B > 0
ha sido sumamente utilizada en estudios vinculados con la relación entre la
eficacia con que una persona puede desempeñar una tarea y la cantidad de
experiencia que el individuo posee. Suponga que Q es la experiencia que Ud.
acumula en su trabajo. Según la formulación presentada, Ud. puede acumular un
máximo de experiencia de B (que es la cota superior de Q); asimismo, conforme
va obteniendo mayores conocimientos, el ritmo al que obtiene experiencia se
reduce proporcionalmente con la constante k.

La dinámica de Q puede ser expresada como

Q '+ k Q = kB
La solución complementaria de esta EDO es Q = Ce u mientras que la solución
particular es constante e igual a Q =B. Así, dado un nivel de experiencia inicial,
Q0, la trayectoria de Q responde a

Q, =ÍQ0 - B ] e k' +B
Ya que k es un número positivo y B es una constante (por lo que su límite
cuando t —>°° es B) Qt es una trayectoria convergente. Del mismo modo, dado
que por construcción B > Qg, Q: es una trayectoria creciente (ver Figura No. 4.5).

4.2.3.1 Aplicación: dinámica de la producción per cápita

En un estudio sobre la evolución de la población (mano de obra) de una econo­
mía como factor productivo, el premio nobel Trygve Haavelmo 6 propone las
siguientes ecuaciones:

6 Ver Trygve Haavelmo, A Study in the Theory o f Economic Evolution. North-Holland,

Amsterdam, 1964, citado en: Peter J. Hammond, y Knut Sydsaeter, op. cit., pp. 624.
132 Apuntes de estudio

Figura No. 4.5:

Q La curva de aprendizaje
B
La c u r v a d e a p r e n d iz a je describe la
relación existente entre la eficacia de
un individuo en ejercer alguna labor y

ft la cantidad de conocimiento adquiri­

do. Esta relación es descrita por una
t
EDO donde cierta cantidad varía ne­
gativa y proporcionalmente con la dis­
tancia existente entre su valor presente y una cota superior, B, de modo que
cuanto mayor sea esta diferencia, mayor el ritmo de acumulación de la men­
cionada cantidad.

— =a-/3— a,fi> 0
N X
X = AN " a > 0, a A 1 ...(**)

donde N representa el tamaño poblacional y X denota la producción total de

esta economía. La primera ecuación sostiene que la tasa de crecimiento poblacional
es una función inversa del recíproco del producto medio de esta economía; en
otras palabras, un producto medio elevado (reducido no) puede sostener nive­
les altos de crecimiento poblacional. La segunda ecuación es una función de
producción, siendo a el parámetro que determina la naturaleza de los rendimien­
tos marginales de la mano de obra. Si se define a x = X/N al producto por
poblador (producto medio), luego de (**) se desprende

...(♦)

Si reemplazamos lo hallado en ( ♦ ) en (*), se tiene

(a -l)x

Multiplicando ambos miembros de esta ecuación por (a—l)x se tiene la siguien­

te EDO
Ecuaciones diferenciales 1 133

x '- a ( a - l)x = - / 3 ( a - 1)

De esta forma la trayectoria del producto por poblador viene dada por

Como se discutió esta trayectoria es convergente si 0 < a < 1, es decir, si la

función de producción presenta rendimientos marginales decrecientes, lo que
es intuitivamente lógico7.

Con el propósito de analizar los resultados del modelo de Haavelmo, considere

las siguientes ecuaciones adicionales a (*) y (**)

Asimismo, bajo la condición de convergencia o la presencia de rendimientos

marginales decrecientes, el valor de estado estacionario d e x e s definido, al igual
que el correspondiente a N y X. Los límites de estas variables vienen dados por

lim ^ = lim jf- = lim -f = 0

Como se aprecia en la Figura No. 4.6, inicialmente existe cierta cantidad poblacional
que varía a una tasa aproximadamente constante. Dados los rendimientos decre­
cientes, la producción crece a un ritmo menor al número de habitantes, lo que
conlleva a una permanente caída en el producto medio. Este deterioro desacelera
el crecimiento poblacional, que a su vez desacelera el crecimiento de la producción
hasta llegar a una situación de estado estacionario o de largo plazo. En ella, el
producto medio, x, permanece constante al igual que N y, por ende, X y sus tasas
de variación se toman nulas. Note que el nivel de vida de esta economía, aproxima­
do por la cuantía dada por el producto medio, depende fundamentalmente de dos
parámetros, a y P, asociados estrechamente con el ritmo de crecimiento poblacional.

7 Vease Michael J. Stutzer, “Chaotic Dynamics and Bifurcation in a Macro M odel”, Federal
Reserve Bank of Minneapolis, S ta ff Report No. 55, 1980.
134 Apuntes de estudio

Figura No. 4.6:

El modelo de Haavelmo

Bajo rendimientos decrecientes, la producción de un país crece a un ritmo menor que su

población lo que, a su vez, reduce el producto medio. Esta disminución desacelera el
crecimiento poblacional con repercusiones en la producción y así sucesivamente. La
economía converge a un estado estacionario donde el producto medio permanece cons­
tante.

Estas cifras vienen dadas exógenamente y están asociadas con la capacidad de la

economía de absorber la nueva entrada de mano de obra así como a un nivel
exógeno del resto de factores productivos como, por ejemplo, capital. Así, un
incremento en el stock de capital de esta economía es consistente con una reduc­
ción considerable en [3. Por su parte, políticas de control de natalidad buscan
reducir a. Ambos fenómenos aseguran un mayor nivel de vida futuro.

4.3 Ecuaciones diferenciales ordinarias lineales de segundo

orden

En esta sección se estudia lo concerniente a EDO lineales de segundo orden. La

materia por tratar es fundamental para el entendimiento tanto de EDO lineales de
orden superior, sistemas de EDO y una mayor gama de modelos dinámicos utilizados
en economía. En general, una EDO de segundo orden puede ser expresada como

y '+ a l ( t) y '+ a 0( t)y = b (t) ...( 8 )

donde a r aQy b son funciones continuas. Como se mencionó, y ' denota la

segunda derivada de y con respecto a t. En adelante trataremos la solución
analítica del caso donde los coeficientes de ( 8 ) son constantes. Soluciones
generales de ( 8 ) podrán ser halladas o aproximadas con los métodos desarrolla­
dos en el siguiente capítulo.