Unidad 2 - ECUACIONES DIFERENCIALES ORDINARIAS DE 1ER ORDEN

Descargar como doc, pdf o txt
Descargar como doc, pdf o txt
Está en la página 1de 53

1

UNIVERSIDAD PRIVADA DE SANTA CRUZ DE LA SIERRA

C A L C U L O III

UNIDAD 2

ECUACIONES DIFERENCIALES ORDINARIAS DE PRIMER ORDEN

1.1.- CONCEPTOS BASICOS Y ALGUNAS IDEAS

1.2.- CONSIDERACIONES GEOMETRICAS, ISOCLINICAS. METODOS DE


EULER

1.3.- ECUACIONES DIFERENCIALES DE VARIABLES SEPARADAS Y


VARIABLES SEPARABLES

1.4.- ECUACIONES DIFERENCIALES EXACTAS

1.5.- FACTORES INTEGRANTES

1.6.- ECUACIONES HOMOGENEAS

1.7.- ECUACIONES DIFERENCIALES LINEALES DE PRIMER ORDEN

1.8.- METODOS DE VARIACION DE LOS PARAMETROS

1.9.- ECUACIONES DE LAGRANGE

1.10.- ECUACIONES DE CLAIRAUT

1.11.- ECUACION DE BERNOULLI

1.12.- APLICACIONES, MECANICA, ETC.


2

1.1.- CONCEPTOS BASICOS Y ALGUNAS IDEAS

Familias de Soluciones

Si tenemos la ecuación diferencial de 1er orden de la forma:

dy
 f  x
dx

su anti derivada es : y   f  x  dx  C

La solución es y  F  x   C cuya derivada es idéntica a f  x  en algún intervalo.

dy
Observamos que se demuestra que cualquier otra función satisface  f  x ,
dx
pero difiere de la primera función por una constante para toda x en el intervalo.

Este concepto nos da a entender que las soluciones no ocurren aisladas, sino
como una familia de soluciones con un parámetro llamada constante arbitraria C .
pero si queremos obtener un miembro de estas soluciones es necesario dar
valores a las variables y, x y obtener el valor del parámetro C .

Ahora si la ecuación diferencial de 1er orden es de la forma:

dy
 f  x, y 
dx

El problema es encontrar soluciones como funciones C 2  2 y    x  que


satisfagan la ecuación cuando se sustituya por la variable independiente y , por lo
general es más difícil, si no imposible pero si existen aparecen como una familia
de soluciones con un parámetro.
3

Ejemplo 1

Dada la ecuación diferencial

dy
 8sen 4 x
dx

Hallar la familia de soluciones y un miembro de dicha familia.

 cos 4 x 
y   8sen4 xdx  8    C  2 cos 4 x  C donde
 4 

y  2 cos 4 x  C familia de soluciones

Para encontrar un miembro de dicha familia hacemos y  6 cuando x  0


Estamos obligados para elegir c  8 entonces:

y  2 cos 4 x  8

Ejemplo 2

Dada la función y  f  x   ce su derivada es y  f   x   2ce


2x 2x

dy
Decimos que la ecuación diferencial  2y
dx

Tiene una familia de soluciones y  f  x   ce


2x

y si buscamos una solución que satisfaga

dy
 2 y ; si y  4 cuando x  0 vemos 4  ce0  c  4
dx

la solución es : y  f  x   4e
2x
4

Ejemplo 3

La ecuación diferencial de segundo orden:

y  f   x   12 x 2

12 x 3
Integrando y   f   x  dx   12 x 2 dx   C1  4 x 3  C1 ; y  f   x   4 x  C1
3

4 x4
Integrando nuevamente: y   f   x  dx    4 x 3  C1  dx   C1 x  C2  x 4  C1 x  C2
4

y  x 4  C1 x  C2

Observamos que existen dos constantes arbitrarias de modo que tenemos una
familia de soluciones con dos parámetros.

Para encontrar los parámetros que son las constantes C1 , C2 para un mismo valor
de x  0 hacemos y  0   1 ; y  0   2 sustituyendo en :

y  f   x   4 x 3  C1

y  x 4  C1 x  C2 tenemos

y  0   0  C1  2  C1  2 Reemplazando en y  x 4  C1 x  C2 tenemos:

y  x 4  2 x  C2 Luego hacemos x  0 con y  0   1 obtenemos:

y  0   0  0  C2  1  C2  1 Reemplazando en y  x 4  2 x  C2 tenemos:

y  x4  2x  1
5

1.2.- CONSIDERACIONES GEOMETRICAS, ISOCLINICAS.

Consideraciones Geométricas

Vemos que una ecuación diferencial de primer orden tiene una familia de
soluciones y manera de entender la naturaleza de estas soluciones es trazar las
graficas representativas de esta familia

Supongamos que   x, y,c   0 representa una familia de curvas.

Una trayectoria ortogonal a esta familia es una curva que corta cada miembro
de la familia en ángulo recto.

Ejemplo Conceptual

Una familia de curvas x 2  y 2  a 2 admite como trayectoria ortogonal toda recta


que pase por el origen, observamos en la figura:

Para hallar trayectorias ortogonales, podemos utilizar los métodos de ecuaciones


diferenciales.

Sabemos que las pendientes de dos curvas se cortan en un ángulo recto, en el


punto de intersección, son inversas y signo contrario.

Por lo tanto para hallar las trayectorias ortogonales de la familia   x, y,c   0 ,


primero determinamos la ecuación diferencial que satisface a la familia
dy dx
(eliminando c en el proceso) y luego reemplazando por  en esa
dx dy
ecuación, la solución de la nueva ecuación diferencial da las trayectorias
ortogonales.
6

Ejemplo 1

Hallar las trayectorias ortogonales de la familia de parábolas

y 2  cx
dy
Sacamos su diferencial 2 ydy  cdx es decir 2y c o sea reemplazando
dx
dy dy
en la familia de parábolas y 2  2 xy es decir y  2x como solución de
dx dx
diferencial de las parábolas.

dx dy 2x
Para las trayectorias ortogonales escribimos y  2 x o sea  es
dy dx y
decir 2 xdx  ydy  0 sacamos su solución 2  xdx   ydy  0 es decir

x2 y 2
2   0c tenemos 2 x 2  y 2  2c dividiendo la ecuación en ambos
2 2
x2 y 2
miembros por 2c tenemos   1 observamos que las trayectorias
c 2c
ortogonales son una familia de elipses confocales de acuerdo a la figura:
7

Ejemplo 2

Trace la gráfica de varios miembros de la familia de soluciones para la ecuación:

dy
 8sen 4 x
dx

De acuerdo al ejemplo anterior vimos:

 cos 4 x 
y   8sen4 xdx  8    C  2 cos 4 x  C
 4 

y  2 cos 4 x  C es la familia de soluciones

Al representar la gráfica de las soluciones de C  2,1, 0, 1

Obtenemos la figura:

Vemos por cada punto del plano pasa exclusivamente un miembro de la


familia de soluciones
8

Ejemplo 3

Trace la gráfica de varios miembros de la familia de soluciones de la ecuación


diferencial de segundo orden:

d2y
 12 x 2
dx 2

Como vimos en el ejemplo anterior la familia de soluciones es:

12 x 3
Integrando y   f   x  dx   12 x 2 dx   C1  4 x 3  C1 ; y  f   x   4 x  C1
3

4 x4
Integrando nuevamente: y   f   x  dx    4 x 3  C1  dx   C1 x  C2  x 4  C1 x  C2
4

y  x 4  C1 x  C2

Trazamos la gráfica de las soluciones a las parejas de constantes:

C1  2, C2  1 ;

C1  0, C2  0 ; C1  0, C2  1 tenemos la figura:
9

Se puede apreciar claramente, hay al menos dos soluciones que pasan por los
puntos  0,1 y por uno cercano  1/ 2, 0  , pero no tienen la misma pendiente en su
punto de intersección.

Si especificamos una solución particular para la ecuación diferencial de 2do


orden decimos que en el punto donde pasa una solución tiene una pendiente
en dicho punto.

Vemos entonces por cada punto en el plano pasa exclusivamente un


miembro de la familia de soluciones que tiene una pendiente dada.

Las isoclinas de una ecuación

En los ejemplos anteriores conocimos algunas propiedades geométricas de las


familias de soluciones por métodos analíticos, ahora veremos métodos
geométricos para encontrar curvas solución.

Consideremos la ecuación de 1er orden:

