UNIVERSIDAD NACIONAL DE SAN CRISTÓBAL DE

HUAMANGA
FACULTAD DE INGENIERÍA DE MINAS, GEOLOGÍA Y CIVIL
DEPARTAMENTO ACADÉMICO DE MATEMÁTICA Y FÍSICA

MÓDULOS
Para la EP de Ingeniería de Sistemas
Semestre 2019-I

VLADIMIR ACORI FLORES

Abril de 2019
Práctica 1

Introducción a las ED

1.1. ED con Matlab

1. Averiguar si las siguientes funciones o expresiones son solución de la ecuación diferencial dada.
1
a) f (x) = 2 e−2x + ex es solución de y0 + 2y = ex para x ∈ R.
3
1 2 2
b) g(x) = + x es solución de y00 − 2 y + = 0 para x ∈ h0, +∞i.
x x x
c) xy − y = c, donde c constante arbitraria, es solución de (2x − 3y) y0 + y = 0.
2 3

Ejercicio:
√
a) y = 2x 1 − x2 es solución de y · y0 − 4x = −8x3 para x > 1.
b) y = e−x cos (x/2) es solución de y00 + 2y0 = e−x cos (x/2).
c) y = a e−x +b ex +c e−2x +d e2x donde a, b, c, d constantes es solución de

y(4) − 5y00 + 3y = 0.

d) y = (3 cos (2x) + sen (2x)) ex es solución particular de

y00 − 2y0 + 5y = 0.

e) y = x e5x cos (2x) es solución de

y(4) − 20y000 + 158y00 − 580y0 + 841y = 0.

cos (5 ln x) sen (5 ln x)
f) y = 20 · −3· es solución de
x x

x3 y000 + 2x2 y00 + 20xy0 − 78y = 0.

√
g) y = 1 + c 1 − x2 es solución de 1 − x2 y0 + xy = x.


Ecuaciones Diferenciales

3 − 2x3


dy y y
h) x3 + y3 = 3cxy para c constante es solución de = .
dx x (2y3 − x3 )
1 2x
2. La función y = + es solución particular de la edo
x 2 − x2

x2 y0 − x2 y2 + xy + 1 = 0

a) ¿en el intervalo I = h−4, −1i?

b) ¿en el intervalo I = h0, 2i?
c) ¿en el intervalo I = h−1, 0i?
d) ¿en el intervalo I = h3, 6i?
e) Dé el intervalo más grande en el cual está definida la solución.

3. Comprobar que la expresión dada sea solución de la ecuación diferencial dada. Determine al menos
una solución explícita en cada caso. Trace las gráficas de las soluciones explícitas y describa el
intervalo I de definición de cada solución.
 
dy 2y − 1
a) = (y − 1) (1 − 2y), ln = x.
dx y−1
b) 2xydx + x2 − y dy = 0, −2x2 y + y2 = 1.

2
x2 + c
4. La función y = + 1 tal que c ∈ R es solución de la edo
4

y0 = 2x y − 1
p

Grafique las curvas solución para c = −1, 0, 1, 3. Que puede decir de y = 1 ¿es solución de la edo?

5. Determine una región del plano xy para la cual la ecuación diferencial dada tenga una solución
única que pase por un punto (x0 , y0 ) en la región.
dy √
a) = xy
dx
dy
b) −y = x
dx
c) 1 + y3 y0 = x2

d) (y − x) y0 = y + x

6. Considere el problema de valor inicial

dy p
= 1 − y2
dx
para y(2) = 0 y y(π/2) = 1. ¿Tiene solución única?

UNSCH 2 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

7. Encuentre las trayectorias ortogonales para la familia de curvas dada. Trace la gráfica de algunas
curvas de la familia y de las trayectorias ortogonales.

a) y = c e−x
b) x2 + y2 = 2cx
c) y = −x − 1 + c ex
1
d) y =
x+c
e) cos y = c e−x Rta. sen y = k e−x

8. Determine un miembro de la familia de trayectorias ortogonales a la familia de curvas dada, que

pase por el punto indicado.

a) 3xy2 = 2 + 3cx, (1, 1).

