UNIVERSIDAD DE CONCEPCION
FACULTAD DE CIENCIAS
FISICAS Y MATEMATICAS
DEPARTAMENTO DE INGENIERIA MATEMATICA
Ecuaciones Diferenciales ordinarias (521218)
Juan Molina Salazar
Capı́tulo 2
EDO de Orden Superior
Definición 1 (Orden Superior)
Una EDO de orden dos es lineal, si ella puede escribirse como
a2 (t)y 00 (t) + a1 (t)y 0 (t) + a0 (t)y(t) = f (t) (1)
Diremos que la EDO anterior está normalizada si se escribe como
y 00 (t) + b(t)y 0 (t) + c(t)y(t) = h(t) (2)
a1 (t)
sobre {t ∈ R : a(t) > 0} o sobre {t ∈ R : a(t) < 0}; obviamente se tiene que b(t) = ,
a2 (t)
a0 (t) f (t)
c(t) = y h(t) = .
a2 (t a2 (t
Ejemplo 1 (EDO normalizada de orden superior)
La EDO
(t − 1) x00 (t) − et x0 (t) + t x(t) = t5
definida para para t > 1 o bien para t < 1, es equivalente a la EDO normalizada
et 0 t t5
x00 (t) − x (t) + =
t−1 t−1 t−1
definida sobre {t ∈ R : t > 1}, o bien sobre {t ∈ R : t < 1}.
Definición 2 (EDO homogénea)
Si en la ecuación (1) el termino conocido f es identicamento nulo, f ≡ 0, decimos que la EDO
es HOMOGENEA, es decir,
a2 (t)y 00 (t) + a1 (t)y 0 (t) + a0 (t)y(t) = 0.
Ejemplo 2 (EDO homogénea y normal)
y 00 (t) + (t − 2)y 0 (t) + 5y(t) = 0.
1
Definición 3 (Operadores Diferenciales Lineales)
Para J1 ⊂ R, J2 ⊂ R consideremos los espacio vectorial
C 1 (J1 , R) = {y : J → R : existe y 0 (t) e y 0 (t) es continua}.
C 0 (J2 , R) = {y : J → R : y(t) es continua}
Entonces definamos la aplicación lineal
dy
D : C 1 (J, R) → C 0 (J2 , R) definida porDy = (3)
dt
En las aplicaciones, usualmente se escribe Dy(t) o y 0 (t) en lugar de Dy; se debe estar atento,
muchas veces se quiere indicar la función y 0 y se escribe y 0 (t).
De manera similar, de manera inductiva para n ∈ N se escribe
C n (J1 , R) = {y ∈ C n (J1 , R) : existe y (n) (t) e y (n) (t) es continua}.
y ponemos
d(n) y
Dn : C n (J, R) → C 0 (J2 , R) definida porDn y = (4)
dtn
que también resulta un operador lineal.
Observación 1
(i) Observe que la EDO
y 00 (t) − y 0 (t) − 6y(t) = cos(t) (5)
se puede escribir como (D2 − D − 6)(y) = cos(t).
(ii) La EDO homogénea asociada a la EDO anterior es
y 00 (t) − y 0 (t) − 6y(t) = 0 (6)
o de modo equivalente
(D2 − D − 6)(y) = 0 (7)
(iii) Observe que toda solución de la EDO homogénea anterior es un elemento del Ker del
operador D2 − D − 6 .
Observación 2 (General)
En general note que toda solución de la EDO lineal homogénea de segundo orden, normal y
de coeficientes variables
y 00 (t) + p(t)y 0 (t) + q(t)y(t) = 0 (8)
es un elelemento que pertenece al Kernel del operador
T = D2 + p(t)D + q(t)
2
Ejemplo 3
Consideremos la EDO a coeficientes constantes
y 00 (t) − y 0 (t) − 6y(t) = 0 (9)
Entonces ya sabemos que determinar todas las soluciones de (9) es equivalente a determinar el
Kernel del operador T = D2 − D − 6.
Se puede comprobar fácilmente que las funciones y1 e y2 definidas por
y1 (t) = e3t , y2 (t) = e−2t
pertenecen al Ker(T ). Mas aún, se puede ver que toda combinación lineal c1 y1 + c2 y2 ( donde
c1 y c2 son constantes reales) también pertenece al Ker(T ) (esto último es obvio, puesto que el
Ker(T ) es un espacio vectorial).
C 1 (J1 , R). De otra parte, observe que {y1 (t) = e3t , y2 (t) = e−2t } es un conjunto linealmente
independiente en el espacio vectorial C 2 (J1 , R). En efecto, sean a, b ∈ R tales que
a y1 + b y2 = θ
Esto es, para todo t > 0, a e3t + b e−2t = 0.
Derivando la relación anterior, se obtiene 3a e3t − 2b e−2t = 0,
y evaluando en t = 0 , sigue el sistema
a+b = 0
3a − 2b = 0.
Resolviendo, se obtiene a = b = 0; esto es, el conjunto {y1 (t) = e3t , y2 (t) = e−2t } es linealmente
independiente en C 2 (J1 , R).
A este punto cabe preguntarse, ¿Existen otras soluciones de la EDO considerada? La respuesta
es NO, y se obtiene del siguiente Teorema:
Teorema 1 (de la Dimensión)
Consideremos la EDO lineal homogénea
y 00 (t) + p(t)y 0 (t) + q(t)y(t) = 0 (10)
donde p(t) y q(t) son dos funciones continuas dadas, entonces el conjunto
S = {u ∈ C 2 (J1 , R) : u00 (t) + p(t)u0 (t) + q(t)u(t) = 0}
es un subespacio vectorial de C 2 (J1 , R), además
Dim(S) = 2.
Observación 3
Observe que el conjunto S del Teorema anterior corresponde al Kernel del operador
T = D2 + p(t)D + q(t); y sabemos que el Kernel de toda transformación lineal es un sube-
spacio vectorial del espacio vectorial del dominio del operador. Por tanto, lo que el Teorema
afirma es que la dimensión del Kernel del operador lineal en cuestión, es dos (el orden la EDO
considerada).
3
Ejemplo 4 Determine todas las soluciones de la EDO
y 00 (t) − y 0 (t) − 6y(t) = 0.
Del Ejemplo en la EDO (9) y del Teorema de la Dimensión sigue que si z(t) es solución de la
EDO dada, entonces existen constantes reales c1 y c2 de modo que
z(t) = c1 e3t + c2 e−2t
Pregunta ¿Se puede determinar los valores de las constantes c1 y c2 ?
Respuesta: Sı́. La respuesta afirmativa seguirá del Teorema de Existenca y Unicidad, que
veremos la próxima clase.
Lunes 27 Marzo 2023
Teorema 2 (Caso general)
Consideremos la EDO lineal homogénea
y (n) (t) + an−1 (t) y (n−1) (t) + an−2 (t) y (n−2) (t) + · · · + a1 (t)y 0 (t) + a0 (t)y(t) = 0 (11)
donde para cada j ∈ {0, 1, · · · , n − 1} las funciones aj (t) son continuas y conocidas. Si escribi-
mos
S = {u ∈ C n (J1 , R) : u es solución de (11)},
entonces resulta que Dim(S) = n.
Ejemplo 5 Consideremos la EDO
y 000 (t) − 7y 0 (t) + 6y(t) = 0 (12)
Resulta que las funciones y1 (t) = et , y2 (t) = e2t e y3 (t) = e−3t son soluciones de la EDO
dada (Ejercicio: compruebe la afirmación anterior). Además, las funciones {y1 , y2 , y3 } son l.i.
en C 3 (J1 , R). Por el Teorema de la dimensión sigue que toda solución u(t) de (12), es de la
forma
u(t) = c1 et + c2 e2t + c3 e−3t
donde c1 , c2 y c3 son constantes arbitrarias.
Teorema 3 (de Existencia y Unicidad de soluciones)
Sean a1 (t), a0 (t) y f (t) funciones continuas y conocidas definidas en un intervalo abierto J,
t0 ∈ J y b0 , b1 dos números reales conocidos. Entonces el Problema con Valor Inicial
00
y (t) + a1 (t)y 0 (t) + a0 (t)y(t) = f (t)
y(t0 ) = b0
y 0 (t0 ) = b1
admite única solución.
4
Ejemplo 6 Considere el PVI
00
y (x) + y(x) = x
y(0) = π
y 0 (0) = −π
Mas adelante en el curso veremos que toda solución u(x) de la EDO y 00 (x) + y(x) = x es de la
forma
u(x) = c1 cos(x) + c2 sen (x) + x.
Ahora, derivando la expresión anterior sigue u0 (x) = −c1 sen (x) + c2 cos(x). Evaluando en
x = 0, se obtiene el sistema lineal de dos por dos:
u(0) = c1 + 0 = π
u0 (0) = c2 + 1 = −π,
Por tanto, la única solución al PVI dado, es
u(x) = π cos(x) + (−π − 1)sen (x) + x.
Definición 4 (Wronskiano)
El Wronskiano de dos funciones derivables f (t) y g(t), se define como
f (t) g(t)
W [f, g](t) := 0 = f (t)g 0 (t) − f 0 (t)g(t).
f (t) g 0 (t)
Observación 4 Note que para cada t ∈ Dom(f ) ∩ Dom(g), el Wronskiano de f y g en t
define una función derivable. Si no hay ambiguedad en las funciones f y g en que se evalua,
simplemente escribiremos W (t).
Ejemplo 7
Consideremos f (t) = e2t , g(t) = t. Entonces
2t
e t
W [e , t](t) := 2t = e2t (1 − 2t).
2t
2e 1
Observe que en este caso
d dh i
W [e2t , t] (t) = W (t)
dt dt
d 2t
= e (1 − 2t)
dt
= −2 e2t (1 + 2t).
Definición 5 (W de orden tres)
De manera inductiva, se define el Wronskiano de las funciones derivables y1 , y2 e y3 , como
y1 (t) y2 (t) y3 (t)
W [y1 , y2 , y3 ](t) := y10 (t) y20 (t) y30 (t)
y100 (t) y200 (t) y300 (t)
5
Teorema 4 (Independecia Lineal)
Sean y1 , y2 dos funciones de C 2 (J, R). Entones si para t0 ∈ J
W [y1 , y2 ](t0 ) 6= 0,
entonces {y1 , y2 } es l.i. en C 2 (J, R).
Demostración
De la combinación lineal nula, c1 y1 + c2 y2 = θ, obtenemos el sistema (cuyas incógnitas son c1
y c2 )
c1 y1 (t) + c2 y2 (t) = 0
c2 y10 (t) + c2 y20 (t) = 0
De la hipótesis sigue que el sistema anterior tiene solución única, a saber, c1 = c2 = 0; lo
anterior dice que {y1 , y2 } es l.i. en C 2 (J, R).
Ejemplo 8 El wronskiano
2t
e t
W [e , t](t) := 2t = e2t (1 − 2t).
2t
2e 1
Entonces podemos asegurar que el conjunto {e2t , t} es l.i. en C 2 (J, R), donde J es cualquier
intervalo contenido en R con la condición que el número (1/2) 6∈ J.
Mas aun, tenemos la siguiente e importante
Observación 5 Si dos funciones f y g son l.i en C 2 (J, R) para algún intervalo J de R, no se
puede concluir que W [f, g](t) 6= 0.
En efecto, basta considerar las funciones f y g como sigue
( (
0 si t < 0, 0 si t < 0,
f (t) = 2 , y g(t) =
t si t ≥ 0. t2 si t ≥ 0.
Es fácil ver que {f, g} es l.i. en C(R, R); sin embargo se puede ver que W [f, g](t) = 0 (para lo
anterior haga el análisis primero para t > 0).
Teorema 5 (Teorema de Abel)
Sean y1 e y2 dos soluciones de la EDO lineal homogénea
y 00 (x) + p(x) y 0 (x) + q(x) y(x) = 0
para x ∈]a, b[. Entonces
h Z i
W (y1 , y2 )(x) = c exp − (p(x) dx
donde c es una constante real.
Corolario 6 El wronskiano de dos soluciones y1 e y2 de la la EDO lineal homogénea
y 00 (x) + p(x) y 0 (x) + q(x) y(x) = 0
es cero, o nunca es cero.
6
Aplicación
(i) Las funciones y1 (x) = ex e y2 (x) = cos(x) no pueden ser solución de
y 00 (x) + p(x) y 0 (x) + q(x) y(x) = 0
en el intervalo ] − π, π[ pues W (ex , cos(x))(x) = −ex sen (x) + cos(x) se anula en (π/4).
(ii) ...
Teorema 7 (Raı́ces repetidas)
Sea L un operador diferencial lineal de orden n.
Si la función y ∈ Ker(L) (esto es, si y(t) es solución de la EDO homogénea L(y) = 0), entonces
la función z(t) = t y(t) ∈ Ker(L2 ) (esto es L2 (z) = 0).
α1 = −4, De manera inductiva, t z(t) = t2 y(t) ∈ Ker(L3 ), etc.
Ejemplo 9 Sea L el operador L = (D − a), note que y(t) = eat ∈ Ker(L), entonces
z(t) = t eat
es solución de L2 (y) = 0.
De concreto, si a = 2 entonces del hecho que e2t es solución de y 0 (t) + 2y(t) = 0, sigue que
z(t) = t e2t es solución de
y 00 (t) − 4y 0 (t) + 4y(t) = 0 (13)
puesto que L = (D − 2)2 conduce a y 00 (t) − 4y 0 (t) + 4y(t) = 0
Note que La EDO homogénea asociada al operador (D−2)3 es y (3) (t)−6y 00 (t)+12y 0 (t)−8y(t) = 0
y una base para el espacio solución, es B = {e2t , te2t , t2 e2t }. Ası́, toda solución u(t) de
y (3) (t) − 6y 00 (t) + 12y 0 (t) − 8y(t) = 0, se escribe como
u(t) = c1 e2t + c2 t e2t + c3 t2 e2t ,
con c1 , c2 y c3 constantes reales arbitrarias.
Objetivo
Nuestro objetivo es resolver EDO lineal no homogénea
y (n) (t) + an−1 (t) y (n−1) (t) + an−2 (t) y (n−2) (t) + · · · + a1 (t)y 0 (t) + a0 (t)y(t) = f (t) (14)
donde para cada j ∈ {0, 1, · · · , n} las funciones aj y f son funciones continuas y conocidas.
Observe que el operador diferencial lineal anterior se puede escribir como
L = Dn + an−1 (t) D(n−1) + · · · + a1 (t) D + a0 (t), (15)
entonces la EDO (14) se puede escribir como
L(y) = f.
El siguiente Teorema nos da la forma que tiene toda solución de (14).
7
Teorema 8 (Forma de la solución general)
Si y = y(t) es una solución de la EDO (14), entonces
y(t) = yp (t) + yh (t) (16)
donde yp (t) es una solución cualquiera de L(y) = f e yh (t) es una solución arbitraria del Ker(L)
Observación 6 Que yh (t) sea una solución arbitraria del Ker(L), significa que si el Ker(L)
tiene base (o sistema fundamental) dado por {y1 , · · · , yn } entonces
n
X
yh (t) = cj yj (t)
i=1
donde los cj son constantes reales arbitrarias.
Ejemplo 10 Consideremos la EDO
x2 y 00 (x) − x y 0 (x) + y (x) = cos ( ln (x)) , x > 0. (17)
En este caso una solución particular yp para (17) es yp (t) = − 21 sen ( ln (x)). Además, la EDO
homogénea asociada, tiene sistema fundamental dado por B = {y1 (x) = x, y2 (x) = x ln (x)},
esto es, yh se escribe como yh (x) = c1 x + c2 x ln (x).
En consecuencia, la solución general de (17), es
1
y (x) = − sen ( ln (x)) + A x + B x ln (x)
2
siendo A, B ∈ R constantes arbitrarias.
Primero resolveremos el caso homogéneo a coeficientes constantes, esto es, la EDO ho-
mogénea
y (n) (t) + dn−1 y (n−1) (t) + dn−2 y (n−2) (t) + · · · + d1 y 0 (t) + d0 (t)y(t) = 0 (18)
donde los coeficientes dj son todas constantes reales.
d
Si ponemos D = entonces la EDO (18) en forma de operadores se escribe
dt
Dn + dn−1 D + · · · + d1 D + d0 D0 = 0
En este caso, si buscamos soluciones del tipo y(t) = eαt , al reemplazar en (18) se obtiene
eαt αn + dn−1 αn−1 + · · · + d1 α + d0 = 0;
ası́, u(t) = eαt es solución de la EDO homogénea (18) si, y sólo si α es raı́z del polinomio
p(α) = αn + dn−1 αn−1 + · · · + d1 α + d0 = 0 (19)
Como los dj son todas constantes reales, del Teorema Fundamental del Algebra sigue que el
polinomio (19) tiene exactamente n raı́ces en el campo complejo C.
8
Para resolver lo anterior y visto que
p(α) = q(α)(α2 + c1 α + c0 )
donde q es un polinomio de coeficientes reales constantes de grado n − 2, podemos reducir el
estudio a las raı́ces de
r(α) = α2 + c1 α + c0 (20)
Teorema 9 (Tres casos)
Para resolver la EDO homogénea a coeficientes reales y constantes
(D2 + c1 D + c0 )y = 0, (21)
relativo a las raı́ces del polinomio r(α) = α2 + c1 α + c0 , existen tres casos distintos:
(i) Caso α1 y α2 reales y distintos. En este caso las dos soluciones l. i. buscadas (SISTEMA
FUNDAMENTAL son: y1 (t) = eα1 t e y2 (t) = eα2 t .
(ii) Caso α1 = α2 . En este caso del Teorema de las Raı́ces Repetidas, las soluciones l.i.
buscadas son y1 (t) = eα1 t e y2 (t) = t eα1 t
(iii) Caso α1 y α2 complejos conjugados, a saber, si α1 = a + ib, entonces las soluciones l.i.
buscadas son: y1 (t) = eat cos(bt), y2 (t) = eat sen (bt).
Ejemplo 11
(i) Caso (i), y 00 (t) + y 0 (t) − 12y(t) = 0. Aquı́, el polinomio es
p(α) = α2 + 4α + 12 con raı́ces α1 = −4, α2 = 3. Solución general y(t) = c1 e3t + c2 e−4t .
(ii) Caso (ii), y 00 (t) − 4y 0 (t) + 4y(t) = 0. Aquı́, el polinomio es
p(α) = α2 − 4α + 4 con raı́ces α1 = α2 = 2. Solución general y(t) = c1 e2t + c2 te2t .
(iii) Caso (iii), y 00 (t) − 4y 0 (t) + 13y(t) = 0. Aquı́, el polinomio es
p(α) = α2 − 4α + 13 con raı́ces α1 = 2 + 3i, α2 = 2 − 3i.
Solución general y(t) = c1 e2t cos(3t) + c2 e2t sen (3t).
Ahora, en el caso de L a coeficientes constantes, exponemos un método para determinar
una solución particular de L(y) = f .
Método de los Aniquiladores
Definición 6 Decimos que el operador diferencial Lineal a coeficientes constantes L ANIQUILA
a la función f , si L(f ) = 0.
Ejemplo 12 D aniquila a toda función constante f (x) = c.
D2 aniquila a la función f (x) = x.
(D − a) aniquila a eat
(D − a)2 aniquila a t eat
D2 − 2a D + (a2 + b2 ) aniquila a eat cos(bt) y también aniquila a eat sen (bt).
9
Observación 7 (i) Si L1 y L2 son operadores a coeficientes constantes de modo que L1 (f ) =
0 y L2 (g) = 0 entonces L1 ·L2 (c1 f +c2 g) = 0 donde c1 y c2 son constantes reales arbitrarias.
(ii) Existen muchas funciones que no son aniquilables, por ejemplo, las funciones f (x) =
(1/x), f (x) = (sen )−1 (x), f (x) = ln (x), etc.
El Método
Sea L un operador diferencial lineal a coeficientes constantes entonces el método de los
Aniquiladores para resolver la EDO lineal L(y) = h, en el caso que h sea una función aniquil-
able, se desarrolla como sigue:
(i) Aniquilar la función h a través de un operador diferencial lineal L1 , esto es, L1 (h) = 0.
(ii) Se determina una propuesta de solución particular, yp (x), para la EDO considerada
L(y) = h. Esta propuesta tiene la forma de la solución homogénea de de la EDO
L1 L(y) = 0.
(ii) LLevar esta propuesta a la EDO L(y) = h para reducir la solución propuesta a una
mı́nima expresión.
Ejemplo 13 Determinar la solución general de de
y 00 (t) − y 0 (t) − 12y(t) = e4t (22)
Solución
Sabemos que la solución general z(x) de la EDO lineal dada L(y) = e4t (aquı́ el operador L
es L = D2 − D − 12) es del tipo
z(x) = yp (x) + yh (x)
donde yp (x) es una solución cualquiera de L(y) = e4t e yh (x) es una solución arbitraria del
Ker(L).
(i) Si f (x) = e4x entonces para determinar la yp (x) de este caso, primero aniquilamos f por
D − 4.
De otra parte, observando que el operador diferencial lineal asociado a la EDO es L = (D −
4)(D + 3), y del hecho que (D − 4)e4x = 0, se sigue que
(D − 4)(D − 4)(D + 3) = 0. Entonces el método de aniquiladores sugiere buscar soluciones
particulares del tipo:
yp (x) = c1 e4x + c2 xe4x + c3 e−3x (23)
donde c1 , c2 y c3 son constantes reales arbitrarias.
Sin embargo de la condición
h i
(D − 4)(D + 3)yp (x) = (D − 4)(D + 3) c1 e4x + c2 xe4x + c3 e−3x = e4x
vemos que nuestra propuesta yp se reduce a yp (x) = c xe4x ; y entonces el problema se reduce a
determinar el valor de la constante c. Para ello, puesto que yp (x) = c xe4x debe ser solución de
que
(D − 4)(D + 3)[c xe2x ] = e4x .
10
Por tanto, debe tenerse h i
e2x 4c − 2c − 6c = e2x ;
haciendo el cálculo correspondiente sigue que c = (1/7).
Ası́, la solución particular buscada es yp (x) = (1/7) xe4x .
De otra parte, como L = (D − 4)(D + 3), vemos que todo elemento yh del Ker(L) es del tipo
yh (x) = c1 e4x + c2 e−3x .
Finalmente, la solución general a la EDO propuesta es:
y(x) = (1/7) xe4x + c1 e4x + c2 e−3x
donde c1 y c2 son constantes reales arbitrarias.
Ejercicio 1 Determine la única solución del PVI
00
y (t) − y 0 (t) − 12y(t) = e4x
y(0) = 1
y 0 (0) = −1
Solución:
· · · La única solución al PVI dado resulta ser
h i
y(t) = (13/49) + (1/7) e4t + (36/49) e−3t
Variación de Parámetros
Ahora queremos resolver
y 00 (t) + a1 (t)y 0 (t) + a0 (t)y(t) = h(t) (24)
en que los coeficientes a1 y a0 son, en general, funciones continuas no constantes. Note que en
este caso el operador asociado a la EDO dada es
L = D2 + a1 (t) D + a0 (t))
y obviamente el método de los aniquiladires no funciona.
También queremos resolver el caso en que la EDO teniendo coeficientes constante, el termino
forzante sea no aniquilable, por ejemplo
1
y 00 (t) − 3y 0 (t) + 2y(t) =
t
En ambos casos, supondremos que ya hemos determinado un sistema fundamental para la EDO
homogénea asociada L(y) = 0.
Ası́, supongamos que {y1 , y2 } es una base para Ker(L).
El método de Variación de Parámetros consiste en buscar una solución particular, yp (t),
para L(y) = h(t) como
yp (t) = c1 (t)y1 (t) + c1 (t)y2 (t) (25)
donde ahora los coeficientes c1 (t) y c2 (t) son funciones incógnitas a determinar.
11
Para determinara las incógnitas c1 (t) y c2 (t), sabemos que se debe cumplir que
h i
L c1 (t) y1 (t) + c2 (t) y2 (t) = h(t), (26)
esto es,
yp00 (t) + a1 (t)yp0 (t) + a0 (t)yp (t) = h(t) (27)
Ahora reemplazando nuestra propuesta yp en (25 ) en la EDO (27 ), esto es,
d2 h i dh i h i
c 1 (t)y 1 (t) + c 1 (t)y 2 (t) +a 1 (t) c 1 (t)y 1 (t) + c 1 (t)y 2 (t) +a 0 (t) c 1 (t)y1 (t) + c 1 (t)y 2 (t) = h(t)
dt2 dt
(28)
Antes de reemplazar nuestra propuesta en la EDO (27 ), notemos que
2 h
X i
yp0 (t) = c0i (t)yi (t) + ci (t)yi0 (t) (29)
i=1
y en la igualdad anterior imponemos que
2 h
X i
c0i (t)yi (t) = c01 (t)y1 (t) + c02 (t)y2 (t) = 0; (30)
i=1
por tanto, la derivada de la solución propuesta, yp (t) se reduce a
2 h
X i
yp0 (t) = ci (t)yi0 (t) , (31)
i=1
por tanto,
2 h
X i
yp00 (t) = c0i (t)yi0 (t) + ci (t)yi00 (t) (32)
i=1
y reemplazado este valor en (27 ) se obtiene que
2 h
X i
(c0i yi0 + ci yi00 ) + (a1 (t)ci (t)yi0 (t) + a0 (t)ci (t)yi (t)) .
i=1
Ahora recordando que L(y1 ) = 0 y L(y2 ) = 0, en la igualdad anterior se obtiene
c01 (t)y10 (t) + c02 (t)y20 (t) = h(t).
Ası́, para determinar las incógnitas c1 (t) y c2 (t) debemos resolver el sistema de dos por dos
c01 (t)y1 (t) + c02 (t)y2 (t) = 0
c01 (t)y10 (t) + c02 (t)y20 (t) = h(t)
y0 y
Ejemplo 14 Para x > 0, resolver y 00 − + 2 = x, sabiendo que las funciones y1 e y2 ,
x x
definidas por y1 (x) = x e y2 (x) = x ln (x) son soluciones del problema homogéneo asociado.
12
Desarrollo:
Puesto que {y1 (x) = x, y2 (x) = x ln (x)} es l.i. en el espacio de las funciones dos veces
derivables (Ejercicio), sigue que la solución general del problema propuesto es
y(x) = A x + B x ln (x) + yp (x) (A y B constantes arbitrarias)
con yp (x) = c1 (x) x + c2 (x) x ln (x) solución particular del problema.
Aquı́ las funciones c1 y c2 (que son las incógnitas) verifican:
(
c01 (x)x + c02 (x)x ln (x) = 0
c01 (x) · 1 + c02 (x)( ln (x) + 1) = x,
resolviendo el sistema anterior se obtiene
x2
c1 (x) = 1 − 2 ln (x)
4
c2 (x) = 1 x2 .
2
Es decir, la yp viene dada por
x2 1 x3
yp (x) = 1 − 2 ln (x) x + x2 x ln (x) = .
4 2 4
Finalmente, la solución general del problema, es:
x3
y(x) = A x + B x ln (x) + .
4
donde A y B son constantes arbitrarias,
Abril, 25 de 2023.
JMS//jms
13