dy
y   f  x, y  donde f  x, y   c es una familia isoclinas
dx

Si asignamos a cada punto  a, b  en la función obtenemos alguna dirección con


pendiente f  a, b   c donde c es un parámetro, lo que obtenemos es un campo
de direcciones de la ecuación diferencial llamada familia de isoclinas, es en tal
sentido para cualquier solución debe tener una gráfica y en cada punto de la
misma tendrá la dirección que se requiere.

Una manera de visualizar es dibujar una pequeña marca en varios puntos de


la solución para indicar la dirección asociada en esos puntos que es la
pendiente en cada punto. Esto se puede hacer sistemáticamente dibujando
primero las curvas llamadas isoclinas que son curvas que indican la dirección
indicada por la ecuación diferencial que es fija y que las curvas solución
tienen en cada punto la pendiente que indican las marcas.

La isoclina nula f  x, y   0 proporciona las líneas en las que pueden estar


situados los máximos y mínimos de las curvas solución
10

Para trazar las curvas solución con mayor exactitud se hallan también el lugar
geométrico de los puntos de inflexión con las segundas derivadas.

Para poder hallar los mismos es necesario hallar la segunda derivada de la

dy
ecuación y   f  x, y  tenemos:
dx

d 2 y f f f f f f
y    y   f  x, y  igualando a cero  f  x, y   0 que es
dx 2
x y x y x y

precisamente el lugar geométrico de los puntos de inflexión.

Para averiguar las direcciones de las concavidades de las curvas solución se


utilizan las segundas derivadas, despejando de las mismas el valor de y .

Dicha función divide al plano xoy en dos partes y observamos la que está situada
por encima de la función con y>0 nos indican que las curvas solución tienen sus
concavides para arriba y la que está situada por debajo de la función con y<0
nos indican que las curvas solución tienen sus concavides para abajo.

Ejemplo 1

Trace gráficamente las curvas solución, considerando las isoclinas de la siguiente


ecuación diferencial:

dy
y
dx

La familia de isoclinas es f  x, y   y  c

Si le damos valores a la función x  0 es decir f  0, y   y  c tenemos isoclinas


que son líneas rectas en cada punto y la dirección impuesta por la ecuación
diferencial dada es el número c
11

Obtenemos la figura:

Por ejemplo se ve en la figura en cada punto de la recta y  1 la ecuación


dy
diferencial  1 determina una dirección con pendiente 1 es decir dicha
dx
pendiente determina la marca solución que tiene la misma dirección a lo largo de
la isoclina y  1 .

El figura se ha dibujado varias curvas solución indicando las direcciones asociadas


con cada isoclina mediante pequeñas marcas.

Para obtener gráficamente las curvas solución se realiza de la siguiente forma:

Empezamos en una marca dirección de cualquier punto en el plano y nos


movemos a lo largo de las marcas dirección obtenemos una curva solución
cuya dirección sea siempre el de las marcas de dirección.

Analíticamente estudiaremos más adelante varios métodos para encontrar la


familia de soluciones, por el momento nos conformaremos en ilustrar.
12

Ejemplo 2

Use el método de isoclinas para bosquejar algunas curvas solución para la


siguiente ecuación diferencial:

dy
 x2  y 2
dx

Observamos en la figura que las isoclinas son círculos x 2  y 2  c con c  

Cuando c  1 la isoclina tiene radio 1 ; cuando c  4 el radio es 2

En los círculos que representan las isoclinas se ha marcado cada una con la
dirección apropiada bosquejando varias curvas que representan soluciones para
la ecuación dada.
13

Ejemplo 3

Trazar aproximadamente las curvas solución de la ecuación diferencial

y  2 x  y

Solución

Para obtener las ecuaciones de las isoclinas utilizamos f  x, y   c es decir para


nuestro ejemplo 2 x  y  c despejamos y  2 x  c damos valores a c  0,1, 2 , 3,...

Las isoclinas que son rectas paralelas como muestra la figura

Si hacemos f  x, y   0 con c  0 tenemos y  2 x que es la isoclina que se


encuentran los puntos de inflexión de las curvas solución, observamos en la figura
que la curva solución pasa de la región de decrecimiento a la región de
crecimiento y viceversa, por lo tanto los puntos son puntos mínimos.

Hallemos ahora veremos el lugar geométrico con la segunda derivada de


y  2 x  y es decir y  2  y  0 reemplazando y por su valor tenemos
y  2  y  2   2 x  y   0 despejamos y obtenemos y  2 x  2 que es la recta
isoclina cuando c  2 pero también es una curva solución.
14

La isoclina y  2 x es la que se encuentran los puntos mínimos de las curvas


solución y está situada por encima de la isoclina y  2 x  2

También observamos que la isoclina y  2 x  2 divide al plano xoy en dos partes y


observamos la que está situada por encima de la recta con y>0 nos indican que
las curvas solución tienen sus concavides para arriba y la que está situada por
debajo de la recta con y<0 nos indican que las curvas solución tienen sus
concavides para abajo.

En nuestro ejemplo la isoclina y  2 x  2 determina el lugar geométrico de


separación donde existen los puntos de inflexión, por debajo no existen puntos
extremos porque las curvas solución no se cortan.

Dicho análisis nos permite trazar aproximadamente las curvas solución de la


ecuación dada.

METODO DE EULER

El método de Euler consiste en sustituir por una quebrada de segmentos rectos la


curva solución buscada en la ecuación diferencial y  f  x, y  que pasa por el
punto M 0  x0 , y0  .

Se puede hallar en el segmento  x0 ,b


la solución aproximada de la ecuación
diferencial que satisface la condición inicial yx  x0  y0

Dividamos el segmento  x0 ,xn  b  en n partes, no necesariamente iguales por


los puntos x0 <x1<x2 <x3 <..............<xn  b .

M
2
M
1 M
3
M
0

x x x x x
0 1 2 3
15

Tracemos por el punto inicial M 0  x0 , y0  de la curva solución una recta M 0 M 1


cuya pendiente es y  f  x0 , y0  hasta el punto M 1  x1 , y1  de intersección con la
recta x  x1

Las ordenada del punto M n  xn , yn  se determina por la fórmula


yn  yn 1  f  xn 1 , yn 1   xn  xn 1 

Las ordenadas de los puntos M 1 ,M 2 ,M 3 de acuerdo a la fórmula son:

y1  y0  f  x0 , y0   x1  x0  para M 1 , y2  y1  f  x1 , y1   x2  x1  para M 2 y

y3  y2  f  x2 , y2   x3  x2  para M 3

Los valores aproximados de la curva solución de los puntos x1 ,x2 ,x3 ..............xn son
y1 , y2 , y3 ,................., yn

Realizando la construcción se obtiene una quebrada que se llama quebrada de


Euler que representa la curva solución que pasa por el punto M 0  x0 , y0  de
acuerdo a la figura.

Generalmente se divide el segmento  x0 ,b  en partes iguales y se utiliza la


notación h  xk  xk 1 y se llama intervalo de variación del argumento.

Ejemplo

Aplicando el método de las quebradas de Euler, hallar a) En el segmento  0,1 la


solución aproximada de la ecuación diferencial y  2 x  y que satisface a la
condición inicial yx 0  1 . y b) Dividir el segmento  0,1 en 10 partes iguales y
comparar los valores de la solución aproximada en los puntos de división con los
valores respectivos de la solución exacta y  2 x  2  e  x

Solución

Los valores de las ordenadas yk en los puntos M k  xk , yk  se calcula con la


fórmula yk  yk 1  f  xk 1 , yk 1   xk  xk 1  para k  1, 2 , 3 ,4 ,5 ,6 , 7 , 8 , 9 ,10
16

La función y  2 x  y  y  f  xk 1 , yk 1   2 xk 1  yk 1 y las diferencias


x1  x0  x2  x1  ........  x10  x9  0.1 porque se ha dividido el segmento dado en 10
partes iguales.

Por lo tanto yk  yk 1   2 xk 1  yk 1   0.1 para k  1, 2 , 3 ,4 , 5 ,6 , 7 , 8 ,9 ,10

De acuerdo a la condición inicial yx 0  1 Es decir x0  0; y0  1

Tenemos para k  1 tenemos y1  y0   2 x0  y0   0.1 o sea

y1  1   2  0    1   0.1  0.9

De acuerdo a la condición inicial yx 0.1  0.9 Es decir x1  0.1; y1  0.9

