Tema 7 Ecuaciones Diferenciales Ordinarias
Tema 7 Ecuaciones Diferenciales Ordinarias
Tema 7 Ecuaciones Diferenciales Ordinarias
Índice
1 Conceptos generales 1
1 Conceptos generales
Una ecuación diferencial es una ecuación en la que la incógnita es una función y que
contiene la incógnita y, la variable x de la que depende la incógnita y las derivadas y (n) de
la incógnita. Es decir, es una expresión de la forma
Se llama orden de la ecuación diferencial al orden de la derivada más alta que aparece
en la ecuación. Ası́, por ejemplo, la ecuación y 00 + y 0 = y es una ecuación de segundo orden,
mientras que la ecuación y 0 + y 2 = 0 es de primer orden.
Una solución de una ecuación diferencial es una función y(x) que al sustituirla en la
ecuación diferencial da una identidad. La solución de una ecuación diferencial no es única.
De hecho el lector ya conoce las ecuaciones diferenciales más sencillas, las que son de la forma
y 0 = f (x). La solución de una tal ecuación no es más que la primitiva de f (x), que ya conoce
el lector está definida salvo una constante. Es decir, sonR soluciones de la ecuación diferencial
anterior todas las funciones y(x) de la forma y(x) = f (x)dx + c, siendo c una constante
arbitraria. Se llama solución general de una ecuación diferencial a una solución de la
ecuación diferencial que depende de constantes y tal que cualquier solución de la ecuación
diferencial se obtiene de esta solución general dándole valores apropiados a las constantes.
Todas las ecuaciones diferenciales que veremos en este curso tienen una solución general que
depende de tantas constantes como sea el orden de la ecuación diferencial, es decir: una
ecuación diferencial de orden n tendrá una solución general que dependerá de n constantes.
Una manera de fijar una solución particular (es decir, de fijar los valores de las constantes
arbitrarias) de una ecuación diferencial es imponer condiciones (a menudo vienen impuestas
estas condiciones por el problema fı́sico que ha originado la ecuación). Las condiciones que
1
C?? Ecuaciones diferenciales 2
estudiaremos son de dos tipos: las condiciones iniciales, que fijan el valor de la función y
sus derivadas hasta el orden n − 1 si la ecuación diferencial es de orden n en un punto dado
x0 , y las condiciones de contorno, que fijan los valores de la función en los extremos de
un intervalo.
Ası́, por ejemplo, la ecuación diferencial y 00 = 0 tiene como solución general y = ax + b.
Si fijamos las condiciones iniciales y(0) = 1, y 0 (0) = 2, entonces la única solución verificando
estas condiciones es y(x) = 2x + 1. Si fijamos las condiciones de contorno y(0) = 1, y(4) = 5,
entonces la única solución verificando estas condiciones es y = x + 1.
Una solución elemental de una ecuación diferencial es una solución y(x) que es una
combinación de las funciones elementales que hemos estudiado en este curso. No siempre (se
podrı́a decir que casi nunca) existe una solución elemental de una ecuación diferencial, pero
hay ocasiones en que es posible, al menos, dar una solución implı́cita usando funciones
elementales. Incluso hay ocasiones en que la solución en modo implı́cito es más cómoda que
la solución explı́cita. Ası́, por ejemplo, la ecuación diferencial 3y 2 y 0 + xy 0 + y = 0 tiene como
solución implı́cita y 3 + xy = C, que es más manejable que la solución explı́cita
1 √ 1
−2 3 x (27C + 729C 2 + 108x3 ) 3
y= √ 1 + 1 .
(27C + 729C 2 + 108x3 ) 3 3 (2) 3
arctg y = x − ln |x + 1| + C, y y = tg (x − ln |x + 1| + C) .
que es una ecuación de variables separables, que se resuelve como indicamos en el caso
anterior.
Veamos un ejemplo. La ecuación (x2 +y 2 ) dx+2 x y dy = 0 es del tipo A dx+B dy, donde
A y B son polinomios homogéneos de grado 2, luego se trata de una ecuación diferencial de
primer orden homogénea, luego conviene hacer el cambio v = y/x, es decir, y = v x, con lo
que
dy = v dx + x dv, x2 + y 2 = x2 + v 2 x2 , 2 x y = 2 x2 v,
y, sustituyendo en la ecuación,
(x2 + v 2 x2 ) dx + 2 x2 v 2 dx + 2 x3 v dv = 0, (1 + 3 v 2 ) x2 dx + 2 x3 v dv = 0,
(ρ y)0 = ρ y 0 + ρ P y, (∗)
ρ y 0 + ρ P y = ρ y 0 + ρ0 y, de donde ρ0 y = ρ P y, i.e. ρ0 = ρ P,
de donde Z R
P dx
ln ρ = P dx + c, ρ=C e ,
y, sustituyendo en (**),
Z
1 1
e y = ex ex dx + C = e2x + C, de donde y = ex + C e−x .
x
2 2
y 0 = f (x, y),
, (1)
y(a) = y0
y queremos encontrar una solución aproximada en el intervalo [a, b]. Para ello comenzamos
dividiendo el intervalo [a, b] en N subintervalos iguales de longitud (paso) h = (b − a)/N , de
modo que tenemos
Para valores de h pequeños, la función y(x), en el intervalo [a, x1 ] se puede aproximar bien
por la función cuya gráfica es la recta tangente en a, y(x) = y(a) + y 0 (a) (x − a) = y(a) +
f (a, y0 ) (x − a), donde hemos sustituido y 0 (a) por la expresión que tiene en nuestro problema
de valores iniciales. Por lo tanto, aproximadamente,
A partir de este valor aproximado de y(x1 ), podemos usar el mismo argumento para calcular
el valor aproximado de y(x2 ):
Esta es la fórmula del método de Euler o método de la tangente, que aproxima la solución
de (1) por (2).
C?? Ecuaciones diferenciales 7
Método del desarrollo en serie de Taylor. Uno de los métodos más efectivos para
resolver ecuaciones diferenciales es suponer que la(s) solucion(es) de la ecuación se pueden
expresar como una serie de potencias, es decir, que sus desarrollos de Taylor son una buena
aproximación de la función. Eso lo usan con frecuencia los fı́sicos en sus desarrollos teóricos
(por ejemplo, en la deducción de la ecuación de ondas).
Veamos como resolver una ecuación de la forma y 0 (x) = f (x, y(x)), con la condición
inicial y(0) = 0, usando el método del desarrollo de Taylor .
En primer lugar consideramosP una aproximación y(n, x) hasta el orden n de la función
incógnita y(x), es decir y(n, x) = ni=0 ai xi . De la condición y(0) = 0 se deduce a0 = 0 .
A continuación , en la expresión f (x, y(x)), cambiamos y(x) por su aproximación de orden
n, y(n, x), y hacemos el desarrollo de Taylor g(n, x) de f (x, y(n, x)) alrededor de x = 0, hasta
el orden n − 1 .
Después igualamos la derivada yp(n, x) de y(n, x) (que será un polinomio de grado n − 1)
con g(n, x) (que será otro polinomio de grado n−1) , lo que exige la igualdad de los coeficientes
de los polinomios yp(n, x) y g(n, x), lo que da un sistema de ecuaciones cuya solución nos
dará los coeficientes a1 , ... , an de la solución aproximada y(n, x) .
La solución y(n, x) se aproximará rápidamente a la real si los coeficientes a1 , ... , an son
pequeños (menores que 1) , mientras que la aproximación será lenta (habrá que tomar un
valor grande de n para tener una buena aproximación) si a1 , ... , an son grandes . Además,
en este caso, aún con n grandes y(n, x) se aproximará a la solución real solo para pequeños
valores de x .
Con Mathematica, todo los anterior se puede expresar de la siguiente manera
Xnn
y[n , x ] := a[i]xxi /.{xx → x, nn → n};
i=0
a[0] = 0;
nn−1
X
yp[n , x ] := (i + 1)a[i + 1]xxi /.{xx → x, nn → n};
i=0
f s[n , x ] := f [xx, y[nn, xx]]/.{xx → x, nn → n};
g[n , x , i ] := D[f s[nn, xx], {xx, ii}]/.{xx → x, nn → n, ii → i};
nn
X 1
(∗g[n , x ] := g[n, 0, i]xi ∗);
i!
i=0
1
coef = Solve[T able[(i + 1)a[i + 1] == (g[n, x, i]/.x → 0), {i, 0, n − 1}], T able[a[i +
i!
1], {i, 0, n − 1}]]
Las soluciones de esta última ecuación dan los coeficientes que permiten conocer la
solución aproximada y[n, x] .
Veamos dos ejemplos :
Ej.1) Resolver aproximadamente, por el método de Taylor, la ecuación diferencial y 0 (x) =
xy − sen x, y(0) = 0.
Aplicamos el método anterior usando la función f (x, y) = xy − sen x. Con Mathematica,
el trabajo se puede hacer como sigue:
In[]= Clear[”@”];
f [x, z] := xxzz − Sin[xx]/.{xx → x, zz → z}
C?? Ecuaciones diferenciales 10
nn
X
y[n , x ] := a[i]xxi /.{xx → x, nn → n};
i=0
a[0] = 0;
nn−1
X
yp[n , x ] := (i + 1)a[i + 1]xxi /.{xx → x, nn → n};
i=0
f s[n , x ] := f [xx, y[nn, xx]]/.{xx → x, nn → n};
g[n , x , i ] := D[f s[nn, xx], {xx, ii}]/.{xx → x, nn → n, ii → i};
nn
X 1
(∗g[n , x ] := g[n, 0, i]xi ∗);
i!
i=0
1
coef [n ] := Solve[T able[(i + 1)a[i + 1] == (g[n, x, i]/.x → 0), {i, 0, n − 1}], T able[a[i +
i!
1], {i, 0, n − 1}]]
coef ic = T able[coef [n], {n, 5, 15, 5}]
1 1
Out[]= {{{a[1] → 0, a[2] → − , a[3] → 0, a[4] → − , a[5] → 0}}, {{a[1] → 0, a[2] →
2 12
1 1 11 19
− , a[3] → 0, a[4] → − , a[5] → 0, a[6] → − , a[7] → 0, a[8] → − , a[9] →
2 12 720 10080
137 1 1
0, a[10] → − }}, {{a[1] → 0, a[2] → − , a[3] → 0, a[4] → − , a[5] → 0, a[6] →
725760 2 12
11 19 137 3767
− , a[7] → 0, a[8] → − , a[9] → 0, a[10] → − , a[11] → 0, a[12] → − , a[13] →
720 10080 725760 239500800
97943
0, a[14] → − , a[15] → 0}}}
87178291200
In[]= asol[n , i ] := coef ic[[n/5]][[1]][[i]][[2]];
Xn
ysol[n , x ] := asol[n, i] xi ;
i=1
Las soluciones aproximadas de órdenes 5, 10 y 15 respectivamente, son
In[]= T able[ysol[n, x], n, 5, 15, 5]
x2 x4 x2 x4 11 x6 19 x8 137 x1 0 x2 x4 11 x6 19 x8
Out[]= {− − , − − − − − , − − − − −
2 12 2 12 720 10080 725760 2 12 720 10080
1
137 x 0 3767 x 21 7943 x 41
− −9 }
725760 239500800 87178291200
Por otra parte, la solución exacta es
In[]= solex = DSolve[{w0 [x] == f [x, w[x]], w[0] == 0}, w, x]
Out[]= r
1 − 1 + x2 π i +x 1 1+i x
{{w → F unction[{x}, − E 2 2 i Erf √ + 2 Erf i √ − Erf i √ ]}}
2 2 2 2 2
Las siguientes gráficas comparan la solución exacta y las aproximadas de órdenes 5, 10
y 15
In[] = Show[Graphics[{T ext[”solex”, {2.11, solex[[1]][[1]][[2]][2.31]}],
T ext[”n = 5”, {2.31, ysol[5, 2.31]}],
T ext[”n = 10”, {2.51, ysol[10, 2.31]}], T ext[”n = 15”, {2.51, ysol[15, 2.31]}]}],
P lot[{solex[[1]][[1]][[2]][x], ysol[5, x], ysol[10, x], ysol[15, x]},
{x, −2, 2.3}, P lotStyle → {{Hue[0]}, {Hue[0.2]}, {Hue[0.4]}, {Hue[0.8]}},
P lotRange → All, DisplayF unction → Identity], P lotRange → All, Axes → T rue]
C?? Ecuaciones diferenciales 11
3.0.1 Ejercicios
1. Resuelve las siguientes ecuaciones diferenciales:
d2 y
a) dy = 4xdx, b) = 6x, c) y 2 + (y 0 )2 = 1.
dx2
y0
2. Resuelve la ecuación diferencial y = 2.
2
3. Encuentra la solución general de la ecuación diferencial xydy − (1+y 1+x2
)dx
= 0, y la
Rsolución
1
particular que pasa por el punto (1, 3). (Ayuda: para resolver la integral
1 a bx+c
x(1+x2 )
dx, usar que x(1+x 2 ) = x + 1+x2 , donde a, b y c son constantes a determinar)
y dy
a) + 1 = 0, (1, 2); b) x2 dy + y 2 dx = 0, (1, 1);
x dx
c) s(1 + r2 )ds + r(1 + s2 )dr = 0, (s = −2, r = 2).
p dr
a) xydy − 1 − y 2 dx = 0, (1, 0); b) t = r2 , (1, 1);
dt
dy
c) − xy 2 + x = 0, (1, 1); d) sin(x) cos(x)dx + cos2 (x)dy = 0, (1, 2).
dx