NOTAS HISTÓRICAS

El campo de las ecuaciones diferenciales se originó en el estudio del cálculo por Isaac
Newton (1642-1727) y Gottfried Wilhelm Leibniz (1646-1716) en el siglo XVII.

Aunque Newton trabajó relativamente poco con ecuaciones diferenciales en sí, su

desarrollo del cálculo y la dilucidación de los principios básicos de la mecánica
constituyeron la base para sus aplicaciones en el siglo XVIII, más notablemente por
Euler. Newton clasificó las ecuaciones diferenciales de primer orden en las formas
dy/dx = f(x), dy/dx - f{y) y dy/dx = f{x, y). Para la última ecuación desarrolló un
método de solución mediante series infinitas cuando f(x, y) es un polinomio en x y y.
La activa investigación de Newton en la matemática cesó a principios del decenio de
1690, excepto por la solución de "problemas desafiantes" ocasionales y la revisión y
publicación de resultados obtenidos mucho tiempo antes.

Leibniz nació en Leipzig y concluyó su doctorado en filosofía a los 20 años de edad

en la Universidad de Altdorf.

Leibniz llegó a los resultados fundamentales del cálculo de manera independiente,

aunque un poco después que Newton, pero fue el primero en publicarlos, en 1684.
Leibniz estaba muy consciente de la conveniencia de una buena notación matemática
y la notación actual para la derivada (dy/dx) y el símbolo de la integral se deben a él.
Descubrió el método de separación de variables en 1691, la reducción de ecuaciones
homogéneas a ecuaciones separables en el mismo año, y el procedimiento para
resolver ecuaciones lineales de primer orden en 1694.

Los hermanos Jakob (1654-1705) y Johann (1667-1748) Bernoulli, de Basilea,

contribuyeron mucho al desarrollo de métodos para resolver ecuaciones diferenciales
y a la expansión de su gama de aplicaciones.

Con la ayuda del cálculo resolvieron varios problemas de mecánica al formularlos

como ecuaciones diferenciales. Por ejemplo, Jakob Bernoulli resolvió la ecuación
diferencial y´= [a3/(b2 y - a3 )] 1/2 en 1690 y en el mismo artículo fue el primero en usar
el término "integral" en su sentido moderno. En 1694, Johann Bernoulli pudo resolver
la ecuación dy/dx = y/ax.

Daniel Bernoulli (1700-1782), hijo de Johann, emigró en su juventud a San

Petersburgo. Se interesaba principalmente en las ecuaciones diferenciales parciales y
sus aplicaciones. Por ejemplo, es su nombre el que se asocia con la ecuación de
Bernoulli en mecánica de fluidos. También fue el primero en encontrar las funciones
que un siglo más tarde se conocerían como las funciones de Bessel.

[Link]
El matemático más grande del siglo XVIII, Leonhard Euler (1707-1783), identificó la
condición para la exactitud de las ecuaciones diferenciales de primer orden en 1734-
1735, desarrolló la teoría de los factores de integración en el mismo artículo, y dio la
solución general de las ecuaciones lineales homogéneas con coeficientes constantes.
En 1743. En 1750-1751 extendió estos resultados a las ecuaciones no homogéneas.
Aproximadamente a partir de 1750, Euler hizo uso frecuente de series de potencias
para resolver ecuaciones diferenciales. También propuso un procedimiento numérico.

Joseph-Louis Lagrange (1736-1813) demostró en 1762-1765 que la solución general

de una ecuación diferencial lineal homogénea de n-ésimo orden es una combinación
lineal den soluciones independientes. Más tarde, en 1774-1775, realizó un desarrollo
completo del método de variación de parámetros. Lagrange es conocido también por
su trabajo fundamental en ecuaciones diferenciales parciales y cálculo de variaciones.

Pierre-Simon de Laplace (1749-1827). La ecuación de Laplace es fundamental en

muchas ramas de la física matemática, y Laplace la estudió extensamente en relación
con la atracción gravitacional. La transformada de Laplace también se llama así en su
honor, aunque su utilidad para resolver ecuaciones diferenciales no fue reconocida
sino hasta mucho tiempo después.

Hacia finales del siglo XVIII se habían desarrollado muchos métodos elementales
para resolver ecuaciones diferenciales ordinarias. En el siglo xix, el interés se volvió
más hacia la investigación de cuestiones teóricas de existencia y unicidad y el
desarrollo de métodos menos elementales, como los basados en desarrollos en series
de potencias.

Otra característica de las ecuaciones diferenciales en el siglo xx fue la creación de

métodos geométricos o topológicos, en especial para ecuaciones no lineales. El
objetivo es comprender al menos el comportamiento cualitativo de las soluciones
desde un punto de vista geométrico así como analítico.

Aunque constituyen un tema antiguo del que se sabe mucho, en el amanecer del siglo
xxi las ecuaciones diferenciales siguen siendo una fuente inagotable de problemas
fascinantes y decisivos por resolver.

[Link]
2.2 DEFINICIÓN Y CLASIFICACIÓN DE LAS ECUACIONES DIFERENCIALES

Una ecuación diferencial es aquella que contiene derivadas de una función de una o más
variables. La ecuación se denomina ordinaria cuando las, derivadas que aparecen en ella,
corresponden a una función en una variable independiente y(x), por ejemplo: las
ecuaciones:

dy d y d y
a 4 5, 10, 5

Son ecuaciones diferenciales ordinarias. En el ejemplo a) aparece la primera derivada y la

función y(x); en el ejemplo b) aparecen la primera y la segunda derivada de la función y(x),
además la función y(x). El ejemplo c) es similar al b), sólo que la segunda derivada de y(x)
aparece elevada al cuadrado.

La ecuación diferencial se denomina parcial, cuando las derivadas que aparecen,

corresponden a una función de dos o más variables independientes. Así por ejemplo, las
ecuaciones:

∂z ∂z ∂ u ∂u
a 2 5 ,
∂ ∂ ∂ ∂y ∂y

Son ecuaciones diferenciales parciales. En el ejemplo a) aparecen la primeras derivadas

parciales en x y de la función z (x, y) y en el ejemplo b), aparecen la segunda derivada
parcial cruzada de la función u (x, y) con respecto a x e y, elevada al cuadrado; y la primera
derivada parcial de la ecuación u (x, y) con respecto a y elevado al cubo.

El orden de una ecuación diferencial lo determina el número de la derivada más alta que
aparece en la ecuación y el grado es el exponente de la derivada que determina el orden.
Para los ejemplos anteriores, tenemos

4 5
!
"
La ecuación diferencial es ordinaria de primer orden y primer
grado EDO (1,1)

10
#! !
"# "
La ecuación diferencial es ordinaria de segundo orden y
primer grado EDO (2,1)

$ % 5 0
#! !
"# "
La ecuación diferencial es ordinaria de segundo orden y
segundo grado EDO (2,2)

[Link]
 2 5
&' &"
&! &!
La ecuación diferencial es parcial de primer orden y primer
EDP (1,1)

$ % $ %
&# ( &'
&" ! &!
La ecuación diferencial es ordinaria de segundo orden y
segundo grado EDO (2,2)

En el curso de Matemáticas II, se resolvieron ecuaciones diferenciales ordinarias cuando se

trató el tema de la integración, ya que toda ecuación que contenga derivadas, es una
ecuación diferencial.

Ejemplo 2.1

Hallar la solución de la ecuación # 10 En otras palabras, se quiere encontrar una

)# !
)"
función en una variable independiente y(x) tal que su segunda derivada sea la constante 10.
Es claro que en este punto del curso el estudiante puede resolver este ejercicio al integrar
dos veces. Veamos,

La ecuación anterior se puede reescribir en la forma

dy
10
dx

Para obtener el diferencial del lado izquierdo integramos miembro a miembro, así:

dy
+ + 10
dx

Cuya solución es $)-% 10

),

La constante de integración c, se añade por aquello de la integración indefinida.

La última ecuación es otra ecuación diferencial, que se puede reescribir como:

.10 /

Integrando de nuevo miembro a miembro obtenemos:

+ +.10 /

Y así:

[Link]
 5 /

Es la solución general de la ecuación dada.

En general representamos una ecuación diferencial de orden n como:

! )# , )0 ,
,
"
F(x, y, , "#
…… "0
.) =0 una forma de escribir una Ecuación diferencial

Otra forma es

d1 y d15/ y
2. , , , … . .,
1 15/
y(n) = f(x, y, y´,….., y (n-1))

Ejemplo
6
1. 3y´´ +5y= 6x por tanto y´´=2x- y
6
2. 3y´´+ 5(y´)3+6y3 =0 por tanto y´´= - (y´)3-2y3
/ 6
7
3. 3(y´´)2 +5(y´)3+6y´´´=0 por tanto, y´´´= - (y´´)2- (y´)3

Ecuaciones diferenciales según su linealidad

Decimos que la ecuación diferencial de orden n
F(x, y, y´,…..y(n) )=0
Es una ecuación diferencial de orden n si y, y´, y(n) son linealmente independiente.
Cualquier ecuación diferencial que no sea de la forma:
an(x)y(n) + an-1(x)yn-1+…..+ a2(x) y´´+ a1(x)y´+ ao (x)y= 0
con an(x) ≠ 0
recibe el nombre de Ecuación diferencial no lineal

[Link]
Ejemplos
1. 3x2y´´´-5xy´´+6y´-7y =0 Ecuación Diferencial Lineal de orden 3
2. 4y´´-5y´+6y = 0 Ecuación Diferencial Lineal
3. (cos x)y´´- 3y=0 Ecuación Diferencial Lineal
4. -2y´´ +4y´-7y2= 0 No es una Ecuación Diferencial Lineal
5. y´´ +3y´+5ex+y= 0 No es una Ecuación Diferencial Lineal
6. y´´ +3y´+5 ex+ln y= 0 Ecuación Diferencial Lineal
7. y´= sen(x+y) No es una Ecuación Diferencial Lineal
Solución de una Ecuación Diferencial
Decimos que la función Φ definida en el Intervalo I es solución de la ecuación F(x, y,
y´…..y(n)) = 0 , si al reemplazar Φ, Φ´…….. Φ(n) en la ecuación diferencial, esta se reduce a
la identidad, es decir que
F(x, Φ, Φ´, Φ(n)) = 0
Observación
La función Φ está definida en un intervalo I
La función Φ es n veces diferenciable
Ejemplo.
La función y = xex es solución de la ecuación diferencial y´´-2y´+y = 0 en el intervalo
I= (-∞,∞)
Entonces
y´= ex + x ex y´= ex(1+x)
y´´ = ex+ 1 ex + x ex y´´ = ex(2+x)
Reemplazando en la ecuación diferencial tenemos:
y´´- 2y´+y = ex (x+2) – 2 ex(x+1) + x ex
= ex((x+2)-2(x+1)+x)
= ex(0)
=0

[Link]
2.3 SOLUCIÓN DE LAS ECUACIONES DIFERENCIALES ORDINARIAS

Como se vio en el ejemplo anterior, la solución general de una ecuación diferencial es una
función que contiene tantas constantes arbitrarias independientes como el orden de la
ecuación, es decir, que si la ecuación diferencial es de segundo orden, su solución general
debe tener dos constantes de integración independientes, como en el ejemplo c1 y c2. La
solución general se dice que es una familia de curvas solución de la ecuación, ya que al
darle valores a las constantes, se obtienen distintas curvas.

Como se muestra en el diagrama de bloques, generalmente se define una de las curvas

solución que se conoce como la solución particular de la ecuación y es aquella que cumple
condiciones iniciales o de frontera. Las condiciones iniciales son aquellas que se dan para
la curva y sus derivadas en un mismo valor de la variable independiente y las condiciones
de frontera en cualquier otro valor de la variable independiente por donde pasa la curva
solución.

Ecuación Técnica de Solución Solución

Diferencial Solución General Particular

Condiciones Iniciales
o de Frontera

Solución General y Particular de una Ecuación Diferencial

Ejemplo 2.2

Para la ecuación diferencial del ejemplo 2.1, halle la solución particular que satisface las
condiciones iniciales y (0) = 1 y y´(0) = 1.

Las condiciones iniciales indican que de la familia de curvas solución, nos interesa aquella
función que pasa por el punto (0,1) y cuya derivada en x= 0, vale 1. Entonces, hacemos en
la ecuación y = 5x2 + c1x + c1, y = 1 y x= 0 para obtener.

1 5.0 / .0

Y así 8 1

Como .0 1, entonces en la ecuación " 10 0 y 1, para

! !
/ , hacemos "
obtener 1 10 .0 / , y así / 1. Luego la solución particular es entonces:

. 5 1

[Link]
Cuya gráfica corresponde a la parábola que se muestra en la figura:

2,5

1,5

0,5

0
-1 -0,5 0 0,5 1

Ejemplo 2.3

Demostrar que . , 9' , 9 5' es la solución general de la ecuación diferencial

: 0, en donde
#!

"# / y son constantes arbitrarias.

Para llevar a cabo la demostración, simplemente debemos derivar dos veces la función y
reemplazar en la ecuación diferencial para verificar una igualdad.

Entonces, si . /9
;"
9 5;" , se sigue que

/ :9 :9 5;" , y / : 9 ;" : 9 5;"

! ;"
#!

" "#

Que reemplazando en la ecuación diferencial, se tiene:

/ : 9 ;" : 9 5;" : . / 9 ;" 9 5;" 0

2.4 ECUACIONES DIFERENCIALES ORDINARIAS DE PRIMER ORDEN Y

PRIMER GRADO

[Link]
A las ecuaciones diferenciales que se han resuelto hasta ahora se les ha aplicado técnicas
de integración directas. Sin embargo, no siempre es posible encontrar la solución aplicando
anti-derivadas. En lo que sigue, se analizan métodos para resolver ecuaciones diferenciales
de primer orden y primer grado cuyas soluciones no pueden obtenerse aplicando andi-
derivadas.

Una ecuación diferencial de primer orden y primer grado se puede escribir en una de las
dos formas siguientes:

1) <. , =. , 0

2. ,
!
"
2)

Observe que si en la forma ii) f solo depende de x o es una constante, entonces la ecuación
diferencial se puede resolver aplicando las técnicas de integración directas ya vistas.

2.4.1 Ecuaciones diferenciales de variables separables.

Si las funciones M y N en la forma i) son funciones de x, y de y, respectivamente, entonces

la ecuación diferencial se puede reescribir en la forma M (x) dx + N (y) dy = 0, y se dice
que es una ecuación diferencial de variables separables, la cual, se puede resolver
integrando término a término, así:

+ <. + =. 8

Una vez se encuentren las integrales para M y N, no es necesario introducir más constantes
de integración ya que la constante C representa la única constante de integración por ser la
ecuación diferencial de primer orden y primer grado.

Ejemplo 2.6

3) Halle la solución general de la ecuación diferencial . . 0 .

Reacomodando la ecuación por factorización, tenemos

. 1 . 1 0

[Link]
Dividiendo todos los términos por el factor . 1 . 1 , se llega a la ecuación
0, que es una ecuación de variables separables.
" !
"# >/ # 5/
!

8
" !
Integrando término a término, tenemos: ? ? ! # 5/
" # >/

Recodando las técnicas de integración por sustitución, se puede encontrar las dos
integrantes anteriores.

Haciendo 1, en la primera integral, entonces 2 que sustituyendo

obtenemos

1 1
+ ln | |
2 2

La misma técnica se aplica a la segunda integral y se llega a la solución,

1 1
ln. 1 ln. 1 8
2 2

Multiplicando por 2 todos los términos, llegamos a ln. 1 ln. 1 28 C.

La constante 2C, es otra constante que se ha reemplazado por A.

Ahora aplicamos las propiedades de los logaritmos ln D ln ln.D y obtenemos

ln. 1 . 1 C

Para obtener el argumento de la función logaritmo, aplicamos la propiedad e;! x, y así,

9 ;! " 9F G
# 5/"! # 5/

Luego,

. 1 . 1 G

De nuevo, la constante 9 " , se ha reemplazado por otra constante B. de la última ecuación se

puede despejar y, para obtener la solución general

G
. H 1
1

Ejemplo 2.7

[Link]
Halle la solución de la ecuación diferencial que describe el hecho de que cuando los
factores ambientales imponen un límite superior sobre su tamaño, la población crece a un
ritmo que es conjuntamente proporcional a su tamaño actual y a la diferencia entre su límite
superior y su tamaño actual.

Solución

Sea P (t) el tamaño de la población en el tiempo t, y L el límite superior impuesto a la

población por el medio ambiente. Entonces,

I
J
B ritmo de crecimiento de la población será

Y la diferencia entre el límite superior y la población K L

Como “conjuntamente proporcional” significa “proporcional al producto”, entonces la

ecuación diferencial que modela el problema es,

L
:L.K L
M

Siendo k, la constante de proporcionalidad

: M o en la forma general
I
Ajustando las variables, la ecuación queda I.N5I

: M 0, integrando término a término se tiene

I
I.N5I

L
+: M + 8
L.K L

La primera integral es inmediata, la segunda se puede resolver aplicando fracciones

parciales así:

, , que se puede reescribir como 1 C.K L GL

/ / O
I.N5I I N5I

Haciendo L K, encontramos G N y L 0, hallamos que C

/ /
N
de esta forma, la
solución de la ecuación diferencial será:

:M ? PNI Q R 8, para obtener :M ln L ln.K L 8

/ / / /
N.N5I N N

Multiplicando por L la última ecuación, se obtiene :KM S L ln.K L 8K, ahora

aplicamos las propiedades de los logaritmos y simplificamos.

[Link]
 K L
:KM S C
L
K L
S C :KM
L
K L
9 ;J ⟹ K L L9 ";J ⟹ K L L9 ";J ⟹ K L.1 9 /5;J
L

Finalmente, la solución general es

K
L.M
1 9 F5;NJ

Ejemplo 2.8

El valor de reventa de cierta maquinaria industrial decrece durante el periodo de 10 años a

un ritmo que depende de la antigüedad de la maquinaria, cuando la maquinara tiene t años,
el ritmo al que cambia su valor es 220 (t-10) dólares por año. Exprese el valor de la
maquinaria como una función de su antigüedad y del valor inicial. Si la maquinaria valía
originalmente US $12.000 =¿Cuánto valdrá cuando tenga 10 años?

Solución

tSea v(T) el valor de la maquinaria al cabo de t años, entonces la ecuación diferencial que
modela el problema es,

U
220.M 10
M

De donde se obtiene V.M ? 220.M 10 M 110M 2200M 8

Como V.0 12.000 entonces V.M 110M 2200M 12.000 y el valor de la

maquinaria después de 10 años es

V.10 110.10 2200.10 12.000 WX$1.000

WX$30.000 Después de 4 años. ¿Cuánto valdrá la maquinaría cuando tenga 8 años?

Para cierto bien las ecuaciones de oferta y de demanda son las siguientes:

[: R 2 I 25

[Link]
 X: R 3 ] 5

Supongamos que si el mercado no está en equilibrio ^ _ ] `, entonces el precio cambia

:^ _ ] `. Sustituya _
_
J
en razón proporcional al exceso de demanda sobre la oferta
y ] , y resuelva la ecuación diferencial resultante para R.M . Pruebe que no importa cuál sea
el precio inicial, el mercado se aproxima eventualmente al equilibrio en R 17.

Suponga que la tasa de crecimiento proporcional ´.M , .M de la población de la tierra es

mil millones. Considerando a 1930 como M 0, determine la población .M de la tierra

una constante. La población en 1930 era de 2 mil millones de habitantes y en 1975 fue de 4

en el instante M. De acuerdo con este modelo, ¿Cuál debió ser la población en 1960?

El ritmo al que se propaga una epidemia en una comunidad es conjuntamente proporcional

a la cantidad de residentes que ha sido infectada y al número de residentes propensos a la
enfermedad que no ha sido infectado. Exprese el número de residentes que ha sido
infectado como una función del tiempo.

2.7 APLICACIONES DE LAS ECUACIONES DIFERENCIALES EN LOS

MODELOS ECONÓMICOS

1. Macro Modelo de Domar

S (t) es el ahorro, l (t) es la inversión, y (t) la renta, todas las funciones del tiempo.

X.M D .1

S.S c
!
J

X.S K.M

d > 0 c > 0

Hallar la solución particular _ .K , K_ .M X_ .M SI f cuando M 0

Realizar un gráfico.

2. Modelos de Deuda de Domar

[ es una deuda nacional y _ .M la renta nacional

[Link]
 d .M
g
J

c
!
J

.0 f

[.0 [f

d > 0 c > 0

Analizar la relación entre la deuda nacional y la renta nacional en M → ∞.

3. Un segundo modelo de deuda de Domar

D es la deuda nacional y y (t) la renta nacional.

c .M
!
J

.0 f

[.0 [f

d > 0 c > 0

Analizar la relación entre la deuda nacional y la renta nacional en cuanto M → ∞.

4. Modelo de ajuste del precio de Evans

.M df d/ L.M Demanda

j.M cf c/ L.M Oferta

k. j
_
J

df > 0 d/ > 0 c/ > 0 k > 0

Hallar R.M si R.0 Lf . Tenga en cuenta que LJ

lm 5no
nm 5pm
es el precio de equilibrio.

[Link]
5. Modelo de Demanda y Oferta

La demanda y la oferta (por unidad de tiempo) de un artículo están dadas respectivamente

por x y y, siendo p el precio unitario:

DR R

decrece a una tasa proporcional a dicho exceso. Encuentre R.M si R.0 Lf . Demostrar
Supóngase que el precio cambia de tal manera que el exceso de la demanda sobre la oferta

que el precio unitario tiende al valor de equilibrio dado por Rq n5p

r5

Ejemplo 2.10

La relación entre el costo de manufactura por artículo, C, y el número de clases de

artículos fabricados, x, es tal que la tasa de incremento del costo de manufactura, a medida
que aumenta el número de clases, es igual a la razón del costo por artículo más el número
de clases, dividido todo entre el número de clases de artículos que se manufacturan.

productos fabricados si 8 8f cuando 1.

Obtener la relación entre el costo de fabricación por artículo y el número de clases de

Solución

s s>"
" "
La ecuación diferencial que modela el problema es que se puede reescribir en la
1, y así, 2. 1.
s s
" "
forma

Aplicando la fórmula, tenemos: ln :

' s
? .'>/ 5' ? "
. Luego la solución
general es

8. ln : siendo : la constante de integración.

Solución particular

Como 8 f Cuando 1 entonces, al reemplazar en la solución general tenemos:

8f .t ln : .u

De donde se obtiene para : 9 so . Reemplazando este valor de : en la solución general,

tenemos la solución particular, 8. ln 9 vw . f ln

. .8f ln

[Link]
TEOREMA DE EXISTENCIA Y UNICIDAD
Sea R una región rectangular del plano xy definida por a≤ x ≤b y c≤y≤d que contiene al
xy
punto (xo, yo) en su interior. Si f(x,y) y x! son continuas en R, entonces existe un
intervalo I0 : x0 – h < x < x0 + h , h>0 contenido en a≤ x ≤b , y una función Única y(x)
definida en I0 que es una solución del problema de valor inicial
y´= f(x,y)
y(x0)= y0
Si y´= y
Cuál será la solución del problema de valor inicial y´=y , y(0)=1
x0=0 y0= 1 donde (x0, y0)= (0,1)
Un problema de valor inicial tiene más de una solución siempre que la Ecuación
Diferencial satisfaga algunas condiciones, si las satisface es única

Ejemplos:
Determine si el Teorema de Existencia y Unicidad garantiza la unicidad de la solución en
los siguientes Problemas de valores iniciales
a) y´= y , y(0)=0

1
xy
x!
Solución: f (x,y)= y ; (x0, y0)= (0,0) ;

xy
x!
tanto f(x,y) como son continuas en R2 ,como (0,0) ∈ R2, entonces el Teorema
Fundamental de Existencia y Unicidad garantiza la unicidad de una única solución del
problema de valor inicial y´= y , y(0)=0 en el intervalo I= R= (-∞,∞)
b) y´=y1/2 , y(0)=0
xy /
x! √!
Solución: f (x,y)= y1/2 ; (x0, y0)= (0,0); =

xy
x!
Donde es continua en R ={(x,y) ∈ R2/ y>0}, como (0,0) ∉R entonces el Teorema de
Existencia y Unicidad no garantiza la unicidad de la solución (pueden ser varias soluciones)
c) y´=y1/2 , y(1)=0
xy /
x! √!
Solución: f (x,y)= y1/2 ; (x0, y0)= (1,0); =

[Link]
 xy
x!
Donde es continua en R ={(x,y) ∈ R2/ y>0}, como (1,0) ∉R entonces el Teorema de
Existencia y Unicidad no garantiza la unicidad de la solución .

EJERCICIOS
a) y´= y1/2 , y(0)=1
b) y´= sen(x+y) , y(π)=0
c) y´= ln (x/y) , y(1)=1

ECUACIÓN DIFERENCIAL DE PRIMER ORDEN

1. ECUACIÓN DIFERENCIAL LINEAL DE PRIMER ORDEN

Dado el operador diferencial lineal L, de C/ .I en C.I , definido por L.y y´

P .x y , donde y pertenece a C/ .I , y, P es una función continua en C.I .

Estamos interesados en hallar la solución de L.y Q.x , donde Q es una función

solución de y´ 1 . .y 0 , en el segundo solucionaré la ecuación y´ P .x y

continua en (I). Para ello presentaré tres teoremas, en el primero hallaré la

0, y , finalmente buscaremos la solución de y´ P.x y Q.x .

ecuación diferencial y´.x y .x 0 , para todo x real, y, además satisface la

Teorema 1. Si b es un real dado, existe una y solo una función f que satisface la

condición y.0 b. Esta función viene dada por f.x be- , para todo x real.

condiciones de teorema o sea y´.x y .x 0, y, y .Q b. En efecto como

Demostración: comprobemos que la función f así definida satisface las

f.x be- , luego f´.x be- , entonces f´.x f.x be- be- 0 para todo x
real, además f.0 bef b. 1 b.

Entonces f.x be- satisface la ecuación diferencial y la condición dada,

probemos que esta función es única, sea g una función que verifica g´.x g.x
0 , y, g.0 b; definimos: h.x g.x e para todo x real, luego h es diferenciable
5-

en R.

Entonces

h´.x g´.x e5- g.x e5-

[Link]
 e5- [g´.x g.x ]

e5- . 0

0 Para todo x real

Luego h es una función constante, y, h.x h.0 g.0 e5f b para todo x real,
luego h.x g.x e 5-
b, o sea g.x be para todo x real, y era a donde
-

deseábamos llegar.

un real dado. Existe una y solo una función f que satisface la ecuación y´.x
Teorema 2. Sea P una función continua en un intervalo I. Sea a un punto de 1, y, b

P.x y.x 0, para todo x en I, además y.a b. Esta función viene dada por
f.x be ?‰ P.t dt.
5†.- -
dondeA.x