Tenemos para k  2 tenemos y2  y1   2 x1  y1   0.1 o sea

y2  0.9   2  0.1   0.9    0.1  0.79

De acuerdo a la condición inicial yx 0.2  0.79 Es decir x2  0.2; y2  0.79

Tenemos para k  3 tenemos y3  y2   2 x2  y2   0.1 o sea

y3  0.79   2  0.2    0.79    0.1   0.671

De acuerdo a la condición inicial yx 0.3  0.671 Es decir x3  0.3; y3  0.671

Tenemos para k  4 tenemos y4  y3   2 x3  y3   0.1 o sea

y4  0.671   2  0.3    0.671   0.1   0.5439

De acuerdo a la condición inicial yx 0.4  0.5439 Es decir x4  0.4; y4  0.5439

Tenemos para k  5 tenemos y5  y4   2 x4  y4   0.1 o sea

y4  0.5439   2  0.4    0.5439    0.1   0.40951


17

De acuerdo a la condición inicial yx 0.5  0.4095 Es decir x5  0.5; y5  0.4095

Tenemos para k  6 tenemos y6  y5   2 x5  y5   0.1 o sea

y6  0.4095   2  0.5    0.4095    0.1   0.26855

De acuerdo a la condición inicial yx 0.6  0.26855 Es decir x6  0.6; y6  0.26855

Tenemos para k  7 tenemos y7  y6   2 x6  y6   0.1 o sea

y7  0.26855   2  0.6    0.26855    0.1   0.121695

De acuerdo a la condición inicial yx 0.7  0.121695 Es decir x7  0.7; y7  0.12695

Tenemos para k  8 tenemos y8  y7   2 x7  y7   0.1 o sea

y8  0.121695   2  0.7    0.121695    0.1   0.0304745

De acuerdo a la condición inicial yx 0.8  0.0304745 Es decir


x8  0.8; y8  0.0304745

Tenemos para k  9 tenemos y9  y8   2 x8  y8   0.1 o sea

y9  0.0304745   2  0.8    0.0304745    0.1   0.13257295

De acuerdo a la condición inicial yx 0.9  0.13257295 Es decir


x9  0.9; y9  0.13257295

Tenemos para k  10 tenemos y10  y9   2 x9  y9   0.1 o sea

y10  0.13257295   2  0.9    0.13257295    0.1  0.060684345


18

Ahora con la solución exacta dada en el problema y  2 x  2  e  x

De acuerdo a la condición inicial yx 0  1 Es decir x0  0; y0  1

x
Tenemos y1  2 x0  2  e 0 y1  2  0   2  e  0  1

De acuerdo a la condición inicial yx 0.1  1 Es decir x1  0.1; y1  1

x
Tenemos y2  2 x1  2  e 1 y2  2  0.1  2  e  0.1  0.89516

De acuerdo a la condición inicial yx 0.2  0.89516 Es decir x2  0.2; y2  0.89516

x
Tenemos y3  2 x2  2  e 2 y3  2  0.2   2  e  0.2   0.78169

Así sucesivamente tenemos y4  0.6592; y5  0.5297; y6  0.3935

y7  0.2512; y8  0.1034; y9  0.0493; y10  0.2066

1.3.- ECUACIONES DIFERENCIALES DE VARIABLES SEPARADAS Y


VARIABLES SEPARABLES

Forma P  x, y  dx  Q  x, y  dy  0

La expresión es una ecuación con variables separables :

P  x, y  dx  Q  x, y  dy  0

dy P  x, y 
es equivalente a la ecuación diferencial 
dx Q  x, y 

No se hace ninguna distinción entre ecuaciones que contienen derivadas y


aquellas que contienen diferenciales.
19

Ejemplo 1

Resolver  1  x   1  y  dx  xydy  0
2 2

Solución

Hacemos  1  x   1  y  dx  xydy
2 2
luego

 1  x  dx 
2
ydy
observamos que la variables están separadas
x  1 y2 
1 
  x  dx
x   ydy resolviendo las integrales en ambos miembros
1 1 y 
2

1 2 1
ln x  x  ln  1  y 2   C que es una solución en forma implícita
2 2

podemos expresar en forma explícita

ln  1  y 2   2 ln x  x 2  2C sacamos antilogaritmos

 1 y   e 2ln x  x 2  2 C 2
2
 e2ln x e x e 2C el valor de la constante llamamos D  e 2C

 y   De 2ln x e x2 solución explícita


2
obtenemos 1  y 2  De 2ln x e x
20

dy
Forma  f  x g  y ; g  y  0
dx

Si la ecuación diferencial de primer orden es de la forma:

dy
 f  x g  y ; g  y  0
dx

decimos que las variables de la ecuación son separables, donde:

dy dy
 f  x  dx integrando encontramos la solución   f  x  dx  C
g  y g  y 

Ejemplo 2

Demostrar la ecuación:

dy
 x2  1  y 2 
1

dx

que es separable y encontrar la ecuación solución

Solución

Hacemos:

x 2 dx   1  y 2  dy observamos que la variables están separadas

Integramos

 x dx    1  y  dy
2 2

x3 y3 x3 y3
 y C   y C  x3  3 y  y 3  C como solución
3 3 3 3
21

Ejemplo 3

Resolver la ecuación diferencial

y 2 dx  2 y 2 dy  dy

Solución

Hacemos  y dx   2 y  1 dy dividimos a ambos miembros por y 2 tenemos:


2 2

 2y 2
 1  1 
dy   2  2  dy integramos   dx  
 2 y  1 dy   2  1  dy
2

 dx 
y2  y  y2   y 2 
1
x  2 y   C multiplicamos por y a ambos miembros
y

 xy  2 y 2  1  C donde x  2 y  y 1  C que es la solución

Ejemplo 4

Hallar la curva que pase por el punto  0 ,2  de modo que su pendiente de la
tangente en cualquiera de sus puntos sea igual a la ordenada del punto
aumentada en 3 unidades.

Solución

Basándonos en y  f  x, y 

Podemos obtener la ecuación diferencial, sabemos que su pendiente es la primera


derivada del 1er miembro de la ecuación diferencial y el segundo miembro es la
función de acuerdo al enunciado del problema la ordenada es y sumada en tres
unidades.

dy
Es decir y   f  x, y   y  3
dx
22

dy dy
Separando las variables  y  3  dy   y  3  dx   dx
dx  y  3
dy
e integrando   y  3    dx  Ln y  3  x  C pero como tenemos el punto  0 ,2 

que representa la condición inicial en la ecuación o sea yx 0  2 reemplazando


en Ln y  3  x  C dichos valores obtenemos el valor de la constante, es decir :

Ln 2  3  0  C  C  Ln 1  0 donde Ln y  3  x sacando antilogaritmos

y  3  ex  y  ex  3

Ejercicios

Resolver las ecuaciones diferenciales dadas para los siguientes ejercicios:

1) x 2 dy  cos 2 ydx  0

2) e dx  x  2  e  dy  0
y 2 y

3) e 2 x  y dx  2e x  y dy  0

4) Un depósito cilíndrico de volumen V0 está lleno de aire atmosférico, que se


k
P  V0 
comprime de un modo adiabático (Ecuación de Poisson  ) sin
P0  V 
intercambio de calor con el medio que le rodea, hasta que su volumen se hace
igual a V1 . Calcular el trabajo invertido durante la compresión.
23

1.4.- ECUACIONES DIFERENCIALES EXACTAS

Dada la Ecuación Diferencial

P  x, y  dx  Q  x, y  dy  0

P Q
Que satisface la condición:  entonces decimos que el 1er miembro de
y x
la ecuación es una diferencial exacta, por lo tanto existe una función f  x, y  tal
f f
que se cumpla df  x, y   P  x, y  dx  Q  x, y  dy donde P  x, y   ; Q  x, y  
x y

Vemos que cualquier función diferenciable satisface a la ecuación f  x, y   C

que es una función implícita.

Ejemplo 1

Resolver la ecuación diferencial

 x  2 y  2  dx   2 x  y  3 dy

Sabemos que P  x, y   x  2 y  2 ; Q  x, y   2 x  y  3
P Q
Hacemos  es decir 2  2 cumple con la ecuación y es una diferencial
y x
exacta.

Encontremos la solución: f  x, y  ?

f f
Sabemos P  x, y   ; Q  x, y  
x y

x2
f  x, y     x  2 y  2  dx   2 yx  2 x  C  y  es decir
2