b) y = c sen x, (0, 2).
1
c) y = , (2, 3).
ln (cx)
d) x + y = c ey , (0, 5). Rta. y = 3 e−x −x + 2
e)

1. Buscar la solución graficando el campo de direcciones de la ecuación diferencial:

dy
a) = e−x −2y
dx
dy
b) = x2 − y2
dx
c)

2. Supóngase que la población P de bacterias en un cultivo al tiempo t, cambia a una razón directa-
mente proporcional a P2 − P. Si inicialmente hay 1000 bacterias y después de 5 horas la población
se redujo a 100 bacterias, determine:

a) La población como función del tiempo.

b) La población después de un tiempo bien grande.

UNSCH 3 Lic. Vladimir Acori Flores

Práctica 2

ED de segundo orden y de orden superior

2.1. PVI y PVF

Cada familia de soluciones es la solución general de la ED en el intervalo indicado. Determine un miem-
bro de la familia que sea solución del pvi.

a) y = c1 ex +c2 e−x , x ∈ R; y00 + y0 = 0, y(0) = 0, y0 (0) = 1

b) y = c1 e4x +c2 e−x , x ∈ R; y00 − 3y0 − 4y = 0, y(0) = 1, y0 (0) = 2

c) y = c1 x + c2 x ln x, x > 0; x2 y00 − xy0 + y = 0, y(1) = 3, y0 (1) = −1

d) y = c1 + c2 cos x + c3 sen x, x ∈ R; y000 + y0 = 0, y(π) = 0, y0 (π) = 2, y00 (π) = −1

Cada familia biparamétrica es una solución de la ED dada, en el intervalo h−∞, +∞i. Determine si se
puede encontrar un miembro de la familia que fatisfaga las condiciones de frontera.

a) y = c1 ex +c2 e−x ; y00 + y0 = 0, y(0) = 1 ,y(1) = 1

b) y = c1 ex cos x + c2 ex sen x, y00 − 2y0 + 2y = 0,

i) y(0) = 1, y(π/2) = 1
ii) y(0) = 0, y(π) = 0

c) y = c1 x2 + c2 x4 + 3, x2 y00 − 5xy0 + 8y = 24,

i) y(−1) = 0, y(1) = 4
ii) y(0) = 3, y(1) = 0

2.2. Wronskiano
Compruebe si las siguientes funciones son li en el intervalo h−∞, +∞i, usando el Wronskiano de las
funciones.

a) y1 = x + 2, y2 = |x| + 2
Ecuaciones Diferenciales

b) y1 = ex/2 , y2 = x ex/2
c) y1 = sen x cos x, y2 = sen 2x
d) y1 = 0, y2 = x, y3 = ex
e) y1 = cos 2x, y2 = 1, y3 = cos2 x
f) y1 = ex , y2 = e−x , y3 = senh x
g) y1 = 1, y2 = senh2 x, y3 = cosh2 x
Compruebe que las funciones dadas forman un CFS de la ED en el intervalo indicado. Forme la solución
general.
a) y00 − 4y0 = 0, y1 = cosh 2x, y2 = senh 2x, h−∞, +∞i
b) x2 y00 − 6xy0 + 12y = 0, y1 = x3 , y2 = x4 , h0, +∞i
c) x2 y00 + xy0 + y = 0, y1 = cos(ln x), y2 = sen(ln x), h0, +∞i
d) x3 y000 + 6x2 y00 + 4xy0 − 4y = 0, y1 = x, y2 = x−2 , y3 = x−2 ln x h0, +∞i
e) y4 + y00 = 0, y1 = 1, y2 = x, y3 = cos x, y(4) = sen x h0, +∞i

2.3. EDL Homogéneas con coeficientes constantes