Demostración: veamos que f.x be5†.- donde A.x ?‰ P.t dt Verifica que
-

f.x ´ p.x f.x 0. Para todo x en I, y, f (a) = b.

En efecto como

f.x be5†.- , y, f´.x bA´.x e5†.-

bP.x e5†.-

Entonces

f´.x P.x f.x bP.x e5†.- P.x be5†.-

0. Para todo x en I.

Como

A.a ?‰ P.t dt 0 , luego y .a be5†.‰ b.

‰

satisface las condiciones dadas, sea esta g entonces g´.x P.x g.x 0 para x
Probemos ahora que esta función es única, para ello tomemos una función que

en I y g.a b; definimos:

h.x g.x e†.‹ para x en l

Derivando h obtenemos:

[Link]
h´ .x g´.x e†.‹ A´.x g.x e†.-

e†.‹ [g´.x P .x g .x ]

0 para x en I.

Luego h es una función constante en I.

h .x h.a g.a o†.‰ b, o sea g.x o†.‰ b , para todo x en I, entonces

g.x be 5†.-
, para todo x en I, lo cual verifica que la solución es única.

Teorema 3. Sean P y Q, funciones continuas en un intervalo I. Sea a un punto de I,

y, b un número real dado.

y´.x P.x y.x Q.x , y, la condición y .a b.

Entonces existe una y solo una función f que satisface la ecuación diferencial

Esta función es f.x be5†.- e5†.- . ?‰ Q.t e†.• dt.

-

Donde A.x ?‰ P.t dt

-

Demostración: es fácil verificar que f así definida satisface las condiciones del
teorema, probemos que esta función es única, para ello tomemos una función g
que satisface las condiciones dadas.

Definimos h.x g.x e†.- para x en I, derivando h tenemos:

h´.x g´.x e†.- A´.x g.x e†.-

e†.- [g´.x A´.x g.x ]

e†.- [g´.x P.x g.x ]

e†.- Q.x

Por el primer teorema fundamental del Cálculo tenemos:

h.x h.a ?‰ e†.• Q.t dt. Para x en I.

-

Luego
-
h.x h.a + e†.• Q.t dt
‰

[Link]
 y g.x e†.- g.a e†.‰ ?‰ e†.• Q.t dt
-

g.x be5†.- e5†.- . ?‰ Q.t e†.• dt para x en I.

-

Ejemplo: hallar la solución de la ecuación diferencial lineal xy´ 2y x 6 , para x

mayor que cero, y, y (1) = 1.

Como x es diferente de cero, entonces la ecuación la podemos escribir así:

y´ x Ž , luego P.x , y, Q.x xŽ

,
- -

Aplicando el teorema 3, la solución es:

-
y.x y .1 e 5†.-
e 5†.-
+ t Ž e†.• dt
/

Donde:
-
2 -
A.x + [ dt] + 2t 5/ dt 2 ln x
/ t /

Luego:

y.x 1 e e ?/ t Ž e5 •• • dt
•• - •• - - #

x x ?/ t Ž . t 5 dt
-

x x ?/ t 5 dt
-

x x ‘ “
-’ /

para todo x mayor que cero.

-# -”

2. ECUACIONES DE VARIABLES SEPARABLES DE PRIMER ORDEN

Una ecuación diferencial de la forma y´ f.x, y se denomina de variables

separables si f.x, y Q.x R.y , o sea se puede expresar como un producto de

y´ Q.x R.y , donde Q y R son funciones dadas, luego la ecuación la podemos

dos factores, uno dependiente de x únicamente y el otro solamente de y; es decir

[Link]
escribir así: P.y y´ Q.x , si R(y) ≠, y, P.y [R.y ]5/; la solución a dicha
ecuación viene caracterizada por la siguiente proposición.

Teorema 4. Sea y y.x una solución cualquiera de la ecuación diferencial

P.y y´ Q .x (1)

Tal que Y es un elemento de c/ .I . Sean P y PoY elementos de c.I , definimos

h´.x P.x , para x en I. entonces y verifica

h.y ?/ Q.x dx c
-
(2)

Para un cierto valor c. Además si y satisface la ecuación (2) entonces y es

solución de (1).

Demostración: supongamos que y Y.x para x en I es solución de P .y ´

Q .x

Luego P [ Y .x ] y´.x Q .x para x en I.

Como h´.x P .x , entonces h´[Y .x ] Y´.x Q .x

Utilizando la regla de la cadena, tenemos que:

[.h0Y .x ]´ Q.x

Luego

h[Y .x ] + Q.x dx c

Para un cierto valor de C. como y y .x entonces

h .y + Q.x dx c

si derivamos obtenemos [h´.y ]y´ Q.x

o sea P.y y´ Q .x para x en I.

Lo cual prueba que si y satisface (2), entonces y es solución de (1).

Es importante notar que como y y.x para x en I, luego dy Y´.x dx

[Link]
Entonces P.y y´ P[Y .x ]Y´.x dx Q .x , e integrando con respecto a x tenemos:

+ P.y y´ + P.y dy + Q .x dx c

),
)-
si colocamos y´ de la forma y lo consideramos como un cociente obtenemos la
relación P.y dy Q .x dx , ponemos los signos de integración a ambos miembros
de esa ecuación y sumamos una constante obtenemos (2), lo cual muestra la
eficacia de la notación de Leibniz.

Consideramos la ecuación de primer orden y´ f .x, y (1)

La cual verifica la propiedad de homogeneidad, es decir

f .tx, ty f .x, y (2), para todo x, y, y, con t φ 0.

Aplicando esta condición para x

/
•
a (1), tenemos:

y´ f $1, -%
,
(3)

Sea u , luego y ux derivando tenemos y´ u´x u

,
-

sustituyendo en (3), tenemos: u´x u f .1, u

o sea x f .1, u u.
)(
)-

La cual es una ecuación de variables separables de primer orden.

3. ALGUNOS PROBLEMAS CUYO MODELO MATEMÁTICO ES UNA

EDUCACIÓN DIFERENCIAL DE PRIMER ORDEN

3.1 CRECIMIENTO O DESINTEGRACIÓN

Sea y f .t la cantidad de sustancia existente de un cuerpo en un tiempo t,

entonces la ecuación diferencial y´ ky (1)

Indica que la intensidad o velocidad de variación con relación a la variable t, de

una cantidad y, es proporcional a la cantidad existente del elemento en ese mismo
instante; si la constante de proporcionalidad k, y, y son positivas, entonces y´ es

[Link]
positiva y f (t) es creciente, en este caso decimos que el problema es de
crecimiento. Pero, si k es negativa, y, y es positiva, entonces y´ será negativa, lo
cual implica que f (t) es decreciente, y el problema es de desintegración.

El valor de K depende del elemento particular que se está tratando; la solución a la

ecuación (1) es f .t f .0 e™• en base al teorema 3; donde f .0 es la cantidad
inicial del elemento.

Si el proceso es de desintegración, luego k es menor que cero entonces f.t

f.0 e™• , en ningún tiempo t se anula, puesto que e™• es siempre positiva, por tanto
no se puede hablar de “tiempo total de vida” de una sustancia radiactiva. Pero es
posible determinar el tiempo necesario para que se desintegre una fracción parcial
de la muestra.

Generalmente se elige la mitad de la sustancia, es decir f .t f.0 , luego t se

/

denomina vida media de la sustancia y tenemos que:

1
f.t f.0 e™• f.0
2

Luego e™•
/
, tomando logaritmo natural, tenemos:

1
ln e™• ln
2
1
kt ln
2

T Pln Q
/ / ••
™ ™
con k menor que cero.

La gráfica de f.t f.0 e™• con k menor que cero es:

[Link]
Ejemplo: la vida media del radio es de 1.700 años, qué porcentaje de radio es
probable que quede al cabo de 50 años, y al cabo de 2.000 años.

Solución: sea y f.t la cantidad de radio existente en un tiempo t, como y´

ky con k mayor que cero, luego la ecuación diferencial que determinará el
problema es:

y´ ky 0 , cuysolución es f.t f .0 e5™•

hallemos k, utilizando la vida media del radio la cual es de 1.700 años, entonces:

f.1700 f.0 f.0 e5/›ff™ , entonces e5/›ff™

/ /

tomando logaritmo natural tenemos:

ln $ % ln e5/›ff™ 1700 k, luego k

/ ••
/›ff

Luego la solución a la ecuación diferencial dada es:

f .t f .0 e5mžoo • ,
ϥ #

Y f.50 f .0 exp. $ . 50% f .0 exp. $ %

5 •• ••
/›ff Ž

Es la cantidad existente en 50 años, o sea esta es la cantidad que no se ha

desintegrado, su porcentaje es:

ln 2
f .0 exp. $ % . 100 10
34 Ÿ 9,79818
f .0 ln 2
exp $ %
34

Para calcular el porcentaje que no se ha desintegrado al cabo de 2.000 años,

hallamos la cantidad de radio que no se ha desintegrado a esta fecha, la cual es:

ln 2 ln2
f .2000 f.0 exp. . 2000 f.0 exp. 20
1700 17

Y el porcentaje correspondiente que no se ha desintegrado es:

[Link]
 20 ln 2
100 f.0 exp. $ % 100
17 44,2432
f .0 20 ln2
exp. $
17 %

3.2 ENFRIAMIENTO

Se calienta agua hasta su punto de ebullición 100º, se retira del fuego y se

mantiene en un recinto que se encuentra a una temperatura constante de 60º. Al
final de 3 minutos la temperatura del agua es de 90º. Calcular la temperatura del
agua al final de 6 minutos y el tiempo en la cual la temperatura del agua será de
75º.

Solución: sea y f .t la temperatura de un cuerpo en un tiempo t ; la velocidad

temperatura y la del medio ambiente (Ley de Enfriamiento de Newton). Sea M.u

de variación de y con respecto al tiempo t es proporcional a la diferencia entre su

la temperatura del medio ambiente. Luego la ecuación diferencial que sirve de

modelo a este problema es:

y´ k ^y M.u ` , o, y´ ky k M .u

Donde k es una constante positiva, esta ecuación es diferencial lineal de primer

orden, luego su solución viene dada por el teorema 3, es decir:
•
f .t f .a e 5™•
e 5™•
+ k M.u e™- dx
‰

donde f .a es la temperatura del cuerpo en un tiempo t a para nuestro problema

tenemos f .0 100, la cual es la condición inicial dada, además f.3 90, y,
M.u 60.

Y la ecuación diferencial que define el problema es:

y´ ky 60 k, su solución es:

f .t 100 e5™• e5™• ?f 60 ke™- dx

•

100 e5™• e5™• ¥60 e™• 60¦

100 e5™• e5™• 60

40e5™• 60

[Link]
Hallamos k utilizando la condición de frontera f .3 90º

f .3 90 40e5 ™
60

30 40e5 ™

e
Ž ™
, tomando logaritmo natural, tenemos:

ln ln e 3k, luego k ln
Ž ™ / Ž

Entonces la solución a la ecuación diferencial es

1 4
f.t 40 exp ln t 60
3 3

Hallemos la temperatura al final de 6 minutos, calculamos

f.6 40 exp $ ln % 6 60
/ Ž

40 e5
¨
••
’ 60

40 e5 •• ª
m©
60

60 82,5º C
Ž6

El instante en el cual la temperatura es de 75º, lo podemos calcular así:

75 f.t 40 exp $ ln % t 60
/ Ž

40 e5’ ln
«
60
Ž

15 40 e
«
ln 40 $ %
5 Ž ’
’
Ž

« «

$ % , luego $ %
/6 ’ ’
Žf Ž ¬ ¬

Aplicando logaritmo natural, hallamos que t 10.2 seg.

[Link]
3.3 APLICACIONES A LA MECÁNICA

Una gota de aceite, cuy masa es de 0,2 grs, cae partiendo del reposo. Cuando su
velocidad es de 40 centímetros por segundo, la fuerza debida a la resistencia del
aíre es de 160 dinas.

Suponiendo que dicha fuerza es proporcional a la velocidad instantánea:

a) Hállense la velocidad y la distancia recorrida en función del tiempo.

b) Determinar la velocidad límite

Solución: la gota afectada por dos fuerzas una debida a su peso W, la cual obliga
a la gota a ir hacia abajo y otra que se opone a esta y es la resistencia del aíre R,
la cual es proporcional a la velocidad instantánea. Consideraremos positivo hacia
abajo y negativo hacia arriba. Entonces según la segunda ley de Newton tenemos:

F W kv

Donde k es la constante de proporcionalidad y es positiva, luego la ecuación

diferencial correspondiente es:

mv´ mg kv

Donde m es la masa, v la velocidad, g la gravedad, esta ecuación se puede

transformar en:

kv
v´ g
m

La cual es una ecuación diferencial lineal de primer orden, su solución es:

•
v .t v .0 e5† .• e 5† .• + ge† .- dx
f

Donde
•
k k
A .t + dx t
f m m

Luego

[Link]
 ™• ™• • ™-
v .t v .0 e5² e5² + ge ² dx
f

™• ™• gm ™• gm
v .0 e5² e5² ³ e² ´
k k
™• gm ™•
v .0 e5² ³1 e5 ² ´
k

Si v .0 0, como en nuestro caso entonces:

gm ™•
v .t ³1 e5 ² ´
k

De acuerdo a la condición dada que la resistencia al movimiento es de 160 dinas

proporcionar a la velocidad instantánea, entonces podemos hallar el valor de k, ya

cuando la velocidad es de 40 centímetros por segundo y dicha fuerza es

que 160 40 k, luego k 4. Reemplazando este valor en la solución final tenemos

que v .t 49 [1 e5 f• ].

Para hallar la distancia recorrida en función del tiempo podemos emplear la

condición que v .t luego:
)µ
)•

ds
49 49 e5 f•
dt

Esta es una ecuación de variables separables, entonces

µ • •
+ ds + 49 dt + 49 e5 f•
dt
f f f

49 5
s .t 49t [e f•
1]
20

La velocidad límite se tiene cuanto t se hace tan grande que podemos decir que
t → ∞, entonces en este caso e5 f• tiende a cero, luego a velocidad límite es 49
centímetros por segundo.

Otro ejemplo de este mismo aspecto es:

fuerza atractiva en 0 que varía directamente con su distancia a 0. Cuando t 0 la

Una partícula se mueve en línea recta hacia un punto fijo 0, bajo la acción de una

[Link]
partícula dista 4 cms de 0 y se mueve hacia 0 con una velocidad de 6 centímetros
por segundo, y una aceleración de 16 centímetros por segundo cuadrado.

a) Hállense su posición y velocidad de partícula en función del tiempo.

b) Obténgase la amplitud, el periodo y la frecuencia del movimiento.
c) Determínese la velocidad y aceleración máxima

Solución:

supongamos que la partícula inicia su movimiento en un punto A, tomemos 0A

Tomemos un punto fijo 0 como origen de un sistema de coordenadas,

como sentido positivo. Sea P la posición de la partícula en un instante t. Como la

magnitud de la fuerza de atracción F hacia 0 es proporcional a su distancia a este
punto tendremos:

F k. x

Con K mayor que cero.

Supongamos x mayor o igual que cero, luego la magnitud de la fuerza es kx.

Como la fuerza F está dirigida hacia 0 (en sentido negativo) tendremos que F
kx con x mayor o igual a cero.

Si x es menor que cero, luego la magnitud de la fuerza es – kx, como la fuerza F

está dirigida hacia la derecha (en sentido positivo) luego F kx con x menor que
cero; luego la ecuación diferencial que describe este fenómeno físico es:

d x
m kx
dt

Para nuestro caso la velocidad y la aceleración está dirigida hacia Q, cuando dista
4 cms, o sea para t 0 tendremos: x .0 4, v.0 6, a.0 16, cada una
con sus respectivas unidades del sistema CGS, ahora podemos calcular el
cociente ² ya que 16 ² . 4, luego ² 4.
™ ™ ™

Y la ecuación diferencial original se transforma en:

[Link]
 d x
4x
dt

Para resolver esta ecuación de segundo orden la transformamos en una de primer

orden, mediante el siguiente cambio de variable, como v .t
)-
)•
, de modo que

d x dv dv dx dv
. v.
dt dt dx dt dt

Luego la ecuación diferencial de segundo orden se convierte en:

dv
v. 4x
dx

La cual es una ecuación diferencial de variables separables, su solución está dada

por:
¸ -
+ v dv + 4x dx
57 Ž

v .t 36
Entonces:
2x 32
2 2

v .t 4.25 x

v.t ±2 º25 x

Como v .t
)-
)•
, reemplazando en la última ecuación tenemos una ecuación de
variables separables.

dx
± 2 º25 x
dt

Su solución es:
-
dx •
+ + 2 dt
Ž √25 x f

arc sen 6 arc sen 2t

- Ž
6
(1)

[Link]
Definimos u arc sen , luego sen u y, cos u
Ž
6 6 6
, luego la ecuación (1) se
transforma en:
x
arc sen 2t u
5

Aplicando la función seno tenemos:

x 5 sen .5 2t

5 cos 2t sen u sen 2 cos u

5 cos 2t sen 2t
Ž
6 6

Luego

x 4 cos 2t 3 sen 2t

v 8 sen 2t 6 sen 2t

Podemos definir el signo, tanto para x como para v, usando la condición dada de
que v .0 6, lo cual implica que:

x .t 4 cos 2t 3 sen 2t

v .t 8 sen 2t 6 cos 2t

a .t 16 cos 2t 12 sen 2t

Para calcular la amplitud, periodo y frecuencia del movimiento recurrimos a la

fórmula:

x .t 5 sen .u 2t

Luego

v .t 10 cos.5 2t

a .t 20 sen .u 2t

[Link]
 5, el periodo 1π y la frecuencia
/
½
Es fácil ver que la amplitud , cada uno con
sus respectivas unidades.

La velocidad máxima es de 10 centímetros por segundo.

La aceleración máxima es de 20 centímetros por segundo cuadrado.

3.4 SOLUCIONES QUÍMICAS

A un tanque que contiene 400 litros de agua fresca, sele incorpora salmuera que
/
¬
contiene de kilogramo por litro de sal a razón de 8 litros por minuto y la mezcla
se mantiene uniforme por agitación, abandona el tanque por un orifico a razón de
4 litros por minuto.

Encontrar:

a) La cantidad de sal presente en el tanque cuando este contenga 500 litros de

salmuera
b) La concentración de sal en el tanque al final de una hora

Solución: sea y f .t la cantidad de sal existente en el tanque en un tiempo t,

luego la rapidez de cambio de la cantidad de sal con respecto al tiempo depende
de la sal que entra menos la cantidad de sal que sale.
/
¬
Como entran 8 litros por minuto y cada litro contiene de kilogramo de sal,
entonces la sal que entra en el tanque en 1 minuto es 1 kilogramo; como entran 12
litros de salmuera por minuto y de la mezcla salen 8 litros por minuto, luego el
número de litros del tanque aumenta 4 litros por minuto de manera que en un
tiempo t el volumen de salmuera en el tanque es de 4 4t, la cual contiene Žff> ŽJ
,

kilogramos de sal por litro, y de esta mezcla salen 8 litros por minuto. Luego en
cada minuto la cantidad de sal que sale del tanue es $Žff> Ž% 8 kilogramos.
,

Entonces el fenómeno está definido por la ecuación diferencial:

2y
´ 1
100 M

O sea:

´ 1 con 2 .0 0
,
/ff> J

[Link]
La cual es una ecuación homogénea lineal de primer orden y su solución está
dada por el teorema 3:
¿
2 .M 2 .0 9 5F .J
9 5F .J
+ ¾ . 9 F ."
f

Donde
¿
2
C .M + 2 ln.100 M ln 100
f 100

100 M
2 ln
100

Luego

100 M J
100
2 .M exp À ln Á . + 9 Àln Á
100 f 100

100 J
100
+
100 M f 100

1 J
+ .100
100 M f

1 M
À10Ž M 10 M Á
100 M 3

como el volumen del tanque está dado por 400 4t , para un tiempo t, luego los
Deseamos hallar la cantidad de sal existente cuando el tanque contenga 500 litros,

500 litros se tendrán para t 25 minutos y ahora podemos hallar f .25 .

2 .25 kilogramos de sal

76

Para hallar la concentración de sal al cabo de una hora, hallamos el valor del
0,072 kilogramos de sal por litro.
 .7f Ž6
à .7f 7 6
cociente

3.5 CIRCUITOS ELÉCTRICOS

[Link]
Los circuitos eléctricos están descritos en su comportamiento por la Ley de
Kirchhoff, en realidad, la teoría eléctrica está regida por un conjunto de ecuaciones
conocidas en la teoría electromagnética como ecuaciones de Maxwell, que se
salen de nuestro alcance, y para nuestro estudio es suficiente por ahora la Ley de
Kirchhoff, la cual describe una ecuación diferencial lineal de primer orden.

Estudiaremos circuitos eléctricos en serie tan elementales como:

Donde sus elementos son: E, es la fuerza electromotriz que produce un voltaje el

cual “origina” una corriente eléctrica que recorre el circuito; R es la resistencia que

intensidad de corriente, y C el condensador que es un elemento que almacena

utiliza dicha energía; L es la inductancia, que se opone a la variación de la

energía.

Para el circuito de la figura 3 su modelo matemático es:

t
K Åu Æ
M

Para el circuito de la figura 4 su modelo matemático es:

¾ ¾
Å Æ
M 8

Donde I es la intensidad de la corriente que circula por el circuito en un tiempo t, y,

Q es la carga la cual se relaciona con la intensidad mediante la siguiente
condición;

¾
u
M

Ejemplo: un generador cuya fem es 100 voltios, se conecta en serie con una

cuando t 0. Hallar la intensidad de la corriente en función del tiempo.

resistencia de 10 ohmios y una inducción de 2 henrios. Si el interruptor se cierra

Solución: su esquema gráfico del circuito está dado por la figura 3, su ecuación
diferencial al caso es:

[Link]
 dI
2 101 100
M

O sea

dI
51 50
M

Utilizando el teorema 3 tenemos:

J
I .t I .0 e 5F .¿
e 5F .¿
+ 50 9 F ."
f

Donde
C .M ?f 5 5M , e, u .0 0
J

Luego
u .M 50 9 56J ?f 9 6"
J

10 .1 9 56J

Su gráfica es:

Ejemplo: una fem decreciente, Æ 200 9 56J , se conecta en serie con una
resistencia de 20 ohmios y un condensador de 0.01 faradios. Supongamos que la

tiempo t, hallar la carga máxima y el tiempo necesario para alcanzarla.

carga es cero en un tiempo cero; hallar la carga y la intensidad de cualquier

Solución: su esquema gráfico está dado por la figura 4, y su ecuación diferencial

es:

dQ
20 100¾ 1009 56J
M

[Link]
 dQ
O sea
5¾ 109 56J
M

Su solución es:
J
Q .t ¾ .0 9 5F .J 9 5F .J . + 10 9 56" 9 F."
f

Donde
C .M ?f 5 5M , y, ¾ .0 0
¿

Luego
¾ .M 9 56J ?f 10
J

10M 9 56J

La intensidad de la corriente en un instante cualquiera es:

u .M 10 9 56J 50M 9 56J

Ç
J

0
Ç
J
Para determinar el máximo de carga, hacemos:

O sea 109 56J 50M 9 56J 0, de la cual deducimos que M

/
6
segundos, y el valor
de la carga máxima es ¾ 9 5/ 0,74 coulomba.
/ /f
6 6

3.6 SALIDA DE UN LÍQUIDO POR UN ORIFICIO

[Link]
Supongamos que el recipiente es el de la figura 6, donde C . es el área de la
sección recta del depósito a una altura y G es el área del orificio con bordes
perfectamente pulimentados. Si no existiera rozamiento no hay pérdida de energía
potencial toda se convertiría en energía cinética las cuales serían las únicas que
aparecerían en este fenómeno físico, luego:

1
ÈÉ ÈV
2

Donde es la altura de caída, luego la velocidad de salida sería de .2 É

m
# metros
por segundo.

Como el área de orificio es G entonces el producto G .2 É # representa el

m

número de metros cúbicos por segundo de líquido que sale por el orificio.

Sea c el coeficiente de rozamiento, es decir, el chorro no es completo por el orifico,

luego la velocidad de descarga es:

G.2 É
m
# metros cúbicos por segundo.

Sea U . el volumen del líquido que está en el depósito a una altura y entonces:

U . ?f C . C .
! Ã
!
, y,

Luego

V /
G .2 É
M

Este resultado es negativo porque el volumen decrece con el tiempo. Usando la

regla de la cadena tenemos:

V V
. C .
M M M

entonces combinando la dos últimas ecuaciones tenemos la ecuación que

describe nuestro suceso.

/
C . G .2 É
M

La cual es una ecuación de variables separables de primer orden.

[Link]
Ejemplo: un embudo tiene la forma de un cono circular recto, cuyo vértice está

tiempo T, encontrar el tiempo requerido para vaciarse completamente.

hacia abajo y está lleno de agua. Si la mitad del volumen se desocupa en un

Solución: el gráfico del recipiente y sus elementos son: H la altura del cono, R el
radio mayor, r el radio de la sección recta C .
2 .M la función que determina la altura del
a una altura y del vértice, como la
altura depende del tiempo, sea

sale en un tiempo T, podemos hallar la altura del líquido que queda en función de
líquido en el recipiente, de acuerdo a la condición dada que la mitad del volumen

H.

Tenemos que U .MËMDS V .MËMDS

½ Ì # Í / ½ Ì # Í
7
, y,

Î Ï
De acuerdo a la semejanza de triángulos, tenemos que: Ì Í

Î#Ï Ï’
Ì# Í Í’
Luego, como la razón entre los dos volúmenes es de 1;2,

Entonces

1 ℎ
2 Ñ

Luego

ℎ 2ℎ

de donde

Ñ
2 .Ò ℎ
√2

Como

[Link]
 Ó Å
C . Ó Ô
Ñ

Luego la ecuación diferencial C . 8 .2 É

m
!
#
J
se transforma en:

Ó Å . /
G .2 É
Ñ M

cuya solución viene dada por:

Ó Å !
3 J
+ + G M
Ñ º2É Í f

luego:

2Ó Å 5 5
³ ´ GM
5Ñ º2É Ñ

Para 0 tenemos el tiempo total, entonces

5
2Ó Å Ñ
2 GM ÒËMDS
5Ñ º2É

Luego

2Ó Å √Ó
MJwJpÕ
5 Gº2É

El cual aparece en función de H, R, C y B, y lo importante es dejarlo en función

únicamente de T, pero la condición dada de que 2 .Ò
Í
√
, luego podemos
reemplazar esta condición en la solución a la ecuación diferencial y tenemos:

1
2Ó Å Ñ
2³ 1 1´ GM
5 G º2É √32
©

2Ó Å √Ñ √32
©

Ò ©
5 G º2É √32 1

Luego

[Link]
 MJwJpÕ 2,28 Ò

3.7 TRAYECTORIAS ORTOGONALES

Supóngase que se tiene una sola familia infinita de curvas. Considérese una
segunda familia compuesta de todas las curvas que intersectan a todas las curvas
de la familia dada en ángulos rectos. Cuando dos familias están relacionadas en la
forma que se ha descrito, se dice que las curvas de una familia son las
trayectorias ortogonales de la otra.

Si se diferencia la ecuación de la familia de curvas dadas y se elimina el

parámetro, se obtiene la ecuación diferencial de la familia dada, una ecuación que
de cualquiera de las curvas, en función de las coordenadas x, y
!
"
da la pendiente
de un punto de la curva.

Ahora bien, la pendiente de la trayectoria ortogonal en el punto (x, y) debe ser la

recíproca negativa de la pendiente de la curva dada, con objeto que se cumpla la
condición de perpendicularidad.

La ecuación diferencial de las trayectorias ortogonales se obtiene entonces

!
"
escribiendo , e igualando a la recíproca negativa del valor que tiene en la
ecuación diferencial de la familia dada.

Ejemplo: hallar las trayectorias ortogonales de las hipérbolas xy c.

! !
" "
Solución: la ecuación diferencial de la familia dada en , luego la ecuación
! "
" !
diferencial de las trayectorias ortogonales es: , la cual es una ecuación de
variables separables, cuya solución es , donde C es una constante.

3.8 CURVAS DETERMINADAS DE PROPIEDADES GEOMÉTRICAS

A menudo se tiene una familia de curvas caracterizada por una propiedad

geométrica indicada en función de las coordenadas de un punto en una de las
curvas y de la primera derivada de una de las coordenadas con respecto a la otra.
Tal enunciado será una ecuación diferencial de primer orden cuya solución
general representa la familia de curvas, cada una de las cuales posee la propiedad
dada. Si se da una condición adicional, tal como que la curva contiene un punto

[Link]
fijo dado, se puede usar dicha condición para determinar la constante en la
solución general.

al eje x tiene una longitud k.

Ejemplo: hallar la familia de curvas de tal manera que la proyección de la normal

en un punto P .x,y , luego para nuestro caso tenemos que DB mide k, luego
Solución: dada la figura 8, en ella observamos la curva, su tangente y su normal

MDÉ
!
"

Para P en cualquiera de los cuatro cuadrantes ´ ya sea u obtuso o agudo, la

longitud de DB es siempre " , o, . Luego la ecuación diferencial para la
! !
"
familia pedida es: ± :
!
"

La cual es una ecuación de variables separables cuya solución es:

± 2:

Donde C es una constante.

ECUACIONES LINEALES DE PRIMER ORDEN

Son ecuaciones diferenciales de la forma

a/ . af . É. donde a/ . ≠ 0
!
"

Observación:

Cuando g. 0 se dice que la ecuación lineal homogénea si g. ≠ 0 entonces al

ecuación diferencial es no homogénea.

Ejemplos:

dy
a 3 j9 . ln.
dx

[Link]
 a/ . 3

af . j9 D .

É. ln.

¿Una ecuación diferencial separable puede ser también una ecuación diferencial lineal?

ECUACIÓN DIFERENCIAL LINEAL NO HOMOGÉNEA

La forma estándar de la ecuación diferencial lineal es

L. ¾.
),
)-
(3)

Donde

L.
po ."
pm ."
y

g.x
¾.
D/ .

La solución de una ecuación diferencial lineal es

v R

Donde v es la solución general de la ecuación diferencial lineal homogénea asociada (3),

es decir, es la solución de )- L. 0
),

Y _ es una solución particular.

L. 0 es una ecuación diferencial separable, en consecuencia

),
Observemos que )-

dy
L.
dx

Entonces

1dy L.

[Link]
 ? L.
,
? !
por lo tanto

ln. + L. ln 8

En consecuencia

v. 9 ? 5I." ">N1 s

v. 9 5 ? I." "

Encontrar una solución particular. Utilizaremos el método de variación de parámetros.

Sabemos que v _ Queremos una _ que sea linealmente independiente de v.

Linealmente independiente _ no es múltiplo de v.

En consecuencia buscamos una _ que no sea múltiplo de v.

_. Ù. v.

Vamos a determinar . Como I es solución entonces tenemos que I satisface la

ecuación no homogénea.

´I L. _ ¾. (4)

Observemos que

´I . . . v. ´

´I . ´. v . ´v .

Reemplazando en (4) tenemos

´. v . v ´. L. . v. ¾.

Factorizando . se tiene ´. v. . [ ´v . L. v. ] ¾.

Recordemos que v es solución de la ecuación diferencial homogénea

L. 0

[Link]
En consecuencia

´v . L. v. 0

Lo cual implica que

´. v. ¾.

Despejando nos queda

¾.
´.
v.

Por el teorema fundamental del cálculo tenemos que ahora, integrando tenemos.

¾.
. +
v.

Por lo tanto

_. . v.

¾.
_. 9 5 ? _." +
"
9 5 ? _." "

¾.
_. 9 5 ? _." +
"
9 5 ? _." "

. v. _.

9 5 ? _." "
9 5 ? _." "
+ ¾. ? I." "

.
? Ç." q ? Ú.Û ÜÛ >v
q ? Ý.Û ÜÛ
Ecuación General

Ejemplo: ejercicio

Resolver la ecuación diferencial lineal paso a paso

3 .j9 ln
!
"
(1)

La forma estándar es:

[Link]
 $ j9 % ln
! / /
"
(2)

Encontramos v. la ecuación D.L.H. asociada es

1
j9 0
3

$ j9 %
! /
"
O equivalente

En consecuencia

1 1
j9
3

1 1
+ + j9
3

1
ln| | . cos j9 K
3

V j9
V Ëj

u Ëj + cos

u . Ëj j9 ln 8

Por lo tanto
/
v. 9 .5" Þßµ ">]q1 "

Ahora la solución particular es

/
I. . 9 .5" Þßµ ">]q1 "

Sustituimos en (2)

[Link]
 / 1 / 1
. .9 .5" Þßµ ">]q1 "
j9 . 9 .5" Þßµ ">]q1 "
ln| |
3 3

1 /
. . cos j9 cos 9 .5" Þßµ ">]q1 "
3
/ 1 / 1
´. 9 .5" Þßµ ">]q1 "
j9 . 9 .5" Þßµ ">]q1 "
ln| |
3 3

u(x)P . cos j9 cos 9 ’.5" Þßµ ">]q1 "

m
j9 . 9 ’.5" Þßµ ">]q1 "
m
/ /

´. 9 ’.5" Þßµ ">]q1 " Q ln| |

m
/

Ahora integrando tenemos:

1
ln| |
+ + 3
/
9 . Ëj j9

1
ln| |
+ 3
/
9 . Ëj j9

EL MÉTODO DEL FACTOR INTEGRANTE

• Método de factor integrante para resolver ecuaciones diferenciales

La ecuación diferencial lineal estándar es

´. L . ¾ .

Multipliquemos por 9 ? I." "

toda la ecuación diferencial.

Es decir:

9 ? I." "
´. 9 ? I." "
L. 9 ? I." "
¾.

[Link]
 ^9 ? I." "
` ^9 ? I." "
` 9 ? I." "

9 ? I." "
+ L. 9 ? I." "

9 ? I." "
L. 9 ? I." "

En consecuencia:

^9 ? I." "
` 9 ? I." "
¾.

