Tema 4 Grado en Ingeniería Mecánica

ECUACIONES DIFERENCIALES DE PRIMER ORDEN

CONOCIMIENTOS PREVIOS

Para poder seguir adecuadamente este tema, se requiere que el alumno repase y ponga al día sus
conocimientos en los siguientes contenidos:
 Cálculo diferencial de funciones de una variable.
 Cálculo integral de funciones de una variable.
 Derivación de la función compuesta e implícita.
 Representación de curvas planas.
 Dibujo de curvas con Matlab.

OBJETIVOS ESPECÍFICOS

Los objetivos específicos de este tema son:

1. Poder definir los conceptos básicos del tema.

2. Dada una función de una variable, saber comprobar que es solución de una ecuación
diferencial ordinaria dada.
3. Saber comprobar si una familia uniparamétrica de curvas es solución general de una e.d.o.
de primer orden dada. Saber obtener una solución particular cumpliendo una condición
inicial. Conocida su solución general, poder encontrar la ecuación correspondiente.
4. Saber si un problema de valor inicial tiene solución y ésta es única.
5. Saber obtener las trayectorias oblicuas a una familia dada.
6. Poder reconocer las ecuaciones de los tipos separables, homogéneas, exactas y lineales,
manejando las técnicas de búsqueda de solución analítica de esas ecuaciones.
7. Realizar cambios de variable en e.d.o. de primer orden.
8. Saber bosquejar un campo de direcciones de una e.d.o. de primer orden y utilizarlo, junto
con las isoclinas para aproximar gráficamente una solución.
9. Poder analizar las características relativas a monotonía y concavidad de las curvas de una
familia uniparamétrica, conociendo únicamente la ecuación diferencial que verifican.
10. Entender los métodos de Euler y Euler mejorado para la aproximación de soluciones de
problemas de valor inicial y saber utilizarlos con la ayuda de calculadora u ordenador.
11. Saber utilizar los polinomios de Taylor para aproximar soluciones de problemas de valor
inicial.
12. Poder escribir la e.d.o. de primer orden que modela un proceso y resolverla analítica o
numéricamente.
2 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

ECUACIONES DIFERENCIALES. DEFINICIONES BÁSICAS

1 Definiciones

Llamaremos ecuación diferencial (ED) a una ecuación diferencial ordinaria (EDO) o a una ecuación
en derivadas parciales (EDP), de acuerdo con las definiciones siguientes:

Ecuación diferencial ordinaria (EDO). Es toda ecuación que liga la variable independiente
 