Determine la solución general de la ED.
a) y00 − 3y0 + 2y = 0
b) y00 − 10y0 + 25y = 0
c) 2y00 − 3y0 + 4y = 0
d) y000 + 3y00 + 4y = 0
e) y000 + 3y00 − 4y0 − 12y = 0
f) y(4) − 2y00 + y = 0
g) y(5) + 5y(4) − 2y000 − 10y00 + y0 + 5y = 0
h) 6.11y000 + 8.59y00 + 7.93y0 + 0778y = 0
i) 3.15y(4) − 5.34y00 + 6.33y0 − 2.03y = 0
Resuelva las siguientes ED sujetas a las condiciones iniciales.
a) y00 − 10y0 + 25y = 0, y(0) = 1, y0 (1) = 0
b) y00 − 10y0 + 25y = 0, y(0) = 1, y0 (1) = 0
c) 2y(4) + 3y000 − 16y00 + 15y0 − 4y = 0, y(0) = −2, y0 (0) = 6, y00 (0) = 3, y000 (0) = 1/2
d) y(4) − 3y000 + 3y00 − y0 = 0, y(0) = 0, y0 (0) = 0, y00 (0) = 1, y000 (0) = 1

UNSCH 5 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

2.4. EDL no Homogéneas con coeficientes constantes

Resuelva las ED por coeficientes indeterminados.

a) y00 + y0 − 6y = 2x

b) y00 + 3y = −48x2 e3x

c) y00 + 2y0 + y = sen 2 + 3 cos 2x

d) y000 − 3y00 + 3y0 − y = x − 4 ex

e) y(4) + 2y00 + y = (x − 1)2

f) y00 − 4y0 + 8y = 2x2 + 3 e2x cos 2x + 10x2 − x − 1 e2x sen 2x

g) y(4) + 2y00 + y = 2 cos x − 3x sen x

Resuelva el pvi

a) 2y00 + 3y0 − 2y = 14x2 − 4x − 11, y(0) = 0, y0 (0) = 0

b) y00 + 4y0 + 5y = 35 e−4x , y(0) = −3, y0 (0) = 1

c) y000 − 2y00 + y0 = 2 − 24 ex +40 e5x , y(0) = 0.5, y0 (0) = 2.5, y00 (0) = −4.5

d) y00 − 2y0 + 2y = 2x − 2, y(0) = 0, y(π) = π

Resuelva las ED por variación de parámetros.

a) y00 + y = sen x

b) y00 + y = sec x tan x

ex
c) y00 − 2y0 + y =
1 + x2
√
d) 2y00 + 2y0 + y = 4 x
√
e) 4y00 − 4y0 + y = ex/2 1 − x2

f) y000 + 4y0 = sec 2x

Resuelva las ED por variación de parámetros, sujeta a las condiciones iniciales.

a) 4y00 − y = x ex/2 ,y(0) = 1, y0 (0) = 0

b) y00 − 4y0 + 4y = 12x2 − 6x e2x , y(0) = 1, y0 (0) = 0

UNSCH 6 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

2.5. ED de Euler-Cauchy
Resuelva la ED respectiva.

a) 4x2 y00 + y = 0

b) x2 y00 + 5xy0 + 3y = 0

c) x3 y000 + xy0 − y = 0

d) x4 y(4) + 6x3 y(3) + 9x2 y00 + 3xy0 + y = 0

e) (x + 2)2 y00 + (x + 2)y0 + y = 0

Resuelva por el método de variación de parámetros.

a) 2x2 y00 + 5xy0 + y = x2 − x

b) x2 y00 − 2xy0 + 2y = x4 ex

c) x3 y000 − 3x2 y00 + 6xy0 − 6y = 3 + ln x3

Resuelva los pvi en el intervalo h0, +∞i

a) 4x2 y00 + y = 0, y(−1) = 2, y0 (1) = 4

b) xy00 + y0 = x, y(1) = 1, y0 (1) − 0.5 =

c) x2 y00 − 5xy0 + 8y = 8x6 , y(1/2) = 0, y0 (1/2) = 0

UNSCH 7 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

Ecuaciones Diferenciales de orden superior con Matlab.

Ejemplo 1. Resolver la ED y00 − 10y0 + 25y = 0.

Solución: Resolviendo con Matlab, se tiene ( aquí, Dy = y0 , D2y = y00 , D3y = y000 , D4y = y(4) , · · · )

>> y = dsolve(’D2y-10*Dy+25*y=0’,’x’)

y =

C3*exp(5*x) + C4*x*exp(5*x)