Por el teorema fundamental del cálculo

+ ^9 ? I." "
` + 9 ? I." "
¾. 8

9 ? I." "
+ 9 ? I." "
¾. 8

Por lo tanto

? 9 ? I." " ¾. 8
.
9 ? I." "

El término 9 ? I." "

se conoce como el factor integrante para las ecuaciones lineales de
primera orden.

Ejercicio

1. Resolver por medio del método del factor integrante el P.V.I.

.0 4
!
"
(1)

Como la ecuación dada está en su forma estándar con L. 1 . Entonces el factor

integrante es:

9? "
9"

[Link]
Ahora multiplicamos la ecuación (1) por este término:

9" 9" 9"

.9 " 9"

Sea

V 9"
V 9

u 9" + 9"

u 9" 9"

Integrando tenemos:

+ .9 " + 9"

9" 9" 9" 8

9" 9" 8
9" 9"
1 8 9 5"

.0 0 1 8 9 5f 4

8 4 1

8 5

Luego:

1 59 5"

Donde

[Link]
 ´ 1 59 5"

1 59 5" 1 59 5"

ECUACIONES DIFERENCIALES EXACTA

Sea 2. , entonces el diferencial de es

à2 à2
à à

Si 2. , 8, Entonces

à2 à2
0
à à

Definición: una expresión diferencial <. , =. ,

en una región Å del plano
es una diferencial exacta

2. , . Una ecuación diferencial de la forma

, si corresponde a la diferencial exacta de alguna función

<. , =. , 0

Recibe el nombre de ecuación diferencial exacta, cuando el lado izquierdo corresponde a

una diferencial exacta.

Ejemplo: 3 2 0 Es una ecuación diferencial exacta. En este caso

<. , 3 y <. , 2

3 2
xy xy
x" x!

2. ,

3 2

En consecuencia es una ecuación diferencial exacta.

[Link]
Demostración (idea de demostración):

A partir del cálculo de varias variables se tiene que las derivadas cruzadas:

à 2 à 2
à à à à

Por otro lado sabemos que:

<. ,
xy
x"
Lo cual implica que:

à< à à2 à 2 à 2

à à à à à à à

à à2 à à=
.=. ,
à à à à

CRITERIO PARA UNA ECUACIÓN DIFERENCIAL EXACTA

Sean <. , y =. , con derivadas parciales continuas en una región Å rectangular

definida por D < < y < < , entonces, una condición necesarias y suficiente para
que <. , =. , sea una diferencial exacta es que

à< à=
à à

Ejemplo: la ecuación diferencial lineal

3 2 0
Es una ecuación diferencial exacta. En efecto:

<. , 3 y =. , 2

6 6
xâ xã
x! x"

à< à=
à à

Ejercicio:

[Link]
Determinemos si las siguientes ecuaciones diferenciales son exactas:

1. . Ëj j9 .1 0

2. .9 !
cos. .2 9 !
cos 2 0

Solución

1. <. , cos j9 ; = . , .1

à<
2
à

à=
à
à=
2 2
à

2. <. , 9 !
Ëj. ; =. , 2 9 !
Ëj 2

29 .cos. j9
xâ !
x!

29 cos. j9
xâ !
x!

29 cos. . j9
xã !
x"

Método de solución

Una vez determinada la exactitud de la ecuación diferencial, utilizamos el hecho de que

<. , =. ,
xy xy
x" x!
y

Para resolver la ecuación diferencial encontrar la primera integral. Sabemos que;

à2
<. ,
à

Integrando con respecto de tenemos que:

[Link]
 à2
+ + <. , å
à

Por el Teorema Fundamental del cálculo

2. , + <. , å.

Por otro lado

.? <. , å.
xy x
x! x!
Es decir que

à
=. , + <. , å´.
à

Despejando å´.

à
å´. =. , + <. ,
à

Por lo tanto

à
å. + ³=. , + <. , ´
à

Una vez encontrado :. s e reemplaza en la expresión que define 2, obteniendo la

solución o primer integral 2. , 8, donde 8 es una constante.

La solución o primera integral 2. , 8 donde 8 es una constante.

Repetir el método de solución pero empezando con =. ,

xy
x!

Ejemplo: resolver la ecuación diferencial.

. cos j9 .1 0

Como sabemos esta ecuación diferencial es exacta, entonces

à2
<. , cos j9
à

Siguiendo el método

[Link]
 2. , + <. , ; 2. , +. Ëj j9 :.

j9
2. , å.
2 2

:´. =. ,
xy "#!
x!
puesto que

=. , .1 Entonces

:´. Por la propiedad conmutativa

:´.

Por lo tanto

:. +
2

En consecuencia

j9
2. ,
2 2 2

Lo cual implica que la familia de curvas solución está dada por :

. 1 j9
2. , 8
2

Ejercicio

Resolver $æçè èéêë .ìè % íì .çìæçè ì îïð ìè çè íè ñ

à2
9 !
cos.
à

2. , +.9 !
Ëj . å.

2. , +9 !
+ Ëj . å.

[Link]
 2. , 9 !
j9 . :.

2. , 9 !
j9 . :.

à2
2 9 !
cos. :`. 2 9 !
cos 2
à

:`. 2

:. +2

:.

En consecuencia 2. , 9 !
j9 .

Lo cual implica que la familia de curvas solución está dada por

2. , 9 !
j9 . , 8

.0 2

Entonces 2 8 8 4

Luego

2. , 9 !
j9 . , 4

Solución del problema de valor inicial

En algunas ocasiones ciertas ecuaciones diferenciales se pueden reducir a ecuaciones

diferenciales exactas (si multiplicamos por un factor integrante las llevamos a exactas).

Ejemplo: la ecuación diferencial

.2 3 20 0

No es exacta en efecto

<. , , = . , 2 3 20

[Link]
 à< à=
4
à à

à< à=
≠
à à

Por lo tanto no es exacta, sin embargo la ecuación diferencial es reducible a exacta.

En efecto, multiplicando la Ecuación Diferencial por tenemos:

[ .2 3 20 ] 0
Ž .2 3 20 0

En consecuencia

ó. ,
< Ž

ó. ,
= .2 3 20

à<ó
4
à

à=ó
4
à

à<ó à=ó
à à

Solución:

1. .ìèç îïð ì ëæô ì íì è .õ ìç íè ñ

Como sabemos esta ecuación diferencial es exacta, entonces

=. , .1
xy
x!
siguiente el método

2. , ? .1 ℎ. entonces

[Link]
 2. , + + ℎ.

2. , ℎ.
2 2
à2
ℎ´. cos j9
à

j9

cos

ℎ. + cos j9

ℎ. +

ℎ.
2

j9
ℎ.
2

j9
2. ,
2 2 2

En consecuencia la familia de curvas solución está dada por

j9
2. , 8
2 2 2

2. .æçè ö îïð.ìè íì .çìæçè ì éêë ìè çè íè ñ

Como sabemos que esta ecuación diferencial es exacta, entonces

2 9 Ëj 2
xy !
x!
Integrando tenemos

2. , +.2 9 !
cos 2 ℎ .

2. , 9 !
j9 . ℎ .

[Link]
 à2
9 !
Ëj . ℎ´. 9 !
Ëj .
à

ℎ´. 0

ℎ. 8

En consecuencia la familia de curva solución está dada por

2. , 9 !
j9 . 8

CRITERIO PARA UNA ECUACIÓN DIFERENCIAL EXACTA

Método de solución

à2
=. ,
à

Integrando con respecto a tenemos

à2
+ + =. ,
à

2. , + =. , :.

Por otro lado

à2 à
+ =. , :.
à à

à2 à
+ =. , :´.
à à

Despejando :´. tenemos

à2 à
:´. + =. ,
à à
Por lo tanto

à2 à
:. +³ + =. , ´
à à

[Link]
ECUACIONES DIFERENCIALES REDUCIBLES A EXACTAS

Supongamos que la ecuación diferencial

<. , = . , 0

No es exacta. Entonces existe en algunos casos un factor integrante.

. , que convierte o transforma la ecuación anterior en una Ecuación Diferencial

exacta.

[< = ]0 o equivalente

< = 0

Dicha ecuación diferencial es exacta si y sólo sí

. < ! . = "

. < . =
x x
x! x"
por la derivada del producto

!= <! "= ="

^<! =" ` "= != (1) Buscamos funciones que dependan de una sola variable

Busquemos soluciones .

Entonces Si depende de , entonces

à
u.g 0
à

Reemplazando en (1) tenemos

^<! =" ` "=

O equivalentemente

<! ="

"
="!

[Link]
 depende de , por lo tanto su derivada también depende de .
âø 5ãù
En el caso de que ÷ ã
solo depende de , es decir ÷.

Entonces la Ecuación Diferencial anterior se escribe como

" ÷.

÷.
'
"
o Solución total. Ecuación Diferencial separable

Observemos que la ecuación diferencial es separable

1
÷.

En consecuencia

K | | + ÷.

La constante no se la asume.

Luego
. 9 ? ú." "

à< à=
Donde

à à
÷.
=

Encuentre el factor integrante cuando .

Entonces busquemos soluciones para .

0
x'
Entonces depende de ; " x"

Reemplazando en (1) tenemos ¥<! =" ¦ "= !<

¥<! =" ¦ !< o equivalente

^<! =" `
!
<

[Link]
 ^<! =" `
÷
<
=" <!
÷
<

Solo depende de . Es decir ÷ ; entonces la ecuación diferencial se escribe como:

! ÷.

Observemos que la ecuación diferencial es separable.

1
÷. !

En consecuencia

ln| | + ÷. !

! 9 ? ú.! ø

Donde

à< à=
à à
÷.
<

Ejercicio

Resolver la Ecuación Diferencial

.2 3 20 0

à= à<
à à 4 3 3
÷.
<

. 9 9 9 •• !
?! ! •• ! ’

[Link]
 . .2 3 20 0
Ž .2 3 6
20 Ž
0

4 4
xâ xã
x! x"
;

Puesto que el factor integrante es entonces:

Ž .2 3 6
20 Ž
0

<. , Ž Ž
=. , 2 3 6
20 Ž

à2
<. , Ž
à

Entonces

2. , + Ž
:.

Ž
2. , :.
2

à2 à Ž
:.
à à 2

à2
2. :´. 2. 3 6
20 Ž

à

:´. 3 6 20 Ž
:. ?.3 6
20 Ž

:. 4
!© 6

En consecuencia
Ž 7
2. , 4 6
2 2

Por lo tanto la familia monoparamétrica de soluciones está dada por

[Link]
 Ž 7
4 6
8
2 2

ECUACIONES DIFERENCIALES HOMOGÉNEAS

Decimos que la función 2. , es homogénea de orden si 2.M , M M r 2. ,

Ejemplo:

a). 2. , es una función homogénea de segundo orden. En efecto,

2.M , M .M .M M . M 2. ,

Ejercicio

Determine si las siguientes funciones son homogéneas

a) 2. ,
"
!

b) 2. , Ž

c) 2. ,
"’
">!

ECUACIÓN DIFERENCIAL HOMOGÉNEA

Decimos que la Ecuación Diferencial

<. , =. , 0

Es una ecuación diferencial homogénea si

<. , =. ,

Son funciones homogéneas del mismo orden; se reducen a separables.

ECUACIONES HOMOGÉNEAS

<. , =. , 0

[Link]
Demuestra que una educación homogénea es reducible y separable.

Es homogénea sí y solo sí < y = son funciones homogéneas de primer orden.

Observación: la ecuación diferencial homogénea

En efecto, supongamos que < y = son funciones homogéneas del mismo orden

Entonces

<.M , M M 1 <. , y =.M , M M 1 =. ,

Sea M entonces < y = se reescriben como

/
"

< $1, % 51
<. ,

= $1, % 51
=. ,

En consecuencia

<. , 1
<.1, û

=. , 1
=.1, û

Reemplazando en la Ecuación Diferencial homogénea