2
 1 f  x, y   x  2 yx  2 x  C  y 
2
24

Entonces

f
 Q  x, y   2 x  C  y    2 x  y  3  C  y   y  3
y

 y2
Integrando C  y    C   y  dy     y  3 dy   3 y  C1
2

Reemplazando en  1 tenemos
x2 y2
f  x, y    2 yx  2 x   3 y  C1
2 2
2
x y2
es decir f  x, y    2 yx  2 x   3 y  C que es la solución
2 2

Ejercicios

En los ejercicios siguientes, resolver la ecuación diferencial dada:

1. x  x  y  4  dx   x  y  2  dy  0
2.  y  3x 2
 4  dx   x  y 2  2 y  dy  0

3.  2 xseny  e x
 
cos y  dx  x cos y  e x seny dy  0
2

1.5.- FACTORES INTEGRANTES

A veces una expresión de la forma

P  x, y  dx  Q  x, y  dy  0 no es una diferencial exacta.

pero existe una función I  x, y  tal que al multiplicar el primer miembro de la


ecuación diferencial :

I  x, y   P  x, y  dx  Q  x, y  dy  resulta una diferencial exacta, si I  0 entonces

coinciden las soluciones :

P  x, y  dx  Q  x, y  dy  0

I  x, y   P  x, y  dx  Q  x, y  dy   0
25

Definición

Por lo tanto si una función no nula I  x, y  convierte a la expresión


P  x, y  dx  Q  x, y  dy  0 en una diferencial exacta se llama FACTOR
INTEGRANTE.

Demostración

Partimos de la ecuación:

I  x, y   P  x, y  dx  Q  x, y  dy   0

es decir  I  x, y  P  x, y  dx  I  x, y  Q  x, y  dy   0

sabemos que dicha expresión para que sea una ecuación diferencial exacta tiene
que cumplirse:

  IP    IQ 
 y como son I ,P,Q funciones de las mismas variables f  x, y 
y x
derivamos parcialmente en cada miembro tenemos:

  P  I   Q  I
I P I Q realizando pasaje de términos
y y x x

 I  I     Q   P   I   I     P   Q 
P Q I   es decir Q P I  
y x  x y  x y  y x 

Obteniendo una ecuación de derivadas parciales que nos permite encontrar


el factor integrante considerando los siguientes casos particulares.
26

1er Caso

Hacemos I  I  x  depende solamente de x en este caso la derivada parcial


I
 0 es cero.
y

 I     P   Q 
La ecuación toma la forma Q 0  I    es decir
x  y x 
   P   Q 
   1
1   I   y x  pero como d  ln I   dI nos queda
 I
I x Q

   P   Q 
  
d  ln I   y ln I   f ( x )dx  I  e 
f ( x )dx
x  Integrando

dx Q

Para que exista un factor integrante no dependiente de y es necesario y


suficiente que el segundo miembro de la ecuación dependa solamente de x

Ejemplo 1

Determinar el factor integrante de la siguiente ecuación diferencial y encontrar la


solución.

 x  y  dx  2 yxdy  0
2

Solución
P Q
Verificamos si es una ecuación diferencial exacta  observamos 2 y  2 y
y x
que no es una diferencial exacta.

Hacemos P  x, y   x  y ;Q  x, y   2 xy
2

   P   Q 
  
Por la fórmula d  ln I   y x  2 y  2 y 4y 2
   
dx Q 2 xy 2 xy x
27

d  ln I  2 2
Es decir   resolviendo d  ln I    dx integrando
dx x x

 2
ln I   d  ln I       dx  2 ln x   ln x por la propiedad de los logaritmos
2

 x

2 2 1
ln I  ln x  ln I x  0 aplicando antilogaritmo tenemos Ix 2  e0  1  I 
x2

1
I Factor Integrante
x2

Comprobando

1  1 y2  2y
 x  y 2  dx  2 yxdy   0 es decir   2  dx  dy  0 verificando si es una
x2  
x x  x
diferencial exacta

P Q 2y
 o sea  2 y   x 2  es una diferencial exacta por lo tanto en la
y x x 2

ecuación

 1 y2  2y 1 y2 2y
  2  dx  dy  0 se separan la variables dx  2 dx  dy  0 o sea
x x  x x x x

1 y 2 dx  2 xydy 1 2 xydy  y 2 dx
dx  0 es decir dx   0 sacando diferenciación en el
x x2 x x2
1  y2 
segundo término del 1er miembro dx  d    0 donde integramos
x  x 

1  y2  y2 y2
 x   x 
dx  d  0 tenemos ln x 
x
 0  ln x 
x
sacando el antilogaritmo

y2
nos queda x  Ce x que es la solución.
28

2do Caso

Hacemos I  I  y  depende solamente de y en este caso la derivada parcial


I
 0 es cero.
x
 I    P   Q 
La ecuación toma la forma 0  P I   es decir
y  y x 
   Q   P 
   1
1   I   x y  pero como d  ln I   dI nos queda
 I
I y P

   Q   P 
  
d  ln I   x y  Integrando ln I   g( y )dy  I  e 
g( y )dy


dy P

Para que exista un factor integrante no dependiente de x es necesario y


suficiente que el segundo miembro de la ecuación dependa solamente de y

Ejemplo 2

Determinar el factor integrante de la siguiente ecuación diferencial y encontrar la


solución.

 
2 xy ln ydx  x 2  y 2 y 2  1 dy  0

Solución

Sabemos P  x, y   2 xy ln y y Q  x, y   x 2  y 2 y 2  1

   Q   P 
  
Por la fórmula d  ln I   x y   2 x  2 x ln y  2 x   2 x ln y   1

dy P 2 xy ln y 2 xy ln y y
29

d  ln I  1 1
Es decir    d  ln I    dy Integrando
dy y y

1
ln I    dy   ln y por la propiedad de los logaritmos ln I  ln y  o  ln Iy  0
y

1 1
sacando antilogaritmo Iy  e0  1  I  es decir I  Factor Integrante
y y

Ahora encontraremos la solución, a la ecuación diferencial dada dividimos por y

2 xy ln ydx


x 2  y 2 y 2  1 dy
0 o sea
 
2 x ln ydx  x 2 y 1  y y 2  1 dy  0 
y y
comprobamos que es una diferencial exacta

P Q 2x 2x
 es decir 
y x y y

x2
Hacemos f  x, y    2 x ln ydx  2 ln y  C  y  derivamos la función con respecto
2
a despejamos y tenemos

f
 x 2 y 1  C   y   Q  x, y   x 2 y 1  y y 2  1 despejamos C   y  nos queda
y

C  y   y y2  1 integrando
C  y    C   y  dy   y y 2  1dy 
y 2
 1 3
 C1
3
reemplazando en f  x, y  tenemos

f  x, y   x 2 ln y 
y 2
 1 3
 C1
que es la solución
3

Por lo general no existen reglas para determinar los factores integrantes, pero en
algunos casos sencillos se emplean expresiones conocidas como factores
integrantes.
30

Ejemplo 3

Resolver x dy   xy  x  1 dx  0
2 2

Solución

P Q
Si hacemos  observamos x  2 x no es una diferencial exacta
y x

Sin embargo notamos x 2 dy  xydx  x  xdy  ydx   xd  xy  es decir

xdy  ydx  d  xy 

Si dividimos a la ecuación original por x , tenemos:

 1
xdy  ydx   x   dx  0 es decir
 x

 1  1
d  xy    x   dx  0 integrando  d  xy     x  x  dx  0 encontramos las
 x
solución:

1 2
xy  x  ln x  C que es la solución
2

1
El factor integrante que hemos utilizado es I  x, y  
x

La clave de esta solución es la observación de xdy  ydx  d  xy  , existen otras que


nos permitirán la resolución de otras ecuaciones que llevan el factor integrante.

Algunas de estas diferenciales son:

1. xdy  ydx  d  xy 
xdy  ydx  y
2. 2
d 
x x
1
3. xdx  ydy  d  x  y 
2 2

2
xdy  ydx  y
4.  d  arctan 
x y
2 2
 x
xdy  ydx 1
5.   ln  x 2  y 2  
x y
2 2
2
31

Observamos en estos ejemplos de diferenciales la facilidad de poder darse cuenta


el factor integrante que se necesita en cada ejercicio particular.

Ejemplo 4

Resolver
xdy  ydx  2 x 2 ln ydy

Solución