Note que 5 es una raíz doble del polinomio característico.

También podemos resolver la ED encontrando las raíces de su ecuación característica (EC): m2 − 10m +
25 = 0

>> R=roots([1 -10 25])

R =

5
5
5x
Vemos que m = 5 es una raíz doble y podemos encontrar dos soluciones de la forma y1 = e5x, y2 = x e
LI, luego el conjunto fundamental de soluciones de la EDL Homogénea es CF = {y1 , y2 } = e5x , x e5x ,
asi la solución general (SG) de la ED es combinación lineal del CF, es decir

y = c1 y1 + c2 y2 = c1 e5x +c2 x · e5x

Si a la ED le imponemos condiciones iniciales o de frontera, se tendría.

Ejemplo 2. Resolver el pvi

y00 − 8y0 + 15y = 0

y(0) = 2
y0 (0) = −1

Solución: Resolviendo mediante Matlab,

>> y=dsolve(’D2y+8*Dy+15*y=0’,’y(0)=2’,’Dy(0)=-1’,’x’)

y =

(9*exp(-3*x))/2 - (5*exp(-5*x))/2

>> pretty(y)

UNSCH 8 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

exp(-3 x) 9 exp(-5 x) 5
----------- - -----------
2 2

Notar que se ha obtenido las constantes 9/2, −5/2.

Luego la solución del PVI es

9 −3x 5 −5x
y= e − e .
2 2
Ejemplo 3. Resolver la ED con condiciones de frontera (CF)

y00 + y0 − 6y = 0
y(1) = 1
y(5) = 3

Solución: Las raíces de su EC son: −3, 2, luego la solución general de la ED es y = A e−3x +B e2x .
Calculamos las constantes A y B de las CF con Matlab.

syms x y A B ;
f=A*exp(-3*x)+B*exp(2*x)-y;
equ1=subs(f,{x,y},{1,1});
equ2=subs(f,{x,y},{5,3});
[A B]=solve(equ1,equ2);
pretty([A B])
vpa([A B],4) %devuelve la expresion con 4 cifras significativas

/ exp(3) (3 exp(12) - exp(20)) exp(-2) (3 exp(12) - 1) \

| - ----------------------------, ----------------------- |
\ exp(20) - 1 exp(20) - 1 /

ans =

[ 20.07, 0.0001362]

También podemos calcular A y B con matrices mediante Matlab,

x = 1, y = 1 ⇒ A e−3 +B e2 = 1
x = 5, y = 3 ⇒ A e−15 +B e10 = 3

>> a=[exp(-3),exp(2);exp(-15),exp(10)];
b=[1;3];
c=inv(a)*b

c =

UNSCH 9 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

2.006532312354813e+01
1.361995106213734e-04

Entonces la solución de la ED con condiciones de frontera es

y = 20.07 e−3x +0.0001362 e2x

Ejemplo 4. Calcular el Wronskiano de

e−3x y senh (x) − cosh (3x)

Solución: Para Matlab podemos proceder de la forma

clear, clc
syms x;
f1 = exp(-3*x);
f2 = sinh(x)-cosh(3*x);
df1 = diff(f1,x,1);
df2 = diff(f2,x,1);
m = [f1 f2;df1 df2];
W = det(m);
pretty(simplify(W))

exp(-2 x) 2 - exp(-4 x) - 3

Luego, el Wronskiano de y1 = e−3x , y2 = senh (x) − cosh (3x) es

W (y1 , y2 ) = 2 e−2x − e−4x −3.

Ejemplo 5. Calcular el Wronskiano de las funciones

ex/2 , x ex/2 y x2 ex/2

Solución: Para Matlab podemos proceder de la forma

clear, clc
syms x;
f1 = exp(x/2);
f2 = x*exp(x/2);
f3 = x^2*exp(x/2);
m = [f1 f2 f3;diff(f1,x,1) diff(f2,x,1) diff(f3,x,1);
diff(f1,x,2) diff(f2,x,2) diff(f3,x,2)];
W = det(m)
%pretty(simplify(W))

W = 2*exp((3*x)/2)

UNSCH 10 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