1
< $1, % 1
=^1, û ` 0

1
Factorizando tenemos

1
¥<^1, û ` =^1, û ` ¦ 0

Como 1
≠ 0, entonces

<^1, û ` =^1, û ` 0 (1)

Realizamos un cambio de variables

!
"
Sea ; es decir,

[Link]
Observemos que

Reemplazando en (1) tenemos

<.1, =.1, [ ] 0

<.1, = .1, =.1, 0

[<.1, =.1, ] =.1, 0

=.1, [< .1, = .1, ]

=.1,
<.1, =.1,

[<.1, ] =.1, ] 1
=[1, ]

É. . ℎ.

Donde

É. ℎ.
[â./,' >ã./,' '] /
ã./,' ü
(2) y

Por lo tanto la ecuación diferencial (2) es separable.

M
/
!
Ejercicio: repita el proceso definiendo , y reduzca la Ecuación Diferencial
homogénea a una Ecuación Diferencial separable.

Ejercicio: determine si la siguiente ecuación diferencial es homogénea, en caso positivo

resolverla por uno de los métodos.

a. .ìç èç íì .ìç ìè íè ñ

<. , Entonces <.M , M .M .M M . M <

=. , . M M, M . M M M . M =. ,

Luego, sea entonces

Entonces

[Link]
 . . [ ] 0

. . . 0

. .

. .1

.1 .1

.1 .1

.1
1

1
1

Integrando a ambos lados tenemos

1
+ +
1

2
+ 1 ln ln 8
1

2 ln|1 | ln ln 8
!
"
Dado que , entonces

2 ln $1 % ln.

ln $1 % ln.

ln ³ $1 % ´

b. .ìç çèç ìè è. õ õ
íì
íè

[Link]
 . 2

. 2 0

<. , 2

<.M , M .M 2.M M . 2 M <. ,

=. , =.M , M M M M M =. ,
!
"
Sea ;

. 2 0

. 2 [ ] 0

¥.1 2 [ ]¦ 0

≠ 0 Entonces

.1 2 0

.1 2 0

.1 3 0

.1 3

1 3

+ +
1 3

1 3

6
1 1
u + ln|1 3 |
6 6

[Link]
 1
ln| | ln|1 3 | ln 8
6
/
ln ý .1 3 7ý ln 8

/
ý .1 3 7ý 8

| 7 .1
3 | 87

|1 3 | 57 7

3 57 7
1

3 57 7
1
7 7
1
3 3
7 7
1
$ %
3 3
7 7
1
3 3

7 7 1
±H
3 3

7
1 3
57

1 3$ %
1 6

[Link]
ECUACIÓN DE BERNOULLI

L. ¾. 1

Método de solución
/51
Para resolverla se hace el cambio de variable , con este cambio la ecuación
diferencial de Bernoulli se reduce a una ecuación diferencial lineal, en efecto.
1
Dividiendo entre la ecuación de Bernoulli se obtiene

51 !
"
L. /51
¾. (3)

Por otro lado

.1 51

O equivalentemente

1
1

1

Reemplazando en (3) tenemos

1
L. þ ¾.
1

o equivalentemente

.1 L. þ .1 ¾.

Lo cual es una ecuación diferencial lineal

Ejemplo:

Resolver la ecuación diferencial de Bernoulli

[Link]
Solución:

Esta ecuación diferencial se reduce a

En este caso 2, por lo tanto /5 5/

Lo cual implica que

Dividiendo la ecuación entre tenemos

5
1 5/

Reemplazando se tiene

/
" "
(3)

Encontramos el factor integrante.

5/
9? " "
9 5 ••|"| 5/

Multiplicando (3) por el factor integrante tenemos

1
5/
1

. 5/
1 ; + . 5/
+

5/
8

Luego

. 8

[Link]
 8

Como

5/

Entonces
5/

Por lo tanto

1
.
.8

ECUACIONES HOMOGÉNEAS

<. , =. , 0 Es homogénea sí y solo sí M y N son funciones homogéneas

de primer orden.

Supongamos que M y N son funciones homogéneas de mismo orden.

Entonces:

<.M , M M 1 <. ,

=. M, M M 1 =. ,

Sea M
/
!
entonces M y N se escriben como

< ,1 51
<. ,

= ,1 51
=. ,
En consecuencia
M(x,y) = 1 < $! , 1%
"

= $ , 1%
1 "
!
N(x,y) =

Reemplazando en la ecuación diferencial homogénea tenemos

[Link]
 1
< ,1 1
= ,1 0

1
³< ,1 = ,1 ´ 0

Puesto que 1
≠ 0 entonces tenemos

< $! , 1% = $ , 1% 0 (1)
" "
!

Ahora realizamos un cambio de variable. Sea , y entonces ,

sustituyendo en (1).

<. , 1 [ ] =. , 1 0

<. , 1 <. , 1 =. , 1 0

[ <. , 1 =. , 1 ] <. , 1

1 [ <. , 1 =. , 1 ]
<. , 1

Luego

É. ℎ.
'
!
esta Ecuación Diferencial es separable.

ECUACIONES DIFERENCIALES REDUCIBLES A SEPARABLES

Una Ecuación Diferencial de la forma " 2.C G 8 , donde A, B, C son

!

constantes, siempre se puede reducir a una Ecuación Diferencia separable.

Demostración;

En efecto, haciendo el cambio de variable C G 8

C G
' !
! "
Obtenemos

Reemplazando en (A) se obtiene

$ C% 2.
/ '
r "
o equivalentemente

[Link]
 2 C; es decir d
'
"

É.
'
"
Ecuación Diferencial separable

Ejemplo: resolver el problema de valores iniciales PVI.

. 2 7 .0 0
!

En este caso 2 En consecuencia

Luego

Reemplazando en el PVI

2 7

.0 2.0 0 0

[Link]
 CAPÍTULO II
EXISTENCIA Y UNIDAD DE LAS SOLUCIONES
PROBLEMAS DE VALORES INICIALES

4. ECUACIÓN DIFERENCIAL LINEAL DE PRIMER ORDEN

orden L definido por:

En la sección anterior hemos visto que si tenemos el operador lineal de primer

K.y ´ L . , dondeP ∈ C .I

Entonces la ecuación L .y Q.x donde Q ∈ C .I , tiene una única solución f, que

verifica la condición 2.D , para a un elemento de I, y b un real dado, esta
función viene dada por: 2 .x 9 5F ." 9 5F ." ?p ¾ .M 9 F .J M
"

Donde
"
C . + L .M M
p

Gráficamente esto significa que para la ecuación K. ¾, y para cada punto

. f , f ∈ u Å, existe 2 que es solución de la ecuación y satisface la condición
2 . f f , o sea se llena la franja dada por u Å.

[Link]
 SOLUCIÓN DE ECUACIONES DIFERENCIALES DE PRIMER,
ORDEN

1° INTRODUCCIÓN:
Con el propósito de que los estudiantes inscritos en los cursos de ecuaciones
diferenciales en la facultad de ingeniería tengan una guía como referencia, en lo
relativo a la solución de las mismas ecuaciones diferenciales, se ha diseñado este
material didáctico que les servirá de gran ayuda a fin de conceptualizar
adecuadamente sobre el proceso.
El documento presenta la solución analítica y numérica de diferentes tipos de
ecuaciones diferenciales de primero, segundo y tercer orden.
El proceso de solución de la ecuación diferencial inicia con la determinación de la
ecuación solución (función) y de sus derivadas intermedias, halladas por cualquiera de
los métodos analíticos posibles vistos en el curso. Una vez encontrada la ecuación de
la función y de sus derivadas intermedias se procede a dibujar sus respectivos lugares
geométricos utilizando el programa de matlab.
El proceso de solución finaliza con la determinación de los diferentes lugares
geométricos de la función y de sus derivadas intermedias, hallados por cualquiera de
los métodos numéricos mediante la simulación de los programas incluidos en los
toolboxes de matlab.
Finalmente el documento presenta las gráficas de los lugares geométricos obtenidos
por los dos métodos, para que el estudiante se pueda dar una idea de la precisión del
método numérico utilizado

2° ECUACIONES DIFERENCIALES DE PRIMER ORDEN

2.1 ECUACIÓN DIFERENCIAL DE PRIMER ORDEN EN DONDE LA DERIVADA

ES FUNCIÓN SOLAMENTE DE LA VARIABLE DEPENDIENTE

íè
í
Desarrollar la ecuación diferencial = - 0.081093 (y - 50), para t = 0, y = 80.
2.1.1 DESARROLLO ANALÍTICO:
íè
Separando variables en la ecuación diferencial, tendremos: .è5 ñ = -
0.081093 dt
Integrando a ambos lados de la ecuación diferencial, resulta:
= ? ñ. ñ õñ dt ; Ln(y – 50) = - 0.081093 t + kint , simplificando,
íè
?.è5 ñ
quedará:
y = 50 + Ko e – 0.081093 t

[Link]
La cual corresponde a la solución general de la ecuación diferencial.
Reemplazando las condiciones iniciales, podremos determinar la constante
de integración.
80 = 50 + Ko e – 0.081093 (0) , luego Ko es igual a 30 y la solución específica o
particular quedará:

y = 50 + 30 e – 0.081093 t

Algunos valores de prueba:

Para t = 10, y(10) = 50 + 30 e – 0.081093 (10) = 63.33°

íè
í
t = 10 = - 0.081093 (63.33 - 50) = - 1.0809

RESULTADOS DE ECUACIONES POR EL MÉTODO ANALÍTICO:

íè
Ecuación diferencial a resolver: í = - 0.081093 (y - 50)
Ecuación solución (Función): y (t) = 50 + 30 e – 0.081093 t

[Link] GRÁFICAS DE LA ECUACIÓN DIFERENCIAL Y DE LA ECUACIÓN

SOLUCIÓN (FUNCIÓN)

A partir de la ecuación diferencial (ecuación de la primera derivada) y de la ecuación para

la función (solución de la ecuación diferencial), podremos dibujar los lugares geométricos
correspondientes, utilizando el software de Matlab, mediante el programa siguiente.
(dibuejemplo4)
%El programa dibuja la función y = 50 + 30.*e- 0.081093.*t
% y su respectiva derivada o ecuación diferencial y´= -0.081093(y – 50)
t = 0 : 0.01 : 15;
Z = 50 + 30.*exp(-0.081093.*t);
figure(1), plot(t,Z), grid on
title('Ecuación solución de la E.D.'),
ylabel('y(t)=50+30exp(-0.081093 t)'), xlabel('Tiempo t'), pause
W = -0.081093.*(Z-50);
figure(2),plot(t,W,'r'), grid on
title('Ecuación diferencial a resolver')
ylabel('y´(t)= -0.081093(y-50)'), xlabel('Tiempo t'), pause
[t.',Z.',W']

[Link]
 Ecuación diferencial a resolver
-0.6

-0.8
X: 10
-1 Y: -1.081

-1.2
y´(t)= -0.081093(y-50)

-1.4

-1.6

-1.8

-2

-2.2

-2.4

-2.6
0 5 10 15
Tiempo t

Ecuación solución de la E.D.

80

75
y(t)=50+30exp(-0.081093 t)

70

X: 10
65 Y: 63.33

60

55
0 5 10 15
Tiempo t

[Link]
2.1.2 DESARROLLO NUMÉRICO:
A continuación se encuentra el programa de matlab por medio del cual
encontraremos la solución numérica de la ecuación diferencial.
(damarquez.m)
%El programa resuelve la ecuación diferencial dy/dt= -0.081093*(y-50)
%paray(0)= 80
function davila
tspan = 0 : 0.01: 15;
y0 =[80];
[T,Y] = ode23tb(@f,tspan,y0);
[T,Y]
Z= -0.081093.*(Y-50);
clf, grid on, hold on
figure(1),plot(T,Z,'r'), title('Ecuación diferencial a resolver')
legend('dy/dt'),xlabel('Tiempo t'), grid on
ylabel('primera derivada de la función dy/dt'), pause
figure(2), plot(T,Y,'-')%Gráfica de y
title('Solución de la E.D. Función y(t)'); %Título del gráfico
legend('y');%Leyenda del gráfico
xlabel('Tiempo t');%Nombre del eje x
ylabel('Función y');%Nombre del eje y
grid on
hold on
function dy = f(t,y)
dy = -0.081093*(y(1)- 50);

ALGUNOS RESULTADOS DE LA SIMULACIÓN

Ecuación diferencial a resolver
-0.6
dy/dt
-0.8
X: 10
-1 Y: -1.081
primera derivada de la función dy/dt

-1.2

-1.4

-1.6

-1.8

-2

-2.2

-2.4

-2.6
0 5 10 15
Tiempo t

[Link]
 Solución de la E.D. Función y(t)
80
y

75

70
Función y

X: 10
65 Y: 63.33

60

55
0 5 10 15
Tiempo t

2.2 ECUACIÓN DIFERENCIAL DE PRIMER ORDEN EN DONDE LA DERIVADA

ES FUNCIÓN SOLAMENTE DE LA VARIABLE INDEPENDIENTE

íè
í
Desarrollar la ecuación diferencial = Cos(5t), para t = 0, y = 2.

2.2.1 DESARROLLO ANALÍTICO:

Separando variables en la ecuación diferencial, tendremos: dy = Cos(5t) dt
Integrando a ambos lados de/ la ecuación diferencial, resulta:
? íè = ? êë. dt y = 6 Sen(5t) + k ;
int

La cual corresponde a la solución general de la ecuación diferencial.

Reemplazando las condiciones iniciales, podremos determinar la constante
de integración.
õ
2 = Sen(5(0)) + kint, luego Kint es igual a 2 y la solución específica o
particular quedará:
õ
y= Sen(5t) + 2

[Link]
Valores de prueba:
/
6
para t = 1.5 , y(1.5) = y = Sen(5(1.5)) + 2 = 2.1875
íè
í
t = 1.5 = Cos(5(1.5)) = 0.34663

íè
RESULTADOS DE ECUACIONES POR EL MÉTODO ANALÍTICO:
í
õ
Ecuación diferencial a resolver: = Cos(5t),
Ecuación solución (Función): y (t) = Sen(5t) + 2

[Link] GRÁFICAS DE LA ECUACIÓN DIFERENCIAL Y DE LA ECUACIÓN

SOLUCIÓN (FUNCIÓN)

A partir de la ecuación diferencial (ecuación de la primera derivada) y de la ecuación para

la función (solución de la ecuación diferencial), podremos dibujar los lugares geométricos
correspondientes, utilizando el software de Matlab, mediante el programa siguiente.
(dibuejemplo5)
%El programa dibuja la función y = (1/5).*Sen(5.*t)+ 2
% y su respectiva derivada o ecuación diferencial y´= Cos(5.*t)
t = 0 : 0.01 : 3;
W = cos(5.*t);
figure(1),plot(t,W,'r'), grid on
title('y´(t)= Cos(5.*t)'),ylabel('Ecuación diferencial'),pause
Z = (1/5).*sin(5.*t)+ 2 ;
figure(2), plot(t,Z), grid on
title('y(t)= (1/5).*Sin(5.*t)+ 2'), ylabel('Función,Solución de
E.D.'),pause
[t.',Z.',W']

[Link]
 y´(t)= Cos(5.*t)
1

0.8

0.6
X: 1.5
Y: 0.3466
0.4
Ecuación diferencial

0.2

-0.2

-0.4

-0.6

-0.8

-1
0 0.5 1 1.5 2 2.5 3

y(t)= (1/5).*Sin(5.*t)+ 2
2.25

2.2

X: 1.5
2.15 Y: 2.188
Función,Solución de E.D.

2.1

2.05

1.95

1.9

1.85

1.8
0 0.5 1 1.5 2 2.5 3

[Link]
2.2.2 DESARROLLO NUMÉRICO:
A continuación se encuentra el programa de matlab por medio del cual
encontraremos la solución numérica de la ecuación diferencial.
(davila.m)
%El programa resuelve la ecuación diferencial y´= Cos(5t),para y(0) =
2.0
function davila
tspan = 0 : 0.01: 3;
y0 =[2];
[T,Y] = ode23tb(@f,tspan,y0);
t = 0 : 0.01 : 3;
Z = [cos(5.*t)];
[T ,Y,Z.']
clf
plot(T,Z),title('Derivada de la función'),xlabel('Tiempo t'),
ylabel('derivada y´'),grid on, pause;
plot(T,Y(:,1),'-')%Gráfica de y
title('Solución de la ecuación diferencial'); %Título del gráfico
legend('y');%Leyenda del gráfico
xlabel('Tiempo t');%Nombre del eje x
ylabel('Función y');%Nombre del eje y
grid on
hold on
function y1 = f(t,y);
y1 = [cos(5.*t)];

Derivada de la función
1

0.8

0.6
X: 1.5
Y: 0.3466
0.4

0.2
derivada y´

-0.2

-0.4

-0.6

-0.8

-1
0 0.5 1 1.5 2 2.5 3
Tiempo t

[Link]
 Solución de la ecuación diferencial
2.25
y
2.2

X: 1.5
2.15 Y: 2.187

2.1
Función y

2.05

1.95

1.9

1.85

1.8
0 0.5 1 1.5 2 2.5 3
Tiempo t
2.3
ECUACIÓN DIFERENCIAL EN DONDE LA DERIVADA ES FUNCIÓN DE LA DE
LA VARIABLE DEPENDIENTE Y DE LA INDEPENDIENTE.

íè
í
Desarrollar la ecuación diferencial = 20 Cos ( 60 t ) – 5y, para t = 0, y = 0.

2.3.1 DESARROLLO ANALÍTICO:

Se puede observar que la ecuación diferencial no es de variables separables,
pero se puede clasificar como una ecuación diferencial lineal de primer
orden considerando a y como función.
La solución de la ecuación diferencial estará conformada por la solución de
la homogénea correspondiente yh(t) y la solución particular. y(t) = yh(t) +
yp(t)
Solución de la homogénea correspondiente:
yh(t) = K e – 5 t

Solución particular de la ecuación diferencial a resolver:

Se asume yp(t) = M Cos(60 t) + N Sen(60 t)
Determinando sus derivada: è´ . = - 60 M Sen ( 60 t ) + 60 N Cos( 60 t )

[Link]
Reemplazando en la ecuación diferencial, tendremos:
- 60 M Sen( 60 t ) + 60 N Cos( 60 t ) + 5 M Cos(60 t) + 5 N Sen(60 t) = 20 Cos (
60 t )

Simplificando:
Sen( 60 t ) (-60M+5N) + Cos( 60 t ) (60N+5M) = 0 Sen( 60 t) + 20 Cos ( 60 t )
Comparando miembro a miembro en la ecuación presentada, resultan dos
ecuaciones con dos incógnitas: -60M+5N = 0 , 60N+5M = 20
Desarrollándolas simultáneamente, tendremos: M = õ = 0.027586 , N = õ
= 0.331034
Por lo tanto, la solución particular quedará:
õ õ
yp(t) = Cos(60 t) + Sen(60 t) ,

√õ
yp(t) = Cos(60 t – 85.23°)

La solución general de la ecuación diferencial quedará.

√õ
y(t) = Cos(60 t – 85.23°) + K e – 5 t

Reemplazando las condiciones iniciales, t = 0, y = 0.

La constante quedará: K = -0.0276, y la solución específica quedará
expresada por:

√õ
y(t) = Cos(60 t – 85.23°) – 0.0276 e – 5 t

y (t) = 0.332182 Cos(60 t – 1.4875r) – 0.0276 e – 5 t

Valores de prueba:
Para t = 0 , y(0) = 0, è´.ñ = 20 Cos(60(0)) – 5(0) = 20

√õ
Para t = 0.1, y(0.1) = Cos(60 (0.1) – 85.23°) – 0.0276 e – 5 (0.1) = - 0.0826984
è´.ñ.õ = 20 Cos (60(0.1)) – 5(-0.0827472) = 19.6171

√õ
Para t = 0.3, y(0.3) = Cos(60 (0.3) – 85.23°) – 0.0276 e – 5 (0.3) = - 0.236508

[Link]
è´.ñ. = 20 Cos (60(0.3)) – 5(-0.236544) = 14.3890
RESULTADOS DE ECUACIONES POR EL MÉTODO ANALÍTICO:

íè
í
Ecuación diferencial a resolver: = 20 Cos ( 60 t ) – 5y

√õ
Ecuación solución (Función): y(t) = Cos(60 t – 85.23°) – 0.0276 e – 5 t

[Link] GRÁFICAS DE LA ECUACIÓN DIFERENCIAL Y DE LA ECUACIÓN

SOLUCIÓN (FUNCIÓN)

A partir de la ecuación diferencial (ecuación de la primera derivada) y de la ecuación para

la función (solución de la ecuación diferencial), podremos dibujar los lugares geométricos
correspondientes, utilizando el software de Matlab, mediante el programa siguiente.
(dibuejemplo6)
%El programa dibuja la función
%y = 0.338122.*cos(60.*t-1.4875)-0.0276.*exp(-5.*t) y su respectiva
%derivada o ecuación diferencial y´= 20.*cos(60.*t)-5.*(y)
t = 0 : 0.001 : 0.5;
W = 20.*cos(60.*t)-5.*(0.338122.*cos(60.*t-1.4875)-0.0276.*exp(-5.*t));
figure(1),plot(t,W,'r'), grid on
title('y´(t)= 20.*cos(60.*t)-5.*(y)')
ylabel('Ecuación diferencial'),pause
Z = (0.338122.*cos(60.*t-1.4875)-0.0276.*exp(-5.*t)) ;
figure(2), plot(t,Z), grid on
title('y(t)= 0.338122.*cos(60.*t-1.4875)-0.0276.*exp(-5.*t)')
ylabel('Función,Solución de E.D.'),pause
[t.',Z.',W']

y´(t)= 20.*cos(60.*t)-5.*(y)
25

20
X: 0.1
15 Y: 19.62

10
Ecuación diferencial

-5

-10

-15

-20
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5

[Link]
 y(t)= 0.338122.*cos(60.*t-1.4875)-0.0276.*exp(-5.*t)
0.4

0.3

0.2
Función,Solución de E.D.

0.1

0 X: 0.1
Y: -0.08388

-0.1

-0.2

-0.3

-0.4
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5

2.3.2 DESARROLLO NUMÉRICO:

A continuación se encuentra el programa de matlab por medio del cual
encontraremos la solución numérica de la ecuación diferencial.
(otrorod.m)
%El programa resuelve la ecuación diferencial y´= 20 Cos(60t)-5y
% para y(o)=0
function Rodolfo
tspan = 0 : 0.001: 0.5;
y0 = [0];
[T,Y,]= ode23tb(@f,tspan,y0);
Z= 20.*cos(60.*T)- (5.*Y);
[T,Y,Z],plot(T,Z,'r'),title('Primera derivada o ecuación diferencial');
grid on, xlabel('Tiempo t'),ylabel('dy/dt'),pause;
plot(T,Y,'-'),title('solucion de la ecuacion diferencial');
grid on, xlabel('Tiempo t'),ylabel('Funcion y');
[Y(:,1)]
function y1 = f(t,y)
y1 = [20.*cos(60.*t)-5.* y];
[t ,y, y1]

[Link]
 Primera derivada o ecuación diferencial
25

20
X: 0.1
Y: 19.62
15

10

5
dy/dt

-5

-10

-15

-20
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Tiempo t

solucion de la ecuacion diferencial

0.4

0.3

0.2

0.1
Funcion y

0 X: 0.1
Y: -0.08278

-0.1

-0.2

-0.3

-0.4
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
Tiempo t

[Link]
5. ECUACIÓN DIFERENCIAL LINEAL DE PRIMER ORDEN N

los teoremas que hallan la solución a las ecuaciones de orden n; hasta el

La teoría general de ecuaciones diferenciales lineales en realidad comienza con

momento únicamente lo hemos hecho para orden uno, para el caso de orden n la
existencia y unicidad de la solución es mucho más compleja de probar, debido a
ello me limitaré a enunciarlo y al final en el apéndice segundo lo demostraré, para

Teorema 5: se L un operador diferencial lineal normal de orden n, definido en un

el lector interesado.

intervalo I, sea u f un punto arbitrario de I; entonces si Df , D/ , … , D15/ , son reales

arbitrarios.

Luego la ecuación K . ℎ con ℎ ∈ 8 .u

Y, . f Df, ´ . f D/ , … , .15/
. f D15/ tiene una única solución.

ecuación diferencial lineal normal de orden n. Ahora problemas que la dimensión

Luego este teorema caracteriza la existencia y unicidad de la solución a una

del espacio solución a dicha ecuación es igual al orden de la ecuación.

Teorema 6: dada la ecuación diferencial lineal normal homogénea de orden n

D1 . .1
D15/ . . 1 … Df . 0 (1)

Definida en I. entonces el espacio solución de esta ecuación es un sub-espacio n-

dimensional de 8 1 .u .

Demostración: para un punto fijo f de u, y las n-etuplas.

.1, 0, 0, … , 0

.0, 1, 0, … , 0

.0, 0, 1, … , 0
.
.
.

.0, 0, 0, … , 1

Por el teorema 5, existen / . , . , … , 1 . , que son soluciones de la

ecuación (1) y satisfacen las condiciones dadas, entonces

[Link]
 /. f 1 ´/ . 0 . 0
.15/
f ... / f

. f 0 ´ . f 1 ... /
.15/ . f 0
.
.
.

1. f 0 ´1 . 0 . 1
.15/
f ... 1 f

O sea . . f , ´. f ,... .15/

. f , para t 1, 2, … . ,

Los cuales son vectores de Å 1 y forman una base de éste.

Probemos que estas soluciones son una base para el espacio solución de (1).

Sean D/ , D , … D1 n-reales tales que

D/ / . D . … D1 1 . 0 para en t.

Luego hallando las 1 derivadas, tenemos

D/ / . D . … D1 1 . 0

D/ ´/ . D ´ . … D1 ´1 . 0
.
.
.
D/ / .15/ . D .15/
. … D1 1
.15/
. 0

Haciendo f, obtenemos

D/ / . f D . f … D1 1 . f 0

D/ ´/ . f D ´ . f … D1 ´1 . f 0
.
.
.
D/ / .15/ . f D .15/
. f … D1 1
.15/
. f 0

Entonces D/ D … . D1 0 , lo cual implica que / , , ... 1 son linealmente

.u , lo único que falta para poder concluir que
. /. , . , … , 1 .
independientes en
es una base del espacio solución de la ecuación (1), es que

[Link]
cualquier solución de dicha ecuación se puede expresar como combinación lineal
de las n-funciones dadas.

Sea y una solución arbitraria de la ecuación (1) y supongamos que . f

/ , ´. f , … , .15/ . f 1

Por el teorema 5, sabemos que esta solución es única.

Definimos ℎ . / /. . … 1 1. para x en I.

Luego ℎ es solución de (1) y satisface las condiciones dadas, pero el teorema 5

/ / … 1 1 , de donde
podemos concluir que / , , … 1 generan el espacio solución de la ecuación (1),
asegura que dicha solución es única, luego

luego este conjunto de funciones es una base a dicho espacio.

Al observar la demostración del teorema 6, vemos que establece un método para

comprobar la independencia lineal de funciones, el cual lo podemos caracterizar
así:

Corolario 1. Sea / , , … , 1 , n-funciones de 8 .u , cada una de las cuales posee

derivadas hasta del orden . 1 inclusive en todo I, y supongamos que existe un
punto f en I tal que los vectores . . f , ´ . f , … , .15/ . f , para t
1, 2, … . Son lineamientos independientes en Å 1 . Entonces / , , … , 1 son
linealmente independientes en .u .

6. EL WRONSKIANO

.u , pero esto también lo podemos caracterizar a través de la idea del

Acabamos de mostrar cuando n-funciones son linealmente independientes en

Wronskiano de . / , , … 1 para n-funciones de 8 .15/ .u , el cual notaremos así:

. / . , . ,…, 1.

Y está definido por el siguiente determinante.

/ . . … 1.
´/ . ´ . … ´1 .
. /. , . ,…, .
1. 0
.
.
/
.15/
. .15/
. … 1
.15/
.

[Link]
En el apéndice uno, demostraré que si D/ , D , … , D1 son n-vectores linealmente
independientes entonces

D// D/ … D/1
D / D … D 1
.
[.D/ , D , … , D1 0
.
.
D1/ D1 … D11

En base a este enunciado tenemos:

/, , … , 1 , n-funciones de 8 15/ .u tal que

. /. , . , … , 1 . ≠ 0 para en u, entonces / , , … , 1 son linealmente
Teorema 7: dadas las funciones

independientes en 8.u .

[9 " , 9 5" ] 9" 9 5"

P Q 2
9" 9 5"
Ejemplo: hallemos el

Luego, 9 " , y 9 5" son linealmente independientes en .Å

Ejemplo: dadas las funciones y | | son linealmente independientes en 8.Å ,

ya que si

/ | | 0 para 1 , y, 1, tenemos que / 0

[ ,| | ] P Q 0 , si
3 3
Pero es mayor o igual a cero.

[ ,| | ] P Q 0 , si
3 3
es menor de cero.

deducir dependencia lineal de un conjunto de funciones en 8.u del hecho de que

Lo cual demuestra que el recíproco al teorema 7 es falso, o sea no podemos

su Wronskiano sea idénticamente cero en u, pero si condicionamos estas

funciones a ser soluciones de una ecuación diferencial lineal normal homogénea,
podemos concluir:

Teorema 8: sean / , , … 1 , n-funciones de 8 1 .u que son soluciones de la

ecuación diferencial lineal homogénea normal de orden n.

D1 ." ! .0
D15/ ." ! .0 m
⋯ Df ." !
0

. /. , . , … , 1 . 0 sobre u. Entonces /, ,…, 1,

linealmente dependientes en 8 1 .u .
Tales que son

[Link]
Demostración: sean 8/ , 8 , … 81 , n-incógnitas tales que 8/ / . 8 . ⋯
81 1 . 0 para en u.
Sea f un punto cualquiera de u, consideremos el sistema de ecuaciones lineales.

8/ / . f 8 . f ⋯ 81 1. f 0

8/ ´/ . f 8 ´ . f ⋯ 81 ´1 . f 0

.
.

⋯
.
8/ /
.15/
. f 8 .15/
. f 81 1
.15/
. f 0

Como . /. , . ,…, 1. 0 para en u.

. / , , … 1 en base al teorema 18 del Apéndice Uno.

Luego el determinante de (1) es cero y el sistema tiene una solución no trivial

Luego la función . ∑1/ / / .

Es solución a la ecuación diferencial y satisface las condiciones . f

0, ´. f 0, … , .15/ . f 0

Como la función . 0 para en u es solución a la ecuación con esas

condiciones, entonces por el teorema de existencia y unicidad de la solución la
ecuación diferencial lineal normal homogénea de orden , tenemos

/ . . ⋯ 1 1. 0, para en u

Luego / , , … , 1 , son n-funciones linealmente dependientes ya que no todas las

/ para t 1,2, … , , no son cero.

el Wronskiano de / , , … , 1 , se anula e n un solo punto de u. Luego podemos

Al observar la demostración vemos que solamente se hace uso del hecho de que

deducir:

homogénea normal de orden en en lineamiento independiente en 8 1 .u , y por

Teorema 9: un conjunto de soluciones de una ecuación diferencial lineal

Wronskiano nunca se hace cero en u.

consiguiente una base para el espacio solución de la ecuación, si y solo sí, su

Ejemplo: el Wronskiano de / , es nulo para todo en un intervalo abierto u,

demostrar que el cociente ' es constante en u.
'#
m

[Link]
 '# ."
'm ."
Solución: para ello tomemos el cociente y derivamos, obtenemos:

´ . /. ´/ . . [ /. , . ]
0
/ . / .

en u, luego '# es constante en u.

'
Para todo
m

el Wronskiano de / ,
´´ D ´ 0 , siendo D y constantes.
Ejemplo: sea soluciones de la ecuación diferencial

´ D 0 y por tanto
. .0 9 5p" .
a) Demostrar que satisface la ecuación de primer orden

b) Suponiendo que / no es idénticamente nula, demostrar que .0 si y solo

'
si '# es constante.
m

Solución:

/. .
. /. , .
´/ . ´ .
a). Como

/. ´ . ´/ . .

Luego

´[ /. , . ] /. , ´ ´. . , ´/ ´.

Entonces

´[ /. , . ] D [ /. , . ]

/. ´ ´ . . ´/ ´ . D /. ´ . D . ´/ .

/. [ ´ ´ . D ´ . ] . [ ´/ ´. D ´/ . ]

/. [ . ] . [ /. ] 0

Que era lo que deseábamos probar.

Aplicando la solución de la ecuación diferencial lineal homogénea de orden uno, a

[Link]
 ´ D 0

tenemos que . .0 9 5F ."

donde C . ?f D M D ,
"

entonces . .0 9 5F."

b). Si / . es diferente de cero para todo en u, y suponemos que

[ / .0 , .0 ] 0, probemos que / / es constante, para ello derivemos dicho
cociente, y obtenemos

. ´ . /. ´/ . .
³ ´
/ . /.

[ /. , . ]

/.

.0 9 5p"
0
/.

Luego el cociente / / es una función constante en u.

Ahora, supongamos que / / es constante para todo elemento de u, derivando

este cociente tenemos:

. ´ . /. ´/ . .
³ ´
/ . /.

[ /. , . ]
0
/.

Luego

[ /. , . ] 0 [ / .0 , .0 ]

7. ECUACIONES LINEALES DE SEGUNDO ORDEN CON COEFICIENTES

CONSTANTES

[Link]
Estudiaremos la solución a la ecuación

´´ D ´ Å . (1)

Donde D, y , son reales y Å ∈ 8 .u . Si Å. 0 para todo en u, la ecuación (1)

solución, para ello, supongamos que D 0, es decir, la ecuación homogénea se

se denomina homogénea, consideremos en principio esta ecuación y hallemos su

transforma en ´´ 0 (2).

Consideremos tres casos:

0, luego la ecuación (2) se transforma en ´´ 0, luego las funciones

1, y
1) Sea
son solución a la ecuación diferencial dada.

Estas funciones son linealmente independientes porque

1
[1, ] 1
0 1

Entonces cualquier solución de ´´ 0, viene dada por

donde / , y, son constantes, en base al teorema 6.

2) Sea mayor que cero, entonces existe un real : tal que : y la ecuación
(2) se transforma en ´´ : , son soluciones a esta ecuación las funciones
Ëj: , y, j9 : , las cuales son linealmente independientes, lo cual se
puede demostrar usando el Wronskiano de estas funciones.

Entonces cualquier solución a ´´ 0, para mayor que cero está dada

por / Ëj: j9 : donde /, y, son constantes y, : √ .

3) Sea menor que cero, entonces existe un real : tal que : , entonces la
ecuación (2) se transforma en ´´ : , son solución a esta ecuación las
funciones 9 ;" , y, 9 5;" , las cuales son linealmente independientes, ya
que al calcular el Wronskiano de estas funciones

Obtenemos 2:, el cual es diferente de cero, implica que cualquier solución de

la ecuación ´´ 0 para menor que cero, es:

[Link]
 /9 9 5;" , donde son constantes y, :
;"
/ y .

funciones en 8 .u tal que 9 # , entonces

0Û

solución de la ecuación ´´ D ´ 0 si y solo sí, satisface la ecuación:

Teorema 10: sean y, es

4 D
´´ 0
4

Demostración: supongamos que:

9
Û
# , es solución a la ecuación ´´ D ´ 0 , entonces la satisface,
9
Û
como # ,

Luego
5p" D 5p"
´ 9 9
2
5p" D 5p" D 5p" D 5p"
´´ ´´ 9 ´ 9 ´ 9 9
2 2 4

Luego

5p" 5p" D 5p" 5p" D 5p" 5p"

´´ 9 D ´ 9 9 D ´ 9 9 9 0
4 2
5p" D D
9 ´´ D ´ D ´ 0
4 2

Luego

D
´´ 0
4

Entonces

4 D
´´ 0
4

Que era lo que queríamos demostrar.

9 # , es solución de ´´
Û
0
Žr5 p#
Ž
Ahora, supongamos que

[Link]
Luego la satisface, al reemplazar tenemos:

p" p" D p" 4 D p"

´´ 9 D ´ 9 9 9 0
4 4
p" D D
9 À ´´ D ´ Á 0
4 4

Luego

´´ D ´ 0

9
Û
De aquí concluimos que # , solución a ´´ D ´ 0

es solución a ´´ 0, como sabemos la solución a esta última

p# 5Žr
Ž
Si
ecuación viene dada por / /. . , donde /, son constantes, y en la
cual

a) /. 1, . , si D 4 0
/. 9 , . 9 , si D 4 es mayor que cero, y : √D 4
;" 5;" /
b)
c) /. cos : , . j9 : , si D 4 es menor que cero, y :
√4 D
/

Es decir, podemos caracterizar este resultado en la siguiente proposición.

Teorema 11: la solución a la ecuación ´´ D ´ 0, viene dada por la

expresión 9 # [ / / . . ] donde / y,
Û
son constantes y, / ,
funciones que dependen del valor de D 4
son

discriminante de la ecuación cuadrática È DÈ 0, de la siguiente forma:

, que denominaremos el

a) Si 0 entonces / . 1, .
es mayor que cero, entonces / . 9 ;" , . 9 5;" donde : √
/
b) Si
c) Si es menor que cero, entonces /. cos : , . j9 : donde
: √
/

Ejemplo: hallar la solución de ´´ 3 ´ 2 0

Solución: le asociamos la ecuación cuadrática È 3È 2 0

Luego

[Link]
 D 3, 2

entonces el discriminante es 1, el cual es mayor que cero, lo cual implica que la

9 # P / 9 #"
’Û m
9 5#" Q
m
/9 9 " , donde
"
solución a la ecuación dada es /, y
son constantes.

En este momento ya estamos en capacidad de hallar la solución a

´´ D ´ Å

Donde D y son constantes y, Å ∈ 8.u . Definimos el operador K de 8 1 .u en 8.u ,

tal que K. ´ D ´ , luego K. Å, la solución a esta ecuación viene
dada por la suma de la solución general a K.
K. Å, como vimos al final del capítulo cero.
, con una solución particular a

De esta última solución nos ocuparemos; un método es el de “variación de

parámetros”, el cual fue utilizado por Lagrange no solo para ecuaciones de orden
dos sino de mayor orden.

El método consiste en hallar funciones M/ y M tales que M/ V/ M V sea solución de

K. Å, donde V/ , V son soluciones linealmente independientes de K. 0.

Como dicha función verifica la ecuación K. 0, entonces derivamos y

reemplazamos, luego:

Sea

M/ V/ M V

´ M´/ V/ M/ V´/ M´ V M V´

´´ M/ V´´/ M V´´ .M´/ V´/ M´ V´ .M´/ V´/ M´ V´ ´

Como

K. ´´ D ´ Å y K.V/ K.V 0

Luego

K. .M´/ V´/ M´ V´ .M´/ V/ M´ V D.M´/ V/ M´ V

Como deseamos hallar M/ y M de tal manera que S. Å, entonces podemos

hacer M´/ V/ M´ V 0 y M´/ V´/ M´ V´ Å que son dos ecuaciones lineales con
las incógnitas M´/ y M´

[Link]
Luego

V Å
M´/
[V/ , V ]

V/ Å
M´
[V/ , V ]

Entonces

V . Å.
M/ . +
[V/ . , V . ]

V/ . Å.
M . +
[V/ . , V . ]

Se puede observar que M/ y M se expresan como integrales indefinidas, cada

una de las cuales queda determinada si le sumamos una constante, es decir M/
/, y M ; podemos definir una función Ï / V/ V y aplicamos el
operador K a esta función, tenemos K. Ï K. K. / V/ V K. ya que K
es un operador lineal, entonces Ï es también solución particular de la ecuación
diferencial no homogénea.

Ejemplo: hallar la solución general a ´´ j9

Solución: para hallar la solución particular a dicha ecuación usando el método de

´´ 0. Como el discriminante es menos cuatro, luego la solución general a la

variación de parámetros, para ello primero hallaremos la solución de la ecuación

ecuación homogénea es / j9 2 cos 2

Sea V/ . j9 2 yV . Ëj2

Hallemos M/ y M de tal manera que M/ V/ M V sea solución de ´ j9

Luego

Ëj2 j9 Ëj j9 j9 j9
M/ . + + +
2 2 2

j9 Ëj Ëj
+ Ëj2 j9 +
2 3 2

[Link]
 M . ? j9 Ëj
]q1 " ]q1" ]q1’ "
y ?

Entonces la solución a la ecuación diferencial dada es:

j9 Ëj Ëj
/ j9 2 Ëj2
3 3 2

ecuación definida por el operador K, en el cual K. ´´ D ´

Otro método por medio del cual podemos calcular un solución particular a una

K. Å, es el llamado el de coeficientes indeterminados, y como su nombre lo

, siendo

dice se trata de dar como solución un “polinomio” en el cual calculamos los

coeficientes, para que este verifique la ecuación, podemos considerar tres casos.

Caso 1: si Å es un polinomio de grado , con ≠ 0, podemos hallar un polinomio

. ∑1; f D; ; que satisface la ecuación K. Å; sustituyendo en la ecuación
podemos hallar los coeficientes Df , D/ , … , D15/ , D1 .

Ejemplo: hallar la solución general a ´´

/9
"

9 5" , donde / y
Solución: la solución general de la ecuación homogénes asociada es
son constantes.

método de coeficientes, ya Å. 1, luego Å es un polinomio de grado

Ahora calculamos una solución particular a la ecuación no homogénea, usando el

1, entonces supongamos que dicha solución es C G, derivando esta solución y

y,

reemplazando en la original tenemos:

C G , de donde concluimos que C 1 y, G 0.

general a la ecuación ´´
Entonces la solución particular es . Luego ya podemos dar la solución
la cual es:

/9 9 5"
"

0, la ecuación ´´ D ´ Å, no se puede satisfacer con un polinomio de

1, siempre que D ≠ 0. Si D 0, su
Si

2, obtenido de dos integraciones

grado ; pero si con uno de grado
solución general es un polinomio de grado
sucesivas.

Caso 2: si Å. R . 9 " , siendo R un polinomio de grado y È es una

. 9 " , la ecuación
diferencial ´´ D ´ Å se transforma en ´´ .2È D ´ .È DÈ
constante; si hacemos el siguiente cambio de variable

R. La cual es una ecuación del caso anterior.

[Link]
Caso 3: si Å. R. 9 " cos o Å. R. 9 " j9 , siendo R un
polinomio; È y constantes. En ambos casos, existe siempre una solución
particular de la forma

. 9 " [¾.
cos Ñ . j9 ] donde ¾ y Ñ son polinomios.

.5.2 Ecuaciones diferenciales lineales de segundo orden no homogéneas

Algunas ecuaciones de segundo orden se pueden resolver por procesos de integración

inmediata o por reducción de orden. Los casos que resolverán aquí son la de forma:

2.
#!

"#
1. Ecuaciones del tipo

La solución general se obtiene integrando dos veces, así

? 2. . [? 2. ]
!
" / /

Ejemplo 1

10
#!

"#
Hallar la solución general de la ecuación diferencial

? 15 5 8/ . ?.5 8/ 8/ 8
! 6 Ž
" Ž

5
. Ž
8/ 8
4

2$ , %
#! !
"# "
2. Ecuaciones del tipo

"
En estas ecuaciones se utiliza el cambio de variable para llevarlas a una ecuación
2. , .
ø

"
diferencial de primer orden de la forma

Ejemplo 2.

Hallar la solución general de la ecuación diferencial . 1 0

#! !
"# "

[Link]
 ! ! #!

" " "#

Haciendo el cambio de variable , entonces , por lo que la ecuación dada se
transforma en

. 1 0
"

"
"5/

0
"
"5/

Cuya solución es

ln ln. 1 ln /

ln . 1 ln /

z . 1 /

z
vm
"5/

Ahora deshacemos el cambio

), vm
" "5/

De donde se obtiene

y.x
v
? "5/
m

Y finalmente la solución general es

y.x / ln. 1 /

Obsérvese que la solución general tiene las dos constantes / y ya que la ecuación
diferencial es de segundo orden.

[Link]
Ejemplo 3.

Hallar la solución general y particular de la ecuación diferencial.

0, con las condiciones iniciales y.1 .1 0

#! / !
"# " "

! #!

" " "#

Haciendo el cambio de variable , entonces , por lo que la ecuación dada se
transforma en

1
0 ⟹ ⟹ ln ln 8/ ⟹ 8/

Deshaciendo el cambio

8/
8/ ⟹ + + 8/ ⟹ . 8
2

La solución general es . 8/ 8

Y la solución particular será I. 1

2. 2$ , %
#! #! !
"# "# "
3. Ecuaciones del tipo

Para reducir el orden de estos dos tipos de ecuaciones, hay que considerar la como la
!
"
nueva variable independiente, por lo que el cambio se completa así:

Al realizar el cambio de variable en estos tipos de ecuaciones, las nuevas ecuaciones

diferenciales de primer orden quedan así:

!
2. y !
2. ,

[Link]
Ejemplo 4.

0.
#!

"#
Hallar la solución general de la ecuación diferencial

Aplicando los cambios

! #!

" ! "#
y tenemos,

0
!

8/ ⟹ º 8/

Deshaciendo el cambio

º 8/ ⟹ ⟹ ?
! ! !
" º! # >sm º! # >sm
?

La integral de la izquierda se realiza con un cambio de variable como que sigue

º8/ tan º8/ j9

+ +
º 8/ º8/ j9 º8/ j9

+ sec ln|sec MD | 8 ⟹ sec tan 8 9"

[Link]
De acuerdo al cambio de variable en

º8/ sec
ºsm >! #
ºsm

º8/ tan
!
ºsm

Entonces

8 9" ⟹ º8/ 8 9" ⟹ º8/ 8 9"

ºsm >! # !
ºsm ºsm

Elevando al cuadrado en ambos lados de la última ecuación

Pº8/ Q .8 9 " ⟹ 8/ 8 9" 28 9 "

⟹ 8/ 89 "
28 9 "

8 9 " 8/
28 9 "
89 "
8/ ⟹ 8 9" 8Ž 9 5"
28 9 "

. 8/ 9 " 8 9 5"

Que finalmente es la solución general.

[Link]
8. EJEMPLOS DE PROBLEMAS FÍSICOS QUE CONDUCEN A ECUACIONES
LINEALES DE SEGUNDO ORDEN CON COEFICIENTES CONSTANTES

8.1 MOVIMIENTO ARMÓNICO SIMPLE

que puede ser el eje , bajo la acción de una fuerza localizada en el origen 0.
Consideremos una partícula que se mueve sin resistencia a lo largo de una recta,

0 en un tiempo M; suponiendo positiva a la derecha y negativa a la izquierda de

Supongamos que la fuerza es proporcional al desplazamiento de la partícula de

0.

será
#"

J#
Puesto que la fuerza es proporcional a la aceleración, la aceleración
proporcional al desplazamiento .

Si la partícula se encuentra a la derecha de 0, el desplazamiento es positivo y la

encuentra a la izquierda de 0, el desplazamiento es negativo y la fuerza (hacia a

fuerza (hacia la izquierda) y la aceleración son negativas; si la partícula se

derecha) y la aceleración son positivas; luego la aceleración y el desplazamiento

son opuestos, y la ecuación diferencial que define el movimiento es:

kx
M

Donde : es la constante de proporcionalidad y es positiva. Su solución en base al

teorema 11 es:

/ j9 M cos M (1)

Donde / y son constantes y √

Sea C º y cos , sen

vm v#
/ F F

DÔ MÉ v#
v
Luego
m

Entonces la ecuación (1) se transforma en Cj9 . M (2)

La ecuación (1) y (2) da las dos formas equivalentes para V

"
J
luego:

[Link]
 V / cos M sen M

V Acos. M

Este movimiento se denomina armónico simple y la partícula vibra entre C y – C, la

½
cual es la amplitud del movimiento. El movimiento es periódico, su periodo es y,

½
la frecuencia es .

La gráfica del desplazamiento contra tiempo y la de velocidad contra tiempo es:

línea recta. Cuando M 0, la aceleración es 9 metros por segundo cuadrado, la

Ejemplo: una partícula se mueve con un movimiento armónico simple en una

velocidad es de 4.5 metros por segundo y el desplazamiento es de 4 metros a la

izquierda del origen.

Encontrar:

a) El desplazamiento cuando ha transcurrido medio periodo

b) El primer tiempo cuando el desplazamiento es cero y cuando es un metro
c) La velocidad máxima

[Link]
 :
#"

J#
Solución: la ecuación diferencial es

Podemos calcular el valor de : ya que 9 4 luego 9 4: ,

#"

J#
cuando
entonces : Ž

Luego la ecuación diferencial se transforma e:

9
M 4

Su solución es

3 3
/ j9 M Ëj M
2 2
3 3 3 3
V / Ëj M j9 M
2 2 2 2

Podemos hallar / y 4,
V 4,5 para M 0 y obtenemos / 3 y, 4
, empleando las condiciones de frontera dadas;

Luego

3 j9 M 4 Ëj M (1)

V Ëj M 6 j9 M (2)

Sea

D º3 . 4 5 , cos y, j9
5Ž
6 6

Entonces

5 j9 $ M % (3)

V Ëj $ M %
/6
(4)

[Link]
 Ž½
Para hallar la respuesta a la parte a) debemos hallar el periodo el cual es , la
½
mitad de este s sustituyendo en la ecuación (1) tenemos

3 j9 4 cos 4 metros
½ ½

Para hallar el tiempo cuando 0 , ó 1 metros, empleamos la ecuación (3),

luego 0 5 j9 $ M %

, como tenemos que j9

5Ž
6
Para ello primero calculamos el valor de y,
cos 6
entonces pertenece al cuarto cuadrante, luego 0,9273,
reemplazando este valor tenemos

3
0 5 j9 M 0,9273
2
Luego

3
M 0,92730
2

Entonces M 0,618 segundos

Cuando 1 , tenemos que hallar el valor de M positivo más pequeño que

satisface:

3
1 5 j9 M 0,9273
2

Luego

3
M 0,9273 DÔ j9 0,2 0,2013
2

M 1,1286 luego M 0,752 segundos

Para la parte c) tenemos que V Ëj $ M %, luego la máxima velocidad es 7,5

/6

metros por segundo.

Ejemplo: un resorte fijo en su extremo superior, soporta un peso en el extremo

inferior, que alarga el resorte 15 centímetros.

[Link]
a) Si el peso se baja 7,5 centímetros adicionales de su punto de equilibrio y se
suelta, encontrar el periodo de vibración y la ecuación del movimiento del peso.

b) Si el peso se hace descender 7,5 centímetros bajo su posición de equilibrio y se

le aplica una velocidad hacia abajo de 30 centímetros por segundo, encontrar
la distancia, bajo la posición de equilibrio, del punto más bajo que alcanza el
peso; la velocidad máxima y el tiempo que se requiere para que el peso
regrese a la posición de equilibrio por primera vez.

Solución: primero planteemos el problema, para hallar la ecuación diferencial que

representa el movimiento vibratorio de un peso colgando del resorte, bajo la
suposición de que no existe resistencia al movimiento.

En principio tenemos un resorte colgado de su extremo superior verticalmente,

Figura 13(a). En la figura 13 (b) el resorte tiene un peso , en el extremo inferior el
cual lo alarga una distancia a metros; el peso se halla en posición de equilibrio en
la figura 13 (b).

En base a la Ley de Hooke, tenemos que la tensión del resorte es proporcional a

su alargamiento. La constante de proporcionalidad se llama constante del resorte.

:D, donde : es la constante del resorte,

En la posición de equilibrio la tensión del resorte es W kilogramos, y el
alargamiento es a metros; luego
entonces : kilogramos por metro.
p

Ahora, se inicia el movimiento del peso ya sea dándole una velocidad vertical
inicial o puede desplazarse el peso verticalmente de su posición de equilibrio y
entonces soltarlo, o bien darle una velocidad vertical inicial. Se toma la posición de
equilibrio como el origen y consideramos positivo por debajo de la posición de
equilibrio y negativo por encima de éste.

La cantidad de desplazamiento del peso en un tiempo M es D , independiente

tensión del resorte es :.D , figura 13(c).

del hecho que el peso se encuentre arriba o abajo de la posición de equilibrio y la

[Link]
El resorte hala entonces el peso hacia arriba con una fuerza : .D y la

:.D : ya que :D
gravedad lo solicita hacia abajo con una fuerza , luego la fuerza resultante que
actúa sobre el peso es

Como la fuerza que actúa sobre un peso es proporcional a la aceleración,

entonces la ecuación que define el movimiento es:

È :
M

Luego el peso ejecuta un movimiento armónico simple con respecto a su posición

de equilibrio.

De acuerdo a las condiciones del problema, tenemos que el peso produce un

0,15å , entonces :
f,/6
alargamiento de 0,15 metros; luego kilogramos por
metro.

Luego la ecuación diferencial se transforma en:

É
0
M 0,15

Su solución es

/ j9 M Ëj M
f f
(1)

V À Ëj M j9 MÁ
f f f
/ (2)

Donde / y
0,075, V 0 para M 0, sustituyendo tenemos que / 0 y,
son constantes que podemos hallar con las condiciones dadas en
a) y estas son
0,075, luego la ecuación del movimiento es 0,075 cos M.
f

Su periodo es 2Ó 0,778 segundos

f

[Link]
Para la parte b) empleamos las condiciones dadas, es decir 0,075, V 0,30
para M 0 y obtenemos que / 0,3 0,075. La distancia bajo la
f
y,
posición de equilibrio del punto más bajo que alcanza el peso es la amplitud del
movimiento; esto es:

3 5É 12
º H
/
10 80É

La velocidad máxima es el valor de V en la ecuación (2); es decir

20É 3 5É 12
H . H
3 /
20 3

Se obtiene cuando el peso pasa por su posición de equilibrio.

0 en la ecuación (1), y resolviendo para el valor positivo más

El tiempo requerido para que el peso alcance su posición de equilibrio se obtiene

pequeño de M:
haciendo

0,15 20É
0 0,3 H j9 H M
É 3

20É
0,075 Ëj H M
3

20É 20É
MDÉ H M H 2,0023
3 4√3

20É
H M Ó MDÉ5/ 2,0023 2,01
3

M 0,2525 segundos
,f/ √
º f

Ejemplo: se a CG figura 14, el alambre, siendo C el punto fijo de suspensión y G

el otro extremo del alambre, al cual se halla unida una masa È. Designemos por
el ángulo formado por al alambre y la vertical C0 en un instante cualquiera.

[Link]
Cuando la masa È se halla en movimiento, actúan sobre ella dos fuerzas, la
tensión Ò de la cuerda y el peso ÈÉ de la masa.

Podemos descomponer al peso ÈÉ en dos componentes una tangente y otra

perpendicular a la trayectoria del movimiento, tenemos que la componente

ÈÉ j9 . Tomemos positivo si È se halla a la derecha de 0 y es negativo si

perpendicular es equilibrada por la tensión, la fuerza tangencial a la trayectoria es

È se halla a la izquierda de 0. Cuando es positiva la fuerza resultante se halla

fuerza resultante es ÈÉ j9 .
dirigida hacia la izquierda y se dirige hacia la derecha si es negativo, luego la

Como la longitud del arco es j S , luego aplicando la Ley de Newton tenemos:

j
È ÈS ÈÉ j9
M M

O sea

j9
#'

J# /
(1)

entre menos 5º y -5º, y es j9

Esta ecuación no se puede resolver exactamente por funciones elementales, pero
podemos hacer la aproximación, para ,
donde aparece expresado en radianes.

Luego la ecuación (1) se convierte en:

É
0
M 1

Su solución es / j9
/
M Ëj
/
M , donde / y son constantes y el

periodo es
½√
√/
.

[Link]
Otro ejemplo de movimiento armónico simple es el de un cuerpo que flota en el
agua. Para analizar este caso resolvamos este problema. Una caja cúbica de 3
metros de lado flota en agua tranquila. Se nota que la caja oscila subiendo y
bajando con un periodo de medio segundo. ¿Cuál es el peso de la caja?

CG8. En la figura 16 el cubo está casi totalmente sumergido en el agua, en esta

Solución: la figura 15 representa el cubo en su posición de equilibrio, dado por

posición hay una fuerza que tiende a empujar la caja hacia arriba. Esta fuerza está
determinada por el Principio de Arquímedes.

Un peso parcial o totalmente sumergido en un líquido, experimenta un empuje

hacia arriba originando un movimiento vibratorio por una fuerza igual al peso del
líquido que desaloja.

Luego en la figura 15, según el principio dado el peso de la parte u sumergida es

igual al peso de volumen de agua que desaloja.

En la figura 16 se muestra la parte u sumergida y otra parte de cubo dentro del

peso de agua desalojada es 3 3 1.000, o sea 9000 , ya que 1.000 kilogramos

agua a una altura , lo cual indica que hay una fuerza no equilibrada, luego el

por metro cúbico es la masa específica del agua.

El valor de la fuerza que tiende a mover el cubo es análoga a la fuerza de

restitución del resorte vibrante. Si designamos por el peso de la caja tendremos
la ecuación diferencial:

. 9000
É M

O sea

88200
0
M

[Link]
 / cos M sen M
¬¬ ff ¬¬ ff
Su solución es

Donde / y, son constantes, el periodo de vibración es:

2Ó 1
√88200 2

Resolviendo esta ecuación tenemos que 559 :Éj

En todos estos problemas hemos considerado que la resistencia es despreciable,

ahora consideraremos en caso en el cual la resistencia no es nula.

8.2 FUERZA DE ATRACCIÓN PROPORCIONAL AL DESPLAZAMIENTO;

RESISTENCIA PROPORCIONAL A LA VELOCIDAD

El movimiento vibratorio con amplitud constante, resulta de una fuerza

proporcional al desplazamiento. Si se toman en cuenta la fricción o fuerzas de
rozamientos que disminuyen la amplitud de las oscilaciones y finalmente detienen
el movimiento, esta fuerza se denomina amortiguadora o resistente; no se conoce
ley exacta que gobierna estas fuerzas puesto que dependen de muchos factores
variables, pero se observa experimentalmente que, si la velocidad es pequeña, la
magnitud de la fuerza amortiguadora es aproximadamente proporcional a la
velocidad instantánea, o sea su magnitud es J donde es la constante de R
"

proporcionalidad.

La fuerza amortiguadora se opone al movimiento en caso de una partícula que se

desplaza sobre la línea recta, sucede que si la partícula se mueve hacia la
J representa la fuerza resistente.
"
derecha, la velocidad es positiva y

J
"
Si la partícula se mueve hacia la izquierda, la velocidad es negativa y
representa la fuerza amortiguadora; similarmente podemos considerar el caso de
movimiento producido por un resorte; luego el modelo matemático que define
estas situaciones es:

È :
M M

O sea

[Link]
 c :
. 0
M È M È

La solución a esta ecuación depende del valor de:

c 4k c 4:È
È È È

1) Si c 4:È 0, entonces la solución es

vJ
x e5
[ / M]

donde / y son constantes, su gráfica es:

1) Positiva
2) Cero Velocidad Inicial
3) Negativa

2) Si 4:È es mayor que cero, la solución es:

vJ
x e5
[ / 9 J
95 J]

√c 4:È , siendo c/ y c
/
Con constantes, su gráfica es:

1) Positiva
2) Cero Velocidad Inicial
3) Negativa

[Link]
3) Si c 4:È es menor que cero, la solución viene dada por:
vJ
x e5
[ / cos M j9 M]

√4:È
/
con , donde / y son constantes, su gráfica es:

1) Positiva
2) Cero Velocidad Inicial
3) Negativa

!
En todos los casos el factor 9 5#" se llama factor de amortiguación.

Ejemplo: una partícula parte del reposo para M 0, con un desplazamiento 5

metros a la derecha del origen y se mueve a lo largo del eje , de acuerdo a la
ecuación diferencial:

1,25 0
M M

Encontrar:

a) El tiempo requerido para que el factor de amortiguación disminuya el 50%

b) El porcentaje de disminución del factor de amortiguación después de un periodo
c) La posición de la partícula después de un periodo y la velocidad en ese instante
d) Hallar la gráfica del desplazamiento contra tiempo y la de velocidad contra
tiempo

Solución:
!
a) El factor de amortiguamiento es 9 5# , para M 0, este es 1, debemos hallar el
tiempo M, en el cual el factor de amortiguamiento es , luego
/

J 1
95
2

[Link]
 M
S 2
2

M 21 S 2 1,39 segundos

calculamos 9 5 0,0432; luego su porcentaje de disminución es 95,7%.

b) Hallemos el factor de amortiguación después de un periodo, para ello

c) Como la solución a la ecuación diferencial es:

!
9 5# [ / j9 M cos M] (1)

Derivando tenemos
!
V 9 5# P$ % cos M . j9 MQ
v# vm
/ (2)

Podemos hallar / y con las condiciones dadas de 5, V 0 , para M 0,

5.
6
reemplazando en (1) y en (2), obtenemos que / y

Luego:

J 5
9 5 ³ j9 M 5 cos M´
2
J 15
V 9 5 ³ sen M´
2

Entonces

.2Ó 9 5½ .5 0,216 metros

O sea la partícula después de un periodo se halla a 0,216 metros a la derecha del

origen.

V.2Ó 9 5½ .0 0 metros por segundo

d) Para realizar las gráficas pedidas expresamos estas funciones de una manera
más sencilla

sea C √ 1 2 25 DÔ MDÉ v m DÔ MDÉ

6 6√6 v /
Ž
y
#

[Link]
 Luego la ecuación (1) se convierte en:

5√5 5 J
9 [ j9 .M ]
2
J 15
V 9 5 ³ sen M´
2

Ejemplo: se coloca un peso de 3 kilogramos en un resorte el cual se estira 153

milímetros. Si se tira del peso 10 centímetros por debajo de su posición de
equilibrio y luego se suelta. Hallar la posición del peso en función del tiempo si
existe una fuerza de amortiguamiento que es 6,12 veces la velocidad.

Solución: la ecuación diferencial que define el movimiento es

3
19,6 6,12
9,8 M M

O sea

20 64 0
M M

Las condiciones iniciales son 0,10 metros, V 0, para M 0. Luego al

reemplazar las condiciones dadas en la solución a la ecuación diferencial:

C9 5ŽJ G9 5/7J

Tenemos

[Link]
 4 5ŽJ 1 5/7J
9 9
9 9

Cuya gráfica se muestra en la figura 22. Es claro que no habrá oscilaciones; se ha

amortiguado tanto al peso que éste solo volverá gradualmente a su posición de
equilibrio, sin pasar por él.

8.3 VIBRACIONES FORZADAS

Hemos visto problemas únicamente de resortes donde se consideran fuerzas de

restitución y amortiguamiento; ahora trataremos casos en que pueden actuar otras
fuerzas externas que dependen del tiempo. Estas fuerzas suelen presentarse, por
ejemplo, cuando el resorte sube y baja de manera prescrita, como un movimiento
periódico, o bien cuando se empuja ligeramente al peso cada vez que alcanza su
posición más baja.

Si designamos por #.M la fuerza externa, la ecuación diferencial del movimiento

del resorte es:

: #.M
É M M

O sea

: #.M
É M M

Ejemplo: se coloca un peso de 3 kilogramos en un resorte al cual la estira 153

milímetros. Si se tira del peso hasta 10 centímetros por debajo de su posición de
equilibrio y luego se suelta.

Hállese la posición del peso en función del tiempo, si existe una fuerza

periódica externa dada por #.M igual a 39,2 cos 8 M.

amortiguadora, la cual es 2,45 veces la velocidad y además actúa una fuerza

Solución: la ecuación diferencial que define el movimiento es:

3
19,6 2,45 39,2 cos 8M
9,8 M M

[Link]
 8 64 128 cos 8M
M M

La cual es una ecuación diferencial lineal no homogénea de orden dos con

coeficientes constantes, entonces su solución está dada por la suma una solución
particular a la ecuación dada y la solución general a la ecuación homogénea
asociada.

La solución a la ecuación homogénea asociada es:

9 5ŽJ .C cos 4√3 M G j9 4 √3 M

Supongamos que D j9 8M cos 8M, es solución particular a la ecuación no

que D 2, 0.
homogénea, reemplazando esta solución en la ecuación diferencial dada tenemos

Luego la solución a ecuación diferencial dada es:

9 5ŽJ ^C cos 4√3 G j9 4√3 M` 2 j9 8M