P Q
Observamos que no cumple  entonces no es exacta
y x
xdy  ydx  y
De acuerdo a la diferencial conocida 2
 d   , observamos que le falta
x x
1
en el denominador x 2 por lo tanto ensayamos el factor integrante I  x, y   2 el
x
xdy  ydx
resultado es  2 ln ydy  0 , verificamos si es exacta, comprobamos
x2

xdy ydx ydx  1  P Q


 2  2 ln ydy  0   2    2 ln y  dy  0 hacemos 
x 2
x x x  y x

1 1
tenemos  2
  2 es un diferencial exacta
x x

 y  xdy  ydx
Para llegar a la solución tomamos d     2 ln ydy
x x2

 y
Integrando ambos miembros  d  x    2 ln ydy es decir

y
 2  y ln y  y   C
x

Ejercicios

En los siguientes ejercicios resolver las ecuaciones diferenciales dadas

1. xdy  ydx   x  y  dx
2 2

2. xdy  ydx   x  y  dy
2 2

3.  2  xy  ydx   2  xy  xdy  0
32

1.6.- ECUACIONES HOMOGENEAS

Una función f  x, y  es homogénea de grado n si cumple

f  tx,ty   t n f  x, y 

Ejemplo

Verificar si la función f  x, y   x  y  xy es homogénea y de que grado es.


2 2

Solución

Hacemos f  tx,ty    tx    ty    tx   ty   t  x  y  xy   t f  x, y 
2 2 2 2 2 2

Vemos que una función homogénea de grado 2

Una función homogénea de grado cero es n  0

Ejemplo

x2  y2
f  x, y   es una función de grado cero porque
x2  y2

 tx    ty 
2 2
x2  y2
f  tx,ty     f  x, y 
 tx    ty  x2  y 2
2 2

dy
Ahora decimos que una ecuación diferencial de la forma y   f  x, y 
dx
es homogénea si la función es homogénea de grado cero.

dy  y
Ecuaciones Homogéneas de la forma  
dx x
33

dy  y
La ecuación homogénea siempre se representa de la forma:  
dx x
y
Introducimos una nueva función incógnita u  despejamos y  ux y sacamos
x
 du 
su diferencial dy   u  x  dx es decir
 dx 

dy du du dy y du
ux x   u      u    u   u o sea x    u  u
dx dx dx dx x dx

Si u  u0 es una raíz de la ecuación   u   u  0

La solución de la ecuación homogénea es u  u0 reemplazando por su variable


y y
original u  tenemos  u0  y  u0 x que es una recta que pasa por el origen
x x
de coordenadas.

Como recomendación, al resolver las ecuaciones diferenciales homogéneas no es


dy  y
necesario reducirlas a la forma     se puede realizar inmediatamente la
dx x
sustitución y  ux .

dy  a x  b1 y  c1 
Ecuaciones de la forma  f 1 
dx  a2 x  b2 y  c2 

dy  a x  b1 y  c1 
Las ecuaciones de la forma  f 1  se reducen a homogéneas
dx  a2 x  b2 y  c2 
trasladando el origen de coordenadas al punto  x0 , y0  de intersección de las
rectas x    x0 ; y    y0

Esto se realiza haciendo una sustitución de las variables

x    x0 ; y    y0 Ecuaciones de Traslación

Este método no se aplica cuando las rectas son paralelas.

a2 b2
Si las rectas x    x0 ; y    y0 son paralelas, sabemos que  
a1 b1

dy  a x  b1 y  c1 
La ecuación  f 1  se puede escribir de la forma
dx  a2 x  b2 y  c2 
dy  a x  b1 y  c1  dy
 f 1   F  a1 x  b1 y  es decir  F  a1 x  b1 y  cuya resolución
dx    a1 x  b1 y   c2  dx
se reduce a una ecuación con variables separables.

Ecuaciones Diferenciales de la forma P  x, y  dx  Q  x, y  dy  0


34

La ecuación diferencial dada será homogénea si P  x, y  ;Q  x, y  son funciones


homogéneas de un mismo grado.

Para reducir la ecuación a homogénea se realiza la sustitución de la variable

y  z siempre y cuando todos los términos de la ecuación son de un mismo


grado , colocando el grado uno a la variable x , el grado  a la variable y y el
dy
grado   1 a la deriva .
dx

Ejemplo

Resolver la ecuación xy   x 2  y 2  y

Solución

Transformamos la ecuación en homogénea haciendo y  1    


y y
x x

y
Es decir u  o bien y  ux derivando tenemos y  xu  u entonces
x
reemplazando los valores de y  ux, y  xu  u nos queda

du
xu   u  1  u 2  u o sea xu  1  u  x  1  u2
2

dx

Hacemos cambio de variables separables

du dx du dx
 Integrando   o sea arcsenu  ln x  ln C1 ;  C1  0 
1u 2 x 1u 2 x

Aplicando la propiedad de logaritmos tenemos arcsenu  ln C1 x pero sabemos

C1 x  C1 x haciendo la notación C1  C obtenemos arcsenu  ln Cx

Aplicando el seno a ambos miembros tenemos u  sen ln Cx reemplazando el valor


y
de u  es decir y  xsen ln Cx que es la solución general.
x

Ejemplo

Resolver la ecuación
35

 x  y  2  dx   x  y  4  dy  0
Solución

Planteamos el sistema de ecuaciones algebraicas lineales:

x y20 x y  2

x  y  4  0  x  y  4

Consideramos el determinante principal

1 1
  2  0 o sea que existe un punto de intersección donde las recta se
1 1
interceptan, resolvemos el sistema de ecuaciones:

2 1 1 2
4 1 2 y 1 4  6
x0    1 y0   3
 2  2

Hacemos x    1; y    3 sacamos su diferenciales dx  d  ;dy  d

Entonces la ecuación  x  y  2  dx   x  y  4  dy  0 toma la forma

   1    3  2  d     1    3  4  d  0 es

decir      d      d  0 es una ecuación homogénea

Hacemos   u y su diferencial d  ud    du obtenemos

   u  d     u   ud   du   0 donde

 d  u d    2 du   2udu   ud   u 2 d  0 o sea

   u   u  u   d   
2 2
  2u  du  0 es decir

 1  2u  u  d    1  u  du  0
2
separamos las variables

d  1  u  du  0 d  1  u  du  0

  1  2u  u 2 
integrando  

 1  2u  u 2 
 1  u  du
En la integral   1  2u  u  2 hacemos un cambio de variables

z  1  2u  u 2 ;dz  2  1  u  du tenemos

 1  u  du  1 u dz 1
  1  2u  u   
2
z 2 1 u
 1
2  z dz  1
2 ln z  12 ln 1  2u  u 2

Nuestra integral resuelta es:


36

d  1  u  du  ln   1 ln 1  2u  u 2  ln C
 

 1  2u  u 2  2

Aplicando la propiedad de logaritmos

ln   12 ln 1  2u  u 2  ln C es decir ln   1  2u  u 2  2  ln C
1

Sacamos antilogaritmos nos queda:

  1  2u  u 2  2  C es decir   1  2u  u   C donde C1  C
1
2 2 2 2

Volvemos a las variables x, y

 y  3  y  3 
2

 x  1 1  2  C1 o sea  x  1  2  y  3   x  1   y  3    C1
2 2 2
 2 
 x  1  x  1   

x 2  2 x  1  2 yx  2 y  6 x  6  y 2  6 y  9  C1 es decir

x 2  2 yx  4 x  y 2  8 y  C1  14 ó x 2  2 xy  y 2  4 x  8 y  C2 donde

C2  C1  14

Ejemplo

Resolver la ecuación

 x  y  1 dx   2 x  2 y  1 dy  0
Solución

Planteamos el sistema de ecuaciones algebraicas lineales

x  y 1  0

2x  2 y  1  0

Es decir
x  y  1
2x  2 y  1

Resolvemos por medio de determinantes

1 1
  2  2  0 En este caso la resolución no es aplicable que en el ejemplo
2 2
anterior.

En este caso podemos hacer la sustitución de variables x  y  z  dy  dz  dx

Para integrar, entonces la ecuación toma la forma


37

 x  y  1 dx   2 x  2 y  1 dy  0   z  1 dx   2 z  1  dz  dx   0 o sea

zdx  dx  2 zdz  2 zdx  dz  dx  0 es decir  2  z  dx   2 z  1 dz  0

2z  1
separando las variables obtenemos dx  dz  0 integramos
2 z