Luego, el Wronskiano de y1 = ex/2 , y2 = x ex/2 y y3 = x2 ex/2 es

W (y1 , y2 , y3 ) = 2 e3x/2

Ejemplo 6. Determinar si las funciones y1 = 1, y2 = 1 − x, y3 = 3x2 + 2x − 4 son LI en h−1, 2i.

Solución: Calculamos el Wronskiano de y1 , y2 , y3 con Matlab,

clear, clc
syms x;
f1 = 1;
f2 = 1-x;
f3 = 3*x^2+2*x-4;
m = [f1 f2 f3;diff(f1,x,1) diff(f2,x,1) diff(f3,x,1);
diff(f1,x,2) diff(f2,x,2) diff(f3,x,2)];
W = det(m)

W = -6

Como W (y1 , y2 , y3 ) = −6 6= 0 las funciones son LI.

Ejemplo 7. Dada la familia de funciones

y1 = e−2x , y2 = e3x , y3 = e5x .

Hallar una EDL Homogénea con el menor orden posible para la cual las funciones {yi } forma parte de
un CFS.

Solución: Como las soluciones son de la forma y = emx , entonces {y1 , y2 , y3 } debe ser un CFS de una
EDL Homogénea con coeficientes constantes de EC

(m + 2)(m − 3)(m − 5) = 0

desarrollando con Matlab se tiene

>> syms m;
>> expand((m+2)*(m-3)*(m-5))

ans = m^3 - 6*m^2 - m + 30

Luego, la ED que tiene por soluciones y1 , y2 , y3 es

y000 − 6y00 − y0 + 30y = 0.

Ejemplo8. Dada la familia de funciones

y1 = sen x, y2 = x.

UNSCH 11 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

Hallar una EDL Homogénea con el menor orden posible para la cual las funciones {yi } forma parte de
un CFS.

Solución 9. {y1 , y2 } es un CFS de la EDL Homogénea, luego resolviendo con Matlab

clear, clc
syms x y dy d2y;
M = [sin(x) x];
M1 = [y,M];
M2 = [dy,diff(M,x,1)];
M3 = [d2y,diff(M,x,2)];
W = det([M1;M2;M3]);
%pretty(simplify(W))

W = d2y*sin(x) + y*sin(x) - d2y*x*cos(x) - dy*x*sin(x)

La ecuación se satisface si W = 0, siendo dy = y0 , d2y = y00 , entonces ED que tiene por soluciones
y1 = sen x y y2 = x es
(sen x − x cos x) · y00 − x sen x · y0 + sen x · y = 0
Ejemplo 10. Obtener la solución particular de la ED

y00 + 3y0 − 4y = 34 sen x

Solución: Calculamos las raíces de la EC: m2 + 3m − 4 = 0 que son m1 = −4 y m2 = 1, luego la solución

de la EDL Homogénea es yh = c1 e−4x +c2 ex .
La propuesta para y p por coeficientes indeterminados es

y p = A sen x + B cos x

Calculamos A y B con Matlab

clear, clc
syms x A B;
yp = A*cos(x)+B*sin(x);
ecu = diff(yp,x,2)+3*diff(yp,x,1)-4*yp-34*sin(x);
ecu1 = subs(ecu,x,0);
ecu2 = subs(diff(ecu),x,0);
[A B] = solve(ecu1,ecu2)

A = -3

B = -5
Luego, y p = −3 cos x − 5 cos x. Entonces la solución general de la ED es:

y = yh + y p
y = c1 e−4x +c2 ex −3 cos x − 5 cos x.

UNSCH 12 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

Ejemplo 11. Obtener la solución particular de la ED

y00 − 6y0 + 9y = 2 e3x

Solución: Calculamos las raíces de la EC: m2 − 6m + 9 = 0 que es m = 3 de multiplicidad dos, luego la

solución de la EDL Homogénea es yh = c1 e3x +c2 x e3x .
La propuesta para y p por coeficientes indeterminados es

y p = Ax2 e3x

Calculamos A con Matlab