(por ejemplo, x ) con alguna derivada de la función incógnita y x , y ' x ,..., y (n x  :

 
F x , y(x ), y (x ), , y (n (x )  0

Ecuación en derivadas parciales (EDP).- Es aquella en la que la función incógnita depende

de más de una variable independiente conteniendo la ecuación derivadas respecto de
diferentes variables.

Orden de la ecuación diferencial.- Es el mayor orden de la derivada que aparece en la

ecuación diferencial

Ejemplo: Las siguientes ecuaciones son ecuaciones diferenciales ordinarias

 y ' 2x  0 , es de orden 1

 y ' 2x  4y '' , es de orden 2

 xy ''' 2xy ''  x 3  1 , es de orden 3

Ejemplo: Las siguientes son ecuaciones diferenciales parciales

u u
 3  u , es de orden 1
x y

 uxx''  uyy''  0 , es de orden 2

 ux''  uyy''  u , es de orden 2

En esta asignatura, trabajaremos únicamente con EDOs.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
3

2 Soluciones

Resolver una ecuación diferencial es encontrar la función que la verifica.

Solución de la EDO.- Es una función y   x  , definida en un intervalo I  a, b  que

tiene al menos n derivadas continuas en I y que satisface la EDO en I .
 
F x , (x ),  (x ), , (n (x )  0

Distinguimos tres tipos de soluciones

 Solución general de la EDO.- Es una familia de funciones dependientes de n

parámetros, de la cual pueda extraerse la EDO por derivación y eliminación de los
parámetros.
 Solución particular de la EDO.- Es cada una de las funciones que se obtiene de la
solución general al dar valores a los parámetros.
 Solución singular de la EDO.- Es una solución que no puede extraerse de la solución
general.

Ejemplo: Consideremos la familia de rectas y  Cx  2C 2 que es la solución general de la

ecuación diferencial y  xy ' 2 y ' . Se puede comprobar que la parábola x 2  8y  0 es una
2

solución singular. Se puede observar que la parábola es tangente en cada uno de sus puntos a una
curva de la familia de rectas, cuando ocurre esto, se dice que la curva es la envolvente de la
familia de curvas. En este caso, x 2  8y  0 es la envolvente a la familia de rectas y  Cx  2C 2

Cualquier curva de la familia sería una solución particular. Así, por ejemplo, si queremos la solución
de la ecuación diferencial que pasa por el punto (2,4), bastará considerar el valor de C=1
obteniéndose la solución particular y  x  2 .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

4 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

3 Problema de valor inicial

Condiciones iniciales para la EDO.- Es un conjunto de n datos que acompañan a la

ecuación, formado por los valores que toma la función incógnita y sus derivadas en un
punto x  xo : y(x 0 )  y 0 , y (x 0 )  y1, y (x 0 )  y2 , , y (n 1(x 0 )  yn 1

Problema de valor inicial de orden n (p.v.i.).- Es una e.d.o de orden n acompañada de

n condiciones iniciales:





 
F x , (x ),  (x ), , (n (x )  0

 y(x )  y 0 , y (x 0 )  y1, y (x 0 )  y2 , , y (n 1 (x 0 )  yn 1

 0

4 Linealidad

EDO lineal.- Es toda ecuación lineal en las variables y, y ',..., y (n , tomando la forma:
an (x )y (n (x )    a1(x )y (x )  a 0 (x )y(x )  b(x )

En una EDO lineal la variable dependiente y sus derivadas sólo aparecen elevadas a grado uno, y
además, cada coeficiente sólo depende de la variable independiente.

5 Modelización

Modelar mediante ecuaciones diferenciales un fenómeno o proceso es escribir en términos de

ecuaciones diferenciales las relaciones de cambio sufridas por unas magnitudes frente a otras en el
desarrollo de ese fenómeno. Una vez establecido un modelo matemático, podrá encontrarse una
solución exacta o aproximada de las variables incógnitas que intervienen en el proceso.

ECUACIONES DIFERENCIALES DE PRIMER ORDEN

6 Conceptos básicos

Se escriben a continuación las distintas formas en las que puede presentarse la ecuación de una
EDO de primer orden.

 La forma implícita de una EDO de primer orden es F x , y, y '  0 .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
5

 Si en la ecuación implícita se puede despejar y ' , se obtiene la forma explícita

y '  f x , y  .

dy
 Teniendo en cuenta que y '  , también se puede escribir la ecuación en forma
dx
diferencial M x , y dx  N x , y dy  0 .

La solución general de la EDO de primer orden es una familia uniparamétrica (dependiente de un

parámetro) de funciones
y = (x,C )

donde C es un parámetro real; para cada valor de C se obtendrá una solución particular.

Ejemplo: Consideremos la ecuación diferencial x 2


 1 y ' 3y  2 . Esta expresión de la
ecuación diferencial sería la forma implícita. Puede pasarse a la forma explícita despejando y’

2  3y 2  3y
y' , siendo f x , y   2 .
x 1
2
x 1

dy
Teniendo en cuenta que y '  se tendrá, la forma diferencial de la ecuación diferencial
dx

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

7 Teorema de existencia y unicidad

Un problema de valor inicial de primer orden es

 
F (x , y, y ) = 0

ó y  = f (x , y )

 

y(x ) = y 0 
y(x ) = y 0
 0  0
y se resuelve (también se dice integra) encontrando la función de la familia y = (x,C ) que
verifica y(x 0 ) = y0 . Geométricamente, la familia solución general es el conjunto de curvas en las
que se cumple que la pendiente en un punto es el valor de la función f en ese punto. De esas
curvas, la que pasa por el punto (x 0 , y 0 ) es la solución del problema de valor inicial.

TEOREMA DE EXISTENCIA Y UNICIDAD.-

Dada la ecuación y  = f (x , y ) , con f definida en D = [a, b ][c, d ] , a < x 0 < b ,
c < y 0 < d , si f y fy son continuas en D , entonces existe una solución única de la
ecuación cumpliendo que y(x 0 ) = y0 ; esta solución única del problema de valor inicial está
definida en un intervalo centrado en x 0 .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

6 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

Ejemplo. ¿Tiene solución el PVI y  = x 2y  y 5 con y 2 = 3 ?

Para esta ecuación se tiene que f x , y  = x 2y  y 5 , que es una función continua en  2 y con
derivada parcial respecto de y, fy' x , y  = x 2  5y 4 , también continua en todo el plano. Por lo
tanto, aplicando el teorema de existencia unidad se puede garantizar que hay una solución y solo
única que pasa por el punto (2,-3).

Ejemplo. ¿Tiene solución el PVI y  = log xy  con y 0 = 2 ?

En este caso la función f x , y  = log(xy ) está definida en el primer cuadrante, siendo continua en
1
dicho cuadrante salvo los ejes coordenados. La función fy' x , y  =
, es continua en todo el plano
y
salvo el eje y=0. Por lo tanto, como el punto (0,2) se encuentra sobre el eje x=0, el teorema no
asegura la existencia de solución única.

Si se tuviera una condición del tipo y x o  = yo con x , y   D  x, y  / x  0, y  0

o o

entonces sí se podría asegurar la existencia de una única solución pasando por dicho punto.

8 Campo de direcciones

La ecuación
y  = f (x , y )

fija, en cada punto del plano donde f está definida, el valor de la pendiente de la curva y = y(x )
solución de la ecuación diferencial.
La terna (x , y, y ) establece la dirección de la recta que pasa por (x , y ) y es tangente a la curva
solución y = y(x ) . Así, la ecuación genera un campo de direcciones; la representación gráfica de
este campo mediante pequeños segmentos situados en una muestra de puntos del plano es lo que
se suele llamar campo de direcciones de la ecuación diferencial.
El campo de direcciones de la ecuación proporciona un “mapa” aproximado de las soluciones, pues
al representar sus tangentes plasma gráficamente cómo van cambiando estas soluciones. Además,
facilita el trazado de una solución particular, al menos de forma aproximada; esta aproximación es

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
7

en ocasiones suficiente, como por ejemplo si únicamente fuera de interés el comportamiento de la

solución en cuanto a monotonía y concavidad.

9 Isoclinas

Para facilitar el trazado del campo de direcciones suelen utilizarse las curvas en las que la
pendiente permanece constante, al ser paralelos sobre sus puntos todos los segmentos que
pintemos. Una isoclina es el lugar geométrico de los puntos del plano donde la pendiente de las
curvas solución es constante, siendo su ecuación f (x , y ) = C . Las pendientes de las soluciones en
los puntos de corte con la isoclina son, por tanto, y   C .

Dando a C un número suficiente de valores próximos podremos dibujar una red de isoclinas que
permita representar el campo de direcciones. Uniendo un número finito de los segmentos del
campo obtenemos una aproximación gráfica, en forma de poligonal, de la solución particular que
se precise. Ver figuras 1 y 2.

Figura 1.- A la izquierda, muestra de isoclinas (hipérbolas) y del campo de direcciones de la ecuación y  = xy ; a la derecha se
superponen algunas soluciones.

Figura 2.- A la izquierda, muestra de isoclinas (parábolas) y del campo de direcciones de la ecuación y  = x 2  y ; a la derecha se
superponen algunas soluciones.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

8 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

10 Métodos elementales de resolución

Veamos a continuación algunos tipos de ecuaciones de primer orden con sus métodos de
resolución.

TIPO 1. ECUACIONES DE VARIABLES SEPARABLES O SEPARADAS

Definición (Ecuación Diferencial Ordinaria de primer orden).- Una EDO de primer orden
es de variables separables si se puede escribir de la forma f1 (x )g1 (y )dx = f2 (x )g 2 (y )dy o

también si se puede escribir de la forma y '  g x  h y 

En el caso de que la ecuación sea de la forma f1 (x )g1 (y )dx = f2 (x )g 2 (y )dy , las funciones que
satisfacen esta ecuación deben cumplir alguna de las siguientes ecuaciones:

f1(x ) g (y )
dx = 2 dy , f2 (x ) = 0 , g1(y) = 0
f2 (x ) g1(y )

RESOLUCIÓN

La solución general se obtendrá por integración de la primera de las ecuaciones anteriores:

f1(x ) g (y)
 f2 (x )
dx =  2 dy
g1(y )

y deberá analizarse si f2 (x ) = 0 y g1 (y ) = 0 proporcionan soluciones singulares.

CAMBIO DE VARIABLE

En las EDO de primer orden es frecuente utilizar un cambio de variable para reducir la dificultad de
la ecuación. Por ejemplo, el cambio u = ax  by  c , reduce la ecuación y  = g (ax  by  c) a
una de variables separables.

Ejemplo: Encontrar la familia de curvas cuya pendiente de la recta tangente en todos sus puntos
2y
está dada por y '  .
x
Es una ecuación diferencial de variables separadas

dy 2y dy dx log y
     log x  k 
dx x 2y x 2

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
9

log y
 log y  log x 2  k   e log x  y  Cx 2 C 0
2
e k

 y  Cx 2 C 

TIPO 2. ECUACIONES HOMOGÉNEAS

Definición.- Una función f (x , y ) es homogénea de grado n en sus argumentos si se

cumple f (tx , ty )  t n f (x , y )

Ejemplo: Por ejemplo f x , y   x 3  xy 2  5y 3 es homogénea de grado 3 ya que


f tx , yt   tx   tx ty   5 ty   t 3 x 3  xy 2  y 3  t 3 f x , y  
3 2 3

Observa que en este caso todos los monomios tienen el mismo grado, que es 3.

Definición (Ecuación diferencial homogénea).- Una ecuación expresada de la forma

y   f (x , y ) es homogénea si f (x , y ) es una función homogénea de grado cero en sus
argumentos.
Si aparece expresada de la forma M (x , y )dx  N (x , y )dy  0 será homogénea si tanto M
como N son homogeneas del mismo orden.

Ejemplo:

x y x y
(a) y   es homogénea. Observamos que f x , y   es homogénea de grado 0.
x y x y
El numerador es homogéneo de grado 1 y el denominador también. Por lo tanto, la
ecuación diferencial es homogénea.

(b) x 2  xy  dx  y 2dy  0 es homogénea ya que M x , y   x 2  xy , N x , y   y 2 son

funciones homogéneas de grado 2 (ambas son suma de monomios de grado 2).
 x3   x3  es homogénea de
(c) y '  sen  2  también es homogénea. f x , y   sen 
 x 2y  y 3 
 x y  y   
3

grado 0.

RESOLUCIÓN

Las ecuaciones diferenciales homogéneas se resuelven reduciéndolas a ecuaciones de variables

y
separables mediante el cambio de variable: z  , o lo que es lo mismo y  zx , con lo cual
x
dy  dz  x  z  dx o y '  z '  x  z .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

10 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

Ejemplo: Encontrar la familia de curvas cuya pendiente de la recta tangente en todos sus puntos
2y
está dada por y '  .
x
Esta ecuación diferencial, que se resolvió como de variables separadas, puede resolverse también
como homogénea. Realizamos el cambio

y
z y '  x z ' z
x
Se tendrá que

dz dx
x z ' z  2z  xz'z    log z  log x  k
z x
Tomando logaritmos

z xC C 0

Deshaciendo el cambio,

y  Cx 2 C 0  y  Cx 2 C 

TIPO 3. ECUACIONES EXACTAS

Definición (EDO exacta).- La EDO, M x , y dx  N x , y dy  0 es una ecuación exacta si

existe una función u x , y  tal que du(x , y ) = M (x , y )dx  N (x , y )dy

Si esto es así, se tendrá que

M (x, y )dx  N (x , y )dy = 0  du(x, y ) = 0

con lo cual la solución general de la ecuación será la familia uniparamétrica u(x , y ) = C

Para conocer si una ecuación es exacta se utilizará el siguiente test de exactitud.

TEOREMA (CRITERIO DE ECUACIÓN EXACTA).-

Supongamos M (x, y ) y N (x, y ) funciones con derivadas parciales continuas. La ecuación
M (x , y ) N (x , y )
M (x, y )dx  N (x , y )dy = 0 es exacta si y sólo si =
y x

RESOLUCIÓN

Si la EDO es exacta, existe una función u(x , y ) tal que

ux (x , y ) = M (x , y ) , uy (x , y ) = N (x , y )

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
11

puesto que la solución general de la ecuación es u(x , y ) = C , debemos hallar u(x , y ) realizando
los siguientes pasos:
1. Integrar M (x, y ) respecto de x

ux (x , y ) = M (x , y )  u (x , y ) =  M (x, y )dx  C (y )

2. Derivar respecto de y la expresión obtenida para u x , y  e igualar a N x , y  para
despejar C ' x , y 

uy (x , y ) = N (x , y )  C (y ) = 

3. Integrar C ' x , y  respecto de y .

Como es obvio, este proceso se puede comenzar igualmente integrando N respecto de y para
luego derivar el resultado respecto de x .

FACTOR INTEGRANTE

En ocasiones es posible convertir ecuaciones que no son exactas en exactas mediante el método
del factor integrante.

Definición (Factor integrante).- La función (x , y ) es factor integrante de

M (x, y )dx  N (x , y )dy = 0 si es exacta la ecuación
(x , y )M (x , y )dx  (x , y )N (x , y )dy = 0

La solución general de esta última es la solución general de la ecuación inicial. Con este nuevo
término surge la pregunta, ¿en qué condiciones puede encontrarse un factor integrante?
Analizamos seguidamente la respuesta.
Sabemos que si (x , y ) es factor integrante de la ecuación M (x, y )dx  N (x , y )dy = 0 , la
ecuación
(x , y )M (x , y )dx  (x , y )N (x , y )dy = 0

es exacta. Por tanto, debe ocurrir que

M N    N M 
 M N     
y x y x  x y 

y entonces

 
y N M   N M
M  N x   M (log )  N (log )  
  x y y x x y

Esta es una ecuación diferencial en derivadas parciales, de la cual sólo vamos a encontrar alguna
solución particular y sólo en algún caso concreto, como son los dos siguientes:
CASO 1.-   (x ) :

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

12 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

El factor integrante depende sólo de x , luego la ecuación anterior queda reducida a

M N

 y x
(log ) 
x N
por lo que, dada una ecuación diferencial, para detectar si estamos en este caso basta ver
que el término de la derecha de la expresión anterior es sólo función de x . En caso
afirmativo, integramos esta expresión para obtener el factor integrante,
M N

y x
 dx
(x )  e N

CASO 2.-   (y ) :

Análogamente si el factor integrante depende sólo de y , luego la ecuación se reduce a la

siguiente
N M

 x y
(log ) 
y M

encontrando por integración la forma del factor (y ) siempre que el segundo término de
la expresión anterior dependa sólo de y . El factor integrante resulta ser
M N

y x
  dy
(y )  e M

Por supuesto estos dos casos no cubren todas las posibles ecuaciones para las que se puede
encontrar un factor integrante. El mismo método que se ha seguido para saber cuándo existe un
factor integrante en función sólo de x o sólo de y , podría aplicarse para factores integrantes en
función de xy , x  y , etc.

Ejemplo: Resolver (4x  3y 3 )dx  3xy 2dy  0 sabiendo que admite un factor integrante
dependiente solo de x, es decir, de la forma    x  .

Comprobamos si es diferencial exacta:

M 

M (x , y)  4x  3y 3   9y 2 

y 


N 
2
N (x , y)  3xy 2   3y 
x 



No es exacta, por lo que sabiendo que admite un factor integrante de la forma  x  se tendrá que
cumplir que la siguiente ecuación es exacta

 x  (4x  3y 3 )dx   x  3xy 2 dy  0

 
M1 x ,y  N 1 x ,y 

es decir,

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
13

M 1 N 1

y x

 x  9y 2   ' x  3xy 2   x  3y 2

 ' x  2
 ' x  3xy 2   x  6y 2  
 x  x

Integrando
2
2  dx
log  x    dx   x   e x  x 2
x
Esta expresión se podría haber obtenido directamente a partir de la expresión deducida en
general:
M N

y x 2
 dx  x dx
(x )  e e  e log x  x 2
2
N

Multiplicando la ecuación diferencial por este factor integrante se tendría que la siguiente ecuación
diferencial es exacta

(4x 2  3x 2y 3 )dx  3x 3y 2dy  0

Y se puede comprobar que su solución es:

x 4  x 3y 3  C

TIPO 4. ECUACIONES LINEALES

Definición.- Una EDO lineal de primer orden es de la forma y   p(x )y = q(x ) donde p(x )
y q(x ) son funciones continuas en la región en la que se integra.
● Si q x   0 la ecuación se llama lineal homogénea.
● Si q x   0 la ecuación se llama completa.

Ejemplo: Son ejemplos de ecuaciones lineales

1 1
 xy ' y  x 2  y ' y  x , observa que p(x )   , q (x )  x .
x x
1
 xy ' y  0  y ' y  0 Ecuación lineal homogénea.
x

MÉTODO 1 DE RESOLUCIÓN

La solución general de la ecuación lineal completa se puede obtener siguiendo un método

específico basado en el siguiente teorema, que se enuncia sin demostración.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

14 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

TEOREMA.-
Si se denota por y p cualquier solución particular de la ecuación completa
y   p(x )y = q(x ) y por y h la solución general de la homogénea asociada y   p(x )y = 0
, entonces la solución general de la ecuación completa es y p  yh .

La consecuencia del teorema anterior es que para resolver una ecuación lineal basta encontrar la
solución general de su homogénea asociada y una solución particular de la completa.
RESOLUCIÓN DE LA HOMOGÉNEA ASOCIADA

No tiene ninguna novedad, pues todas las EDOs lineales homogéneas son de variables separables:

y  = p(x )y  y = Ce g (x ) ó y0

donde

g (x ) =  p(x )dx

La única dificultad puede venir del cálculo de primitivas necesario para la resolución.
SOLUCIÓN PARTICULAR DE LA COMPLETA

Para encontrar una solución particular de la ecuación completa se puede utilizar el Método de
variación de constantes, que consiste en buscar una solución particular convirtiendo en función
indeterminada la constante de la familia uniparamétrica de soluciones de la homogénea asociada:

yh = Ce g (x )  y p = C (x )e g (x )

Para determinar C (x ) se deriva y se impone que se verifique la ecuación diferencial

y p = e g (x ) (C (x )  C (x )p(x ))  y p  p(x )y p = C (x )e g (x ) = q (x )

de aquí se despeja C '(x ) y se obtiene C (x ) por integración.

C (x ) = q(x )eg (x )  C (x )   q(x )e dx

g (x )

Por tanto, la solución general es:

yg = y p  yh  e g (x )  q (x )e g (x )dx  Ce g (x )

Sacando factor común e g (x ) y agrupando la constante de integración con la que se obtendrá al

calcular C (x ) , resulta finalmente,

yg = e   q(x )e
 p (x )dx dx
 p (x )dx

Ejemplo: Resolver la ecuación diferencial y ' 2xy  4x

Solución

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
15

La solución general de la ecuación completa es y x   y p x   yh x  .

Primero se resuelve la ecuación homogénea asociada, que es de variables separadas

dy
y ' 2xy  0   2xy
dx

dy
   2xdx  log y  x 2  log C  y  Ce x
2

y
Buscamos ahora una solución particular de la ecuación completa:
y p (x )  C x e x  
y ' p (x )  e x C ' x   2xC x  
2 2

Sustituyendo en la ecuación inicial

y p ' 2xy p  4x  
e x C ' x   2xC x   2xC x e x  4x 
2 2

e x C ' x   4x  C x    4x e x dx  2e x
2 2 2

y(x )  Cex  2 e x e x  Ce x  2
2 2 2 2

MÉTODO 2: RESOLUCIÓN COMO EXACTA

Escribiendo la ecuación en forma diferencial p(x )y  q(x )dx  dy  0 , se observa que siempre
va a tener un factor integrante de la forma   (x ) . Es decir,
M N

y x

 e
dx p (x )dx
(x )  e N

Multiplicando la ecuación por este factor integrante y resolviendo como exacta se llega a la
expresión

1
(x ) 
y(x )  q(x )(x )dx

También se llega a esta solución teniendo en cuenta que si la EDO es de la forma

y ' p x  y  q x 

Multiplicando por el factor (x )  e 

p (x )dx
se cumple que

e y ' e  p x  y  e  q x   x y  '   x q x 
p(x )dx p( x )dx p (x )dx


 
 e  p ( x )dx y '

Integrando

1
 x  y x     x q x dx  y x    x    x q x dx

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

16 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

Ejemplo: Resolver la ecuación diferencial y ' 2xy  4x

Solución
Sabiendo que admite un factor integrante de forma
px dx
 x   e   e
2 xdx
 ex
2

se tendrá

e x y ' 2xe x y  4xe x

2 2 2

Opción 1: Dado que

e x y ' 2xe x y  4xe x 

2 2 2


  ex2y 
 4xe
x2
 dx
 
'
e y ' 2xe y  e y 
x2 x2 x2




 e x y  2e x  C  y  2  Ce x
2 2 2

Opción 2: La resolvemos directamente como exacta

2xe x2

y  4xe x dx  e x dy  0
2 2

Buscamos u(x,y) de forma que

ux'  2xe x y  4xe x

2 2
(1)
uy'  e x  u x , y   e x y  h(x )
2 2
(2)

Derivando respecto de x la expresión de u obtenida en (2) y aplicando (1)

2xe x y  h ' x   2xe x y  4xe x  h ' x   4xe x

2 2 2 2

 h x   2e x 
2

u x , y   e x y  2e x  C
2 2

O lo que es lo mismo y  2  Ce x
2

MODELADO DE PROCESOS FÍSICOS Y PROBLEMAS GEOMÉTRICOS

11 Ejemplos de aplicación

VACIADO DE TANQUES

Utiliza la herramienta Vaciado de tanques de la página

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
17

https://personales.unican.es/alvareze/CalculoWeb/CalculoII/integral_superficie.htmll

^2
para ver un ejemplo de aplicación de una ecuación diferencial.

LA RAZÓN DE CAMBIO DE UNA VARIABLE ES PROPORCIONAL AL VALOR DE LA VARIABLE EN

CADA INSTANTE

Podemos citar varios procesos que responden a este modelo: el ritmo de cambio de una
población, la velocidad a la que una sustancia A se transforma en otra en una reacción química
simple, los procesos de desintegración radiactiva o el ritmo de cambio de la velocidad de un
móvil debido a fuerzas de rozamiento.
La ecuación diferencial que rige estos procesos es

dy
 ky
dt

La constante k introduce en la ecuación una medida de la rapidez con la que se produce el cambio
en la magnitud estudiada. El valor de k es propio de cada proceso y si no se conoce, se podrá
determinar conociendo dos datos del valor de y , en t  0 y en otro instante t .

 La cantidad inicial, y(0) , acompañará a la ecuación para definir un problema de valor

inicial y se utilizará para encontrar la forma de su solución particular.

 El dato en otro valor del tiempo se impondrá a la solución anterior para determinar el valor
de k y con ello la función y(t ) .

DESINTEGRACIÓN

Dadas las características de los materiales radiactivos, con el paso del tiempo ocurre una
desintegración o decaimiento del material. No quiere decir que el material desaparezca si no que la
configuración de sus átomos cambia y dejan de ser radiactivos.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

18 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

El tiempo que transcurre para que la cantidad del material inicial decaiga a la mitad se llama vida
media del material radiactivo.

CAMBIO DE TEMPERATURA DE UN OBJETO

Experimentalmente se ha comprobado que la temperatura de un cuerpo cambia

proporcionalmente a la diferencia entre su temperatura y la del medio en que se encuentra. Esta
es una formulación básica de la ley de enfriamiento de Newton.
Si para el momento en que se empieza a estudiar el proceso, t  0 , la temperatura del cuerpo es
T (0)  T0 , y la temperatura del medio es Tm , esto se modeliza a través de la ecuación

dT
 k (T  Tm )
dt

La ecuación irá acompañada de la condición inicial T (0)  T0 y se deberá contar con otro dato de
la temperatura para determinar el valor de la constante k al igual que ocurría en la aplicación
anterior.

 Cuando T0  Tm , ocurrirá un enfriamiento y T(t) decrece por lo tanto

dT
 0 , T  Tm  0 por lo tanto, k será negativa.
dt

 Cuando T0  Tm , ocurrirá un calentamiento y T(t) crece por lo tanto

dT
 0 , T  Tm  0 por lo tanto, k será negativa.
dt
TRAYECTORIAS ISÓGONAS

Se llaman trayectorias isógonas a las familias de curvas que se cortan bajo un ángulo dado.
Si todos los cortes entre curvas de dos familias se producen bajo un ángulo de  radianes, quiere
decir que las rectas tangentes en los puntos de corte forman ése ángulo  .

Por lo tanto, si

 y1  tg , es la pendiente de las curvas de la primera familia, e

 y2  tg , es la pendiente de las trayectorias isógonas

se cumple,      , y la relación entre las pendientes de ambas familias es:

tg   tg y   tg
y1  tg   tg       2
1  tg tg 1  y2 tg

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
19

Suponiendo que se conozca la ecuación diferencial de una de éstas familias, la ecuación diferencial
de la otra familia se obtendrá sin más que sustituir la derivada conocida por su expresión en
función de la derivada de la familia buscada.
Por ejemplo, si conocemos la edo y   f (x , y ) , la edo de sus trayectorias isógonas con ángulo 
es:

y   tg
 f (x , y )
1  y tg

TRAYECTORIAS ORGONONALES

En el caso de que las trayectorias sean ortogonales, la relación entre las pendientes es:

1
y1  
y2

Lo cual, aplicado a la EDO y   f (x , y ) , conduce a la EDO de las trayectorias ortogonales,

1
  f (x , y )
y

LÍNEAS DE FLUJO O LÍNEAS DE FUERZA

Recordemos que las líneas de flujo de un campo vectorial, son las trayectorias seguidas por una
partícula cuyo campo de velocidad es el campo vectorial dado. De esta forma los vectores del
campo son tangentes a las líneas de flujo en cada punto.
Llamando,
F(x , y, z )  M (x , y, z)i  N (x, y, z ) j  P(x, y, z )k
dr  dxi  dyj  dzk

la condición para que la curva cuyo vector desplazamiento es dr sea una línea de flujo del campo
F se traduce en el siguiente sistema,

dx dy dz
 
M N P

que establece la proporcionalidad entre las componentes de los vectores F y dr . La resolución de

este sistema conduce a la familia de curvas buscada.

FLUJO DE CALOR EN UNA DIMENSIÓN EN ESTADO ESTACIONARIO

El flujo de calor en una dimensión se gobierna por la ecuación

dT
H  kA
dr

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

20 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

donde H es la cantidad de calor fluyendo a través de un material ( cal / sg ), k es el coeficiente de

conductividad térmica ( cal / (ºC  cm  sg ) ), A es el área perpendicular al flujo del calor ( cm 2 ),
T es la temperatura ( ºC ) y r es la distancia ( cm ).
Ejemplo. Por una tubería de sección transversal circular por la que está circulando vapor, el flujo de
calor es perpendicular al flujo de vapor. Supondremos tanto simetría axial como radial. Si tomamos
una porción de tubería de 1 cm de longitud, el área de la superficie para un radio r es
2r  1  2r . La variable r varía desde el radio interno ( r0 ) hasta el radio externo. Llamaremos
T0 a la temperatura del vapor en r0 , es decir, en la pared interna del tubo. Si la tasa de flujo, H , y
el coeficiente k de conductividad térmica son conocidos, la ecuación anterior modela la
distribución de temperatura hacia el exterior de la pared de esa parte del tubo:

dT H

dr 2rk

Junto con la condición T  T0 en r  r0 .

Ejercicios propuestos

1 c) y1 (x )  2  3(cos 2  cos 2x ) ;
Halla la e.d.o. de primer orden cuya
y2 (x )  log(e 2  e  e x )
solución general es la familia y = x 2  Cx
Solución: 2xydy = (x 2  y 2 )dx
3
Una familia de curvas también puede
venir dada por una expresión implícita. Por
2
Las expresiones ejemplo (x  1)2  (y  2)2 = C 2 es la familia
1) y1 (x ) = C  3 cos 2x de circunferencias centradas en (1, 2) y radio
2) y 2 (x ) = log(C 2  e x ) |C | .
a. Utilizando derivación implícita, prueba que
corresponden a dos familias uniparamétricas de
todas las curvas de la familia anterior
curvas porque cada valor real del parámetro C
1x
proporciona una curva diferente. verifican la relación y (x ) =
a. Comprueba que todas las curvas de la y 2
familia 1) verifican la relación b. b. Encuentra la curva de esa familia que
y '  6 sen 2x  pasa por el punto P (3, 0) .
b. Encuentra la ecuación diferencial de la c. Analiza si es posible encontrar una curva de
familia 2). esa familia pasando por cualquier punto
c. Encuentra la curva de cada familia que pasa P (x 0 , y 0 ) del plano.
por el punto P 1, 2 . Solución: b) (x  1)2  (y  2)2 = 8
d. Lo que has realizado en el apartado
anterior ¿se podría hacer para cualquier
punto P xo , yo  del plano con ambas 4
familias? Para la familia x 3y 3  x 2  C
Solución: b) y   e x y a. Encuentra la relación entre x , y e y  que
verifica cada curva de la familia; esto es

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
21

más sencillo si se utiliza derivación Solución: a) decreciente: x  , 1 ,

implícita.
b. Determina la curva de esa familia que pasa creciente: x  1,  ; cóncava.
por el punto P (1, 4) . b) crecientes y cóncavas si y  1 , crecientes y
c. Analiza si es posible encontrar una curva de convexas si y  1
esa familia pasando por cualquier punto
P (x 0 , y 0 ) del plano.
8
2  3xy 3 Para cada una de las siguientes
Solución: a) y   ecuaciones,
3x 2y 2
1) y '  xy 2) y '  x  y 3) y '  x 2  y
b) x 3y 3  x 2  63 contesta a estas preguntas:
a. Traza a mano una muestra de las isoclinas
5 correspondientes a cada ecuación. Sobre
Comprueba que las familias de curvas cada isoclina traza los segmentos que
definidas por las funciones siguientes son indican la pendiente de esa isoclina.
ortogonales1 : x 2  y 2  ky , x 3  3xy 2  C b. Estudia la monotonía y la concavidad de las
curvas solución y comprueba que
concuerda con el campo de direcciones
6 representado en el apartado anterior.
Dada la ecuación diferencial y  = 2x
a. Busca todas las funciones y = y(x ) que la
Solución: 1) Crecientes en xy  0 , (1º y 3er
verifican. Debes obtener una familia
cuadrantes); decrecientes en xy  0 (2º y
uniparamétrica de curvas (esto es lo que se
4ºcuadrantes); cóncavas en y  0 ; convexas en
llama solución general de la ecuación).
b. Si no fuera tan sencillo encontrar la familia y 0.
de soluciones, una herramienta que facilita 2) Crecientes en y  x  0 ; decrecientes en
el estudio de su comportamiento es el y  x  0 ; cóncavas en y  x  1  0 ;
método de las isoclinas (una isoclina es el convexas en y  x  1  0 .
lugar geométrico de los puntos del plano 3) Crecientes en y  x 2 ; decrecientes en
donde la pendiente de las soluciones es la
y  x 2 ; cóncavas en y  x 2  2x ; convexas en
misma) y el dibujo del campo de
direcciones (que es una muestra de y  x 2  2x .
pequeños segmentos de las rectas
tangentes a las curvas solución). 9
Representa las isoclinas de pendientes Encuentra la solución de cada uno de
1,  0.5, 0, 0.5, 1 y sobre ellas el los siguientes problemas de valor inicial,
comprobando si se verifican o no las
campo de direcciones. condiciones de existencia y unicidad de
c. Representa las soluciones que pasan por los solución:
puntos 0, 0.5 , 0, 0 y 0, 0.5 . a. (1  e x )yy  = e y , y(0) = 0
Solución: a) y  x 2  C b. y log ydx  xdy = 0 , y(1) = 1

c. x 1  y 2 dx  y 1  x 2 dy = 0, y(0)  1
7
Sin encontrar su expresión, analiza 1  ex
crecimiento y concavidad de las curvas solución Solución: a) e y (1  y ) = 1  x  log
2
de cada una de la siguientes ecuaciones:
b) y(x ) = 1 c) y(x ) = 1 ,
a. y '  x  1
1  x 2  1  y2 = 1
y '  y  1
2
b.

1Todas los cortes entre curvas de ambas familias se 10

Halla la solución general de las
producen perpendicularmente, luego las pendientes en
ecuaciones siguientes transformándolas en una
los puntos de corte distan  / 2 radianes; puesto que
de variables separables mediante los cambios
tg (   / 2) =  1 /tg  esto significa que la derivada
de variable que se indican.
de una familia es inversa y opuesta a la de la otra.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

22 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

a. 2x  2y  1 = y (x  y  2) ;
cambio: z  x  y 14
b. y   cos(x  y  5) ; Encuentra la solución general de las
cambio: z  x  y  5 siguientes ecuaciones, utilizando un factor
integrante. En el primer caso, la ecuación
Solución:
admite factor integrante de la forma  = (x )
2x  y
a) x  y  1 = C exp , (C  0) ; y en el segundo de la forma  = (y ) :
3
 x  y  5  a. (1  x 2y )dx  x 2 (y  x )dy = 0
b) x  C   cotg  
 2  b. xdx  ydy  x (xdy  ydx ) = 0
Solución: a) xy 2  2x 2y  2 = Cx
11 b) x 2  2y  1 = C (1  y )2
Encuentra la solución general de las
siguientes ecuaciones:
a. x(2x 2  y 2 )dx  y(x 2  2y 2 )dy = 0 15
Para cada uno de los casos siguientes,
b. (3x  y )dx  2xydy = 0
2 2 encuentra la solución general y la particular si la
ecuación se acompaña de una condición inicial:
 2
2x  x  y dx = x  y dy y ' 2xy  x 3e x , y 0  1
2 2 2 2
c.  a.
 x 2y  xy 2

b. 4x  3y  y (2y  3x )  0
Solución: a) x 4  x 2y 2  y 4 = C
b) x 3  xy 2 = C c) x 2  y 2  x 3y = Cxy c. xy   (2x  3)y = 4x 4
d. ydx  (xy  x  3y )dy = 0
12 e. uv  = v(1  u tg u )  u 2 cos u
Dada la ecuación diferencial
x sin(ax / 2)  aydx  k x dy , se pide: f. r   = 3  r 4e r  r 5e r
a. Determinar a mano el valor de k para que la e 2  e 2
g. (2) = 8Ch 2 Nota : Ch 2 
ecuación diferencial sea exacta. 2
b. Resolver con Matlab la ecuación diferencial Solución
para el valor de K calculado. 1 1 2
a) y g (x ) =  x 2  e x  Ce x ,
2

c. Encontrar con Matlab la curva solución que 4  2 

 a 
pasa por el punto  ,1 1 1 2 9 2
 2  y p (x ) =  x 2  e x  e x ,
4  2  8
Solución: a) k=a b)
b) y 2  3xy  2x 2  C
C 4sen ax / 2 2 cos ax / 2
y   c) yg (x ) = x 3 (2  Ce 2x )
x a 3x a2
e y y 1
13 d) x g (y ) = C 3 ,
Encuentra la solución de las siguientes y y
ecuaciones: e) vg (u) = u(u  C )cos u
f) p (r ) = r 3 (C  e r  e r  re r )
a.  
3x 2  4xy  2x 2  y y ' = 0
 2
b. ydx  x  dy  0 16
 y  Hallar la solución general de las
dy 2  ye xy siguientes ecuaciones diferenciales:
c.  y
dx 2y  xe xy a) y  = ;
x y
y2
Solución: a) x 3  2x 2y  C = 0 2y  x
2 b) y  = ;
2x  y
b) xy  2 log y  C c) 2x  e xy  y 2  C

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
23

 y y y  a. ¿Al cabo de cuánto tiempo la velocidad

c) x sen  y cos dx  x cos   dy = 0 disminuirá hasta 1 cm/sg?
 x x   x 
b. ¿Qué distancia recorrerá el bote hasta
Solución: a) y(x ) = C  e x /y ; detenerse?
y Solución: a) t  49, 4 sg b) s  14, 8 m
b) y  x = C  y  x  ;
3
c) x  sen =C
x
20
17 Se dispone inicialmente de 2 gr de un
2
Se vierte agua con 0.5 kg de sal por elemento radiactivo en proceso de
litro (l) a un tanque a razón de 2 l/min, y la desintegración. Al cabo de 3 sg quedan 1,4 gr.
mezcla homogénea sale al mismo ritmo. Tras 10 ¿Cuánto tardará en descomponerse 0,5 gr más?
min, se para el proceso y se vierte agua limpia a Solución: t  6, 7 sg
razón de 2 l/min, con la nueva mezcla saliendo
a la misma razón. Si inicialmente había 100 litros
de agua pura en el depósito, se pide: 21
Dos recipientes de igual capacidad se
Escribe las ecuaciones que rigen este proceso
llenan de líquidos diferentes. El recipiente A se
para la cantidad (en kg) de sal en el tanque en
mantiene a 10ºC y el B a 80ºC. Una pequeña
función del tiempo (observa que habrá dos
pieza metálica con temperatura inicial de 30ºC
ecuaciones, según que el tiempo sea menor o
se introduce en el recipiente A, donde su
mayor que 10 min.)
temperatura baja a 25ºC durante el primer
Encuentra la cantidad de sal en el tanque a los
minuto; a los 2 minutos (desde el comienzo) se
20 min., resolviendo analíticamente las
saca de A y se introduce en B, donde pasado
ecuaciones anteriores.
otro minuto alcanza los 40ºC. ¿Cuánto tiempo,
Solución: b) S (20)  7.42053 kg ;
contado desde el principio del proceso, tardará
la pieza en llegar a 60ºC?
18 Nota: Experimentalmente se ha comprobado
Un tanque con agua está recibiendo que la temperatura de un cuerpo cambia
tierra arrastrada por una riada a razón de 700 proporcionalmente a la diferencia entre su
dm 3 / sg . El agua arrastra 53 gr / dm 3 de temperatura y la del medio en que se
dT
tierra y el tanque contiene inicialmente 113,2 encuentra:  k (T  Tm ),(k  0)
dt
m 3 de agua limpia. Un desagüe en el tanque
Donde, Tm es la temperatura del medio
permite salir la mezcla a razón de 140 dm 3 / sg
. Solución: La temperatura durante los 2
t
a. Encontrar la ecuación que verifica S (t ) , 3
primeros minutos es T (t )  10  20   y
que es la cantidad de tierra en el tanque,  4 
transcurrido el tiempo t . (Recordar que la
para los siguientes minutos
tasa de cambio de S (t ) será la diferencia t 2
235  32 
entre la tasa de entrada y la tasa de salida). T (t )  80    . La temperatura será
b. Resolver analíticamente la ecuación 4  47 
obtenida en el apartado anterior. de 60ºC para t  4, 8 min .
Solución: b) S (1800)  5.6042  10 gr
7

22
19 Un depósito contiene inicialmente 110
Un bote de masa m , disminuye su gr de sal disueltos en 10 litros de disolución.
movimiento bajo la acción de la fuerza de Durante los 7 primeros minutos se vierte
resistencia del agua, la cual es proporcional a la disolución a razón de 25 cl/min con 12 gr/l de
velocidad del bote. La velocidad inicial del bote sal. Tras esos 7 min se abre un orificio de salida
es 1,5 m/sg, y al cabo de 4 s su velocidad es 1 (el flujo de entrada se mantiene igual y con la
m/sg. misma concentración).
CASO 1: Supón que el flujo de salida es 25
2Tomado
cl/min. Calcula en ese supuesto la cantidad de
de Applied Mathematical Methods for Chemical
Engineers de Loney.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

24 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

sal en el depósito cuando hayan transcurrido 7 x

min. más. Solución: e.d.o.: y    ; 4y 2  x 2  C
4y
CASO 2: Supón que el flujo de salida es 1/6
l/min. Calcula en ese supuesto la cantidad de sal
en el depósito cuando se empiece a desbordar, 26
sabiendo que su capacidad total es 18 litros. Obtener la ecuación de las trayectorias
Solución: a) t  132, 4 min. b) t  211, 7 gr. que cortan a la familia de curvas de ecuación
x 2  y 2  a 2 , bajo un ángulo de 45º .
y x
23 Solución: e.d.o. y   ;
Obtener la ecuación diferencial de las y x
1
  y
circunferencias de radio 1, cuyos centros se
log x 2  y 2  arctg  C
encuentran en la recta y  2x . 2 x
Solución: (y  2x )2 (1  y 2 )  (2y  1)2
27
Obtener la ecuación de las líneas de
24 flujo de los siguientes campos vectoriales:
Obtener la ecuación diferencial de las
parábolas cuyo eje de simetría es paralelo al eje a. F  xi  yj
OY y que son tangentes, simultáneamente, a las b. F  xi  yj
rectas y  0 e y  x . c. F  (1, 2x , 0)
Solución: xy  2  2yy   y  0 Solución: a) xy  C ;
b) y  mx ;
25 y  x 2  C 1
Obtener las trayectorias ortogonales a la c) 
z  C 2
familia de curvas y  Cx 4 . 

Test de autoevaluación
D) Ninguna de las anteriores.
1
La curva que pasa por el punto 0,1 y
cumple que la pendiente de la recta tangente 3
en cualquier punto es el triple que la abscisa La ecuación diferencial de la familia
del punto de contacto es: uniparamétrica de curvas definida por
3 y  C x 2  4 es:
A) y  x2
2 A) (x 2  4)y   xy
3 B) (x 2  4)y   2xy
B) y  x2 1
2 C) (x 2  4)  2xyy 
3 D) Ninguna de las anteriores.
C) y  x2  1
2
D) Ninguna de las anteriores.
4
La ecuación de la curva que pasa por el
2 9x
El cambio de variable x  e t punto 1,1 y tiene pendiente y    , es:
16y
1
transforma la ecuación 3xy   y  en: A) 8y 2  9x 2  25
x
B) 16y 2  9x 2  25
A) y   y  et
C) 16y 2  9x 2  25
B) 3y   y  e t .
D) Ninguna de las anteriores.
C) 3y   y  e t

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
25

5 8
Las trayectorias ortogonales a la familia Las constantes A y B, tales que
x  Cy , son:
2
y p (x )  A sen x  B cos x sea una solución de
A) 2y 2  x 2  K y   y  2 sen x , son :
B) y2  x2  K A) B  0, A  2
C) 2y 2  x 2  K
B) B  2, A  0
D) Ninguna de las anteriores.
A  1, B  1
C)
6 D) Ninguna de las anteriores.
Sea f (x , y ) una función homogénea de
grado 2. Si f (1,1 / 2)  3 , entonces
9
A) f (4, 2)  48 Una solución particular de la ecuación
B) f (4, 2)  16 diferencial y   2xy  4x es de la forma,
f (4, 2)  24 y p  C (x )e x
2
C) A)
D) Ninguna de las anteriores.
y p  C (x )e x
2
B)
y p  Ce x
2
C)
7
Dada la ecuación diferencial D) Ninguna de las anteriores.
ydx  (1  y 2  x )dy  0 , se puede afirmar
que: 10
1 La ecuación de la isoclina y   2 de la
familia de curvas y  tg log(Cx ) es:
A) es un factor integrante.
x2
1 A) x  y2
B) es un factor integrante.
y2 B) 2x  1  y 2
C) No tiene factor integrante. C) 2x  1  y 2
D) Ninguna de las anteriores. D) Ninguna de las anteriores.

1 2 3 4 5 6 7 8 9 10
C B A B C A B C A C

Ejercicios resueltos
CONCEPTOS BÁSICOS

1
Comprobar que una familia de curvas es solución general de una ecuación diferencial y
encontrar una solución particular que cumple distintas condiciones iniciales.
Realizar los ejercicios: 1, 2, 3, 4, 5, 14 del siguiente documento. Pinchando en el botón S podéis
acceder a la solución.

2
Encontrar la ecuación diferencial cuya solución general es la familia de curvas
y  C x  x2 .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

26 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

Solución
Para resolver este ejercicio se debe derivar la familia de curvas y eliminar la constante C del
sistema formado por
y  C x  x 2


y   C  2x


y
Despejando C de la segunda ecuación y sustituyendo en la primera, resulta, y   x .
x

3
Desde la página
https://www.giematic.unican.es/index.php/edos-primer-orden/material-interactivo
realiza el ejercicio 1 del bloque Conceptos básicos, existencia y unicidad.

4
Sea C una constante arbitraria. Determinar si la ecuación xy 2  y 3  C es solución de
 
la ecuación diferencial 2x  3y y ' y  0 .

Solución
Basta derivar implícitamente respecto a x en la expresión de la familia de curvas y ver que
cumple la ecuación. Derivando

 
y 2  2xyy ' 3y 2y '  0  y ' 2xy  3y 2  y 2  0  y ' 2y  3y   y  0

Sí es solución para todo valor de C.

 
5
Encontrar la ecuación diferencial cuaya solución general es y Cx  x 2  1 .

Solución

 
Derivando la ecuación diferencial se tendrá: y ' Cx  x 2  y C  2x   0 .

1  yx 2
Eliminando la constante, C  , se tiene que
yx

1  1  yx 2 
y'  y   2x   0
y  yx 

Operando

1 1  yx 2  2yx 2  y 2x 2  y
y'    0 y'
y  x  x

6
Encontrar todas las soluciones de la ecuación diferencial y '  x y .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
27

Solución

La ecuación y '  x y es de variables separables. Dividiendo por y, en el caso de que y  0 , se

tiene que

y' dy x2
x   xdx  2 y  C
y y 2

Además, como y  0 es solución se tiene que será una solución singular.

7
(a) Se considera la siguiente figura

Indicar, justificando la respuesta, si puede ser una representación del campo de direcciones de
la ecuación diferencial y '  x 2 . Explicar cómo se obtiene el campo de direcciones de una edo
de primer orden.

(b) Sin resolver la ecuación diferencial y '  e x , ¿puede ser la siguiente figura el campo de
2

direcciones de la edo dada? Justifica la respuesta.

Solución a)
No puede ser el campo de direcciones porque las pendientes en cada punto de la gráfica no son
positivas.
Solución b)
No puede ser porque la pendiente de la solución en cada punto es positiva y, según el gráfico
que muestra el campo de direcciones, hay puntos del plano en los que la pendiente es negativa.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

28 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

 
8
Encontrar la solución de la ecuación diferencial x 2y ' y  1 que verifica y 1  1 .

Solución

y' 1
Resolviendo x 2y '  1  y   si y  1, x  0 . Antes de buscar la solución
1y x2
general veamos si y  1, x  0 son soluciones de la ecuación diferencial.

 
Como y  1 es solución de la ecuación diferencial y cumple con la condición dada y 1  1
, se concluye que y  1 es la solución buscada.

Observad que la solución es única aplicando el teorema de existencia y unicidad considerando

1y
y '  f x , y   .
x2

9
Comprobar con Matlab Comprobar que y t   A  Be kt /m  mgt / k es una solución
de la ecuación diferencial my ''  mg  ky ' .

Solución
Se verifica con Matlab que la curva verifica la ecuación diferencial.
syms t y k m g A B
y=A+B*exp(k*t/m)+m*g*t/k;
aux1=m*diff(y,t,2)
aux2=-m*g+k*diff(y,t)
simplify(aux1-aux2)

E_01
Obtén la pendiente de la recta tangente a la curva solución de la EDO
y  cos x 
2

y' en el punto 0, 3 . ¿Existe solución única de la edo que pase por dicho
3x 2  y
punto?
Solución

La curva solución y x  que es solución del problema de valor inicial dado en la pregunta debe
32  cos 0 10
cumplir la ecuación luego y '(0,3)   .
30  3
2
3
Aplicando el teorema de existencia y unicidad de soluciones de ecuaciones diferenciales de
primer orden3, el problema de valor inicial tiene una única solución al ser f y fy' funciones
continuas en un entorno del punto (0,3) siendo

3
Nota: Leer las páginas 5 y 6 de los apuntes

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
29

f x , y  
y 2  cos x 
f 'y x , y  
  
2y 3x 2  y  y 2  cos x  
3x 2  y
3x 
2
2
y

E_02
Representa una muestra del campo de direcciones de la ecuación diferencial del
apartado 2 en el rectángulo [0,3]x[-4, 4] considerando 5x6 puntos. Dibuja en la misma gráfica 6
curvas solución de la ecuación diferencial 2  x  2y  dx  dy  0 dando 6 valores a la
constante de integración entre -10 y 4.
Solución
El código pedido es

%Campo de direcciones de la ecuación diferencial

[X,Y]=meshgrid(linspace(0,3,5),linspace(-4,4,6));
%Se dibujan como vectores en lugar de como segmentos
H=quiver(X,Y,ones(size(X)),-2-X-2*Y);
set(H, 'ShowArrowHead','off');
hold on
%Solución del problema de valor inicial
syms y(x) C1
sol(x,C1)=dsolve(diff(y)==-2-x-2*y)
%Representación de 6 curvas con 6 valores entre 10 y 4 para C1
for k=linspace(-10,4,6)
fplot(sol(x,k),[0 3])
end
hold off

RESOLUCIÓN DE EDOS CON MATLAB

dsolve(ecuación,condiciones)
Devuelve la solución de la ecuación diferencial con las condiciones iniciales o de frontera
especificadas.
syms y(x)
ecuac=diff(y)+3*x==0;condion= y(0) ==4
sol(x)=dsolve(ecuac, y(0) ==4)
fplot(sol(x),[-3,3])

https://personales.unican.es/alvareze/CalculoWeb/CalculoII/apuntes/tema4n.pdf

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

30 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

10
(a) Resolver la EDO y '  1  y 2 y representar algunas curvas de la familia solución.

(b) Resolver la EDO y '  1  y 2 y obtener la solución particular correspondiente a la condición

inicial (0,1).
Solución a)

syms y(x)
ecuacion=diff(y)==1+y^2;
solu=dsolve(ecuacion)

La solución que da Matab es

solu(x) =
tan(C1 + x)
1i
-1i

Para dibujar algunas curvas de la solución general se considera el siguiente código.

syms C1
sol(x,C1)=solu(1)
hold on
for k=1:5
fplot(sol(x,k),[-3,3])
end
hold off

Solución b)
syms y(x)
ecuacion=diff(y)==1+y^2;
condicion=y(0)==1
sol(x)=dsolve(ecuacion,condicion)
fplot(sol(x),[-3,3])

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
31

11
Se considera la ecuación diferencial y   sen y  sen t  . Representar varias curvas
solución.
Solución
syms y(x)
ecuacion=diff(y)==sin(x)*sin(y);
solu=dsolve(ecuacion)
sol(x,C1)=solu(1)
syms C1
hold on
for k=1:5
fplot(sol(x,k),[-3,3])
end
hold off

12
Encontrar una curva cuya gráfica pase por el punto (1,1) y que la pendiente de la recta
tangente en cada punto de la gráfica sea igual al producto de las coordenadas del punto.
Representar dicha curva.
Solución
syms y(x)
ecuacion=diff(y)==x*y;

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

32 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

vinicial=y(1)==1;
sol(x)=dsolve(ecuacion,vinicial)
fplot(sol(x),[-1 1.5])
hold on
plot(1,1,'o')
hold off

13
Escribir el código Matlab para
x
(a) Comprobar si la curva y es solución de la ecuación diferencial
cos x
xy ' y  x 2 sec x tan x .

(b) Representar 10 curvas solución de la ecuación diferencial y '  y  sen x  junto con la
curva solución que pasa por el punto (1,2).
x
(c) Representar el campo de direcciones de la ecuación diferencial y '   en el rectángulo
y
R=[0,3]x[-1,4] y representar la solución de la ecuación diferencial que pasa por el punto
(1,2). ¿Qué representa el campo de direcciones de una ecuación diferencial?
Solución (a)
Basta sustituir en la ecuación diferencial y ver que se cumple. Con el código Matlab siguiente
syms x
y=x/cos(x);
simplify(x*diff(y)-y-x^2*sec(x)*tan(x))

obtenemos que sí es solución al obtener 0 como valor de la última orden Matlab.

Solución (b)
%Resolvemos la ecuación diferencial dada
syms y(x)
sol(x,C1)=dsolve(diff(y,x)==y+sin(x))
hold on
for k=1:10
fplot(sol(x,k),[-1,3])
end
%Representamos la solución particular que pasa por (1,2)

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
33

sol1(x)=dsolve(diff(y,x)==y+sin(x),y(1)==2);
fplot(sol1(x),[-1,3])
plot(1,2,'o')
hold off
exp(-1) exp(x) (cos(1)/2 + sin(1)/2 + 2) -...- cos(x)/2

20

18

16

14

12

10

-1 -0.5 0 0.5 1 1.5 2 2.5 3

x

Solución (c)

[x,y]=meshgrid(0:0.3:3,-1:0.3:4);
[n,m]=size(x);
q=quiver(x,y,ones(n,m),-x./y)
set(q,'ShowArrowHead','off')
hold on
syms y(x)
ecuacion=diff(y,x)==-x/y;condicion=y(1)==2;
sol(x)=dsolve(ecuacion,condicion)
fplot(sol(x),[0,3])
plot(1,2,'o')
hold off

21/2 (5/2 - x 2/2)1/2

2.2

1.8

1.6

1.4

1.2

0.8

0.6

0.4

0.2

0 0.5 1 1.5 2 2.5 3

x

Acercándonos más al punto (1, 2) se aprecia el campo de direcciones

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

34 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

21/2 (5/2 - x 2/2)1/2

2.2

2.1

1.9

1.8

1.7

0 0.5 1 1.5
x

EDOS DE VARIABLES SEPARABLES

14
Encontrar todas las soluciones de las siguientes ecuaciones diferenciales, por separación
de variables:

a) dy  e 3x 2ydx b) (1  y 2 )dx  xdy


c) 2x 2y dy  1  x 2 dx  d) dy  y tg x dx  0

Solución a)
2 
e 2ydy  e 3 x dx  2y  log  e 3x  C 
 3 

Solución b)

dx dy
(1  y 2 )dx  xdy   , x  0, 1  y2  0
x 1  y2
Integrando se obtiene la solución general:

log x  arctg y  log C  y  tg log(Cx )

x  0 , es solución singular.
Solución c)
Separamos las variables

1  x2 1 
2 y dy  dx   2  1 dx , x 0
x 2
x 

1  1
 2 ydy    2  1dx  y 2    x  C
x  x

x  0 , es solución singular.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
35

Solución d)
Procedemos de forma análoga a los casos anteriores:

dy dy sen x
y
  tg x dx , y  0, cos x  0   y
  tg x dx  
cos x
dx

log y  log C  log cos x  log C cos x  y  C cos x

15
Resolver la ecuación diferencial y(xy  1)dx  x (1  xy  x 2y 2 )dy  0 , utilizando el
cambio de variable xy  z .

Solución

z xdz  zdx
y  dy 
x x2
Sustituyendo en la e.d.o. el valor de dy y escribiendo la ecuación en las nuevas variable z, x
se tiene,
 z   2
z x  1 dx  x 1  x z  x 2 z  xdz  zdx  0
 x   
x   x x 2  x2

z xdz  zdx
(z  1)dx  x (1  z  z 2 ) 0
x x2

z (z  1)dx  (1  z  z 2 )(xdz  zdx )  0

(z 2  z  z  z 2  z 3 )dx  x (1  z  z 2 )dz  0


z 3dx  x 1  z  z 2 dz  0 
que es de variables separables:

dx 1  z  z2
 dz , x  0, z  0
x z3
Integrando miembro a miembro, se llega a

dx 1  z  z2 1 1 1 1 1
   dz  log x    z    dz   2   log z  C
x z3 3
z 2
z  2z z

y, deshaciendo el cambio de variable se obtiene la solución general,

2x 2y 2 log y  2xy  1  Cx 2y 2

A esta solución añadiremos las soluciones singulares: x  0, y  0.

Se puede comprobar que la solución general obtenida es correcta, derivando implícitamente

respecto de x la expresión,

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

36 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

2 1
2 log y   2 2 C
xy x y

Derivamos
2y  2(y  xy ) 2(xy 2  x 2yy ) 2y  2(y  xy ) 2(y  xy )
  0    0
y x 2y 2 x 4y 4 y x 2y 2 x 3y 3

dividiendo entre 2 y multiplicando por x 3y 3 , se tiene

x 3y 2y   xy 2  x 2yy   y  xy   0  x (x 2y 2  xy  1)y   y(xy  1)  0

dy
cambiando y   , se obtiene finalmente,
dx
y(xy  1)dx  x (1  xy  x 2y 2 )dy  0

que es la ecuación diferencial del enunciado.

E_01 y'
Calcular todas las soluciones de la ecuación diferencial x y  0.
2
Solución

y'
La ecuación es de variables separables, para resolver la ecuación  x y  0 se tiene
2

dy
 xdx y0
2 y

Por lo tanto, la solución general es

x2
y  C y0
2
Como y=0 es solución de la EDO, se trata de una solución singular. Como se puede ver en el
gráfico esta solución es la envolvente de la familia de curvas general.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
37

EDOS HOMOGÉNEAS

16
¿Es homogénea la función f x , y   3 x 2  7xy ?
 7xy  y 2 
¿Y la función g x , y   cos  2  ?
 x  yx 

Solución
Sí, la primera es de orden 2/3, la segunda de orden O.

x y   2
 
17
Determinar si la siguiente función f x , y  sen     e2y /x  6 es
 y  x2
homogénea y, en caso afirmativo, indicar su grado.
Solución


Es homogenea de orden 0 ya que f tx , ty  t 0 f x , y  

18
Resolver las ecuaciones diferenciales homogéneas siguientes:

dy y x
a) ydx  x  y dy  0 b)   c) (x 2  2y 2 )dx  xydy  0
dx x y
Solución a)

Comprobamos que la ecuación diferencial ydx  x  y dy  0 es homogénea. Despejamos

y  , así:

dy y ty
y'    f (x , y ) , siendo f (tx , ty )   f (x , y )
dx x y t x ty

Como la función f (x, y) es homogénea de grado cero en sus argumentos, la e.d.o. es

homogénea. Se resuelve haciendo el cambio:

y
u  y  ux  dy  udx  xdu
x
Sustituyendo en la ecuación diferencial inicial y separando variables, tenemos:
u xdx  (x  ux )(u dx  x du )  0

(u x  x u  u 2x )dx  (x  u x )xdu   u 2xdx  x 2 (1  u )du

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

38 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

dx 1u dx 1u 1 
   1 du

x

u2
du, x  0, y0 
x
  u2
du    u 2 
u 


1
 log x    log u  log C
u
Deshaciendo el cambio:

x y y
 log x    log  log C  log   log x  log (e x /y )  log K
y x x

y Ke x /y
log  log  y  Ke x /y
x x

Soluciones singulares: x  0, y  0.

Solución b)
La segunda ecuación también es homogénea, ya que:
dy y x ty tx
y'    f (x , y ) , siendo f (tx , ty )    f (x, y )
dx x y tx ty

y
Efectuamos el cambio: u  y  ux  dy  udx  xdu .
x
Sustituimos en la ecuación:

xdu  udx 1 1 1
u  xdu  udx  udx  dx  udu   dx
dx u u x
Integrando:

dx u2
 udu  
x

2
  log x  log C

Deshaciendo el cambio quedará:

y2 y2
  log x  log C   log x  K
2x2 2x2

Solución c)

Es una ecuación homogénea, por tanto, el cambio es:

y  zx  dy  zdx  xdz

Sustituyendo en la edo, se obtiene una ecuación de variables separables

dx zdz
(1  z 2 )dx  zxdz  0   2
x z 1
Y, resolviendo,

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
39

1
log Cx  log(z 2  1)  Cx  z 2  1
2
Finalmente, deshaciendo el cambio se obtiene la solución general y 2  x 2  Cx 4 .

Además son soluciones singulares x  0 , y  x , y  x .

19
Encuentra la solución general y las soluciones singulares de la ecuación diferencial
3x 2

 y 2 y '  2xy .

Solución
Como es una EDO homogénea, consideramos el cambio y  zx  y '  z ' x  z , se tendrá

3  z  z ' x  z   2z  3  z  x z '  2z  z 3  z  
2 2 2

3  z  dz  dx
2

 2

 3  z x z '  z  z 3

z z  1
2
x
si z  0, 1

Teniendo en cuenta que

3  z   A  B  C  A z  1z  1  Bz z  1  Cz z  1  3  z
2 2


z z  1
2
z z 1 z 1 A  3, B  1, C  1


Integrando la expresión

3  z  dz  dx
2

 z z  1
2  x

3 1 1 dx
 z
dz  
z 1
dz  
z 1
dz   x
se obtiene como solución general

log z 2  1  3 log z  log x

z2 1 y 2 
  1  C y  y 3  C y 2  x 2
 
3
 Cx   x 2 
z3   x2

Las soluciones singulares son y  x , y  x , que se obtienen de considerar las soluciones de

la ecuación diferencial no incluidas en la solución general a partir de z  0, z  1, z  1 .

 
20
Resolver la ecuación diferencial 4x 2  3y 2 dx  2xydy  0

Solución
La ecuación diferencial es homogénea

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

40 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

4x 2  3y 2
y'
2xy
y
haciendo el cambio z   y '  z ' x  z se tiene
x

4  3z 2
z 'x  z  
2z
Operando
4  3z 2 4  3z 2  2z 2 4  5z 2
z 'x   z  
2z 2z 2z
2
4  5z 2z dx
z 'x    dz  
2z 4  5z 2 x
Resolviendo
1
    C
1/5 1 C
log 4  5z 2   log x  C  4  5z 2  x e  4  5z 2 
5 x5
Deshaciendo el cambio
 y 2 C
4  5     4x 5  5x 3y 2  C C  
 x  x 5

Nota: También se puede resolver utilizando un factor integrante.

E_01 x 2y  xy 2 x 5
Determinar si la función f x , y   x 2  y 4  x 2y 2   es
x y3
homogénea justificando la respuesta.
Solución
Comprobamos que la función es homogénea de orden 2
t 3x 2y  t 3xy 2 t 5x 5
f tx , ty   t 2x 2  t 4y 4  t 4x 2y 2   
tx t 3y 3
x 2y  xy 2 x5
 t 2x 2  t 2 y 4  x 2y 2  t 2  t2  t 2 f x , y 
x y 3

EDOS EXACTAS. FACTOR INTEGRANTE

21
Resolver las ecuaciones diferenciales siguientes exactas:

a) (sen x sen y  xe y )dy  (e y  cos x cos y )dx b) (x 2  y )dx  (y 2  x )dy  0

Solución a)
Comprobamos si es diferencial exacta:

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
41

M 

M (x , y )  e y  cos x cos y   e y  cos x sen y 

y 


N 
N (x , y )  sen x sen y  xe y   e  cos x sen y 
y

x 


Buscamos la solución:

  
f (x , y )   e y  cos x cos y dx   xe y  sen x cos y  (y ) 
Derivando esta solución e igualando a N (x , y ) , se obtiene (y) :

f 
 sen x sen y  xe y   (y )
y   (y )  C

N (x , y )  sen x sen y  xe y


La solución es:

xe y  sen x cos y  C

Solución b)

Es exacta, ya que las funciones M (x , y )  x 2  y , N (x , y )  y 2  x son continuas en  2 y

además verifican:

 M (x , y )  N (x , y )
  1 .
y x

Para obtener la función f (x , y ) que nos resuelve el problema, planteamos el sistema de

ecuaciones diferenciales:

 f 

 

 1) M (x , y )  x 2  y  
 x 

 


  f 


 2) N (x , y )  y  x 
2

 y 

 

Integramos ahora bien la ecuación (1) respecto a x , o bien la ecuación (2) respecto a y . En
este caso, integramos la ecuación (1) así:

f x3
 x
dx   (x 2  y )dx  f (x , y ) 
3
 y x  g(y )

La función g(y ) es una función de y que tenemos que determinar, actúa como una constante
de integración en el proceso de integrar respecto a x la función f (x , y ) , que depende de dos
x
variables x e y . La expresión de f (x , y ) que hemos obtenido al integrar la ecuación (1) la
sustituimos ahora en la ecuación (2), resultando:

  x 3 
y2  x    y x  g (y )  x  g '(y )
 y  3 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

42 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

Identificando ambos miembros de la igualdad, tenemos:

y3
y 2  g '(y )   g '(y )dy   y 2dy  g(y ) 
3
 C1

Sustituyendo en la función f(x,y) nos queda como solución de la ecuación diferencial

x3 y3
f (x, y )  y x  C
3 3

22
Encuentra la solución general de la ecuación diferencial
e x
 
 y dx  2  x  ye y dy  0 
Encontrar la solución que pase por el punto (0,1).
Solución
La ecuación diferencial es exacta
M N
M x , y   e x  y N x , y   2  x  ye y 1
y x

Existe entonces una función u x , y  de forma que du  Mdx  Ndy . Para calcular la función

u
 e x  y  u x, y   e x  yx  h y 
x
u
 2  x  ye y (1)
y

Derivando u respecto de la variable y, y teniendo en cuenta (1)

u
 2  x  ye y  x  h ' y 
y

h ' y   2  ye y  h y   2y  ye y  e y  C

Por lo tanto, la solución general de la ecuación diferencial es

u x, y   C
e x  yx  2y  ye y  e y  C

Si se quiere que pase por el punto (0,1) , el valor de C es 2.

     
23
Resolver la siguiente edo sen xy  xy cos xy  x 2 cos xy y '  0 .

Solución
Se trata de una ecuación diferencial exacta ya que

sen xy   xy cos xy dx  x 2 cos xy dy  0

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
43


y
sen xy   xy cos xy   x cos xy   x cos xy   x 2ysen xy 

x
     
x 2 cos xy   2x cos xy  x 2ysen xy


Calculamos f de forma que

I  fx'  sen xy   xy cos xy 

II  fy'  x 2 cos xy   f x , y   xsen xy   h x 

Sustituyendo en (I)

sen xy   xy cos xy   h ' x   sen xy   xy cos xy   h ' x   0  h x   C

     
La función es f x , y  xsen xy y la solución es xsen xy  C siendo C un número real.

24
Determina si la siguiente ecuación y ' 2y  3t admite un factor integrante de la
forma eat b
Solución
Se multiplica la ecuación por el posible factor integrante

eat b 2y  3t dt  e

at b
dy  0
 N t ,y 
M t ,y 

M N
  2eat b  aeat b  a  2
y t


La función  t  e 2t b es un factor integrante para cualquier valor de b.

25
Dada la siguiente ecuación diferencial,

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

se pide encontrar la solución general sabiendo que admite un factor integrante de la forma
 x  .

Solución
Se puede comprobar que la ecuación diferencial no es exacta. Como se dice en el enunciado
que admite un factor integrante del tipo  x  se deberá cumplir que la siguiente ecuación
tendrá que ser exacta:



 


 x  y 3  2e x y dx   x  e x  3y 2 dy  0
M (x ,y ) N x ,y 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

44 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

En consecuencia, se tendrá que cumplir N x'  M y'

 
 ' x  e x  3y 2   x e x   x  3y 2  2e x  
 
 ' x  3y 2  e x   x  e x  3y 2  
 ' x 
 x    ' x    1  log  x   x   x   e x
 x 

Como la ecuación es exacta:

e y
x 3
 
 2e 2 x y dx  e 2x  3e x y 2 dy  0 
Se busca u de forma que

 
du  e x y 3  2e 2x y dx  e 2x  3e x y 2 dy  
Se cumplirá

u 'x  e x y 3  2e 2x y 1
u 'y  e 2x
 3e y x 2
2
De la igualdad (1) se tendrá

u  e x y 3  e 2x y  h y 

Aplicando (2)

3e x y 2  e 2x  h ' y   e 2x  3e x y 2  h ' y   0  h y   A

La función u buscada es

u x , y   e x y 3  e 2x y  A

Y la solución de la ecuación diferencial es u x , y   C

e x y 3  e 2x y  C

26 1
Resolver la ecuación diferencial xy ' 2y  x 2  x  1 , y 1  sabiendo que admite
2
un factor integrante de la forma  x  .

Solución
La ecuación


2y  x  x  1 dx 

2
x dy  0
N x ,y 
M (x ,y )

no es exacta. Para obtener el factor integrante que dependa de x, se tendrá que cumplir que

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
45

 
 x  2y  x 2  x  1 dx   x  x dy  0
 
M 1 (x ,y ) N 1 x ,y 

sea exacta. Por lo tanto,

2 x    ' x  x   x    ' x  x   x  

 ' x  1
   log  x   log x   x   x
 x  x

El factor integrante será entonces  x   x . Multiplicando la ecuación diferencial dada por

este factor se tendrá que la siguiente ecuación diferencial es exacta

2yx  x 3

 x 3  x dx  x 2dy  0

Obteniendo la solución general de dicha ecuación diferencial se tendrá la solución pedida

x4 x3 x2 1
x 2y    
4 3 2 12

27
Determinar un factor integrante de la forma   x m y n para la ecuación diferencial
3y  2x y dx  4x  3x y dy  0
3 3 4 2

Solución
Al multiplicar la ecuación diferencial por el factor integrante se tendrá que la siguiente ecuación
será exacta

  
x my n 3y  2x 3y 3 dx  x my n 4x  3x 4y 2 dy  0 
   
3x my n 1  2x m 3y n 3 dx  4x m 1y n  3x m 4y n 2 dy  0
    
 M   N 

Por lo tanto, se tiene que cumplir M y'  N x' , es decir,

3x m n  1 y n  2x m 3 n  3 y n 2  4 m  1 x my n  3 m  4 x m3y n 2

3 n  1 x my n  2 n  3 x m 3y n 2  4 m  1 x my n  3 m  4 x m 3yn 2

3 n  1  4 m  1 3n  4m  7
2 n  3  3 m  4 2n  3m  18

6n  8m  14
6n  9m  54 m  40 3n  7  160  153  n  51

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

46 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

El factor integrante es  x, y   x 40y 51 .

28
De la ecuación diferencial M x , y dx  N x , y  dy  0 se sabe que admite un factor
integrante que depende solo de y . Deducir la expresión de dicho factor integrante. Aplicarlo a

 
la siguiente ecuación diferencial xdx  x 2y  4y dy  0 y calcular la solución cumpliendo

y 4   0 .

Solución

Si la ecuación diferencial M x , y dx  N x , y  dy  0 , admite un factor integrante

dependiente de y se tendrá que la siguiente ecuación es exacta

 y  M x , y  dx   y  N x , y dy  0
 
M1 x ,y  N 1 x ,y 

Por lo tanto, se cumplirá

M 1 x , y  N 1 x , y 

y x
es decir,

 ' y  M x , y    y  M y' x , y    x  N x' x , y 

Se cumplirá que

 ' y  M x, y    y  M y' x, y   N x' x , y 

 

 ' y  M y' x , y   N x' x , y 


 y  M x , y 

Integrando

N x' x , y   M y' x , y 
log  y    dy
M x , y 

N x' x ,y My' x ,y 
 dy
 y   e
M x ,y 

En el caso de la ecuación diferencial dada, será

2 xy

 y   e  e
dy 2 ydy
 ey
2
x

Multiplicando la ecuación diferencial por este factor integrante se tendrá,

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

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
47

Buscamos f de forma que df  e y x dx  ey x 2y  4y dy  0 . Para ello se tendrá que  

2 2

cumplir

x2
fx'  e y x  f x , y   e y  h y 
2 2

  
fy'  e y x 2y  4y  e y x 2y  4y  e y 2y  x2
 h ' y   h ' y   4ye y  h y   2e y
2 2 2 2 2

x2
La función es, por lo tanto, f x , y   e y  2e y . La solución de la ecuación diferencial será
2 2


ey x 2  4  C 
2

Para hallar la curva solución para y 4  0 , se calcula el valor de C que cumple
e 0 16  4  C  C  20 . La curva solución que pasa por el punto (4,0) es


e y x 2  4  20 . 
2

29
Dada la ecuación diferencial x  y  y  x  y '  0 , ¿admite un factor integrante de la
1
forma ? Nota: No se pide resolver la ecuacion diferencial.
x  y2
2

Solución
Veamos que la ecuación diferencial siguiente es exacta
1 1
2 2
x  y dx  y  x dy  0
x y x  y2
2

Como se cumple

M x , y  
1 x 2

 y 2  x  y  2y x 2  2xy  y 2
x  y   M y' x, y   
x 2  y2
x  x 
2 2
2
 y2 2
 y2

N x , y  
1 
 x 2  y 2  y  x  2x  x 2  2xy  y 2
y  x   N x' x, y   
x 2  y2
x  x 
2 2
2
 y2 2
 y2

Es factor integrante ya que se cumple M y x , y  N x x , y

'
  '
 

EDOS LINEALES

30
Resolver las ecuaciones diferenciales siguientes lineales:

a) xy   3y  x 4 b) (1  x 2 )dy  2xydx  cotg xdx

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

48 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

Solución a)

a) xy   3y  x 4

3
Dividimos entre x , suponiendo x  0 : y   y  x 3 (forma canónica)
x
Método del factor integrante:
3 1
 x dx log 1
(x )  e e x3

x3
Multiplicando por este factor resulta la ecuación diferencial

1 
'
1 3
y   4 y  1   3 y   1
x 3
x  x 

Integrando

1
y  dx  y  x x  C   x  Cx 3
3 4

x3
x  0 , no es solución.

Método de la solución particular más la homogénea: y  yh  y p

Homogénea:

3 dy 3
y  y  0   dx  log y  log Cx 3  yh  Cx 3
x y x
Particular:

y p  x 3C (x )  3x 2C (x )
3
yp  x C (x ) 

Sustituyendo en la ecuación completa se tiene,

x 4C (x )  3x 3C (x )  3x 3C (x )  x 4  C (x )  1  C (x )  x

Por lo tanto,

y  yh  yp  Cx 3  x 4

Solución b)

Resolvemos la ecuación lineal (1  x 2 )dy  2xydx  cotg xdx

Método del factor integrante:

Dividimos entre 1  x 2  dx  0 para obtener la forma canónica de la ecuación:

2x cotg x
y  y
1x 2
1  x2

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
49

2x
 1x 2 dx
(x )  e  e log(1x )  1  x 2
2

y
1
(x )
 q(x )(x )dx   1 1x  cotg xdx   1 1x log sen x  C 
2 2

31
Resolver la ecuación diferencial y ' 2xy  4x

Solución

La solución general de la ecuación completa es y x   y p x   yh x  .

Primero se resuelve la ecuación homogénea, que es de variables separadas

dy
y ' 2xy  0   2xy
dx

dy
   2xdx  log y  x 2  log C  y  Ce x
2

Buscamos una solución particular de la ecuación completa:

y p (x )  C x e x  
y ' p (x )  e x C ' x   2xC x  
2 2

Sustituyendo en la ecuación inicial

y p ' 2xy p  4x  
e x C ' x   2xC x   2xC x e x  4x 
2 2

e x C ' x   4x  C x    4x e x dx  2e x
2 2 2

y(x )  Ce x  2 e x e x  Ce x  2
2 2 2 2

32
Obtener la curva solución de la ecuación diferencial

1  x  dx
dy2
 2xy 
x
1x 2

que pasa por el punto (1,2).

Solución
Se trata de una ecuación diferencial lineal
dy 2x x
 y
dx 1  x
 
2 2
1  x2

dy 2x
Se resuelve la ecuación homogénea:  y  0 que es de variables separadas
dx 1  x 2

dy
y

2x
1  x2
dx  log y   log 1  x 2  yH    C
1  x2

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

50 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

C x 
Se resuelve la ecuación particular de la completa considerando: y p  . Sustituyendo en
1  x2
la ecuación completa

1  x  y
2 '
p
 2x y p' 
x
1  x2

C ' x  C x  2x 2xC x  x x
    C ' x  
1x
1  x  1  x  1  x  1  x2
2 2 2 2
2 2 2

 C x  
1
2

log 1  x 2 

La solución de la ecuación diferencial será: y  yH  y p 

C

log 1  x 2  
1  x2 2 1  x2  
Para encontrar la curva que pasa por el punto (1, 2) se sustituye x=1, y=2 y se obtiene el valor
de C asociado a dicha curva

C log 1  1 C log 2 log 2

2    2 C  4
11 2 1  1 2 4 2

La curva será y 
8  log 2

log 1  x 2  

2 1  x2  
2 1  x2 
33
Dada la ecuación diferencial y '  cos x  y . Se pide

a) Encontrar la solución general

b) Justifica si se verifica el teorema de existencia
c) Representar la isoclina de pendiente cero.
Solución a)
Es una ecuación lineal. La solución general del homogeneo es

dy
y ' y  0   dx  log y  x  yH  ce x
y

Para encontrar la solución particular de la completa se considera y p  C x e x . Como se tiene

que cumplir

C ' x e x  C x e x  cos x   C x e x  C ' x e x  cos x  

C x    ex cos x  dx 
ex
2

cos x   sen x  

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
51

cos x  senx
Por lo tanto, y p  C x e x  . La solución general es
2

cos x  senx
y  yH  y p  ce x 
2

Nota: Al ser una ecuación lineal y ' p x  y  q x  también se puede resolver como exacta
px dx
utilizando el factor integrante  x   e   ex

Solución b)

Como tanto f x , y   cos x  y como fy' x , y   1 son funciones continuas en todo el

plano, existe solución única de la ecuación cumpliendo y x o   0 para cualquier xo . Es decir,
para cada punto del plano pasa una solución y sólo una.
Solución c)

La isoclina de pendiente 0 es y '  h x , y   0  cos x  y  0

34
Encuentra las soluciones de la siguiente ecuación diferencial:

y ' 3y  e 3x  0, y 1  5

Solución
La solución general de la ecuación homogénea y ' 3y  0 , es

yH x   Ce 3x C  

Como solución particular se considera y p x   C x e 3x , se cumplirá

C ' x  C x e 3x

 3C x e 3x  e 3x  C x   x y p x   C x e 3x  xe 3x

Por tanto, la solución general es yG x   x  C e 3x , si se impone que pase por el punto (1,
5) es
y 1  1  C e 3  5  C  5e 3  1

 
La solución es y x   x  5e 3  1 e 3x .

35
Resolver la siguiente ecuación diferencial usando el método de variación de las
constantes y '  5y  cos x . 
Solución

La ecuación diferencial es lineal, luego la solución general es yG  yGH  y p .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

52 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

 La solución general del homogéneo es

dy
y '  5y   5 dx  log y  5x  C  y  e 5xeC  y  Ce 5x
y
 Para calcular la solución particular

y p  C x e 5x y p'  C ' x e 5x  5C x e 5x

Sustituyendo en la ecuación diferencial:

C ' x e 5x  5C x e 5x  5C x e 5x  cos x   C ' x   cos x e 5x

Integrando por partes

e5x
C x  
26

5 cos x   sen x  
5 1
Por lo tanto y p   cos x   sen x 
26 26
5 1
La solución de la ecuación diferencial es yG  Ce 5x  cos x   sen x 
26 26
36
Resolver el siguiente problema de valor inicial xy ' 4y  x 6e x , y 1  2

Solución
e trata de una ecuación diferencial de primer orden lineal. y ' p(x )y  q (x )

4 4
y ' y  x 5e x p x    , q x   x 5e x
x x
4
 x dx 1
Consideramos como factor integrante:  x   e  e 4 log x  .
x4
Por lo tanto

1 4
y ' 5 y  xe x
x 4
x
y se tiene

1
 x  y   xe xdx y  xe x  e x  C
x4
Luego la solución general es

y  x 5e x  x 4e x  Cx 4

para y 1  2  C  2 . La solución particular es la curva y  x 5e x  x 4e x  2x 4 .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
53

37 dy
Encontrar la solución de la ecuación diferencial x  y  x 2senx .
dx
Solución
Es una ecuación diferencial lineal

dy 1 1
 y  xsenx  p x    q x   xsenx
dx x x
Consideramos el factor integrante
1
 p x dx   dx 1
 x   e   e x  e  log x 
x
Multiplicando la ecuación diferencial por este factor integrante, se tiene

1 
'
1 1
y ' 2 y  senx   y   senx
x x  x 

Integrando

1
y   cos x  C  y  x cos x  Cx
x

E_01
Dado el PVI,

dy
x  y  x 2 cos 2x  , y x o   1
dx

a) ¿Qué significa que exista una única solución de este problema? ¿Qué debería verificarse
para asegurarlo?
dy
b) Encontrar la solución general de la ecuación diferencial x  y  x 2 cos 2x 
dx
Solución
La ecuación es lineal. Para encontrar la solución general buscamos la solución general de la
ecuación homogénea asociada y la solución particular.
dy
- Solución general de la homogénea: x y  0
dx

dy dy dx
x y  0    log y  log x  C  y  x eC  y  x C C  
dx y x

dy
- Solución particular de la completa x  y  x 2 cos 2x 
dx

Consideramos y p  x C x  y p'  C x   x C ' x  . Se debe cumplir

xy p'  y p  x 2 cos 2x

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

54 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

Por lo tanto, se cumple

 
x C x   x C ' x   xC x   x 2 cos 2x

sen2x
x 2C ' x   x 2 cos 2x  C ' x   cos 2x  C x  
2
x sen 2x 
La solución particular es y p 
2
x sen 2x 
La solución general de la ecuación es: yG  Cx  .
2

APLICACIONES

38
Un isótopo radiactivo plutonio 241 se desintegra de acuerdo a la ecuación diferencial

dQ
 0.0525Q
dt
donde Q t  es la cantidad de material en el instante t medido en miligramos y t en años.

a) Determinar la vida media del plutonio 241, es decir, el tiempo necesario para que la
cantidad de la sustancia se reduzca a la mitad.
b) Si actualmente tenemos 50 mg de plutonio, ¿cuánto quedará en 10 años?
Solución
La ecuación diferencial es de variables separadas
dQ
 0.0525Q
dt
dQ
 0.0525dt  log Q  0.0525 t  C  Q t   e
0.0525 t
C
Q

Considerando que Q 0  Qo el valor de la constante será
0.0525 0
Qo  e C  C  Qo

Por lo tanto, Q t   Qo e
 0.0525 t

Para obtener la vida media del plutonio, suponiendo se deberá calcular el valor de t de forma que
Qo Qo 1
Q t  
 0.0525 t 0.0525 t
 Qo e   e 
2 2 2
log 2
  0, 0525 t   log 2  t  13, 2028 años
0, 0525
Si actualmente tenemos 50 mg. de plutonio, en 10 años se tendrá

Q 10   50 e
0,0525 10
 29, 5778 mg .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
55

Código Matlab
syms C Q(x)
sol(C,t)=dsolve(diff(Q,x)==-0.0525*Q,Q(0)==C)
tiempo=solve(sol(C,x)-C/2)
sol50(t)=sol(50,t)
sol50(10)

39
Sabemos que un material radiactivo se desintegra proporcionalmente a la cantidad
existente en cada momento. En una prueba realizada con 60 mg de este material se observó
que después de 3 horas quedaba un 80% de esta masa. Hallar
a) La función que expresa la cantidad de masa restante cuando ha transcurrido un tiempo
t
b) ¿Qué cantidad queda cuanto t=5 horas?
c) ¿En qué instante la cantidad de material es de ¼ de la cantidad inicial?
Solución
Se denota por Q la cantidad de material en miligramos de la sustancia en el instante t. La
ecuación diferencial que describe la desintegración es
dQ
 kQ , Q 0  60
dt

Resolviendo la ecuación diferencial quedará:

log Q  kt  A  Q t   Ce kt

Para obtener la constante C se considera la condición inicial: 60  Ce k 0

 48  1 4 
Q 3  60e 3k  48  60e 3k  3k  log    k  log  
 60  3  5 
Como

1 4
log  t
La solución entonces es Q t   60e 3  5 
1. .

5 4
log 
2. Cuando t=5 se tendrá Q 5  60e
 5 
3
 41.36 mg

3. El valor de t pedido será el que verifique

3 log 4 
4
1
log  t 1 1 4
 log    log   t  t 
 5 
15  60e 3
 18.6377 horas
 4  3  5  log 4   log 5

40
Un termómetro que está a una temperatura de 18oC se coloca en una sala que se
encuentra a una temperatura de 68oC . Si después de 5 minutos la temperatura del
termómetro es de 28oC , encontrar

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

56 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

a) Tiempo que tardará en tener la temperatura de 37oC

b) La temperatura después de 20 minutos
Solución

Considerando T t  la temperatura del termómetro en el instante t medido en minutos, la

ecuación diferencial que modeliza este problema es la siguiente

dT
 k T  68 T 0  18
dt
Donde k será positiva y el termómetro aumentará la temperatura. Resolviendo la ecuación se
tiene

dT
k  log T  68  kt  C  T  68  e ktC C  0
T  68
T  68  e ktC C     T t   68  e ktC

Teniendo en cuenta que T 0  18  18  68  C  C  50 . Es decir

T t   68  50e kt . Para obtener la constante k que depende del medio utilizaremos que
T 5  28 , por lo tanto

1  40 
28  68  50e kt  50e 5k  40  k  log    0.0446
5  50 

Luego
1  40 
log t
T t   68  50e 5  50 

Para obtener el tiempo que tardará en alcanzar una temperatura de 37º

 31 
 40   40  log  
1
log t
1
log t
1  40   31   50 
 50 
log   t  log    t 
 50 
37  68  50e 5
 50e 5
 31   10.71 min .
5  50   50  1  40 

log  
5  50 

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
57

La temperatura pasados 20 minutos es

1  40 
log 20
T 20  68  50e
 50 
5
 47, 52º

Código Matlab
syms T(t) K
solu(K,t)=dsolve(diff(T,t)==K*(68-T),T(0)==18)
%Cálculo de la constante
K1=solve(solu(K,5)-28)
T(t)=solu(t,K1)
hold on
fplot(T(t),[0,100])
%Temperatura al cabo de 20 minutos
temperatura=T(20)
plot(20,temperatura,'or')
%Calculo del tiemp0
tiempo=solve(T(t)-37);
double(tiempo)
plot(tiempo,37,'*r')
T(20)
hold off

41
Determinar la familia de curvas ortogonal a la familia x 2  3y 2  C .

Solución
La ecuación diferencial que tiene por solución la familia dada es 2x  6yy '  0 , es decir,

x
y' 
3y
Encontramos la ecuación diferencial de la familia ortogonal como solución general de la
3y
ecuación diferencial: y '   . Resolviendo
x

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

58 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

dy 3 3 C
  dx  log y  log x C  y
y x x3

42
Las curvas equipotenciales de un determinado campo electrostáticio se pueden
considerar elipses x 2  2cx  2y 2  0 . Calcular las líneas de fuerza, es decir, las trayectorias
ortogonales a las curvas equipotenciales.
Solución
Hay que calcular la
familia ortogonal a la dada. Derivando implícitamente
2C  2x 2y 2  x 2 2y 2  x 2
2x  2C  4yy '  0 y'  donde se ha sustituido 2C 
4y 4yx x
al despejar en la expresión x 2  2cx  2y 2  0 .

La familia de curvas ortogonales es la solución de la ecuación diferencial

4yx 4yx
y'  2
2y  x
2 2
x  2y 2

que es una ecuación homogenea. Para resolverla se hace el cambio

y
z  y  zx y '  z 'x z
x

4zx 2 4z 1  2z 2 dx
z 'x  z    dz 
x  2z x
2 2 2
1  2z 2 3z  2z 3
x

1
3
2
 
log z  log 3  2z 2  log x  C
3

3x 
2
2
 2y 2
Deshaciendo el cambio, la familia de curvas será: C 
y

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
59

43
Dada la familia F de curvas 4y  a x 2  1  Ce 2y  0 , se pide:

a) Obtener a mano la ecuación diferencial cuya solución general es la familia de curvas

dada.
b) Obtener con Matlab la familia de trayectorias ortogonales a la familia F.
c) Representar con Matlab la curva de la familia F que pasa por el punto a, 0 junto
con la curva ortogonal en dicho punto.

Solución
La ecuación diferencial se obtiene primero derivando implícitamente:
4y ' 2a x  C 2e 2y y '  0

y posteriormente eliminando la constante, teniendo en cuenta que Ce 2y  4y  a x 2  1

 
4y ' 2a x  4y  a x 2  1 2y '  0

Despejando, la ecuación diferencial cuya solución es la familia dada es

ax
y'
4y  ax 2  3

La familia ortogonal buscada es la solución de la ecuación diferencial

4y  ax 2  3
y'
ax

Código Matlab
a=3 %Por ejemplo
syms y(x) C1
sola(C1,x)=dsolve(diff(y,x)==(-4*y-a*x^2+3)/(a*x))
valorC=solve(sola(C1,a))
hold on
fplot(sola(valorC,x))
%Representación de la curva F

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

60 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

valorC=solve(4*0+a*a^2+1+C1*exp(2*0))
syms z
funcion(x,z)= 4*z+a*x^2+1+valorC*exp(2*z)
fplot(funcion(x,z))
plot(a,0,'o')
axis equal
hold off

Nota: Este código es general para cualquier valor de a, se puede simplificar sustituyendo desde
el principio el valor de a por el que correspondiera.

44 y 3 
Calcula las líneas de flujo del campo F x , y    , xy 2  .
 3 

Solución

y3
Llamando M x , y   M x , y   xy 2 , se cumplirá que las líneas de flujo son curvas
3
y(x) que tienen como vector tangente el campo, es decir,

N x , y  xy 2 3x
y'  y' 
M x , y  y /3
3
y

Resolviendo la EDO que se obtiene será

y 2 3x 2
ydy  3xdx    C  y 2  3x 2  C
2 2
Y

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
61

45
Un tanque de 10L de capacidad que está inicialmente lleno de agua pura, recibe una
disolución salada con una concentración de sal de 0,3kg/L a una velocidad de 2L/min. La
solución dentro del tanque se mantiene agitada y fluye hacia el exterior a una velocidad de
2L/min.

a) Determinar el tiempo necesario para que la concentración de sal en el tanque sea de

0,2kg/L.

b) ¿Cuál es la máxima cantidad de sal que se puede acumular en el tanque?

Solución
a) Sea t la variable independiente (tiempo en minutos).
Sea y (t ) la cantidad de sal en el tanque en el tiempo t .

La variación de la cantidad de sal en el tanque en el intervalo t, t  t  es:

y(t )  sal entrante  sal saliente

donde,

sal entrante  0, 3 kg  2 L  t min  0, 6t kg
 L min

 y(t ) kg L
sal saliente  2  t min  0, 2y(t )t kg
 10 L min
Por lo tanto,
 y (t )  0, 6 t  0, 2y (t ) t

Dividiendo entre t y tomando límites cuando t  0 , se obtiene la ecuación diferencial

que modela el proceso, es decir que expresa la variación de la cantidad de sal en el tanque en
cada instante t (tasa de acumulación).
Puesto que en el instante inicial no hay sal, el problema de valor inicial es:
y (t )  0, 6  0, 2y(t )




y(0)  0

Esta ecuación se puede resolver por separación de variables:

dy dy
 0, 6  0, 2y   dt
dt 0, 6  0, 2y
o como lineal:

(t )  e 0,2t  y(t )  e 0,2t  0, 6e 0,2t


dt  Ce 0,2t  3

Imponiendo la condición inicial se obtiene C  3 , por tanto y(t )  3 1  e 0,2t  kg.

kg
La concentración de 0, 5 , implica una cantidad de sal en el tanque de 2 kg, por lo tanto
L
resolviendo 2  3 1  e 0,2t  , obtenemos t  5,5min.

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

62 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

c) Cuando t   , y (t )  3 kg, que es la cantidad máxima de sal que se puede

acumular en el tanque.

46
En una población de 10000 habitantes, cincuenta de ellos tienen una enfermedad
contagiosa. La velocidad a la que se propaga la enfermedad es proporcional al producto de
personas contagiadas por las no contagiadas todavía, con una constante de proporcionalidad
0,2. Si al cabo de tres días, se observa que son 250 las personas contagiadas, calcular el número
de enfermos pasados 12 días (Modelo de Verhulst).
Nota: En este caso se está considerando un modelo simple de dinámica de poblaciones, existen
otros modelos en el que se tienen en cuenta la existencia de circunstancias que limitan el
crecimiento de la población.
Solución
syms P(t) k
sol(t,k)=dsolve(diff(P,t)==k*P*(10000-P),P(0)==50)
valork=solve(sol(3,k)-250)
P1(t)=sol(t,valork)
double(P1(12))

47
Una bola de naftalina pierde masa por evaporación con una rapidez proporcional a su
área superficial. Si la mitad de la masa se pierde en 100 días,
a. ¿cuánto tiempo se necesita para que el radio disminuya a la mitad de su valor inicial?,
b. ¿cuánto tiempo pasará hasta que la bola desaparezca por completo?

4
Indicación: Área de la superficie de la esfera es 4R2 , su volumen es R 3 . Suponer la
3
densidad constante.
Solución
a) Plantear la ecuación diferencial teniendo en cuenta que la masa de la bola de la
4
naftalina, al ser una esfera, es m  V   R 3 siendo  la densidad.
3
b) Resolver la ecuación diferencial.
c) Calcular el radio en el instante t=100 teniendo en cuenta que la mitad de la masa se
pierde en 100 días.
d) A partir del dato anterior calcular el valor de la constante de proporcionalidad de la
ecuación diferencial
e) Calcular el tiempo para que el radio se reduzca a la mitad
f) Calcular el tiempo para que la bola desaparezca

syms r0 r(t) K d
sol(t,K,d,r0)=dsolve(diff(r,t)==-K/d,r(0)==r0)
valorK=solve(sol(100,K,d,r0)==r0/2^(1/3),K)
r(t,d,r0)= sol(t,valorK,d,r0)
tiempoMitad=double(solve(r(t,d,r0)-r0/2,t))

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
63

tiempoDes=double(solve(r(t,d,r0)))

a) 242.366 días b) 484.732 días.

48
Calcular la familia ortogonal de curvas a la siguiente cos y  Cex siendo C  

Solución
En primer lugar obtenemos la edo asociada a la familia de curvas
y ' sen y  Cex
x
Teniendo en cuenta que cos y  Ce la edo es entonces
cos y
y ' sen y   cos y  y ' 
seny
sen y 
La edo de la familia ortogonal es y '  
cos y 
cos y
dy  x  log seny  x  C  sen y   Cex C  
seny

E_01
Calcular con Matlab la solución de la ecuación diferencial 2  x  2y  dx  dy  0
que pasa por el punto (1,1). Encuentra una curva ortogonal a la obtenida como solución que
pase por el mismo punto
Solución
Para resolver el problema con Matlab, basta escribir el siguiente código

clear all
syms y(x)
ecuacion=diff(y)==-2-x-2*y;
sol(x)=dsolve(ecuacion, y(1)==1)
%Solución (9*exp(-2*x)*exp(2))/4 - x/2 - ¾

Para calcular una curva ortogonal a la dada habrá que resolver el problema siguiente,
1
y' y 1  1
2  x  2y

A la hora de corregir esta prueba, se ha considerado la máxima puntuación si se escribía el

código Matlab o se explicaba cómo obtener la curva ortogonal aunque no se diera el restultado
ya que Matlab no da ninguna solución cuando se escribe el siguiente código,
syms y(x)
ecuac=diff(y)==1/(2+x+2*y);
sol(x)=dsolve(ecuac, y(1)==1)

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

64 T4 ECUACIONES DIFERENCIALES DE PRIMER ORDEN

No obstante, se puede obtener la solución de forma analítica. Para ello se debe hacer el cambio
z ' 1
de variable z  1  x  2y , siendo z '  1  2y '  y '  . En consecuencia, la ecuación
2
diferencial a resolver será de variables separables
z ' 1 1 2 2z
  z '  1  z ' 
2 z z z
z  2  
  dz  dx  z  2 log z  2  x  C
dz  dx  1 
2z 
 2  z 
Deshaciendo el cambio

2  x  2y   2 log 4  x  2y   x  C
Si se quiere que la solución pase por el (1,1) , se cumplirá

5  2 log 7  1  C  C  4  2 log 7

En consecuencia, la curva ortogonal será entonces:

2  x  2y   2 log 4  x  2y   x  4  2 log 7
La gráfica muestra en la misma figura tanto la curva solución de la ecuación diferencial dada
como su curva ortogonal en el punto (1,1).

E_02
Una familia de curvas en cualquier punto del plano x , y  tiene como pendiente
f x   4  2x . Determinar la curva de esta familia que pasa por el punto 0, 0 y calcular
también la curva ortogonal que pasa por dicho punto.
Solución

La familia de curvas verifica y '  4  2x , integrando y  4x  x 2  C . De toda la familia de

curvas, la que pasa por el (0,0) es la curva y  4x  x 2 .

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
65

1
La curva ortogonal en dicho punto deberá verificar y '  . Integrando
4  2x

1
y log 4  2x   C
2

1
Si debe pasar por el (0,0), la curva será 0  log 4   C  C   log 2 . La curva
2
1
ortogonal es y  log 4  2x   log 2
2

Puedes ver más ejercicios resueltos sobre ecuaciones diferenciales de primer orden en
la página de Giematic UC
https://www.giematic.unican.es/index.php/edos-primer-orden/material-interactivo

Prof. Elena E. Alvarez Saiz Universidad de Cantabria