2z  1 2z 1
 dx   2 z
dz  0 es decir  dx   2  z dz   2  z dz  0 hacemos

2z 1
 dx   z  2 dz   z  2dz  0 resolviendo utilizando las formulas

x x b 1 1
 ax  b dx  a  a 2
ln ax  b y  ax  b dx  a ln ax  b tenemos

x  2  z  2 ln z  2   ln z  2  C o sea x  2 z  3 ln z  2  C

Volviendo a las variables originales x, y obtenemos x  2  x  y   3 ln x  y  2  C

Ejemplo

Resolver la ecuación

x 2
y 2  1 dy  2 xy 3 dx  0

Solución

Hacemos la sustitución y  z  dy   z 1dz donde consideramos por ahora 


es un número arbitrario que se determinará.

Sustituyendo en  x y  1 dy  2 xy dx  0 tenemos  x z  1  z dz  2 xz dx  0
2 2 3 2 2  1 3

O sea  x z  z  1   dz  2 xz 3 dx  0
2 3 1

Para saber si la ecuación diferencial es homogénea todos los grados de


todos los términos tienen que ser iguales.

Observamos la suma de los exponentes de las dos variables x,z en todos los
términos.

Vemos que el grado del término x 2 z 3 1 es 2  3  1  3  1 ahora el grado del

término z 1 es   1 y el grado de xz 3 es 1  3 .

Para que sean iguales es necesario que se cumpla 1  3    1 donde   1

 1 1 1
Entonces sabemos que y  z  y  z  es decir y 
z z

La ecuación diferencial x z
2 3 1
 z 1   dz  2 xz 3 dx  0 toma la forma
38

1
z 2
 x 2 z 4  dz  2 xz 3 dx  0 o bien multiplicando por
z4
tenemos

z 2
 x 2  dz  2 xzdx  0

Resolvamos la ecuación diferencial homogénea haciendo

z
u  z  ux  dz  xdu  udx es decir  u 2 x 2  x 2   xdu  udx   2 xuxdx  0 o sea
x

dividiendo por x 2 tenemos u 2


 1  xdu  udx   2udx  0 hacemos cambio de

variables separables tenemos u 2 xdu  u 3 dx  xdu  udx  2udx  0 o

u2  1
u 2
 1 xdu   u 3  u  2u  dx  0 nos queda
u u
3
1
du  dx  0 integramos
x

u2  1 1
 u 3  u du   xdx  0 Resolvemos utilizando las fórmulas de las tablas de
integrales

1.7.- ECUACIONES DIFERENCIALES LINEALES DE PRIMER ORDEN

La ecuación lineal de primer orden y de grado n es:

n n 1 n2
 dy   dy   dy  dy
c0  x     c1  x     c2  x     ........  cn 1  x   cn  x  y  f  x 
 dx   dx   dx  dx

La ecuación lineal más general de primer orden es:

dy
c0  x   c1  x  y  f  x 
dx

Dividiendo a esta ecuación por c0  x 


dy c1  x  f  x dy
 y es decir  P  x  y  Q  x  cuando c0  x   0
dx c0  x  c0  x  dx

Observamos que la ecuación posee un FACTOR INTEGRANTE I  x  que es


función de x solamente.

dy
Es decir si multiplicamos por I  x  a la ecuación  P  x y  Q  x
dx

dy
Tenemos I  x  I  x  P  x  y  I  x  Q  x  (1)
dx

dy dI
Si derivamos d  I  x  y  x    I  y  I  x  Q  x  (2) observamos coinciden
dx dx

ambas ecuaciones (1) y (2)


39

dI
Por comparación de (1) y (2) hacemos  I  x  P  x  dicha ecuación es una
dx
ecuación diferencial de variables separables, resolvemos:

dI dI
 P  x  dx integrando  I  x    P  x  dx tenemos
I  x

lnI  x    P  x  dx sacando el antilogaritmo tenemos

I  x   e
P  x  dx
obtenemos el Factor Integrante.

Observamos que la ecuación podemos integrar fácilmente

dy dI
d  I  x  y  x    I  y  I  x Q  x es decir
dx dx

 d  I  x  y  x     I  x  Q  x  dx resolviendo

I  x  y  x    I  x  Q  x  dx reemplazando el valor de I  x   e 
P  x  dx
en la ecuación
anterior tenemos

e y  x    e
P  x  dx P  x  dx
Q  x  dx  C estamos suponiendo que dicha integral está
resuelta del segundo miembro, despejando y  x 

y  x  e   P x  dx
 P  x  dx   P  x  dx
 e Q  x  dx  Ce que es la solución en forma explicita

Ejemplo 1

dy
Resolver:  2 xy  2 x 3
dx

Solución

dy
De acuerdo a la ecuación  P  x y  Q  x hacemos P  x   2 x y
dx
Q  x   2 x 3 , también hay que recordar la fórmula del factor integrante I  x   e  P x  dx

2x2
 P  x  dx   2 xdx   x 2 o sea I  x   e x como factor integrante.
2
Es decir
2

dy 2 dy
Multiplicamos I  x   e x a
2 2 2
 2 xy  2 x 3 tenemos e x  2 xye x  2 x 3e x
dx dx

Que es lo mismo 
d ex y
2

  2x e
3 x2 hacemos  2

 2
d e x y  2 x 3e x dx
dx

de 
y   2 x 3e x dx resolviendo por integración por partes el
x2 2
integrando
segundo miembro.
2
ex
hacemos u  x ; dv  xe dx es decir du  2 xdx; v   xe dx 
2 x2 x2

2
40

resolvemos

 x 2e x 
2 2

  ex
e y   2 x e dx  2  udv  2 uv   vdu   2 
x2 3 x2
 2 xdx   x 2 e x  2 e x xdx
2 2

 2 2 

volvemos a integrar
2
2e x
x2 2 x2
 e x  x 2  1  C es decir
2
e yx e 
2

C
y   x 2  1  2 que es la solución
ex

Ejercicios

Resolver las siguientes ecuaciones diferenciales de 1er orden

dy
1.  y   3 x 2  4 x  3 e x
dx
dy xy
2.  2 x
dx x  1
dy 3y x
3.  
dx 2 x  3 2x  3

1.8.- METODOS DE VARIACION DE LOS PARAMETROS

ECUACIONES DIFERENCIALES DE LA FORMA f  y, y    0; f  x, y   0

S i en estas ecuaciones se puede despejar y podemos resolver con variables


separables.

FORMA f  y, y   0

Caso a)

FORMA DONDE SE PUEDE DESPEJAR y de f  y, y   0

Despejamos la variable y tenemos y    y 

Hacemos y  p o sea y    p  diferenciando esta ecuación dy     p  dp

dy
Sabemos que y   p  dy  pdx tenemos pdx     p  dp despejamos
dx

   p  dp    p  dp    p  dp
dx 
p
integrando  dx   p
 C es decir x
p
C

que es la solución.

Resumen

Obtenemos la solución general de la ecuación en forma paramétrica


41

   p  dp
x C
p

y    p

Ejemplo

2 3
 dy   dy  a,b  ctes
Resolver la ecuación diferencial y  a    b  
 dx   dx 

Solución

dy
Hacemos p o sea y  ap 2  bp 3 aplicamos el diferencial tenemos
dx

dy   2ap  3bp 2  dp reemplazando dy  pdx es decir pdx   2ap  3bp 2  dp

3bp 2
dx  2adp  3bpdp integrando x  2ap  C
2

La solución general es

y  ap 2  bp 3

Caso b)

FORMA DONDE NO SE PUEDE DESPEJAR y de f  y, y   0

Observamos que no se puede despejar la variable y ni y o si se puede despejar


es con dificultad.

Podemos trabajar para solucionar este problema en forma paramétrica mediante


dy
el parámetro t es decir hacemos y    t  ; p    t  ; p 
dx

Entonces dy  pdx    t  dx sacando diferencial a y    t  tenemos dy     t  dt

Reemplazando el valor de dy    t  dx en la anterior ecuación nos queda

 t   t 
  t  dx     t  dt es decir dx  dt integrando x   dt  C
  t   t

que es la solución.

Resumen

Obtenemos la solución general de la ecuación en forma paramétrica

 t 
x dt  C y   t
  t

Ejemplo
Resolver la ecuación diferencial y 3   y  3  1
2 2
42

Solución

dy
Hacemos y  cos 3 t y y  p  sen 3t pero como y   p  sen 3t
dx