clear, clc
syms x A B;
yp = (A*x^2)*exp(3*x);
ecu = diff(yp,x,2)-6*diff(yp,x,1)+9*yp-2*exp(3*x);
ecu1 = subs(ecu,x,0);
%ecu2 = subs(diff(ecu),x,0);
[A ] = solve(ecu1)

A = 1

Luego, y p = x2 e3x . Entonces la solución general de la ED es

y = yh + y p
y = c1 e3x +c2 x e3x +x2 e3x .

Ejemplo 12. Obtener la solución particular de la ED

y00 + y = x2 − x − sen x.

Solución: Las raíces de la EC m2 + 1 = 0 son m = i y m = −i, luego yh = c1 cos x + c2 sen x. La solución

particular es
y p = Ax2 + Bx +C + Dx cos x + Ex sen x

clear, clc
syms x A B C D E;
yp = A*x^2+B*x+C +D*x*cos(x)+E*x*sin(x);
ecu = diff(yp,x,2)+yp-x^2+x +sin(x);
ecu1 = subs(ecu,x,0);
ecu2 = subs(diff(ecu,1),x,0);
ecu3 = subs(diff(ecu,2),x,0);
ecu4 = subs(diff(ecu,3),x,0);
ecu5 = subs(diff(ecu,4),x,0);
[A B C D E] = solve(ecu1,ecu2,ecu3,ecu4,ecu5)

A = 1

UNSCH 13 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

B = -1

C = -2

D = 1/2

E = 0
Luego, y p = x2 − x − 2 + 21 x cos x. Entonces la solución general de la ED es

y = yh + y p
1
y = c1 cos x + c2 sen x + x2 − x − 2 + x cos x.
2
Ejemplo 13. Obtener la solución particular de la ED

e3x
y00 − 6y0 + 9y = .
x2
Solución: Las raíces de la EC m2 − 6m + 9 = 0 son: m = 3 de multiplicidad dos, luego y1 = e3x , y2 = x e3x
soluciones LI de la EDL Homogénea. Entonces la solución de la EDL Homogénea es

yh = c1 e3x +c2 x e3x

La solución particular por variación de parámetros es

y p = v1 y1 + v2 y2

Calculando con Matlab y p se tiene

clear,clc
syms x ;
y1= exp(3*x);
y2= x*exp(3*x);
W= det([y1,y2;diff(y1,x,1),diff(y2,x,1)]); %Wronskiano de y_1, y_2
w= [0;exp(3*x)/x^2];
a1= [y1;diff(y1,x,1)];
a2= [y2;diff(y2,x,1)];
w1= [w,a2];
w2= [a1,w];
W1= det(w1);
W2= det(w2);
v1= int(W1/W,x)
v2= int(W2/W,x)
yp= simplify(v1*y1+v2*y2)

yp = -exp(3*x)*(log(x) + 1)

UNSCH 14 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

Luego, y p = − (1 + ln x) e3x . Por lo tanto, la solución general de la ED es:

y = yh + y p
y = c1 e3x +c2 x e3x − (1 + ln x) e3x .

Ejemplo 14. Resolver la ED

2x2 y00 + 5xy0 + y = x2 − x.
Solución: ED de Euler Cauchy no Homogénea de segundo orden.

La EC de la ecuación Homogénea asociada es 2m(m − 1) + 5m + 1 = 0 con raíces m = −1/2 y m = −1.

Luego las soluciones LI son y1 = x−1/2 y y2 = x−1 , entonces la solución de la ecuación Homogénea
asociada es
y = c1 x−1/2 + c2 x−1
La solución particular y p por VP con Matlab es

clear,clc
syms x ;
y1= x^(-1/2); %sol de la Homogenea
y2= x^(-1); %Sol de la Homogénea
W= det([y1,y2;diff(y1,x,1),diff(y2,x,1)]);
w= [0;(x^2-x)/(2*x^2)]; %cambiar para cada ejercicio
a1= [y1;diff(y1,x,1)];
a2= [y2;diff(y2,x,1)];
w1= [w,a2];
w2= [a1,w];
W1= det(w1);
W2= det(w2);
v1= int(W1/W,x)
v2= int(W2/W,x)
yp= simplify(v1*y1+v2*y2)