Podemos hallar el valor de C y G, empleando las condiciones iniciales, las cuales

son 0,1 V 0 para M 0, de donde tenemos C 0,1 G 1,33√3 y por
ende

0,1 9 5ŽJ ^cos 4√3 M 13,3 √3 j9 4 √3 M` 2 j9 8M

Cuya gráfica es la siguiente:

En la última ecuación observamos que los términos que contienen 9 5ŽJ se vuelven
despreciables M es grande; estos términos se llaman “Términos Transitorios” y sólo
son significativos cuando M es cercano a cero.

Algunas veces, cuanto estos términos transitorios de la solución son significativos,

reciben el nombre de “Solución Transitoria”.

El término 2 j9 8M permanece aún cuando los términos transitorios sean

despreciables y se denomina “Término de Estado Estacionario o Solución de

[Link]
Estado Estacionario”, porque indica el comportamiento del sistema cuando se han
estabilizado las condiciones de funcionamiento. Se ve que la solución de estado
estacionario (curva de trazos) es periódica y que tiene la misma frecuencia y
periodo que la fuerza externa aplicada.

8.4 EL FENÓMENO DE RESONANCIA

Cuando la frecuencia de una fuerza externa periódica aplicada a un sistema se

relaciona de una manera sencilla, que se describirá, con la frecuencia natural del
sistema, puede aparecer el fenómeno de resonancia, que aumenta la magnitud de
las oscilaciones hasta valores tan grandes que el sistema llega a disgregarse. Es
por ello que una compañía de soldados que marche al paso sobre un puente
puede provocar la caída del mismo, aún cuando éste podría soportar mucho más
soldados si no caminaran al paso. En forma similar puede ocurrir que una nota
musical con una frecuencia característica adecuada rompa un cristal.

El fenómeno de resonancia es funesto en los casos mecánicos mientras que en

los sistemas eléctricos sus efectos son útiles. La radio, la televisión, el radar y las
comunicaciones serían virtualmente imposibles sin la resonancia eléctrica; en
estos casos la corriente y, en consecuencia, la energía eléctrica generada puede
aumentarse hasta grandes intensidades que necesitan estos campos.

El que podamos sintonizar nuestro radio a la frecuencia de la estación

radiotransmisora para conseguir una recepción clara se debe a la resonancia
eléctrica.

Un ejemplo de resonancia mecánica es:

Ejemplo: dado un resorte al cual se le coloca un peso de 3 kilogramos que lo

estira 153 milímetros. Si se tira el peso hasta 10 centímetros por debajo de su

movimiento si se le aplica una fuerza externa dada por 4,9 cos 8M.
posición de equilibrio y luego se suelta. Hallar la ecuación que describe el

Solución: la ecuación diferencial es:

3
19,6 4,9 cos 8M
9,8 M

O sea

64 16 cos 8M
M

[Link]
La solución general a la ecuación homogénea es

C cos 8M G j9 8M

Ï M .D cos 8M j9 8M
sustituyéndola en la ecuación diferencial, obtenemos D 0 y 1. Luego la
Supongamos que la solución particular es

solución general es:

C cos 8M G j9 8M M j9 8M

0,1 V 0 para M 0
M j9 8M
Como las condiciones iniciales son de donde
obtneemos que

Ya que C G 0. La gráfica de esta última función está comprendida entre las

rectas M y M , tal como aparece en la figura 24.

En esta gráfica podemos observar que la magnitud de las oscilaciones aumentará

ilimitadamente, naturalmente el resorte se romperá al cabo de poco tiempo.

resonancia se produjo porque la frecuencia de la 2 fuerza externa aplicada era

Se observará que en este ejemplo se despreció el amortiguamiento y que la

igual a la frecuencia natural del sistema no amortiguado. Esta es un principio

general.

En el caso de que haya amortiguamiento, las oscilaciones no aumentan

ilimitadamente pero, no obstante, pueden hacerse muy grandes; en este caso, la
resonancia se produce cuando la frecuencia de la fuerza externa aplicada es
ligeramente menor que la frecuencia natural del sistema.

[Link]
8.5 CIRCUITOS ELÉCTRICOS

En el capítulo anterior se estudiaron circuitos eléctricos ÅK y Å8 que conducía a

ecuaciones diferenciales lineales, ahora estamos interesados en circuitos ÅK8, tal
como en la figura 25.

Entonces, según la Ley de Kirchhoff,

$
K Åu Æ .Ò
Ç
J s
(1)

Como u
Ç
J
, la ecuación (1) se transforma en:

¾ ¾ ¾
K Å Æ .Ò
M M 8

Si comparamos esta expresión con la ecuación general de las vibraciones

forzadas del numeral 8.3. Observamos una sorprendente analogía entre las
cantidades mecánicas y las eléctricas.

La carga ¾ corresponde a la posición .

La inducción K corresponde a la masa È, o sea

La resistencia Å corresponde a la constante de amortiguamiento

La inversa de la capacidad v , corresponde a la constante del resorte :

/

La fuerza electromotriz Æ .M corresponde a la fuerza externa aplicada # .M

La intensidad de la corriente u corresponde a la velocidad V

Ç "
J J

Debido a la notable analogía que existe entre estas cantidades mecánicas y

eléctricas, la cual se verifica en los casos más complicados, la mayoría de las
conclusiones que se dedujeron para los sistemas mecánicos se aplican a los
eléctricos, e inversamente. Es más, en la industria se utiliza con frecuencia esta
analogía para estudiar sistemas mecánicos donde la construcción de estos
resultaría muy cara o complicada, y cuando las consecuencias quizá pudieran ser
muy peligrosas.

[Link]
Ejemplo: se conectan en serie un inductor de 0,5 henrios, una resistencia de 6

alterno dado por 24 j9 10M, para M mayor o igual a cero.

ohmios, un condensador de 0,02 faradios, un generador que tiene un voltaje

Hallar la carga y la intensidad de la corriente en un tiempo M cualquier, si la carga

existente en el condensador es nula para un tiempo M 0.

Solución: de acuerdo a la Ley de Kirchhoff, tenemos:

¾ ¾
0,5 6 50 ¾ 24 j9 10M
M M

O sea

¾ ¾
12 100 ¾ 48 j9 10M
M M

La solución a la ecuación homogénea es

¾ 9 57J .C cos 8M G j9 8M

Supongamos que C j9 10M cos 10M , es solución particular de la ecuación no

homogénea dada, sustituyéndola en esta tenemos D 0 y 6
, luego la
solución general a la ecuación no homogénea es

2
¾ 9 57J .C cos 8M G j9 8M cos 10M
5

Con las condiciones iniciales ¾ 0, u 0, para M 0

Ç
J

Obtenemos que C 6
y, G
/f
, luego la solución buscada es

¾ 9 57J .0,4 cos 8M 0,3 j9 8M 0,4 cos 10M

Derivando tenemos:

u 59 57J j9 8M 4 j9 10M

Nótese que el término con el factor 9 57J es la solución transitoria que pronto se
hace despreciable. La solución estacionaria la constituye RD ¾ 0,4 cos 10M , y
para u 4 j9 10M ; las cuales se conservan después de que el término transitorio
ha desaparecido virtualmente.

[Link]
8.6 FLEXIÓN DE VIGAS

considerarse como compuesta de fibras tales como #´´#, todas originalmente de

En la figura 26 se muestra una sección de una viga flexionada. Puede

longitud j. La superficie L¾, que contiene fibras cuyas longitudes no se alteran

cuando se flexiona la viga, se llama superficie neutra.

La curva de una de estas fibras G ´G se denomina curva elástica de la viga. Las

cortan cuando la viga se flexiona. Supongamos que la fibra #´#

fibras debajo de la superficie neutra se alargan y las que están por encima se a

abajo de la superficie neutra, se alarga una distancia 9 por la

fuerza X C, donde X, es el esfuerzo por unidad de área transversal y C es el área
A una distancia

transversal de la fibra.

Sea Å la longitud de G8 (radio de curvatura) y sea ¾¾´ perpendicular a G´G y G8.

Por la Ley de Hooke, tenemos que el esfuerzo X por unidad de área es

proporcional al alargamiento por la longitud de la fibra #´# ; esto es
q
]

X Æ
q
]
(1)

Donde la constante de proporcionalidad Æ, es el módulo de Young, o modulo de

elasticidad. Además se tiene que:
q %
] Ì
(2)

Entonces de (1) y (2) tenemos:

Æ
X
Å

[Link]
El momento de la fuerza X C con respecto al ¾¾´ es

Æ
X C C
Å

Integrando para toda la sección transversal de la viga se tiene el momento

flexionante M:

Æ
< + C
Å

Pero

C u

Donde u es el momento de inercia del área transversal de la viga con respecto al

&$
eje ¾¾´ luego <
Ì

Ahora, se toma el eje horizontal a través de algún punto de la fibra G´G, tomada
como el origen y el eje positivo hacia arriba.

La fórmula para radio de curvatura es Å pero para la flexión reducida ´

./>!´#
! ´´
es pequeña y ´ puede despreciarse en comparación con la unidad, de manera
que una aproximación bastante precisa del radio de curvatura es Å
/
!´´
.

Reemplazando esta aproximación en la última ecuación tenemos:

< Æu ´´

Ejemplo: una viga horizontal de 21 metros de longitud está apoyada en sus

extremos. Hallar la ecuación de su curva elástica y su máxima deformación
vertical cuando tiene una carga uniformemente distribuida de kilogramos.

Solución: tomamos 0 como origen en el extremo izquierdo como en la figura 27.

Luego la ecuación diferencial asociada es

1
Æu 1
2

Luego

[Link]
 1 1 1
Æu 1 Ž
1
6 24 3

5 1Ž
p"
24Æu

8.7 EL CABLE SUSPENDIDO

suspendido entre dos soportes C y G como se muestra en la figura 28. El cable se

Supongamos que se tiene un cable de peso uniformemente distribuido ,

combará y tendrá su punto más bajo U.

Deseamos determinar la curva formada por e l cable suspendido. La curva se

denomina catenaria. Tomando ejes de coordenadas como se muestra en figura
29.

Tomando U en eje vertical, sea j lo longitud del cable desde U a un punto L de

coordenadas . , . Entonces la porción del cable desde U hasta L está sujeto a
fuerzas, una de ellas es la fuerza gravitacional j actuando hacia abajo a través
del centro de gravedad de la porción de cable desde U hasta L, otra fuerza es la
tensión Ò/ actuando tangencialmente en L, y finalmente la tensión Ò actuando
horizontalmente en U. La tensión Ò/ es variable y Ò es constante.

Suponemos que el sistema está en equilibrio, entonces

Ò/ j9 j 0 (1)

Ò/ Ëj Ò 0 (2)

[Link]
Dividiendo (1) entre (2) tenemos:

j
MDÉ
Ò

Entonces

j

Ò

Como X es la longitud del arco de la curva entonces derivando la última fórmula

tenemos

j
.
Ò
]
"
Si es constante, entonces la solución a la ecuación diferencial es

2Ò

Donde es una constante. En este caso el cable toma la forma de una parábola.

1 $ %
] ] ]
" " "
Si no es constante, como

Luego

j
H1
Ò

Ecuación que podemos resolver transformándola en una ecuación de primer orden

si definimos ℎ
!
"
, luego

ℎ
º1 ℎ
Ò

La cual es una ecuación de variables separables, al integrar tenemos

"
ℎ ºℎ 1 89 ¿#

[Link]
elevando al cuadrado tenemos

1 " "
R 9 9
5
¿ ¿
2
# #

Integrando tenemos

Ò " "
9 ¿# 9 ¿# [
5
2

Donde [ es una constante, si utilizamos funciones hiperbólicas, tendremos

Ò
cos ℎ [
Ò

La gráfica de esta curva recibe el nombre de catenaria.

[Link]
3° ECUACIONES DIFERENCIALES DE SEGUNDO ORDEN
3.1 ECUACIÓN DIFERENCIAL DE SEGUNDO ORDEN EN DONDE LA
SEGUNDA DERIVADA ES FUNCIÓN DE UNA SOLA CONSTANTE

'ç ö
, para y(0) = 0, y.f
Þ²•. ´
' ç µ().#
Desarrollar la ecuación diferencial = -9.8 = 20
3.1.1 DESARROLLO ANALÍTICO.
La ecuación diferencial de segundo orden se puede convertir en dos ecuaciones
diferenciales de primer orden.
'ö )* 'ç ö
Hacemos Z = ' , entonces: )• = ' ç
Remplazando la nueva variable en la ecuación diferencial de 2do orden, tendremos:
)*
)•
= - 9.8
'ö
'
Desarrollando la ecuación diferencial, quedará: Z = - 9.8 t + ko =
Remplazando la condición inicial, y.f
´
= 20, 20 = -9.8 (0) + Ko , luego Ko = 20, y la
'ö
'
ecuación para la primera derivada quedará: Z = - 9.8 t + 20 = o

'ö
'
quedará: Z = - 9.8 t + 20 = o
La ecuación anterior es una ecuación diferencial de primer orden, por lo tanto, procedemos
a desarrollarla.
),
)•
= - 9.8 t + 20 , separando variables: dy = (-9.8 t + 20) dt

•#
Integrando la ecuación: y = -9.8( ) + 20 t + Ko = - 4.9 t2 + 20 t + Ko

Remplazando la condición inicial y(0) = 0 , 0 = - 4.9 (0)2 + 20 (0) + Ko , luego, Ko = 0

La solución quedará:
y(t) = - 4.9 t2 + 20 t cmt.
Valores de prueba:
y.´´ y.´
Þ²•. Þ²•.
µ().# µ().
Para t = 2 seg = -9.8 , = 0.4 , y(2) = 20.4 cmt

[Link]
[Link] GRÁFICAS DE LA ECUACIÓN DIFERENCIAL, DE LAS ECUACIONES
PARA DERIVADAS INTERMEDIAS Y DE LA ECUACIÓN SOLUCIÓN
(FUNCIÓN)

A partir de la ecuación diferencial (ecuación de la segunda derivada), de las ecuaciones

para las derivadas intermedias y de la ecuación para la función (solución de la ecuación
diferencial), podremos dibujar los lugares geométricos correspondientes, utilizando el
software de Matlab, mediante el programa siguiente. (Grafejemplo1)
%El programa dibuja la ecuación diferencial y´´= - 9.8,
%la ecuación de la primera derivada y´ = - 9.8 t + 20 y la ecuación
% de la solución (Función) y = - 4.9 t^2 + 20 t
t = 0 : 0.001 : 5;
V = (t./t)-10.8;
figure(1),plot(t,V,'r'), grid on
title('y´´(t)= - 9.8')
ylabel('Ecuacion diferencial'),pause
W = -9.8.*t + 20;
figure(2),plot(t,W), grid on
title('y´(t)= -9.8.*t + 20')
ylabel('Primera derivada'),pause
Z = -4.9.*t.^2 + 20.*t ;
figure(3), plot(t,Z), grid on
title('y(t)= - 4.9.*t^2 + 20.*t')
ylabel('Función,Solución de E.D.'),pause
[t.',Z.',W.',V.']

[Link]
 y´(t)= -9.8.*t + 20
20

15

10

5 X: 2
Y: 0.4
Primera derivada

-5

-10

-15

-20

-25

-30
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

y´´(t)= - 9.8
-8.8

-9

-9.2

-9.4
Ecuacion diferencial

-9.6 X: 2
Y: -9.8
-9.8

-10

-10.2

-10.4

-10.6

-10.8
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

[Link]
 y(t)= - 4.9.*t2 + 20.*t
25

20
X: 2
Y: 20.4
15

10
Función,Solución de E.D.

-5

-10

-15

-20

-25
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

3.1.2 DESARROLLO NUMÉRICO.

A continuación se encuentra el programa de matlab por medio del cual encontraremos
la solución numérica de la ecuación diferencial.

(romarquez)
%Ejemplo.m soluciona la ecuac dif, d2y/dt2 = -9.8,
%para y(0)=0,y´(0)= 20
% Notaciones: dy(2)=d2y/dt2 , y(2)= y´ , dy(1)= dy/dt , y(1)= y
function Davila
tspan = [0:0.001:5];%Vector de la variable independiente
y0 = [0, 20]; %Condiciones iniciales
[T,Y] = ode45(@luis,[0:0.001: 5],[0 20]);%esta sentencia se puede
cambiar
% por ¨ode45(@luis,tspan,y0)
[T,Y]% imprime los valores
Z = [T./T]-10.8;
figure(4),plot(T,Z,'r'),grid on,title('Aceleración d2y/dt2') , pause
figure(1),plot(T,Y(:,2),'-'),title('Velocidad dy/dt'),xlabel('Tiempo
t')
,ylabel('dy/dt'),grid on, pause;
figure(2),plot(T,Y(:,1),'-'),title('Posición y')
xlabel('Tiempo t'),ylabel('Función y '),grid on,pause;
figure(3),plot(T,Y(:,1),'-',T,Y(:,2),'r'), grid on %Gráfica de y y y´
% El signo - y la r ,para que queden de distinto color
title('Velocidad y posición'); %Título del gráfico

[Link]
legend('y' , 'dy/dt');%Leyenda del gráfico
xlabel('Tiempo t');%Nombre del eje x
ylabel('y y dy/dt');%Nombre del eje y
hold on
function dy = luis(t,y)%Función que depende de y , t
dy = zeros(2,1); %Vector columna con el número de ecuaciones a resolver
dy(1) = y(2); % Ecuación donde decimos que y(2)es la derivada de y(1)
dy(2) = -9.8;%Ecuación diferencial a resolver

Algunos resultados del proceso numérico:

Tiempo Función y´ y´´
2.0000 20.4000 0.4000 -9.8000

4.0000 1.6000 -19.2000 -9.8000

Aceleración d2y/dt2
-8.8

-9

-9.2

-9.4
Ecuación diferencial

-9.6 X: 2
Y: -9.8
-9.8

-10

-10.2

-10.4

-10.6

-10.8
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

[Link]
 Velocidad dy/dt
20

15

10

5 X: 2
Primera derivada dy/dt

Y: 0.4
0

-5

-10

-15

-20

-25

-30
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Tiempo t

Posición y
25

20
X: 2
Y: 20.4
15
Solución de la E.D. (Función) y

10

-5

-10

-15

-20

-25
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Tiempo t

[Link]
 Velocidad y posición
30
y
dy/dt
20

10
X: 2
y y dy/dt

Y: 0.4
0

-10

-20

-30
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Tiempo t

2. ECUACIÓN DIFERENCIAL DE SEGUNDO ORDEN EN DONDE LA SEGUNDA DERIVADA ES

FUNCIÓN DE SOLO LA VARIABLE INDEPENDIENTE
'ç ö Þ²•.
' ç µ().#
Desarrollar la ecuación diferencial = 200 Sen(5t-152.46°) ,

para y(0) = 0.709 y.f

´
= -18.49

3.2.1 DESARROLLO ANALÍTICO.

La ecuación diferencial de segundo orden se puede convertir en dos ecuaciones diferenciales

de primer orden.
'ö )* 'ç ö
' )• ' ç
Hacemos Z = , entonces: =

Remplazando la nueva variable en la ecuación diferencial de 2do orden, tendremos:

)*
)•
= 200 Sen (5t - 152.46°)

'ö
'
Desarrollando la ecuación diferencial, quedará: Z = - 40 Cos(5t - 152.46°) + ko =

[Link]
Remplazando la condición inicial, y.f
´
= - 18.49, -18.49 = - 40 Cos(5(0) - 152.46°) + Ko ,

luego Ko = 0.4213, y la ecuación para la primera derivada quedará:

'ö
'
Z = - 40 Cos(5 t - 152.46°) + 0.4213 = o

La ecuación anterior es una ecuación diferencial de primer orden, por lo tanto, procedemos a

desarrollarla.

'ö
'
= - 40 Cos(5 t - 152.46°) - 53.958, separando variables:

dy = (- 40 Cos(5 t - 152.46°) - 53.958 ) dt , integrando la ecuación, la solución general

quedará:

y = - 8 Sen( 5 t – 152.46°) - 53.958 t + Ko

Remplazando la condición inicial y(0) = 0.709 ,

0.709 = - 8 Sen( 5(0) – 152.46°) - 53.958 (0) + Ko , arroja Ko = 3.3915

Por lo tanto, la solución quedará

[Link]
Valores de prueba:

y./ y./
´´ Þ²•. ´ Þ²•.
µ().# µ().
Para t = 1 seg , = 143.81 , = - 26.164 , y(1) = - 56.31 cmt

y.´´ y.´
Þ²•. Þ²•.
µ().# µ().
Para t = 2 seg , = 174.0702 , = - 73.658 , y(2) = - 111.486 cmt

RESULTADOS DE ECUACIONES POR EL MÉTODO ANALÍTICO:

[Link] GRÁFICAS DE LA ECUACIÓN DIFERENCIAL, DE LAS ECUACIONES

PARA DERIVADAS INTERMEDIAS Y DE LA ECUACIÓN SOLUCIÓN
(FUNCIÓN)
A partir de la ecuación diferencial (ecuación de la segunda derivada), de las
ecuaciones para las derivadas intermedias y de la ecuación para la función
(solución de la ecuación diferencial), podremos dibujar los lugares
geométricos correspondientes, utilizando el software de Matlab, mediante el
programa siguiente. (Grafejemplo2)
%El programa dibuja la ecuación diferencial y´´= 200.*sin(5t-152.46°),
%la ecuación de la primera derivada y´ = 200 Sen(5t-152.46°) y la
ecuación
% de la solucion(Función) y = = -8Sen(5 t –152.46°)-53.958 t+ 3.3915
t = 0 : 0.01 : 3;
V = 200.*sin(5.*t-2.6609);
figure(1),plot(t,V,'r'), grid on
title('y´´(t)= 200.*sin(5.*t-152.46°)')
ylabel('Ecuacion diferencial'),pause
W = -40.*cos(5.*t-2.6609)-53.958;
figure(2),plot(t,W), grid on
title('y´(t)= -40Cos(5t-152.46°)-53.958')
ylabel('Primera derivada'),pause
Z = -8.*sin(5.*t -2.6609)-53.958.*t+3.3915 ;
figure(3), plot(t,Z), grid on
title('y(t)= -8Sen(5t–152.46°)-53.958t+3.3915')
ylabel('Función,Solución de E.D.'),pause
[t.',Z.',W.',V.']

[Link]
 y´´(t)= 200.*sin(5.*t-152.46°)
200

150 X: 2
Y: 174.1

100
Ecuacion diferencial

50

-50

-100

-150

-200
0 0.5 1 1.5 2 2.5 3

y´(t)= -40Cos(5t-152.46°)-53.958
-10

-20

-30

-40
Primera derivada

-50

-60
X: 2
-70 Y: -73.66

-80

-90

-100
0 0.5 1 1.5 2 2.5 3

[Link]
 y(t)= -8Sen(5t–152.46°)-53.958t+3.3915
20

-20
Función,Solución de E.D.

-40

-60

-80

X: 2
-100 Y: -111.5

-120

-140

-160
0 0.5 1 1.5 2 2.5 3

3.2.2 DESARROLLO NUMÉRICO.

A continuación se encuentra el programa de matlab por medio del cual encontraremos
la)solución
#
Þ²•. diferencial.
y.´ f
, numérica de la ecuación
• # µ() .#
) = 200 Sen(5t-152.46°) , para y(0) = 0.709 = -18.49
(Nuevo.m)
%El programa [Link] la ecuación diferencial
%y´´= 200 Sen(5t-2.661),para: yo=7.09,y1 = -18.49
% Notaciones: dy(2)= d2y/dt2 = Z , dy(1)= dy/dt = Y(2) , Función y =
Y(1)
function nuevo
tspan = 0 : 0.01: 3;,
y0 = [7.09 ; -18.49];
[T,Y] = ode45(@luis,[0:0.01: 3],[7.09 -18.49]);
Z= 200.*sin(5.*T-2.661);
figure(1),plot(T,Z,'r'),grid on,title('Ecuación diferencial a
resolver')
ylabel('Segunda derivada y´´(t)'), pause
[' Tiempo ' , ' Función ' , ' y´ ' , ' y´´ ']
[T,Y, Z]
figure(2),plot(T,Y(:,2)), grid on, title('Primera derivada')
ylabel('y´(t) 0 dy/dt'), pause
figure(3),plot(T,Y(:,1)), grid on, title('solución de la E.D.,función
y')

[Link]
ylabel('Función y(t)'), pause
function dy = luis(t,y)%Función que depende de y , t
dy = zeros(2,1); %Vector columna con el número de ecuaciones a
resolver
dy(1) = y(2); % Ecuación donde decimos que y(2)es la derivada de y(1)
dy(2) = [ 200.* sin(5.*t-2.661)];%Ecuación diferencial a resolver

Algunos resultados del proceso numérico:

Tiempo Función y´ y´´

1.0000 -56.3169 -26.1513 143.8320

2.0000 -111.4862 -73.6729 174.0604

Ecuación diferencial a resolver

200

150 X: 2
Y: 174.1

100
Segunda derivada y´´(t)

50

-50

-100

-150

-200
0 0.5 1 1.5 2 2.5 3

[Link]
 Primera derivada
-10

-20

-30

-40
y´(t) 0 dy/dt

-50

-60
X: 2
-70 Y: -73.67

-80

-90

-100
0 0.5 1 1.5 2 2.5 3

solución de la E.D.,función y
20

-20

-40
Función y(t)

-60

-80

X: 2
-100 Y: -111.5

-120

-140

-160
0 0.5 1 1.5 2 2.5 3

[Link]
3.3. ECUACIÓN DIFERENCIAL DE SEGUNDO ORDEN EN DONDE LA
SEGUNDA DERIVADA ES FUNCIÓN DE SOLO LA VARIABLE
DEPENDIENTE O DE LA FUNCIÓN
'ç ö
' ç
Desarrollar el problema de valor inicial: = 2 y , para t = 0, y = 1, y´ = 1

3.3.1 DESARROLLO ANALÍTICO:

), '+ 'ç ö '+ 'ö '+

)• ' ' ç = 'ö ' 'ö
Hacemos z = , entonces = = z , luego los reemplazos serán:

), )# , ),
)• )•# ),
=z , = z

Reemplazando estas expresiones en la ecuación diferencial tendremos:

), ),
), ),
z -2 y = 0 ó z = 2 y, la cual es una ecuación diferencial de primer orden.

)# ,
)•#
'ö
Por lo anterior, la ecuación diferencial de 2do orden - 2 y = 0, quedó reducida a dos ecuaciones diferenciales de primer orden:
),
), '
(1) z = 2y , (2) = z, las cuales habrá que desarrollarlas en el orden numérico presentado.

),
),
Desarrollando la E. D. (1) : z = 2 y , tendremos.

Separando variables la ecuación (1) quedará: z dz = 2 y dy

Integrando indefinidamente a ambos lados de la ecuación: ? z dz = ? 2 y dy + ko

,# 'ö

'
2 2 2

-õ
Desarrollando ambos lados de la ecuación inmediatamente anterior, tendremos: = y + ko ó z = 2 y + k, ó =

ºç ö ç , ó z = º2 y ç k/

º2 y k/
), ),
)• )•
Desarrollando la E. D. (2): = z; = , tendremos:

),
º ,ç >™m
Separando variables, la ecuación (2) quedará: = dt

? dt + ko
),
?º
,ç >™m
Integrando indefinidamente a ambos lados de la ecuación: = ó

),
?º
,ç >™m
= t + ko (A)

Desarrollando por aparte el lado izquierdo de la ecuación (A), tendremos

),
?º 2 2
,ç >™m
= hacemos: 2 y = u , o sea que √2 y = u, √2 dy = du,

[Link]
reemplazando en la E. D.

=
), / )(
?º √
? º(# >™
,ç >™m m
, integral que se puede resolver utilizando la fórmula n° 9 del formulario y que indica lo siguiente:

= ln | x √ a | + c , en donde, u = x , a = ºk/
)-
(9) ? √-# > ‰#

ln| √2 y + º2 y k/ | + C
), / )( /
Por lo tanto, ? ? º(# >™ ç
º ,ç >™ m √ m √
= =

Reemplazando la solución de la integral en la ecuación (A), ésta quedará:

ln| √2 y + º2 y k/ | =
/ ç
√
t + k2, en donde la constante C se incluye dentro de la constante ko y ésta pasa a ser k2

Simplificando la ecuación se puede convertir a:

√2 y + º2 y ç k/ = .√ • > ™#
= ke
√ •

Por lo tanto la solución general implícita es:

Evaluación de las constantes:

Reemplazando las condiciones iniciales, t = 0, y = 1, en la ecuación (3), tendremos:

√2 (1) + º2 .1 ç k/ = ke
√ .f
, por lo tanto º2 k/ = k - √2 (A)

Reemplazando las condiciones iniciales, t = 0, y´ = 1, en la ecuación (2), tendremos:

1 = º2 .1 k/ ; 1 = 2 + k1 , por lo tanto
k1 = - 1

Sí regresamos a la ecuación (A) y reemplazamos el valor de k1, encontraremos que √2 1 = k - √2 , por lo tanto,

Solución específica:

Reemplazando los valores de las constantes encontrados, la solución específica implícita de la ecuación diferencial quedará:

Simplificando la expresión: º2 y 1 = (1 + √2 e√ - √2 y
•

Elevando al cuadrado ambos miembros de la ecuación:

[Link]
(2 y 1 = (1+√2 ) e - 2(1+√2 ) √2 e y + 2 y , que simplificando quedará:
2 √ • √ •

1 = (3+2√2 ) e √ •
- (4+2√2 ) e
√ •
y despejando la función, esta quedará:

= ( )e + ( )e
/>. > √ (#√# « 5√ 5√ • >√ √ •
.Ž> √ (√# « Ž Ž
y=

La función explícita quedará: y = ( ). + ( ).

ç5√ç 5√ç ç>√ç √ç
o

RESULTADOS DE ECUACIONES POR EL MÉTODO ANALÍTICO

[Link] GRÁFICAS DE LA ECUACIÓN DIFERENCIAL, DE LAS ECUACIONES

PARA DERIVADAS INTERMEDIAS Y DE LA ECUACIÓN SOLUCIÓN
(FUNCIÓN)

A partir de la ecuación diferencial (ecuación de la segunda derivada), de las ecuaciones

para las derivadas intermedias y de la ecuación para la función (solución de la ecuación
diferencial), podremos dibujar los lugares geométricos correspondientes, utilizando el
software de Matlab, mediante el programa siguiente. (Grafejemplo3)

%El programa dibuja la ecuación diferencial y´´= 2 y,

%la ecuación de la primera derivada y´ = (2*y^2-1)^(0.5) y la ecuación
% de la solucion(Función)
% y = 0.14644.*exp(-1.4142.*t)+0.85355.*exp(1.4142.*t)
t = 0 : 0.01 : 3;
Z = 0.14644.*exp(-1.4142.*t)+0.85355.*exp(1.4142.*t);
figure(1), plot(t,Z), grid on
title('y(t)= 0.14644.*exp(-1.4142.*t)+0.85355.*exp(1.4142.*t)')
ylabel('Función,Solución de E.D.'),xlabel('Tiempo t'),pause

[Link]
W = (2.*Z.^2-1).^(0.5);
figure(2),plot(t,W), grid on, title('y´(t)= (2*y^2-1)^(0.5)')
ylabel('Ecuación primera derivada, y´(t)'),xlabel('Tiempo t'),pause
V = 2.*Z;
figure(3),plot(t,V,'r'), grid on, title('y´´(t)= 2y')
ylabel('Ecuacion diferencial a resolver, d2y/dt2'),pause
[t.',Z.',W.',V.']

[Link]
 y´(t)= (2*y2-1)(0.5)
90

80

70
Ecuación primera derivada, y´(t)

60

50

40

30 X: 2
Y: 20.41
20

10

0
0 0.5 1 1.5 2 2.5 3
Tiempo t

[Link]
 y(t)= 0.14644.*exp(-1.4142.*t)+0.85355.*exp(1.4142.*t)
60

50
Función,Solución de E.D.

40

30

20 X: 2
Y: 14.45

10

0
0 0.5 1 1.5 2 2.5 3
Tiempo t
3.3.2
DESARROLLO NUMÉRICO.
A continuación se encuentra el programa de matlab por medio del cual encontraremos la
solución numérica de la ecuación diferencial.

)# , Þ²•.
)•#
= 2y µ().#
, para y(0) = 1.0 y.f
´
= 1.0

(probvalorinicial.m)
%El programa resuelve la ecuación diferencial y´´= 2*y,
%yo = 1,y´ = 1
% Notaciones: dy(2)= d2y/dt2 = Z , dy(1)= dy/dt = Y(2),Función y = Y(1)
function nuevo
tspan = 0 : 0.01: 3;,
y0 = [1 ; 1];
[T,Y] = ode45(@luis,[0:0.01: 3],[1 1]);
Z= 2.*Y(:,1);
figure(1),plot(T,Z,'r'),grid on,ylabel('Segunda derivada y´´(t)')
title('Ecuación diferencial a resolver'),xlabel('Tiempo t'), pause
[T,Y, Z]
figure(2),plot(T,Y(:,2)), grid on, title('Primera derivada de la
función')
ylabel('dy/dt = y´(t´'),xlabel('Tiempo t'), pause
figure(3),plot(T,Y(:,1)), grid on, title('Solución de la E.D., función')
ylabel('Función y(t)') , pause
function dy = luis(t,y)%Función que depende de y , t

[Link]
dy = zeros(2,1); %Vector columna con el número de ecuaciones a resolver
dy(1) = y(2); % Ecuación donde decimos que y(2)es la derivada de y(1)
dy(2) = [ 2*y(1)];%Ecuación diferencial a resolver

[Link]
 Ecuación diferencial a resolver
120

100
Segunda derivada y´´(t)

80

60

40 X: 2
Y: 28.9

20

0
0 0.5 1 1.5 2 2.5 3
Tiempo t

Primera derivada de la función

90

80

70

60
dy/dt = y´(t´

50

40

30 X: 2
Y: 20.41
20

10

0
0 0.5 1 1.5 2 2.5 3
Tiempo t

[Link]
 Solución de la E.D., función
60

50

40
Función y(t)

30

20 X: 2
Y: 14.45

10

0
0 0.5 1 1.5 2 2.5 3

Algunos resultados del proceso numérico:

Tiempo Función y´ y´´

1.0000 3.5465 4.9148 7.0930

2.0000 14.449 20.41 28.898

3.4. ECUACIÓN DIFERENCIAL DE SEGUNDO ORDEN EN DONDE LA

SEGUNDA DERIVADA ES FUNCIÓN DE SOLO LA PRIMERA DERIVADA
DE LA FUNCIÓN
'ç ö 'ö Þ²•.
Desarrollar la ecuación diferencial: ' ç
= - 9.8 – 0.035 '
( µ().# ), para y(0) = 0 y.f
´
= 10

3.4.1 DESARROLLO ANALÍTICO.

La ecuación diferencial de segundo orden se puede convertir en dos ecuaciones diferenciales

de primer orden.

[Link]
 'ö )* 'ç ö
' )• ' ç
Hacemos Z = , entonces: =

Remplazando la nueva variable en la ecuación diferencial de 2do orden, tendremos:

)*
)•
= - 9.8 – 0.035 Z = - 0.035(Z + 280)

'ö
'
Desarrollando la ecuación diferencial, quedará: Z = ko e- 0.035 t - 280 =

Remplazando la condición inicial, y.f

´
= 10, 10 = ko e- 0.035 (0) - 280, luego, Ko = 290, y la

'ö
'
ecuación para la primera derivada quedará: Z = 290 e- 0.035 t - 280 = o

La ecuación anterior es una ecuación diferencial de primer orden, por lo tanto, procedemos a

desarrollarla.

),
)•
= 290 e- 0.035 t - 280, separando variables:

dy = (290 e- 0.035 t - 280) dt , integrando la ecuación, la solución general quedará:

y = - 8285.71 e- 0.035 t - 280 t + Ko , remplazando la condición inicial y(0) = 0,

0 = - 8285.71 e- 0.035 (0) - 280 (0) + Ko , arroja Ko = 8285.71. Por lo tanto, la solución quedará

[Link]
Valores de prueba:

Þ²•. Þ²•.
Para t = 1 seg , y./
´´
= - 9.8009 µ().#
, y./
´
= 0.0255 µ().
, y(1) = 4.9835 cmt

Þ²•. Þ²•.
y.´´ = -9.4638 µ().# , y.´ µ().
Para t = 2 seg , = - 9.605 , y(2) = 0.1652 cmt

RESULTADOS DE ECUACIONES POR EL MÉTODO ANALÍTICO

[Link] GRÁFICAS DE LA ECUACIÓN DIFERENCIAL, DE LAS ECUACIONES

PARA DERIVADAS INTERMEDIAS Y DE LA ECUACIÓN SOLUCIÓN
(FUNCIÓN)

A partir de la ecuación diferencial (ecuación de la segunda derivada), de las ecuaciones

para las derivadas intermedias y de la ecuación para la función (solución de la ecuación
diferencial), podremos dibujar los lugares geométricos correspondientes, utilizando el
software de Matlab, mediante el programa siguiente. (Grafejemplo4)
%El programa dibuja la ecuación diferencial y´´= -9.8-0.035 dy/dt,
%la ecuación de la primera derivada y´ = 290exp(-0.035 t)-280,y la
ecuación
% de la solucion(Función), y = -8285.71exp(-0.035t)-280t+8285.71
t = 0 : 0.01 : 3;
Z = -8285.71.*exp(-0.035.*t)- 280.*t+8285.71;
figure(1), plot(t,Z), grid on
title('y(t)= -8285.71 exp(-0.035 t)- 280 t + 8285.71')
ylabel('Ecuación solución de E.D., función y'),pause
W = 290.*exp(-0.035.* t)-280;
figure(2),plot(t,W), grid on
title('y´(t)= 290 exp(-0.035 t)-280 ')
ylabel('Ecuación primera derivada'),pause
V = -9.8-0.035.*W;
figure(3),plot(t,V,'r'), grid on

[Link]
title('y´´(t)= -9.8-0.035 dy/dt')
ylabel('Ecuacion diferencial a resolver'),pause
[t.',Z.',W.',V.']
y´´(t)= -9.8-0.035 dy/dt
-8.8

-9
Ecuacion diferencial a resolver

-9.2

X: 2
-9.4 Y: -9.464

-9.6

-9.8

-10

-10.2
0 0.5 1 1.5 2 2.5 3

y´(t)= 290 exp(-0.035 t)-280

10

5
Ecuación primera derivada

-5
X: 2
Y: -9.606

-10

-15

-20
0 0.5 1 1.5 2 2.5 3

[Link]
 y(t)= -8285.71 exp(-0.035 t)- 280 t + 8285.71
5

X: 2
Y: 0.1652
Ecuación solución de E.D., función y

-5

-10

-15
0 0.5 1 1.5 2 2.5 3

3.4.2 DESARROLLO NUMÉRICO.

A continuación se encuentra el programa de matlab por medio del cual encontraremos la
'ç ö î/ .
numérica de'öla ecuación
ö.´ñ
solución diferencial.
'
ç ð.0.ç '
= - 9.8 – 0.035 ( ) , para y(0) = 0 = 10
(bolalanzada.m)
%El programa resuelve la ecuación diferencial y´´= -9.8 - 0.035 dy/dt,
%yo = 0, dy/dt = 10
% Notaciones: dy(2)= d2y/dt2 = Z , dy(1)= dy/dt = Y(2) , Función y =
Y(1)
function nuevo
tspan = 0 : 0.001: 3;,
y0 = [0 ; 10];
[T,Y] = ode45(@luis,[0:0.001: 3],[0 10]);
Z= -9.8-0.035.*Y(:,2);
figure(1),plot(T,Z),grid on,title('Aceleración de la bola')
ylabel('Ecuación diferencial a resolver y´´(t)'),xlabel('Tiempo t'),
pause,[T,Y,Z]
figure(2) , plot(T,Y(:,2)), grid on, title('Velocidad de la bola')
ylabel('primera derivada y´(t)'),xlabel('Tiempo t'), pause
figure(3) , plot(T,Y(:,1)), grid on, title('Posición de la bola')
ylabel('Ecuación solución de la E.D. y(t)'),xlabel('Tiempo t'), pause
function dy = luis(t,y)%Función que depende de y , t
dy = zeros(2,1); %Vector columna con el número de ecuaciones a resolver
dy(1) = y(2); % Ecuación donde decimos que y(2)es la derivada de y(1)

[Link]
dy(2) = [ -9.8-0.035.*y(2)];%Ecuación diferencial a resolver

Aceleración de la bola
-8.8

-9
Ecuación diferencial a resolver y´´(t)

-9.2

-9.4

-9.6
X: 1
Y: -9.801
-9.8

-10

-10.2
0 0.5 1 1.5 2 2.5 3
Tiempo t

Velocidad de la bola
10

5
X: 1
Y: 0.02557
0
primera derivada y´(t)

-5

-10

-15

-20
0 0.5 1 1.5 2 2.5 3
Tiempo t

[Link]
 Posición de la bola
5
X: 1
Y: 4.984

0
Ecuación solución de la E.D. y(t)

-5

-10

-15
0 0.5 1 1.5 2 2.5 3
Tiempo t

Algunos resultados del proceso numérico:

Tiempo Función y´ y´´

1.0000 4.9837 0.0256 -9.8009

2.0000 0.1655 -9.6058 -9.4638

3.5. ECUACIÓN DIFERENCIAL DE SEGUNDO ORDEN EN DONDE LA SEGUNDA

DERIVADA ES FUNCIÓN DE LA VARIABLE INDEPENDIENTE Y DE LA FUNCIÓN

Desarrollar el problema de valor inicial siguiente:

'ç ö
' ç
= tan(t) + y , para y(0) = 1.0 , ö.ñ
´
= 1.0

3.5.1 DESARROLLO ANALÍTICO

PROCESO ANALÍTICO: MÉTODO GENERAL

[Link]
Ecuación diferencial lineal no-homogénea de coeficientes constantes, en donde el término
independiente es trigonométrico, pero no está dentro de la lista para ser desarrollado por el
método de los coeficientes indeterminados, por lo tanto, se procede a utilizar el método
general.

y(t) = k1 cos(t) + k2 sen(t) + cos(t) Ln[ sec(t) + tan(t)]

Determinando la primera derivada de la función:

îïð .
ö.´ = - k1 sen(t) + k2 cos(t) + sen(t) Ln[ sec(t) + tan(t)] - .ð.î. > 12.
[ sec(t) tan(t) + seî ç . ]

Þßµ.•
y.•
´
= - k1 sen(t) + k2 cos(t) + sen(t) Ln[ sec(t) + tan(t)] – .µ(Þ.• >•‰•.•
sec(t) (tan(t) + sec(t))

y.•
´
= - k1 sen(t) + k2 cos(t) + sen(t) Ln[ sec(t) + tan(t)] – 1

CAPÍTULO III
ECUACIONES LINEALES DE SEGUNDO ORDEN

[Link]
9. SERIES DE POTENCIAS

Hasta el momento hemos visto cómo solucionar la ecuación

´´ D ´ Å

Donde D y son constantes y Å es una función continua en u. Ahora trataremos

de solucionar la ecuación

L. ´´ L/ . ´ Lf . 0

Donde L/ , L , Lf son funciones analíticas, es decir se pueden expresar como una

serie de potencias. La solución a esta ecuación involucra series de potencias, es
decir de la forma
4

. 3 D1 . f
1

1 f

O de una forma más general:

4

. 3 D1 . f
1
. 5.
1 f

Donde 5. es una función en .

Las mismas ideas se pueden aplicar a las ecuaciones de otros órdenes, y también
directamente a ecuaciones no-homogéneas.

de las series de potenciales. Para cada serie de potencias ∑4 1 f D1 .

Antes de empezar este análisis, es necesario repasar las propiedades pertinentes
1
f , existe
un valor Å, para Å un real no negativo fijo denominado radio de convergencia, tal
que cuando | f | < Å , la serie converge y cuando | f | > Å, ella diverge.

converge, es decir ∑4 1 f D1 .
Decimos que la serie converge absolutamente si la serie de valores absolutos
1
f converge para el mismo radio.

La convergencia absoluta siempre implica la convergencia ordinaria.

Las series de potencias son simplemente una generalización de los polinomios, y

muchas veces se pueden manipular de manera análoga. Las series se pueden
sumar, restar, multiplicar, y si el denominado no es igual a cero, dividir. La

[Link]
derivación y la integración de funciones representadas por series de potencias se
permiten término a término dentro del radio de convergencia.

analíticas, es decir, de la forma ∑4

1 f D1 . , con un radio de convergencia Å.
Anteriormente dijimos que los coeficientes de la ecuación eran funciones
1
f

En este caso es fácil verificar que D1

y .0 ."
1 !
, y de hecho esta igualdad se puede
utilizar para construir la serie ya que viene dada en base a la derivada de la
función que es analítica.

Utilizaremos varios criterios para analizar la convergencia de nuestras series de

potencias. Lo más importante son los de la razón, de comparación, y series
geométricas. Se encuentran también los criterios de la raíz y de la integral.

Nuestro interés es una ecuación diferencial lineal homogénea.

.1
L/ . .15/
⋯ L1 . 0

Donde L/ , L , … , L1 , son coeficientes analíticos con radio de convergencia Å,

, , … ,
puede demostrarse que existen soluciones linealmente independientes
/ 1 , cada una de las cuales es analítica con el mismo radio de
convergencia.

Demostraremos este teorema para ecuaciones diferenciales de orden dos.

Teorema 12: sean L/ , L , funciones analíticas con un radio de convergencia Å,

tales como:
4

L/ . 3 1. f
1

1 f

L . 3 1. f
1

1 f

Entonces la ecuación diferencial

´´ L/ . ´ L . 0 (1)

Tiene dos soluciones linealmente independientes / y que son analíticas con el

mismo radio de convergencia.

[Link]
 ∑41 f D1 . f
1

entro de un radio dado. Hallemos los coeficientes D1 , para que se satisfaga

Demostración: sea solución de la ecuación (1), convergente

nuestra condición.

Como es solución de la ecuación, entonces la verifica, luego derivando una

dos veces tenemos:
4 4

´ 3 D1 . f
15/
3. 1 D 1>/ . f
1

1 / 1 f

4 4

´´ 3 . 1 D1 . f
15
3. 2 . 1 D 1> . f
1

1 1 f

Los productos L/ . ´yL. vienen dados por las series de potencias:

4 1

L/ . ´ 3 73.: 1 D ;>/ 15; 8 . f

1

1 f ; f

4 1

L. 3 73 D ; 15; 8 . f
1

1 f ; f

Sustituyendo en la ecuación (1), tenemos:

4 1

3 9. 2 . 1 D1> 3[.: 1 D;>/ 15; D; 15; ]: . f

1
0
1 f ; f

Esto se tiene si
1

. 2 . 1 D1> 3[.: 1 D;>/ 15; D; 15; ]

; f

0, 1, 2, …, esta fórmula expresa D1> en función de los coeficientes

anteriores Df, D/, . . , D1>/ , y los coeficientes de las funciones dadas L/ y L .
Para

/ un punto fijo con / ≠ f , en el intervalo . f Å, f Å y sea M

| / f |, como las series que definen L/ y L convergen absolutamente para
Sea

/ los términos de esas series están acotados, es decir

,

[Link]
 | ; | M
;
≤ </ y | ; | M ; ≤ <

Para ciertos </ y < positivos. Sea < el mayor de </ y M< entonces

| ;| ≤ y | ;| ≤
â â
J< J <m

Luego de la fórmula de recurrencia tenemos:

1
< <
. 2 . 1 |D1> | ≤ 3 =.: 1 |D;>/ | |D; | >
M 15; M 15;>/
; f

1 1
<
93 .: 1 |D;>/ |M ;>/ 3|D;>/ |M ;>/ |Df | |D1>/ |M 1>/ :
M 1>/
1 f ; f

1
<
93 .: 2 |D;>/ |M ;>/ |Df |:
M 1>/
; f

1>/
<
3 .: 1 |D; |M ;
M 1/
; f

Hagamos Cf Df , C/ D/ y definimos C C , …, mediante la fórmula de

recurrencia:

. 2 . 1 C1> ∑1 /
; f .: 1 C; M ; para ≥0
â
J 0 m
(2)

Entonces |D1 | ≤ C1 para todo mayor o igual a cero.

Luego

∑4
1 f D1 . f
1
es menor que ∑4
1 f C1 | f|
1

De (2) tenemos que:

. 2 . 1 C1> ∑1>/
; f .: 1 C; M ;
â
J 0@m
(3)

1
<
. 1 C1>/ 3 .: 1 C; M ;
M1
; f

[Link]
 M 15/ . 1 C1>/ ∑1; f .: 1 C; M ;
â
J 0 m
(4)

De (3) restamos (4), y obtenemos:

. 2 . 1 C1> M 5/ . 1 C1>/ <. 2 C1>/

Luego

. 1 . 2 <M
C1> C1>/
. 2 . 1 M

Hallemos

C1> | f| . 1 . 2 <M | f|
1>
| f| ⟶ D Ë ⟶ ∞
C1>/ | f |1>/ . 2 . 1 M M

Este límite es menor que 1 si | f| < M, luego ∑41 f D1 . f

1

| f | < M. Pero como M | / |

converge si
f /
. f Å, f Å .
y es un punto arbitrario del intervalo

Luego la serie ∑4
1 f D1 . f
1
converge para todo en . f Å, f Å .

Como los coeficientes de ∑4 1 f D1 . f

1
aparecen en función de Df y D/ sea /
la serie de potencias solución tal que Df 1 y D/ 0, entonces:

/. f 1 y ´´/ . f 0

Y una solución con Df 0 y D/ 1 entonces

. f 0 y ´ . f 1

Entonces las soluciones / y serían funciones analíticas linealmente

independientes.

Hallar la solución de ´´ 0

Sea ∑4
1 f D1
1
solución a la ecuación dada, en este caso tomamos f 0

Entonces

[Link]
 4 4

´´ 3 . 1 D1 15
3. 2 . 1 D1> 1

1 1 f

Y
4 4

3 D1 1>/
3 D15/ 1

1 f 1 /

Luego
4 4

´´ 3. 2 . 1 D1 2 "0
3 D15/ 1

1 f 1 /

2D 3[. 2 . 1 D1> D15/ ] 1

0
1 /

Lo cual es cierto si

D 0 y . 2 . 1 D1> D15/ 0

Si 0, 1, 2, …

Luego

D1>
p0 m
1, 2, 3 ….
.1> .1>/
para

De donde podemos deducir que:

D
po
1 1, 2, 3 ….
. .6 .7 …. 15/ 1
para

D
pm
1>/ 1, 2, 3 ….
.Ž .7 .› …. 1 . 1>/
para

D 1> 0 para 0, 1, 2 ….

De esta manera vemos que todos los coeficientes D1 están determinados en base
a Df , y D/ .

Agrupando todos los términos que contienen como factores a Df y D/ , obtenemos:

[Link]
 "?
Df P1 … Q D/ P … Q
"’ "© "¨
. 7 . 6 . . Ž . › . 7 . Ž .

De donde podemos decir que

4

/. 1 3
2 . 3 . 5 . 6 … .3 1 3
1 /

4 1>/
. 3
3 . 4 . 6 . 7 … .3 .3 1
1 /

Veamos que / es convergente, empleando el criterio del cociente; si expresamos

4

/. 1 3 1 .
1 /

Tenemos que:

1 1. 2 . 3 . 5 . 6 … .3 1 .3
1>
ý ý
1. 2 . 3 . 5 . 6 … .3 1 .3 .3 2 .3 3 . 1

→0
|"|’
→∞
. 1> . 1>
si

Luego / converge para todo real.

Dejamos al lector la demostración de la convergencia de y el hecho que / y

son linealmente independientes.
10. ECUACIÓN DE LEGENDRE

Algunas de las ecuaciones diferenciables importantes que surjan en problemas

físicos, son ecuaciones lineales de segundo orden con coeficientes analíticos. Una
de ellas es la ecuación de Legendre.

K. .1 ´´ 2 ´ D.D 1 0

Donde D es una constante. Si esta ecuación la escribimos:

2 D.D 1
´´ ´ 0
1 1

[Link]
Vemos que las funciones /5" # y 0, ya que tienen un
5 " p.p>/
/5" #
son analíticas en
desarrollo en serie de potencias y convergen para | | < 1.

Nuestro interés es hallar una base para las soluciones a esta ecuación, sea
4

3 D1 1

1 f

Solución a la ecuación de Legendre, luego:

4

´ 3 D1 15/

1 /

4 4

2 ´ 3 2 D1 1
3 2 D1 1

1 / 1 f

4 4

´´ 3 . 1 D1 15
3. 2 . 1 D1> 1

1 1 f

3 . 1 D1 1

1 f

Entonces
4

K. 3[. 2 . 1 D1> . 1 D1 2 D1 D.D 1 D1 ] 1

1 f

3[. 2 . 1 D1> .D 1 .D D1 ] 1
0
1 f

Luego los coeficientes de todas las potencias de , deben ser iguales a cero, es
decir

. 2 . 1 D1> .D 1 .D 1 D1 0 para 0, 1, …

[Link]
 D1> D1
.p>1>/ .p51
0, 1, …
.1> .1>/
para

Al desarrollar esta fórmula recurrente, tenemos que

. 1 1 .D 2 1 .D 2 3 … .D 1 .D .D 2 … .D 2 2
D Df
1
.2

. 1 1 .D 2 .D 2 2 … .D 2 .D 1 .D 3 … .D 2 1
D D/
1>/
.2 1

Todos los coeficientes pueden determinarse a partir de Df y D/

Entonces

. Df /. D/ .

donde

.D 1 .D 3 .D 1 D.D 2
/. 1 Ž
⋯
2 4
4
. 1 1 .D 2 1 .D 2 3 … .D 1 .D .D 2 … .D 2 2
1 3 1
.2
1 /

.D 2 .D 1 .D 4 .D 2 .D 1 .D 3
. 6
⋯
.3 .5
4
. 1 1 .D 2 .D 2 2 … .D 2 .D 1 .D 3 … .D 2 1
. 3 1>/
.2 1
1 /

Las series / y son soluciones de la ecuación de Legrende, que corresponden

a las siguientes constantes:

Df 1 D/ 0

Df 0 D/ 1

Respectivamente, estas funciones forman una base para las soluciones, ya que

/ .0 1 .0 0

[Link]
 ´/ .0 0 ´ .0 1

Luego

[ /, ] ≠ 0

Es importante notar que si D es un entero par negativo, entonces / solamente

tiene un número finito de términos diferentes de cero. Más aún, en este caso /
es un polinomio de grado que contiene solamente potencias pares de .

Por ejemplo:

/. 1 para D 0

/. 1 3 para D 2

/. 1 10 para D 4
6 Ž

coeficientes se anula. Ocurre una situación similar cuando D es un entero impar

Y la solución no es un polinomio en este caso, ya que ninguno de sus

no negativo.

Se deja al lector la demostración de que / y son convergentes para | | < 1 .

11. ECUACIONES LINEALES CON PUNTOS SINGULARES REGULARES

Dada la ecuación diferencial:

K. Lf . .1
L/ . .15/
⋯ L1 . 0

Donde Lf , L/ , L , … , L1 son analíticos en algún punto f y vamos a examinar el

caso en el cual Lf . f 0, dicho punto se llama un punto singular de la ecuación
diferencial dada; trataremos de resolver el problema alrededor de una vecindad de
f , entonces a dicho punto lo llamaremos punto singular regular. Ilustremos el
método con un ejemplo.

Hallemos la solución de K. ´´ ´ 0, la cual tiene un punto

singular regular en el origen. Restringimos a > 0, solamente.

Supongamos que la solución a la ecuación es de la forma

[Link]
 4
Î
3 D1 1

1 f

Con Df ≠ 0. Esta idea sencilla es efectiva, veamos. Hacemos formalmente las

operaciones y buscamos las condiciones que deben ser satisfechas por Ô, Df , D/ , …
para que esta función sea solución a la ecuación dada, tenemos que:
4

´ 3. Ô D1 1>Î5/

1 f

´´ 3. Ô . Ô 1 D1 1>Î5

1 f

´´ 3. Ô . Ô 1 D1 1>Î

1 f

4
3 3
´ 3 . Ô D1 1>Î
2 2
1 f

3 D1 1>Î>/

1 f

Reemplazando en la ecuación original tenemos:

4
3
K. 3 =³. Ô . Ô 1 . Ô ´ 1>Î 1>Î>/
D1 >
2
1 f

4
3 3
³Ô.Ô 1 Ô´ Df Î
3 =³. Ô . Ô 1 . Ô ´ D1 D15/ > 1>Î
2 2
1 /

4
1 1
Ô Ô D Î
3 ³. Ô Ô D D15/ ´ 1>Î
0
2 f 2 1
1 /

Luego

Ô $Ô % 0 y . Ô $ Ô %D D15/ 0 1, 2, …
/ /
para

[Link]
Sea

5.Ô Ô $Ô % y 5.Ô . Ô $ Ô % 1, 2, …
/ /
para

Luego K. 0 se verifica si

5.Ô 0 y 5.Ô D1 D15/ 0 para 1, 2, …

Al polinomio 5 se llama polinomio indicial de la ecuación diferencial dada. Además

es el coeficiente de la mínima potencia de que aparece en K. , y vemos que
sus raíces son los únicos valores posibles de Ô para los cuales existen soluciones
de la forma
4
Î
3 D1 1

1 f

En nuestro caso las raíces de 5 son Ô/ 0, Ô

/

Como 5.Ô È D1 D15/ 0 para 1, 2, …

Entonces

D1 1, 2, ….
5p0 m
C.Î>1
para

Si

Ô/ 0, 5.Ô/ 5. ≠0 para 1, 2, ….

Como la otra raíz de 5 es Ô

/
, en forma similar, tenemos

5.Ô 5$ %≠0 1, 2, ….
/
para

Haciendo Df 1 y Ô Ô/ 0 tenemos a / de la forma:

4
. 1 1 1
/. 1 3
5. 5. 1 … 5.1
1 /

Y haciendo Df 1 y Ô Ô
/
obtenemos a , otra solución:

[Link]
 4
/ / . 1 1 1
. 5 5
3
1 3 1
1 /5$ %5$ %…5$ %
2 2 2

Dejamos al lector la prueba de que / y convergen para todo real y que son
linealmente independientes, es decir, forman una base del espacio solución de la
ecuación:

3
´´ ´ 0
2

El siguiente teorema justifica el método utilizado anteriormente.

Teorema 13: consideremos la ecuación:

´´ L/ . ´ L. 0 (1)

Donde L/ y L son funciones analíticas, cuya series de potencia, tienen un radio de

convergencia Å. Sean Ô/ y Ô las raíces del polinomio indicial, tal que Åq .Ô/ ≥
Åq .Ô .

5.Ô Ô.Ô 1 L/ .0 L .0

Para Å, existe una solución / de la forma:

4

/.
Î
1 3 D1 1
.Df 1
1 f

Si Ô/ Ô no es nulo, ni un entero positivo, entonces, existe una segunda solución

para Å , de la forma:
4

. Î
2 3 DD1 1
.DDf 1
1 f

Donde ambas series convergen para Å.

Los coeficientes D1 , D1 pueden obtenerse, mediante sustitución de las soluciones

en la ecuación diferencial.

Demostración: supongamos que existe una solución de la forma:

[Link]
 4

. Î
3 D1 1
.DDf ≠ 0
1 f

Para la ecuación K. ´´ L/ . ´ L. 0

Donde L/ . ∑4
1 f 1
1
y L. ∑4
1 f 1
1
para Å

Entonces
4

´. Î5/
3. Ô D1 1

1 f

´´. Î5
3. Ô . Ô 1 D1 1

1 f

Y entonces los productos L/ . ´ L. vienen dados por:

4 1

L/ . ´ Î
3 73.: Ô D; 15; 8
1

1 f ; f

4 1

L. Î
3 73 D; 15; 8
1

1 f ; f

´´ Î
3. Ô . Ô 1 D1 1

1 f

Luego
4 1

K. Î
3 9. Ô . Ô 1 D1 3[.: Ô D1 15; D; 15; ]: 0
1 f ; f

Luego debe cumplirse

1

. Ô . Ô 1 D1 3[.: Ô D; 15; D; 15; ] 0

; f

Para 0, 1, 2, …, luego

[Link]
 15/

[. Ô . Ô 1 . Ô f f ]D1 3[.: Ô 15; 15; ] D; 0

; f

Y para 0 tenemos que Ô.Ô 1 Ô f f 0 ya que Df ≠ 0

El polinomio 5 de segundo grado, dado por

5.Ô Ô.Ô 1 Ô f f

tomar Ô, son las raíces de 5 , además tenemos:

Se llama polinomio indicial de la ecuación (1), y los únicos valores que puede

5.Ô D1 1 0 para 1, 2, 3, ….

Donde

1 ∑1; f[. Ô 15; 15; ]D; para 1, 2, 3, ….

Notemos que 1 es una combinación lineal de Df , D/ , … , D15/ cuyos coeficientes

involucran las funciones conocidas L/ , L , Ô.

Definimos [/ .Ô .Ô / / Df

[/ .Ô
8/ .Ô
5.Ô 1

Y en general
15/

[1 .Ô 3[.: Ô 15; 15; ] ; .Ô

; f

81 .Ô
g0 .Î
1, 2, 3, ….
C.Î>1
para

Los 81 así definidos, son funciones racionales de Ô y en los únicos puntos

donde no están definidos es en los puntos Ô para los cuales 5 .Ô 0 , para
algún 1, 2, … De éstos, solamente existen dos puntos posibles. Definimos:
4

#. , Ô Df Î Î
3 1 .Ô
1

1 /

Si esta serie es convergente para 0 Å entonces

[Link]
 K[#. , Ô ] Df 5.Ô Î

Hemos llegado a la situación de que si la función

4
Î
3 D1 1

1 f

Es solución de ´´ L/ . ´ L. 0

Entonces Ô debe ser una raíz del polinomio indicial 5 , y además todos los
D1 . ≥ 1 están determinados unívocamente en términos de Df y Ô para ser
las 1 .Ô siempre y cuando 5.Ô 1, 2, 3 … y
recíprocamente, si Ô es una raíz de 5 , y si las 1 .Ô se pueden determinar,
sea distinto de cero para

entonces la función . # . , Ô es solución de la ecuación

´´ L/ . ´ L. 0

Para toda selección de Df , siempre y cuando pueda demostrarse que # . , Ô es

convergente. Ejercicio que dejamos al lector, como también el hecho que las
funciones / y son linealmente independientes.

Luego como las raíces son Ô/ y Ô tales que 5 .Ô/ ≠ 0 5 .Ô ≠ 0 para

todo 1, 2, 3, … luego las soluciones son
4

/.
Î
13 1 .Ô/
1

1 f

. Î
23 1 .Ô
1

1 f
Siempre y cuando Ô/ Ô no sea un entero positivo y tampoco cero, ya que
5.Ô ≠ 0, para 1, 2, 3, ….

Si Ô/ Ô , o Ô/ Ô es un entero positivo, podemos caracterizar las soluciones

mediante el siguiente teorema.

Teorema 14: consideremos la ecuación:

´´ L/ . ´ L. 0

[Link]
Donde L/ y L tienen desarrollos en series de potencias, los cuales son
convergentes para Å. Sean Ô/ , Ô (con Åq .Ô/ ≤ Åq .Ô ) las raíces del polinomio
indicial

5.Ô Ô.Ô 1 L/ .0 Å L .0

Si Ô/ Ô existen dos soluciones /, linealmente independientes para | | ≤ Å,

las cuales tienen la forma de:

/. | |Ô/ [#/ . ]

. | |Ô/ 1[# . ] [SËÉ | |] #/ .

Donde #/ , # pueden desarrollarse en series convergentes de potencias para

| | ≤ Å y #/ .0 ≠ 0.

Si Ô/ Ô /,
independientes y definidas para | | ≤ Å que tienen la forma:
es un entero positivo, existen dos soluciones linealmente

/. | |Ô/ [#/ . ]

. | |Ô [# . ] <[SËÉ | |]#/ .

Donde #/ y # tienen desarrollados en series de potencias convergentes para

| | ≤ Å y además #/ .0 ≠ 0, # .0 ≠ 0, < es una constante; puede suceder que
< sea igual a cero.

Demostración: teníamos en la demostración anterior a la función

4

#. , Ô Df Î Î
3 D1 1

1 /

, K[#. , Ô ] Df 5.Ô Î

Además, 1 .Ô se determina por recurrencia mediante las fórmulas:

f .Ô Df ≠ 0

5.Ô 1 .Ô [1 .Ô
15/

[1 .Ô 3[.: Ô 15; 15; ] ; .Ô

; f

[Link]
para 1, 2, 3, ….

Como 5.Ô/ 0, 5´.Ô/ 0

De donde

à à#
¥K [# . , Ô ]¦ K³ ´. ,Ô
àÔ àÔ

Df [5´.Ô [log ]5. ] Î

Y si tenemos Ô Ô/ Ô , Df 1 , entonces

à#
. [. , Ô/ ]
àÔ

La cual nos daría una solución de la ecuación, siempre y cuando la serie

involucrada converja, entonces
4 4

. Î
1 3 ´ .Ô/ 1 [log ] Î
13 1 .Ô/
1

1 f 1 f

4
Î
1 3 ´ .Ô/ 1 [log ] /.
1 f

Donde / es la solución que ya se había obtenido

4

/.
Î
13 1 .Ô/
1

1 f

Nótese que ´1 .Ô/ existe para para todo 0, 1, 2, … ya que 1 es una función
racional de Ô, cuyo denominador no se anula en Ô Ô/. También f .Ô 1 implica
que ´f .Ô 0, y así que la serie que en multiplica a Î
1 comienza con la
primera potencia de .

Ahora supongamos que Ô/ Ô È , donde È es un entero positivo.

Si Df esta dado,

/ .Ô , .Ô , … , 5/ .Ô

[Link]
Existen todos y tienen valores finitos, pero dado que

5.Ô È .Ô [ .Ô ,

Encontramos dificultad para calcular .Ô ; como

5.Ô .Ô Ô/ .Ô Ô

Luego

5.Ô È .Ô Ô .Ô È Ô

Si [ .Ô también tiene a Ô Ô , como factor, esto implica que se puede cancelar

el mismo factor en 5 .Ô È , y 8 .Ô 1.Ô
2.Ô , …, existen todos. Más aún tenemos W de la forma
existe lo cual implica que

. Î
23 1 .Ô
1
[ f .Ô 1]
1 f

Siempre podremos arreglarla de tal manera que [ .Ô 0, para ello escogemos

8f .Ô Ô Ô

Como 5.Ô 81 .Ô V

Vemos que [1 .Ô es lineal homogénea en 8f .Ô , 8/ .Ô , … , 815/ .Ô

consiguiente [1 .Ô tiene a 81 .Ô Ô Ô , como factor.
y por

Así .Ô existirá en forma de número finito. Definimos

4

E. , Ô Î
3 1 .Ô
1
[ f .Ô Ô Ô]
1 f

Encontramos que

K [E. , Ô ] .Ô Ô 5.Ô Î

Haciendo Ô Ô , obtenemos entonces una solución É dada por

É. E. , Ô

Sin embargo 8f .Ô 8/ .Ô ⋯ 5/ .Ô 0

[Link]
Así, la serie que define a É realmente comienza con la meésima potencia de
entonces É tiene la forma:
y

É. Î
2 È[#. ] Î
1[#. ]

Donde # es una serie de potencias. De donde concluimos que É es un múltiplo

constante de la solución / que ya se había obtenido.

Lo cual nos lleva a obtener la solución asociada con Ô , derivamos a K [E. , Ô ] con
respecto a Ô, obteniendo:

à àE
¥K [E . , Ô ]¦ K³ . ,Ô ´
àÔ àÔ

5.Ô Î
.Ô Ô [5´.Ô .log 5 .Ô ] Î

Y haciendo Ô Ô , hallamos que dada por

àE
. . ,Ô
àÔ

La cual es solución, siempre y cuando las series involucradas sean convergentes,

tiene la forma:
4 4

. Î
2 3 ´1 .Ô 1
.log Î
23 1 .Ô
1

1 f 1 f

Donde f .Ô Ô Ô Dado que

f .Ô / .Ô ⋯ 5/ .Ô 0

O sea
4

. Î
2 3 ´1 .Ô 1
< .log /.
1 f

Donde < es una constante.

El método empleado en esta sección para obtener soluciones alrededor de algún

punto singular, se llama MÉTODO DE PROBENTUS.

Ejercicios interesantes por su aplicación, es hallar la solución es:

[Link]
1) .1 ´´ [ .D 1 ] ´ D 0

En la cual D, y son paramétricos fijos, se denomina la Ecuación

Hipergeométrica de Gauss.

2) ´´ .1 ´ D 0

Donde D es un real, se denomina la Ecuación de Laguerre; si D es un entero

no negativo, una solución a la ecuación es un polinomio.

3) ´´ 2 ´ 2R 0

Donde R es una constante, se llama la Ecuación de Hermite.

4) .1 ´´ ´ R 0

Donde R es una constante, se denomina la Ecuación de Chebyshev.

5) ´´ ´ . R 0

Donde R es un real, se denomina la Ecuación de Bessel y sus soluciones

reciben el nombre de funciones de Bessel.

12. ECUACIÓN DE BESSEL

Nuestro objetivo es solucionar la ecuación

´´ ´ . R 0

solución de la ecuación citada, alrededor de cero, bajo la condición de que R es

Utilizando la técnica presentada anteriormente, o sea hallar una base del espacio

un real, como vimos anteriormente se denomina la ecuación de Bessel de orden R.

Tenemos que el polinomio indicial asociado a la ecuación es

Ô R 0

Cuyas raíce son ±R , de donde el teorema 14 garantiza que la ecuación de Bessel

de orden R posee una solución, válida para todo , de la forma

[Link]
 4

/.
_
3 D1 1
, R 0
1 f

Para hallar los D1 , observamos que

4 4

. R /.
_
3 D15 1 _
3 R D1 1

1 1 f

/ ´.
_
3. R D1 1

1 f

´/ ´. _
3. R . R 1 D1 1

1 f

Al reemplazar en la ecuación original tenemos:

4 4

3[. R . R 1 . R R ] D1 1
3 D15 1
0
1 f 1

o
4

.2R 1 D/ 3[ .2R D1 D15 ] 1

0
1

Y tenemos por tanto

D15
D/ 0, D1 , ≥ 2
.2R

De esto se sigue de inmediato que

D/ D D6 ⋯ 0
Dw
D
2.2R 2
Dw
DŽ
2 ∗ 4 .2R 2 .2R 4

Y en general

[Link]
 Dw
D . 1 1
1
2 ∗ 4 ∗ 6 … .2 .2R 2 .2R 4 … .2R 2
Dw
. 1 1
2 1 ! .R 1 .R 2 … .R

De donde
4
. 1 1
/. Df 3 1>_
2 1 ! .R 1 .R 2 … .R
1 f

En donde Df ≠ 0 es una constante arbitraría.

Por varias razones, resulta conveniente hacer:

1
Df
2RΓ.R 1

En donde denota la tan conocida función gamma definida por

4
Γ.R + 9 5J M _5/ M, R > 0
f

Se puede demostrar que esta integral converge para todo R > 0, diverge hacia
∞ cuando R 0 y tiene los valores

Γ.1 1, Γ .R 1 RΓ.R , R > 0

En particular, Γ. 1 ! Siempre que sea un entero positivo, y por esta

razón a la función gamma se le conoce también como la función factorial
generalizada.

Volvemos ahora al hecho de que si hacemos Df

/
_H._>/
para obtener la primera

Esta solución se conoce como la función de Bessel de orden R de primera clase, y

de las dos soluciones particulares necesarias para resolver la ecuación de Bessel.

se presenta por I_ . Tenemos concretamente

4 1>_
I_ . 3. 1 1
2 1>_ ! .R 1 .R 2 … .R Γ.p 1
1 f

[Link]
Una expresión que puede volverse a escribir en la forma más sencilla como
4
. 1 1 1>_
I_ . 3 $ %
Γ.n 1 Γ.p n 1 2
1 f

En particular, cuando p 0, tenemos

4
. 1 1 1
If . 3 $ %
.n! 2
1 f

Y, en general, cuando p es un entero no negativo k

4
. 1 1 1>;
I; . 3 $ %
n! .n k ! 2
1 f

La gráfica de If , I/ y I aparecen esbozadas en la figura 30.

segunda solución linealmente independiente de I_ .

Para completar el estudio de la ecuación de Bessel, nos falta encontrar una

dependientes del valor de R.

Usamos de nuevo el teorema 14, dividiendo nuestra argumentación en caos

Caso 1: R > 0, 2R no es entero. En tal caso las raíces de la ecuación indicial no

argumento anterior con R en lugar de con R. Obviamente esto nos llevará a una
difieren en un entero, y se puede obtener una segunda solución repitiendo el

serie cuyos coeficientes tendrán la misma forma que anteriormente, y como la

función gamma está definida para valores negativos no enteros, de donde
podemos escribir

[Link]
 4
. 1 1 15_
I5_ . 3 $ %
Γ.n 1 Γ.n p 1 2
1 f

Para > 0. Finalmente observamos que la función está definida incluso cuando R
es de la forma : , : un entero, y de nuevo da lugar a una solución que es
/

linealmente independiente de I_ .

De todo ello se concluye que siempre que R no es un entero, la solución general

de la ecuación de Bessel de orden R es

. / I_ . I5I .

Caso 2: R 0. Aquí la ecuación de Bessel toma la forma

´´ ´ 0

Y su polinomio indicial tiene a cero como raíz repetida. Por tanto, por el teorema
14, podemos encontrar una segunda solución de la forma
4

f. 3 1
1
If . ln
1 /

Con If como anteriormente.

Para hallar los valores de 1, tenemos que

4

f. 3 15
15/
If . ln
1

4
If .
´f . 3 1
15/
I´f . ln
1 /

4
If .
´´. 3 . 1 1
15/
I´´f . ln 2I´f .
1 /

Al reemplazar en la ecuación tenemos:

4

/ 4 3[ 1 15 ] 15/ [ I´´f . I´f . xJf. ] ln 2I´f . 0

1

[Link]
y como

I´´f . I´f . xJf. 0

tenemos
4

/ 4 3[ 1 15 ] 15/
2I´f .
1

Como
4
2
I´f . 3. 1 1 15/
2 1. !
1 /

De donde
4 4
4
4 3[ ] 15/
3. 1 1>/ 15/
/ 1 15
2 1. !
1 1 /

Para facilitar la evaluación de los 1 multiplicamos ahora esta expresión por y

separamos la serie del primer miembro en sus partes par e impar para obtener
4 4

3 [.2 1 1 15/ ] 4 3[.2 ]

1>/ 1
/ 1 1 15
1 / 1

4
4
3. 1 1>/ 1
2 1. !
1

Luego / ⋯ 0

mientras que

4 1 y .2 . 1 > 1
1>/ Ž1
1 15 #0 .1! # ´

De donde

/
#

[Link]
 1 1 1 1
1 1
Ž
2 . 4 2 2Ž .2! 2
.
.

1 1 1
.
. 1 1>/
³1 ⋯ ´
1
2 1 . ! 2

Y de ello se sigue que

4
. 1 1>/ 1 1 1
f. 3 ³1 ⋯ ´$ % If . ln
. ! 2 2
1 /

En el trabajo teórico con funciones de Bessel es práctica común reemplazar f

por una cierta combinación de If y f . La función resultante se conoce como la
función de Bessel de orden cero de segunda clase, y se define por la fórmula
4
2 . 1 1
1 1 1 2
f. 3 ³1 ⋯ ´$ % I. Pln kQ
Ó . ! 2 2 Ó f 2
1 /

En donde k 0,57721566 … , y se conoce como la constante de Euler. La gráfica

de f la podemos observar en la figura 31

[Link]
Caso 3: R
2 0 , y el teorema 14, afirma que la segunda solución de la ecuación de Bessel
, un entero. Esta vez las raíces de la ecuación indicial difieren en

es de la forma
4

1. 3 ;
;>Î
<I1 . ln
; f

En donde < es una constante. Aquí también podemos calcular los ; y < por el
método de coeficientes indeterminados, similar al método empleado por el
teorema, razón por la cual lo omitiré, de donde tenemos

1 !
15/
1 . : ;51 Ñ1 1
1. 3 $ % $ %
2 ! 2 2 ! 2
; f

4
1 . 1 ; [Ñ; Ñ1>; ] ;>1
3 $ % I1 . ln
2 :! . : ! 2
; /

En donde

1 1
Ñ1 1 ⋯
2

Es importante notar que es habitual reemplazar 1 por una combinación lineal de

I1 y 1 , representada por 1 , que es denominada la función de Bessel de orden
de segunda clase. Se define por la fórmula

: 1 !
15/
1 . ;51 Ñ1 1
1. 3 $ % $ %
Ó :! 2 Ó. ! 2
; f

4
1 . 1 ; [Ñ; Ñ;>1 ] ;51 2
3 $ % I . Pln kQ
Ó :! . : ! 2 Ó 1 2
; /

[Link]