dy dy
Despejamos dx   pero sacando su diferencial a y  cos 3 t tenemos
p sen 3t

dy  3 cos 2 tsentdt Reemplazando

dy 3 cos 2 tsentdt 3 cos 2 tdt


en dx     3 cot ag 2tdt es decir
sen 3t sen 3t sen 2t

dx  3 cot ag 2tdt integrando

x  3 cot ag 2tdt  3   cot agt  t   C  3 cot agt  3t  C

Resumen

Obtenemos la solución general de la ecuación en forma paramétrica

x  3 cot agt  3t  C

y  cos 3 t

FORMA f  x, y   0

Caso a)

FORMA DONDE SE PUEDE DESPEJAR x de f  x, y   0

Despejamos la variable x tenemos x    y 

Hacemos y  p es decir x    p  Sacamos el diferencial dx     p  dp

dy dy
Pero y   p despejamos dx  reemplazamos el valor dx o sea
dx p

dy
   p  dp  o dy  p   p  dp integrando y   p   p  dp  C
p

Resumen

Obtenemos la solución general de la ecuación en forma paramétrica


43

y   p   p  dp  C

x    p

Ejemplo

2
dy  dy 
Resolver la ecuación diferencial a  b    x
dx  dx 

Solución

dy
Hacemos y   p (1) reemplazando en la ecuación diferencial dada tenemos
dx

x  ap  bp 2 sacamos el diferencial o sea dx  adp  2bpdp

De (1) despejamos dy  pdx y reemplazamos el valor de dx tenemos


ap 2 bp 3
dy  apdp  2bp 2 dp integrando y  2 C
2 3

Resumen

Obtenemos la solución general de la ecuación en forma paramétrica

ap 2 bp 3
y 2 C
2 3

x  ap  bp 2

Caso b)

FORMA DONDE NO SE PUEDE DESPEJAR x de f  x, y   0

Observamos que no se puede despejar la variable x ni y o si se puede despejar


es con dificultad.

Podemos trabajar para solucionar este problema en forma paramétrica mediante


dx
el parámetro t es decir hacemos x    t  ; p    t  ; p 
dy

Entonces dx  pdy    t  dy sacando diferencial a x    t  tenemos dx     t  dt

Reemplazando el valor de dx    t  dy en la anterior ecuación nos queda

 t    t 
  t  dy     t  dt es decir dy  dt integrando y   dt  C
  t   t

que es la solución.

Resumen
44

Obtenemos la solución general de la ecuación en forma paramétrica

  t 
y dt  C x    t
  t

Ejercicios

Resolver las siguientes ecuaciones diferenciales:

1) y   y  e y
2 

2) x  ln y  seny

3) y  y ln y

1.09.- ECUACIONES DE LAGRANGE

La ecuación de lagrange tiene la forma y    y   x  y 

dy
Hacemos y  p es decir y    p   x  p  pero sabemos y   p  dy  pdx
dx

A la ecuación y    p   x  p  diferenciamos tenemos

dy     p  dp    p  dx     p  dp reemplazamos por dy  pdx nos queda

pdx     p  dp    p  dx     p  dp es decir  p    p   dx     p     p   dp

   p     p      p     p  
dx   dp integrando x   dp  C  r  p,C 
 p    p    p    p  

o sea x  r  p,C 

Resumen

La solución general de la ecuación diferencial de Lagrange es:

x  r  p,C 

y    p   x  p 

Ejemplo

Aplicando el método de Lagrange encontrar la solución de la ecuación diferencial:

y  2 xy  ln y
45

Solución

Hacemos y  p entonces y  2 xp  ln p diferenciamos dy  2 pdx  2 xdp  p 1dp

dy
Sabemos y   p  dy  pdx entonces pdx  2 pdx  2 xdp  p 1 dp donde
dx

pdx   2 x  p 1  dp  0

dividiendo toda la ecuación por pdp nos queda

dx 1 1
 2 xp 1  p 2  2 x  2 es una ecuación diferencial de 1er orden
dp p p

dx 1 1
 2 x   2 Resolviendo la ecuación diferencial de 1er Orden tenemos:
dp p p

2
Planteamos el Factor Integrante I  p   e  p
dp
 e 2 Lnp  eln p  p 2 tenemos
2

dx
p2  2 px  1 sabemos d  p 2 x   p 2 dx  2 pxdp dividimos en ambos miembros
dp
d  p 2 x  dx d  p 2 x 
por dp nos queda p 2
 2 px  1 o sea  1 es decir
dp dp dp

d  p 2 x   dp integrando  d  p x     dp
2
nos queda p 2 x   p  C despejando

C 1
x hallamos x  reemplazando el valor
p2 p

2C
hallado de x nos queda y  ln p  2
p

1.10.- ECUACIONES DE CLAIRAUT

La ecuación de Clairaut es de a forma

y  xy    y 

El método de resolución es el mismo que para la ecuación de Lagrange

La solución general de la ecuación de Clairaut tiene la forma

y  Cx   C 

La ecuación de Clairaut puede tener también una solución singular que se obtiene
eliminando p entre las ecuaciones

y  xp   p 

Diferenciando en ambos miembros tenemos

dy  pdx  xdp     p  dp
46

dy
Sabemos y  p   dy  pdx de las dos ecuaciones igualamos dy nos queda
dx

pdx  pdx  xdp    p  dp  xdp     p  dp  0   x     p   dp  0 es decir

x   p   0

Ejemplo

Aplicando el método de la ecuación de Clairaut resolver la ecuación diferencial


dada

a
y  xy  ;a  cte
2 y

Solución

a a 2
Hacemos y  p obtenemos y  xp  diferenciando dy  pdx  xdp  p dp
2p 2

dy a
Donde y   p  dy  pdx reemplazando pdx  pdx  xdp  p 2 dp es decir
dx 2

a 2  a 
xdp  p dp  0 o  x  p 2  dp  0 analizamos los dos factores de esta
2  2 

ecuación tenemos:

Si dp  0 integrando p   dp  0  C  p  C y la solución general de ecuación es

a
y  Cx 
2C

a 2 a
Si x  p 0 x  eliminando p esta ecuación con la solución general
2 2 p2

a a a a
y  xp  tenemos x   p2   p reemplazando
2p 2p 2
2x 2x

2
 a 
2 x   a
a 2 xp 2  a 2 x 2a a elevando al cuadrado
y  xp       
2p 2p a a a
2 2
2x 2x 2x

a a2
y y2   2ax  y 2  2ax
a nos queda a
2x 2x

1.11.- ECUACION DE BERNOULLI


47

La ecuación de Bernoulli se expresa

dy
 P  x y  Q  x yn ; n  1
dx

Se puede reducir a una ecuación lineal, dividiendo ambos miembros por y n

dy dz dy
obteniendo y n  P  x  y1 n  Q  x  si hacemos z  y1 n y   1  n  y n
dx dx dx

1 dz dy
o sea  y n
1  n dx dx

1 dz
entonces se transforma  P  x  z  Q  x  que es una ecuación lineal
1  n dx

z
una vez determinada z , la sustitución y  da la solución; veremos en un
1 n
ejemplo

Ejemplo 1

Resolver

 1  x  dy
2

dx
 xy  x y 3 3

Dividimos a ambos miembros por  1 x 


2
resulta una ecuación de Bernoulli de
tercer grado con n  3 y dividiendo por y 3
tenemos

dy
3 xy 2 x3
y  
dx  1  x 2   1  x 2 

Haciendo z  y1 n  y13  y 2 y obtenemos su diferencial con respecto a x o sea

dz dy
 2 y 3 observando la ecuación anterior vemos que falta 2 , entonces
dx dx
multiplicamos ambos miembros por 2 y nos queda

dy 2 xy 2 2 x3 dz 2 xz 2 x3
2 y 3    es decir   
dx  1  x 2   1  x2  dx  1  x 2   1  x2 
ecuación Bernoulli

dy
 P  x y  Q  x yn ; n  1
dx

2x 2 x3
por semejanza tenemos P  x    ;Q  x  
 1  x2  1  x2

I  x   e
P  x  dx
Para encontrar el factor integrante utilizamos es decir

2 x
 1 x2 dx
1
 
2  ln 1 x 2 

 ln 1 x 2 
I  x  e e 2 
e por propiedades de logaritmos tenemos
48

  el factor integrante es I  x    1  x 2 
  1  x2 
 ln 1 x 2 1 1
I  x  e multiplicamos