yp = (x*(2*x - 5))/30
x
Luego y p = (2x − 5). Entonces la solución general de la ED es
30
y = yh + y p
x
y = c1 x−1/2 + c2 x−1 + (2x − 5) .
30
Ejemplo 15. Encontrar la solución particular y p por el método de CI y VP.

y000 − 2y00 + y0 + 2y = e3x

Solución: Calculamos las raíces de la EC con Matlab

UNSCH 15 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

>> roots([1 -2 -1 2])

ans =

-1.0000
2.0000
1.0000
las raíces son m = −1, m = 2 y m = 1, luego la solución de la EDL Homogénea es yh = c1 e−x +c2 e2x +c3 ex .

La y p por coeficientes indeterminados es y p = A e3x , calculando A con Matlab

>> syms x A;
yp= A*exp(3*x);
ecu= diff(yp,x,3)-2*diff(yp,x,2)-diff(yp,x,1)+2*yp-exp(3*x);
ecu1= subs(ecu,x,0);
[A]= solve(ecu1)
A =1/8
1 3x
Luego, la solución particular es y p = e .
8
La solución particular y p por VP con Matlab es
clear,clc
syms x ;
y1= exp(-x);
y2= exp(2*x);
y3= exp(x);
W= det([y1,y2,y3;diff(y1,x,1),diff(y2,x,1),diff(y3,x,1);
diff(y1,x,2),diff(y2,x,2),diff(y3,x,2)] );
w= [0;0;exp(3*x)]; %cambiar para cada ejercicio
a1= [y1;diff(y1,x,1);diff(y1,x,2)];
a2= [y2;diff(y2,x,1);diff(y2,x,2)];
a3= [y3;diff(y3,x,1);diff(y3,x,2)];
w1= [w,a2,a3];
w2= [a1,w,a3];
w3= [a1,a2,w];
W1= det(w1);
W2= det(w2);
W3= det(w3);
v1= int(W1/W,x);
v2= int(W2/W,x);
v3= int(W3/W,x);
yp= simplify(v1*y1+v2*y2+v3*y3)

yp = exp(3*x)/8
Note la diferencia entre estos dos métodos.

UNSCH 16 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

Ejercicios:

Ejercicio 1. De las siguientes ED, indicar cuáles son problemas de valor inicial (PVI) y problemas de
valores en la frontera (PVF).

a) x2 y00 − xy0 + y = x, y(0) = 1, y00 (0) = 2 Rta. . . . c) (x − 2)y00 + 3y = x, y(0) = 0, y(1) = 2 Rta. . . .

b) y00 − 3y0 + 4y = 0, y(1) = 1, y0 (1) = 2 Rta. . . . d) x2 y00 − xy0 + 3y = 0, y(0) = 1, y0 (1) = 5 Rta. . . .

Ejercicio 2. La familia biparamétrica y = A ex cos x + B ex sen x es una solución de la ED y00 − 2y0 + 2y = 0

en el intervalo h−∞, +∞i. Determine un miembro de la familia que satisfaga las condiciones

a) y(0) = 1, y(π) = 0 Rta.

b) y(0) = 1, y(π) = −1 Rta.

c) y(0) = 1, y(π/2) = 1 Rta.

d) y(0) = 0, y(π) = 0 Rta.

e) y(0) = 0, y0 (0) = 0 Rta.

Ejercicio 3. Indique si el conjunto de funciones es LI o LD en el intervalo h−∞, +∞i.

a) y1 = x + 2, y2 = |x| + 2 Rta.

b) y1 = ex+2 , y2 = ex−3 Rta.

c) y1 = 0, y2 = x, y3 = ex Rta.

d) y1 = 1 + x, y2 = x, y3 = x2 Rta.

Ejercicio 4. Determine la solución general. Comprobando previamente que las funciones dadas forman
un conjunto fundamental de soluciones, CFS, de la ED en el intervalo dado.

a) y00 + y0 − 12y = 0, y1 = e−3x , y2 = 2 e−3x , h−∞, +∞i Rta.

b) y00 + y0 − 12y = 0, y1 = e4x , y2 = 2 e4x , h−∞, +∞i Rta.