dz 2 xz 2 x3
miembro a miembro la ecuación   
dx  1  x 2   1  x 2  obtenemos
dz 2 xz 2 x3
 1  x2      es decir
1 2 1 2 1
 1  x   1  x
dx  1  x 2   1  x2 
d  1  x 2  z 
1
dz 2 xz 2 x3
 1 x 
2 1

dx  1  x 2  2
  o sea     2x
3
integrando
 1  x2 
2
dx  1  x 
2 2

x3
 d  1  x  z   2 1  x 2 2 dx resolviendo la integral tenemos
 2 1 

 
 1 1  1
 1 x 
2 1
z  2   ln  1  x 2
 C    ln  1  x 2   C
 2  1  x  2 1 x 
2 2


es decir

z  1   1  x 2  ln  1  x 2   C  1  x 2  pero sabemos que z  y 2 reemplazamos

1
y 2  1   1  x 2  ln  1  x 2   C  1  x 2  es decir y 
2

1   1  x  ln  1  x 2   C  1  x 2 
2

sacando raíz cuadrada encontramos la solución

1
y
1   1  x  ln  1  x 2   C  1  x 2 
2

Ejemplo 2

Resolver y  6 y 2  x  1 dx  2 xdy  0

Agrupamos los términos de acuerdo a las potencias de y

6 y 3 dx  y   x  1 dx  2 xdy  0 es decir

dy
2x  y  x  1  6 y 3
dx

Resulta una ecuación de Bernoulli con n  3 si hacemos z  y1 n  y13  y 2

dz  2 y 3dy o sea 2 xdy  y  x  1 dx  6 y dx dividimos


3
Diferenciando tenemos
por  y3 tenemos  x 2 y 3 dy  y 2  x  1 dx  6dx donde xdz  z  x  1 dx  6dx
dz
dividiendo por x o sea dz  z  1  x 1  dx  6 x 1dx  z  1  x 1   6 x 1
dx
49

I  x   e
P  x  dx
Para encontrar el factor integrante utilizamos es decir

 1 x  dx I  x   xe x factor integrante
I  x   e
1

 e
x  ln x 
 e x e ln x  xe x es decir

dz
xe x  ze x  x  1  6e x es decir xe x dz   z  x  1  6 e x dx  0 que es una diferencial
dx

P Q
 o sea e  xe  e  x  1
x x x
exacta es decir
x z

Integrando la ecuación xe dz   z  x  1  6 e dx  0


x x

xe x  dz    z  x  1  6  e x dx  0 es decir xe x  dz   zxe x dx   ze x dx   6e x dx  0

xe x z  z  xe x  e x   ze x  6e x  C  0 dividimos por e x

xz  z  x  1  z  6  Ce  x  0 es decir xz  zx  z  z  6  Ce  x  0 o sea

2 xz  6  Ce  x  0 reemplazando su variable original tenemos z  y 2

y 2 tenemos 2 x  y  6  Ce   0 es decir
2 x
2 xy 2  6  Ce  x  0 multiplicamos por

y 2  6  Ce  x   2 x
2x
donde y  que es la solución
6  Ce  x

Ejercicios

Resolver las siguientes ecuaciones diferenciales, planteando la ecuación de


Bernoulli

1. 6 y dx  x  2 x  y  dy  0 en la variable dependiente x
2 3

2. 2 x dy  y  y  2 x  dx  0
3 2 2

1.12.- APLICACIONES, MECANICA, ETC.

Las ecuaciones diferenciales se aplican tanto en la mecánica, física e ingeniería


como en las ramas de las matemáticas, veremos algunas aplicaciones de
50

ejemplos sencillos de mecánica en fórmulas elementales como las Leyes de


Newton y algunas aplicaciones a problemas geométricos.

En estas aplicaciones es importante determinar primeramente la ecuación


diferencial que describe la aplicación dada y luego hallar la solución particular de
la aplicación particular.

Veremos algunos ejemplos típicos:

Ejemplo 1

Hallar la ecuación de la curva tal que su pendiente en un punto P  x, y  sea igual a


la diferencia de los cuadrados de las distancias de P  x, y  a los puntos  1, 0  ;  4, 0 

Las distancias de P  x, y  a  1, 0  ;  4, 0  designamos d1 y d 2

Planteamos la siguiente ecuación diferencial de acuerdo al enunciado

dy
pero como d12   x  1  y 2 ; d 22   x  4   y 2 tenemos
2 2
 d12  d 22
dx

dy 
  x  1  y 2    x  4   y 2   x 2  2 x  1  y 2  x 2  8 x  16  y 2  6 x  15
2 2

dx    

dy
luego integramos  6 x  15 o sea dy   6 x  15  dx
dx
6 x2
y   dy    6 x  15  dx   15 x  C donde
2

y  3 x 2  15 x  C que es la solución

Ejemplo 2
51

Cada pie de un cable uniforme pesa w libras. Suponiendo que este cable es
flexible y está suspendido de sus extremos, hallar la ecuación de la curva que
forma el cable.

De acuerdo a la figura vemos que las fuerzas sobre un objeto deben estar
balanceadas.

Las fuerzas que actúan sobre el cable son aquellas debidas a la gravedad y a la
tensión del cable que no le permite caer.

Sea A  a,b  el punto más bajo del cable y P  x, y  otro punto con x > a , la
componente horizontal de la tensión del cable que no le permite caer, será una
constante H, luego designamos con L  x  a la fuerza vertical de la gravedad
dirigida hacia abajo.

Si  es ángulo de la tangente P, descomponemos el vector tensión T  x  que es la


tangente del cable en sus componentes vertical y horizontal.

Designando con L  x  , T  x  y H las magnitudes del los vectores, tenemos:

T  x  sen  L  x  ;T  x  cos   H sabemos que la tangente es

sen dy L  x 
tan     la fuerza de gravedad L  x  se debe solamente al peso
cos  dx H

del cable.

Designamos s  x la longitud del cable tenemos L  x   ws  x  reemplazando


dy L  x  ws  x  w
en la ecuación de la tangente   hacemos  k  cte es decir
dx H H H
dy
 ks  x  y teniendo en que la longitud del cable se escribe por esta ecuación
dx
2 2
 dy  ds  dy 
ds  1    dx .   1 
 dx  dx  dx 

De la ecuación de la tangente tenemos dy  ks  x  dx y sacamos nuevamente su


d2y
diferencial tenemos d 2 y  kdsdx 2 es decir  kds reemplazando el valor ds
dx 2
2
d2y  dy 
tenemos 2
 k 1 
dx  dx 

Ejemplo 3
52

Un resorte de peso despreciable está suspendido verticalmente y su extremo libre


está sujetado a una masa m  kg  .
 m 
Si la masa se mueve a una velocidad v0   , cuando el resorte está sin alargar,
 seg 
hallar la velocidad en función del alargamiento x  mts 

Solución

Según la Ley de Hooke, la fuerza del resorte (fuerza opuesta al alargamiento) es


proporcional al alargamiento.

Es decir Ft  ma  W  FR de acuerdo a la figura

Figura

Fuerza de tensión Ft  ma
Fuerza de alargamiento W  mg
Fuerza del resorte FR  kx
dv
Sabemos que la aceleración es a  reemplazando dichos valores en
dt
Ft  ma  W  FR nos queda

dv
Ft  m  mg  kx
dt
dv dv dx dv
En la aceleración hacemos a    v reemplazando en Ft tenemos
dt dx dt dx

dv
Ft  m v  mg  kx es decir mvdv  mgdx  kxdx integrando
dx
v2 x2
m  vdv  mg  dx  k  xdx o sea m  mgx  k  C multiplicando por 2 nos queda
2 2

mv 2  2mgx  kx 2  C

Para cuando la masa m está en la posición x  0 v  v0 tenemos


mv0  2mg  0   k  0   C  C  mv0 2 sustituyendo el valor de C
2

mv 2  2mgx  kx 2  mv02 Fórmula que nos da la velocidad en función del


alargamiento.

Ejercicios
53

1) La población de cierta ciudad, aumenta proporcionalmente a la misma. Si


en 50 años aumenta de 50.000 a 100.000 personas. Cuál será su
población al cabo 80 años.

2) En cultivo de levadura, la cantidad de fermento activo aumenta con una


velocidad proporcional a la cantidad presente. Si esa cantidad se triplica en
una hora, que cantidad habrá al cabo de 4½ horas?

También podría gustarte