c) y00 + y0 − 12y = 0, y1 = e−3x , y2 = e4x , h−∞, +∞i Rta.

d) x2 y00 − 6xy0 + 12y = 0, y1 = x3 , y2 = x4 , h0, +∞i Rta.

e) x2 y00 − 6xy0 + 12y = 2, y1 = x3 , y2 = x4 , h0, +∞i Rta.

Ejercicio 5. Si y p1 = 3 e2x y y p2 = x2 + 3x son soluciones particulares de las ED respectivamente

y00 − 6y0 + 5y = −9 e2x y y00 − 6y0 + 5y = 5x2 + 3x − 16

entonces, determine soluciones particulares para

a) y00 − 6y0 + 5y = 5x2 + 3x − 16 − 9 e2x Rta.

UNSCH 17 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

b) y00 − 6y0 + 5y = −10x2 − 6x + 32 + e2x Rta.

Ejercicio 6. La ED en su forma estándar y00 + p(x)y0 +q(x)y = 0 para p(x) y q(x) continuas en el intervalo
I. Si y1 es una solución conocida de la ED, la segunda solución y2 LI con y1 se obtiene aplicando la
fórmula de reducción Z − R p(x)dx
e
y2 = y1 · dx
y21
Determine la solución general de la ED dada, suponiendo conocida la solución y1 .

a) y00 − 4y0 + 4y = 0, y1 = e2x Rta.

b) x2 y00 − 7xy0 + 16y = 0, y1 = x4 Rta.

c) 1 − x2 y00 + 2xy0 = 0, y1 = 1


Rta.

d) 1 − 2x − x2 y00 + 2(1 + x)y0 − 2y = 0, y1 = x + 1 Rta.

Ejercicio 7. Dada las ED, proponer una solución particular y p por el método de coeficientes indetermi-
nados (MCI).

a) y(4) + y000 = x + x2 e−x Rta.

b) y00 + −6y0 + 9y = x + 1 − 5 e3x Rta.

c) y00 + y = x sen x Rta.

d) y00 + 4y = 2x2 − 3 sen 2x + e2x Rta.

Ejercicio 8. Encontrar la EDL Homogénea de coeficientes constantes dada las raíces de su ecuación
característica.

a) m = 1 mult. 2, m = −2 mult. 2, m = 3 mul. 1 Rta.

b) m = 0 mult. 3, m = i mult. 1 Rta.

Ejercicio 9. Encontrar la ED de Euler- Cauchy dada las raíces de su ecuación característica.

a) m = −1 mult. 1, m = 2 mult. 1 Rta.

b) m = 0 mult. 1, m = 2 + i mult. 1 Rta.

Ejercicio 10. Determinar la solución particular y p por los métodos de CI y VP y obtener la solución
general de las ED.

a) y000 − 2y00 − y0 + 2y = e3x Rta.

b) y00 − y = sec3 x − sec x Rta.

c) x2 y00 − xy0 + y = 4x ln x Rta.

d) y000 − 5y00 + 6y0 = 2 sen x + 8 Rta.

UNSCH 18 Lic. Vladimir Acori Flores

Ecuaciones Diferenciales

Ejercicio 11. Un partícula se mueve a la largo del eje x según la ecuación

x00 + 9x0 + 20x = 0

A partir de un punto a 3 metros a la derecha del origen, la partícula en el tiempo t = 0 segundos se dispara
hacia la izquierda con una velocidad de v = 21m/seg, se pide:

a) x en función del tiempo t.

b) El tiempo en que la partícula pasa por el origen.

c) El desplazamiento máximo negativo.

d) La velocidad máxima negativa.

Ejercicio 12. Una partícula se mueve a lo largo del eje x de acuerdo con la ley

x00 + 4x0 + 13x = 0

Si esa partícula empieza su movimiento en x = 0, con una velocidad inicial de v = 6m/seg, hacia la
izquierda, se pide:

a) x en función de t.

b) Los tiempos en que se producen las paradas.

Rta. a) x = −2 e−2t sen 3t, b) t = nπ/3, n = 0, 1, 2, . . .

UNSCH 19 Lic. Vladimir Acori Flores