Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 67

UNIVERSIDAD NACIONAL DEL CALLAO

FACULTAD DE INGENIERÍA ELÉCTRICA Y ELECTRÓNICA


EP INGENIERÍA ELÉCTRICA

Curso:
Ecuaciones Diferenciales
Proyecto:
Los 500 problemas resueltos
Docente:
Fernández Juan Raymundo
Alumnos:
Alvarado Linares John Dilan
Código:
2223120241

Año:
Perú, Callao 2024
“Año del Bicentenario, de la consolidación de nuestra Independencia, y de la
conmemoración de las heroicas batallas de Junín y Ayacucho”
ANTI-DEMIDOVICH, TOMO 8 (ECUACIONES
DIFERENCIALES)
A.K. Boiarchuk, G.P. Golovach
INTRODUCCIÓN: CONCEPTOS FUNDAMENTALES. CONSTRUCCIÓN
DE ECUACIONES DIFERENCIALES
Escribir las ecuaciones diferenciales correspondientes a las siguientes familias de curvas
Prob.1: x 2−C y 2−2 y =0
Solución: Sea y= y (x)una solución diferenciable con continuidad de la ecuación dada,
donde C es un parámetro que no depende de x. Entonces se debe cumplir la identidad
2 2
F ( x ,C )=x +C y ( x )−2 y ( x )=0 …(1)
Donde x ∈ X ⊂ R, siendo X cierto conjunto. La función F es diferenciable respecto a x.
Hallemos su derivada
∂ F (x , C)
≡2 x +2 y ( x ) y ' (x)C−2 y ' (x)≡0
∂x
De donde
'
y −x ( '
C= y y ≢ 0 ) …(2)
yy '
Sustituyendo (2) en (1), obtenemos la ecuación diferencial buscada

( x 2− y ) y ' −xy=0
Prob.2: Cy−senCx=0
Solución: Procediendo de forma análoga al ejemplo anterior, hallamos
'
C y ( x )−CcosCx ≡0
Si C=0, no está definida la familia de curvas para las cuales se construye la ecuación
diferencial; por tanto, asumimos que C ≠ 0. Del sistema de ecuaciones
'2 2 2 2 2
y =cos Cx C y =sen Cx … (1)
Encontramos
'2
2 1− y
C= 2
, y ≠ 0 … (2)
y
Sustituyendo (2) en (1), resulta

(√ )
'2
'2 2 1− y
1− y =sen 2
x
y
O bien

y ' =cos ( x √ 1− y ' 2


y )
Prob.3: (x−C 1)2−C 2 y 2=1

Solución: Derivando la identidad (x – C 1)2+ C2 y 2 (x) – 1 ≡0dos veces respecto a x ,


obtenemos
2(x – C 1)+C 2 y (x ) y '( x)≡0

1+ ( ( y ( x ) ) + y ( x ) y ' '(x )) C2 ≡0
' 2

Eliminando de las tres identidades las constantes C 1 y C2resulta


2
y y + ( y + yy ' ' ) =0
3 '' '2
Prob.4: y=eC 2

y'
Solución: Derivando respecto a la variable x hallamos y ' ( x )=C e C , de donde C= x
( y ≠ 0).
y
De este modo, la ecuación diferencial correspondiente a la familia dada es
y'
( )x
y
y=e

Prob.5: x−C 1 y 2−C 2 y−C 3=0

Solución: Derivando tres veces la igualdad dada respecto a la variable x , obtenemos


' C1 '
1−2 y y −C2 y =0

2 ( y '2+ y y ' ) C1 +C 2 y ' ' =0


'''
2 ( 3 y ' y ' '+ y ' y ' ' ' ) C 1+C 2 y =0 …(1)

De la última igualdad de (1), hallamos


'''
C2 y
C 1= ( 3 y ' y ' ' + y ' y ' ' ' ≠ 0 ) … (2)
2 (3 y y + y y )
' '' ' ' ''

Sustituyendo (2) en la segunda igualdad de (1), obtenemos


'2 ''' ' '' 2 ''2 ' '' '
y y −3 y y =0 , o bien3 y − y y =0
Pues, y ' ≠ 0, como se deduce de la primera igualdad de (1)
Prob.6: Hallar la ecuación diferencial de todas las circunferencias definidas del plano
Solución: Del curso de geometría analítica se conoce que la ecuación canónica de la
circunferencia con centro en el punto M (C 1 , C 2) y radio R=C3 tiene la forma
2 2
(x−C 1) + ( y −C2 ) −C 3=0 …(1)

Considerando que toda circunferencia se describe mediante dos funciones tres veces
diferenciabas con continuidad y= y (x), a partir de (1) encontramos que

x−C 1+ ( y ( x ) −C2 ) y '(x )≡ 0

1+ y ( x ) + ( y ( x )−C 2 ) y ' ' ( x) ≡0


'2

3 y ' ( x ) y ' ' ( x )+ ( y ( x )−C 2 ) y ' ' '(x )≡ 0

Simplificando y ( x )−C 2 en las dos últimas identidades, finalmente obtenemos

y ' ' ' ( 1+ y ' 2) −3 y ' y =0


' '2

Prob.7: Obtener la ecuación diferencial de todas las circunferencias de radio 1 cuyos


centros se hallan en la recta y=2 x
Solución: De acuerdo con las condiciones del ejercicio, se debe tomar C 3=1 y C 2=2 C1 en la
ecuación (1) del ej.6:
2 2
(x−C 1) + ( y −2C 1 ) −1=0

Esto es una familia uniparamétrica de circunferencias. Derivando la identidad


2 2
(x−C 1) + ( y (x )−2 C 1) −1=0

Respecto a la variable x (considerando que y= y (x) es una función diferenciable con


continuidad), hallamos

(x−C 1) + ( y (x )−2 C 1) y '(x )=0


2

Eliminando C 1 de las dos identidades obtenidas, finalmente resulta

( 2 x− y )2 ( y '2 +1 ) −(2 y ' + 1)2=0


Prob.8: Escribir la ecuación diferencial de la familia de parábolas con ejes paralelos al eje
Oy , y que son tangentes a las rectas y=0 e y=x
Solución: La ecuación de la familia de parábolas con ejes paralelos al eje Oy tiene la forma
2
y=C 1 x +C 2 x +C 3 donde C j son parámetros arbitrarios ( j=1 ,2 , 3). De la condición de
tangencia a la recta y=0 , resulta
'
y =2 C1 x t +C 2=0
2
C 1 x t +C 2 x t +C 3=0

Donde x t es la abscisa del punto de tangencia. De aquí se deduce que


2
C2
C 3= ( C ≠ 0 ) …(1)
4 C1 1

De la condición de tangencia con la recta y=x , obtenemos


'
y ≡ 2C 1 x t + C2=1
xt = yt
2
2 C2
y t =C1 x t +C 2 x t +
4 C1

1 1
lo que nos proporciona C 2= . Sustituyendo el valor deC 2 en (1), hallamos C 3= . De
2 16 C1
2 x 1
este modo, la familia buscada y=C 1 x + + , donde C 1 es un parámetro arbitrario.
2 16 C1
Simplificándolo en las identidades
' 1
y ≡ 2C 1 x +
2
' 2 x 1
y ≡ C1 x + +
2 16 C1

llegamos a la ecuación buscada: x y' 2= y(2 y ' −1)

Prob.9: Hallar la ecuación diferencial de las circunferencias tangentes a las rectas y=0,
x=0 , y que encuentran en el primer y tercer cuadrantes
Solución: Evidentemente, los centros de tales circunferencias pertenecen a la recta y=x ;
por tanto, C 1=C2 (v. ej. 6). Dado que las circunferencias son tangentes a los ejes
coordenados, hallamos C 3 = |C 1|. De esta forma, la familia analizada satisface la ecuación
2 2
(x−C 1) + ( y −C1 ) −C 1=0

A partir de las identidades respecto a x


2 2 2
(x−C 1) + ( y (x )−C 1 ) −C 1 =0

x−C 1+ ( y ( x ) −C1 ) y '(x )=0

se deduce que la ecuación diferencial requerida es

y '2 ( x 2−2 xy ) −2 x y ' + y 2−2 xy =0


Prob.10: Obtener la ecuación diferencial de la siguiente familia de cicloides
x=C (t−sen t)
y=C (1−cos t)
Solución: Derivando las funciones x e y respecto a t , y dividiendo y '(t ) por x ' (t),
obtenemos
dy y ' (t) sen t 1
= = =
dx x ' (t) 1−cos t t
tg
2
t=2 arctg
1
y'
+2 kπ , kϵ Z ( )
Sustituyendo la expresión de t en la igualdad
( t−sen t ) y−x ( 1−cos t ) =0
tras una serie de transformaciones obtenemos la ecuación diferencial buscada
' x+ yy '
y =ctg '2
y (1+ y )
Prob.11: Demostrar que la ecuación diferencial de las curvas de segundo grado tiene la
forma:
' '2 V ' ' ' ' ' IV ' ' '3
9 y y −45 y y y + 40 y
Solución: Consideraremos la ecuación general de la familia de curvas de segundo grado
2 2
a x +2 bxy+ c y +2 dx +2 ey + f =0
donde a ≠ O . Sin pérdida de generalidad se puede suponer que a=1Derivando cinco veces,
obtenemos

x +b ( y + x y ' ) + cy y' +d + e y ' =0 …(1)


1+b ( 2 y ' + x y ' ' )+ c ( y '2 + yy ' ' ) + ey ' '=0 …(2)
b ( 3 y ' ' + x y ' '' ) +c ( 3 y ' y ' '+ yy ' ' ' ) +ey ' ' '=0 …(3)
b ( 4 y ' ' ' + x y IV ) +c ( 3 y ' '2+ 4 y ' y ' ' ' + y y IV ) +e y IV =0 …(4)
b ( 5 y IV + x y V ) +c ( 10 y ' ' y ' ' '+ 5 y ' y IV + y y V ) +e y V =0 …(5)
De las ecuaciones (2), (3) y (4) resulta
IV 2
3 y ' ' y −4 y ' ' '
b=

2 '' ' ' '' 2 IV
3 y ' ' y + 4 y ' y −3 y ' y ' ' y
c=

Donde
2 '' IV '2 ' ' '2 ' ' '2 ' '' ' '4
∆=3 y ' y y −4 y y −6 y y y +9 y
Eliminando la constante e de las ecuaciones (2) y (5), y sustituyendo en la ecuación
obtenida las expresiones (6) de b y c, hallamos finalmente
IV4

9 y ' '2 y V −45 y ' ' y ' ' ' y + 40 y ' ''3=0
que es lo que se pedía demostrar.
y
Prob.12: Demostrar que la familia de curvas diferenciables arctg −lnC √ x + y =0
2 2
2
satisface la ecuación diferencial
Solución: Tomando la diferencial total de la identidad
y (x )
arctg −ln √ x 2 + y 2 ≡lnC
x
Obtendremos

(
d arctg
y (x)
x
−ln √ x 2+ y 2(x ) ≡ 0 )
De donde
xdy− ydx xdx+ ydy 2 2
2 2
− 2 2 =0(x + y ≠ 0)
x +y x +y
O bien
( x + y ) dx−( x− y ) dy=0

Prob.13: Obtener la ecuación diferencial de todas las circunferencias tangentes al eje Ox

Solución. Si en la ecuación de la familia de circunferencias (v. ej. 6) tomamos C 3=|C 2| ,


entonces obtendremos la ecuación de la familia de circunferencias con la propiedad
requerida:
2 2 2
(x−C 1) + ( y −C2 ) −C 2 =0

Hallando C 1 y C 2 a partir de las identidades


x−C 1+ y ' ( y−C 2) =0
'2
1+ y + ( y−C 2 ) y ' '=0

y sustituyendo sus valores en (1), conseguimos la ecuación diferencial buscada

y 2 y ' '2 +2 y ( 1+ y ' 2 ) y ' ' − y '2 (1+ y '2)2=0


Hallar los sistemas de ecuaciones diferenciales cuyas soluciones son las familias de curvas
dadas a continuación:
Prob.14: ax + z=b , y 2 + z 2=b 2
Solución: Sean
x=x , y= y ( x ) , z=z (x)

las ecuaciones paramétricas de la curva, donde y y z son funciones diferenciables con


continuidad. Sustituyéndolas en las ecuaciones dadas obtendremos las siguientes
identidades respecto a x:
2 2 2
ax + z ( x ) ≡ b , y ( x )+ z ( x ) ≡ b …(1)
Derivando estas identidades respecto a x, obtenemos
ax + z ' ( x ) ≡ 0 , y ( x ) y '( x)+ z ( x ) z ' (x)≡0

de dondea=−z ' . Sustituyendo el valor de a en la primera de las identidades (1)


obtenemos que b=z – xz ' . Esta expresión y la segunda identidad de (1) nos conducen a la
igualdad y 2 +2 zz ' x−x 2 z'2=0. De esta forma, el sistema buscado es
' ' 2 ' 2 '2
y y + z z =0 , y +2 xzz −x z =0
Prob.15: x 2+ y 2=z 2−2 bz , y=ax+ b
Solución: Al igual que en el ejemplo anterior, se obtiene
' '
2 x+ 2 y y =2 z z −2 bz '
'
y =a
1 ' '
De donde a= y ' , b= ' (z z −x− y y ). Sustituyendo las expresiones de a y b en las
z
ecuaciones originales, obtenemos:
'
b= y−x y

z ' ( y −x y ' )=z z ' −x− yy '


2 2 2 '
x + y =z −2 z( y−x y )
Las dos últimas ecuaciones diferenciales componen el sistema pedido
Prob.16: Hallar una solución particular de cierta ecuación diferencial cuya solución
general tiene la forma
y=C 1 cosαx+C 2 senαx
Siendo
'
y ( 0 )=1 , y ( 0 )=1
Solución: Derivando y dos veces respecto a x , hallamos fácilmente la ecuación diferencial
correspondiente:
'' 2
y + α y=0
Para encontrar una solución particular de esta ecuación, se deben utilizar las condiciones
iniciales con el objetivo de determinar las constantes C 1 y C 2. Tenemos:
y ( 0 )=¿

y ' ( 0 )=α ¿

De esta manera, C 1=1 y C 1=0. Por consiguiente, la solución particular requerida es


y=cosαx.
Prob.17: Dada la solución general de cierta ecuación diferencial
3
y=C 1 +C 2 lnx+C 3 x

Hallar la solución particular que satisface las condiciones

y ( 1 )=1 , y ' ( 1 )=0 , y ' ' ( 1 )=2


Solución: A partir de las condiciones dadas, resulta
y ( 1 )=¿

y ' ( 1 )= ( Cx +3 C x )¿ =0
2
3
3
x=1

y ' ' ( 1 )=( + 6 C x ) ¿ =2


−C 2
3 x=1
x
−2 −2 2
de donde encontramos C 1= , C 2= ,C 3= . Queda solamente escribir la solución
9 3 9
particular:
−2 2 2 3
y= − lnx+ x
9 3 9
Dejamos a cargo del lector obtener la ecuación diferencial correspondiente.
Prob.18: Supongamos que cierta solución particular satisface el problema de Cauchy
' 2
y =x + y , y ( 0 ) =1
¿Puede esta solución satisfacer también el siguiente problema de Cauchy?
' ' '
y ' =1+2 xy + 2 y , y ( 0 )=1 , y ' ( 0 )=1
Solución: Si la solución particular es dos veces diferenciable con continuidad, entonces del
primer problema de Cauchy obtenemos

y ' ' =1+2 y y ' =1+ 2 y ( x + y 2 )=1+2 xy +2 y 3


así como y '(0)=1. Por consiguiente, la respuesta es afirmativa.
Prob.19: Supongamos que la solución general de cierta ecuación diferencial tiene la forma

y=C 1 x +C 2 e x +C 3 ( x+ x 2 ) ,−∞< x<+ ∞

¿Puede la función y=x +1 ser una solución particular de esta ecuación?


Solución. No puede, pues para ningún valor de las constantes arbitrarias (incluidos ± ∞ ), la
solución particular examinada no se deduce de la solución general.

CAPÍTULO 1: ECUACIONES DIFERENCIALES DE PRIMER ORDEN


1) Ecuaciones de variables separables
Resolver las siguientes ecuaciones diferenciales

Prob.20: √ y 2+ 1 dx=xydy

Solución: Ésta es una ecuación de la forma (1). Dividiendo sus dos miembros por el
producto x √ y 2 +1, obtenemos
dx ydy
− 2 , x≠0
x √ y +1
De donde resulta
dx ydy
∫ =∫ 2 +C
x √ y +1
O bien
ln |x|−√ y 2+1+C , x =0

Prob.21: ( x 2−1 ) y ' +2 x y 2=0 , y ( 0 )=1

Solución: Primero hallamos todas las soluciones de esta ecuación. Tenemos:

( x 2−1 ) dy +2 x y 2 dx=0
de donde, separando las variables x e y , llegamos a la ecuación
dy 2 xdx
2
+ 2 =0
y x −1
Integrando ambos miembros de la ecuación, obtenemos
−1
+ ln| x −1|=C …(1)
2
y
Para obtener todas las soluciones de la ecuación inicial, a la familia (1) de curvas
integrales le agregamos la solución y=0
A partir del conjunto de curvas integrales, separamos la curva que pasa por el punto (0,1).
Tomando x=0 e y=1 en (1), obtenemos C=−1. De esta forma, la función
1
y=
1+ln |x 2−1|
es la solución de la ecuación planteada con la condición y ( 0 )=1

Prob.22: x y' + y = y 2 , y ( 1 )=0.5

Solución: Escribiendo esta ecuación en la forma

xdy + ( y − y 2 ) dx=0
Y separando las variables, obtenemos
dy dx
+ =0
y− y
2
x
Integrando, obtenemos
xy ( 1− y )=C

Nótese que a pesar de que ambos miembros de la ecuación se dividieron por el producto
2
x ( y− y ), las soluciones x=0 , y=0 e y=1no se perdieron. Finalmente, tomando en (2)
1
x=1 , y=0 , 5; encontramos que C= . Así pues, llegamos a que la curva diferenciable
4
4 xy ( 1− y )−1=0

Es la solución del problema planteado

Prob.23: e
−s
( 1+ dsdt )=1
Solución: Escribamos esta ecuación en la forma
ds s
=e −1
dt
Separando las variables s y t, obtenemos
ds
s
=dt
e −1
Integrando la ecuación obtenida, hallamos

| |
s
e −1
ln s
=t+lnC
e

O bien
t
S=−ln ⁡( 1+C e )

Prob.24: y ' =cos ⁡( y −x)


dz dy
Solución: Haciendo z= y −x, obtenemos que − −1 .De esta manera, la ecuación
dx dx
inicial se lleva a la forma
dz
=dx
cosz−1
Después de integrar, hallamos
z
x−ctg =C
2
O bien
y−x
x−ctg =C
2
A estas soluciones hay que agregar las soluciones que perdimos z=2 kπ , kϵ Z
Prob.25: y ' − y =2 x−3
Solución: Primeramente, escribimos esta ecuación en la forma dy =( y +2 x−3)dx .
Cambiando la variable de integración según la fórmula z= y +2 x — 3, obtenemos la
ecuación de variables separables
dz= ( z +2 ) dx
dz
de donde se deduce que (hemos considerado que z ≠ 2), Integrando esta ecuación,
z+2
obtenemos
ln |z +2|=x+ lnC

de donde, finalmente, hallamos y=1−2 x +C e x. Es evidente que la solución z=−2, es decir,


y=1−2 x , pertenece a la familia de curvas integrales obtenida (esta solución puede
incluirse en la familia para C = 0).
Para las ecuaciones siguientes, hallar las soluciones que satisfacen las condiciones
indicadas cuando x →+ ∞:
2 9
Prob.26: x y +cos 2 y =1 , y ( x )= π
4
Solución: Separando las variables, obtenemos
dy dx
= , x ≠ 0 , cosy ≠ 0
2cos y x
2

Después de integrar hallamos


1 1
tgy=C−
2 x
O bien
y=arctg 2C−( 2
x)+2 kπ

Utilizando la condición complementaria dada, encontramos que


9
π= lim y=¿
4 x→+ ∞

( (
¿ lim arctg 2 C−
x→+∞
2
x )
+ 2 kπ =¿ ¿
)
¿ arctg 2 C +2 kπ
π π 1
Como |arctg 2C|< , obtenemos que k =1, arctg 2C= , C= , resultante finalmente
2 4 2
2
y=arctg 1− + 2 π
x ( )
2
3y
Prob.27: 3
dy−2 xdx=0
y −8
Solución: Separamos las variables x e y
2
3y
3
dy=2 xdx , y ≠ 2
y −8
Después integramos los dos miembros de la ecuación y obtenemos
2
3y
∫ y 3−8 dy=2 ∫ xdx +C , ln|y 3 −8|=x2 +C
Suponiendo que C=ln C 1 (C1 >0), tenemos
| y 3 −8|=C1 e x
2

Es evidente que la solución y=2se puede incluir en la familia de curvas integrales si


tomamos C 1=0 . De esta manera, todas las soluciones de la ecuación original se
determinan por la fórmula
| y 3 −8|=C1 e x (C1 ≥ 0)
2

De la fórmula obtenida se deduce que la curva y = 2 es la única que satisface la condición


complementaria planteada en el problema.
Prob.28: Demostrar que toda la curva integral de la ecuación
'
y= 4
y 2+1
t +1 √
tiene dos síntomas horizontales
Solución: Separando las variables x e y en la ecuación diferencial dada y seguidamente
integrando, obtenemos
y x

∫ 3 2 =∫ 3 dt4 … (1)
du
y √ u +1
0 z √ t +10

donde ( x 0 , y 0 ) es un punto arbitrario del plano Oxy . Supongamos que x →+ ∞ en la


+∞
dt
ecuación (1). Entonces en virtud de la convergencia de la integral impropia ∫ 3 4 existe
x √ t +1
0

el límite
y( x)
du
lim ∫ =a
x→+∞ y0 √u2 +1
3

+∞

Puesto que la integral impropia ∫ 3 du es divergente, entonces el límite


x √ t +1
4
0

lim y ( x )= y ( + ∞ ) existe y es finito, siendo además y ( +∞ ) > y , ya que a>0


x→+∞ 0

Supongamos ahora que en la ecuación (1) x →−∞ . Entonces


y( x)
du
lim ∫ 3 2 =b
x→+∞ y √ u +1
0

Donde
−∞
dt
b=∫ 3 4 < 0
+ ∞ √ t +1

Por consiguiente, −∞ < y (−∞ ) < y 0 . De este la fórmula (1) describe la familia de curvas
integrales que poseen dos asíntotas horizontales y (−∞ ) e y (+ ∞)
Prob.29: Investigar el desplazamiento de las curvas integrales de la ecuación
'
y=

ln ⁡(1+ y )
senx
En un entorno del origen de coordenadas. Demostrar que de todo punto de la frontera del
primer cuadrante sale una curva integral contenida en ese cuadrante.
Solución. El dominio D de la función del segundo miembro de la ecuación es
D=¿ k=−∞ ¿ +∞ M k
Donde:
M k = {( x , y ) ∈ R2|2 kπ < x < ( 2 k +1 ) π ,O ≤ y ≤+ ∞ } ⋃
⋃ {( x , y ) ∈ R2|( 2 k +1 ) π < x < ( 2 k +1 ) π ,−1≤ y ≤0 }
Sea 0 ≤ y <+ ∞ , 0< x < π , 0< x 0 < π
Entonces, separando las variables en la ecuación e integrándola, obtenemos
x y

∫ √ sent ∫ ln du
dt
= … (1)
0 x 0 √ ⁡(1+u)
Cuando u → 0 se tiene que ln (1+ u) u. Por consiguiente, la integral en el segundo miembro
de la igualdad (1) converge por el criterio de comparación. Puesto que la función
subintegral es positiva, tenemos
y

∫ lndu >0
0 √ ( 1+u )
Luego
x
dt
∫ √ sent >0
0 x

por tanto, x > x 0. Así pues, para x > x 0 tenemos que y >0 y viceversa, es decir, toda función
que está determinada por la ecuación (1) y que sale del punto ( x 0 ,0 ) se encuentra en el
primer cuadrante. De la desigualdad y '>0 se deduce que y= y ( x )crece al aumentar x. Sin
embargo, de la divergencia de la integral en el segundo miembro de la ecuación (1)
cuando y →+ ∞, y de la convergencia de la integral en el primer miembro de dicha
ecuación cuando x → π−0 , se deduce que y tiende a un límite finito cuando x → π−0 .
Llegamos a una situación análoga si analizamos la igualdad
x y
dt du
∫ √ sent =¿ ∫ ¿
0 y √ ln ( 1+u )
0

0 ≤ y <+∞ , 0< x< π , 0 ≤ y 0 <+∞

Sea −1< y ≤ 0 ,−π < x <0. Entonces, separando las variables en la ecuación inicial e
integrando, resulta
x y
dt du
∫ √−sent =¿ ∫ … (2)¿
0 0 √−ln ( 1+ u )

Dado que y <0 , la integral del segundo miembro de la ecuación (2) es negativa. Esto
significa que x < x 0. Por consiguiente, la curva integral que parte del punto ( x 0 ,0 )( x 0 <0)
tiende hacia abajo y hacia la izquierda. Como se deduce de la ecuación diferencial
considerada, la derivada y ' →+ ∞ cuando x →−π + 0 o cuando y →−1+0 . De esta forma, las
curvas integrales se aproximan asintóticamente a las rectas x=−π o y=−1. Además, las
curvas se aproximan a la recta x=−π por la tangente, mientras que a la recta y=−1 lo
hacen por la normal.

2) Problemas geométricos y físicos que conducen a ecuaciones de


variables separables
Prob.30: Hallar las curvas que poseen la siguiente propiedad, la longitud del segmento del
eje de abscisas seccionado por la tangente y la normal a la curva en un punto cualquier es
igual a a 2.
Solución: Como vemos en la figura 1, el área del triángulo
1
examinado es S= | KN| y , donde y es la ordenada del punto
2
M. Puesto que tgα= y ' (lo que se deduce del significado
2
y
geométrico de la derivada), entonces S= ' , y ’> 0. De esta
2y
2
y
forma tenemos la ecuación diferencial: =a2 y '
2

Considerando que y ≠ 0 y separando las variables, obtenemos


2 dy dx
2
= 2
y a
de donde encontramos
−2 x
= +c
y a2
o bien
2
−2 a
y= 2
Ca +x

2
−y 2
Si y '<0 (fig. 2), entonces S= '
=a . Integrando esta ecuación, obtenemos:
2y
2
2a
y= 2
x−C a
Finalmente, denotando C a2=−C , reunimos ambas soluciones en una:
2
2a
y=
C± x
Prob.31: Hallar las curvas que poseen la siguiente propiedad, la longitud del segmento del
eje de abscisas seccionado por la tangente y la normal a la curva en un punto cualquier es
igual a 2 a.
y
Solución: De la figura 3 vemos que |KL|+|ln|= + yy ' . De esta manera la ecuación
y'
diferencial buscada tiene la forma:
y '
+ y y =2 a
y'
Resolviendo respecto a la derivada, obtenemos

a a2

y '= ± 2 −1
y y
Sea O< y ≤ a. Entonces de la ecuación
diferencial hallamos:
ydy
=dx
a ± √ a2 − y 2
O bien
2 2
d(a − y )
=−2 dx
a ± √ a2 − y 2
Integrando esta igualdad, resulta:

√ a2− y 2−aln ( a ± √ a2− y 2 ) ± x =c


El caso −a ≤ y <0 e y ' <0 se analiza analógicamente.
Prob.32: Hallar las curvas que poseen la siguiente propiedad: si desde un punto
cualquiera de la curva se trazan líneas perpendiculares a los ejes coordenados hasta
cortarlos, entonces el rectángulo obtenido es dividido por la curva en dos regiones cuyas
áreas se encuentran en la relación 1:2.
Solución: De acuerdo con la interpretación geométrica de la integral, tenemos (fig. 4)
x
S2=∫ y ( t ) dt
0

Puesto que S1 + S2=xy entonces teniendo en


cuenta (1) se obtiene:
x
2
S2= xy =∫ y (t ) dt
3 0

Diferenciando esta igualdad respecto a x,


encontramos:
2
( x y'+ y )= y
3
O bien
dy dx
=
y 2x
de donde y=C √ x , o bien x=C y 2. Evidentemente, si intercambiamos el papel de las
variables x e y , el problema tiene la solución y=C x 2.

Prob.33: Hallar las curvas cuyas tangentes en todo punto toman un mismo ángulo con el
radio polar y con el eje polar.
Solución: Utilizando la relación entre los ángulos α y φ (fig. 6), así como la interpretación
geométrica de la derivada, tenemos
'
dx y ( φ ) ρ' senφ+ ρcosφ
= ' = ' =−tgα
dy x ( φ ) ρ cosφ−ρsenφ
1
α = (π −φ)
2
Tras una serie de transformaciones sencillas
obtenemos la ecuación diferencial
ρ' φ
=tg
ρ 2
Integrando la misma, encontramos

| | φ
lnρ=−2 ln cos + ln 2 C
2
O bien
C
ρ=
1+cosφ
Prob.34: Hallar la ecuación en coordenadas polares de una curva tal que la tangente del
ángulo γ comprendido entre el sector de posesión de un punto cualquiera de la curva, y su
tangente en dicho punto es igual al ángulo polar, φ .

Solución: Utilizando la condición tgγ=φ y la relación φ=α +γ (fig.8) podemos escribir


tgφ−tgα
φ=tgγ=tg ( φ−α )=
1+ tgφtgα
'
ρ tgφ+ ρ ρ'
Puesto que tgα= , de la ecuación anterior obtenemos que φ= . Resolviendo esta
'
ρ −ρtgφ ρ
ecuación diferencial, hallamos
C
ρ= ( φ ≠ 0)
φ
Nota: Si en lugar de la figura 8 consideramos la figura 9, obtenemos de un modo
ρ'
análogo la ecuación diferencial φ= , de la cual se deduce que ρ=Cφ.
ρ
Prob.35: La normal MQ a cierta curva corta el eje de abscisas Ox , en el punto Q . Demostrar
que si la abscisa del punto Q es dos veces mayor que la abscisa del punto M, entonces la
curva es una hipérbola equilátera.
Solución: A partir de la figura 10 vemos que
x 0−x
=tgα
y
Utilizando la condición x 0=2 x y el hecho de que y ' =tgα , de (1) resulta la ecuación
diferencial
x= yy '
Separando las variables e integrando,
obtenemos
ydy =xdx
2 2
y −x =C
Del curso de geometría analítica se sabe que la
ecuación canónica de la hipérbola tiene una de
las formas:
2 2
x y
2
− 2 =1
a b
O bien
2 2
x y
2
− 2 =−1
a b
Si a=b , se dice que la hipérbola es equilátera. De este modo, es fácil ver que la segunda
ecuación de (2) determina dos familias de hipérbolas equilátera (para C> 0 y para C< 0).
Nota: Si C=0, a partir de (2) obtenemos un par de rectas y=± x , denominadas
hipérbolas degeneradas.

Prob.36: Demostrar que toda curva cuyas normales pasen por un mismo punto fijo es una
circunferencia.
Solución: Sea M 0 (x 0 y 0) un punto por el cual pasan todas las normales y sea M 1 (x 1 y 1 ) un
punto perteneciente a la curva. En este caso la ecuación de la recta que pasa por esos dos
puntos tiene la forma.
−y
y= ( x −x 0 )+ y 0
y'1 ( x1 )

Como las coordenadas del punto M 1 (x 1 y 1 ) también satisfacen esta ecuación, se tiene
−y
y 1 (x 1)= ( x1− x0 ) + y 0
y'1 ( x1 )

Separando las variables x 1 e y 1 e integrando, obtenemos


( y 1− y 0 ) d ( y 1− y 0 ) + ( x 1−x 0 ) d ( x 1−x 0 ) =0
2 2 2
( y 1− y 0 ) + ( x1 −x 0) =C
De este modo, si C ≠ 0 tenemos la ecuación de una circunferencia de radio C con centro en
el punto M 0.
Prob.37: Consideremos un recipiente de 20 litros de capacidad que contiene inicialmente -
una mezcla de nitrógeno (80%) y oxígeno (20%). En el recipiente entra 0,1 litro de
nitrógeno por segundo (el cual continuamente se mezcla con el contenido) y sale la
misma cantidad de mezcla ¿En cuánto tiempo en el recipiente habrá un 99% de
nitrógeno?
Solución: Sea Q(t ) la cantidad de nitrógeno en el recipiente en un instante t desde que
comenzó el proceso de mezcla. Entonces, 0 , 1 dt litros de mezcla contienen (0 , 1 ⋅Q dt)/20
litros de nitrógeno. De acuerdo con las condiciones de partida, durante un intervalo de
tiempo dt en el recipiente penetran 0 , 1 dt litros de nitrógeno y salen (0 , 1 ⋅Q dt)/20 litros.
Por consiguiente, la cantidad dQ de nitrógeno que queda en el recipiente al cabo del
tiempo dt es igual a 0 , 1 (1−Q/20 ) dt litros. De esta forma, obtenemos la ecuación
diferencial.

dQ=0 , 1 1− ( Q
20 )
dt
O bien
dQ dt
=
20−Q 200
Integrando, hallamos
−0,005t
Q=20−C e
Para determinar la constante C, utilizamos la condición Q ( t ) ¿t =0=16 l. Por tanto, C=4 y la
función.
−0,005 t
Q=20−4 e …(1)
es la solución del problema planteado. Sustituyendo t=T y Q=19 , 8 1 en la expresión (1)
(lo que constituye el 99 % de 20 litros) encontramos que T =200 ln 20 s=599 , 2 s ≈10 min .
Esto quiere decir que, al cabo de 10 min desde el inicio de la mezcla, en el recipiente habrá
un 99 % de nitrógeno.
Prob.38: En un tanque hay 100 l de una disolución que contiene 10kg de sal. En el tanque
continuamente ingresan 5 l de agua por minuto, los cuales se mezclan con la disolución
contenida en el tanque. La mezcla se derrama con la misma velocidad con que ingresa el
agua ¿Qué cantidad de sal quedará en el tanque al cabo de una hora?
Solución: Sea Q(t ) kg la cantidad de sal en el tanque en el instante t después de que la
disolución comenzó a derramarse. Entonces su concentración en la disolución
Q Q
considerada es , y la cantidad de sal que fluye del tanque en el tiempo dt es ⋅5 dt .
100 100
Por consiguiente, tenemos la ecuación diferencial
dQ=−0 , 05 Qdt
Aquí el signo "−¿" significa que la cantidad de sal en el tanque disminuye. Integrando la
ecuación obtenemos que Q=C e−0.005 t . Puesto que parat=0 en el tanque había 10 kg de sal,
entonces C = 10. De esta manera, Q=10 e−0.005 t es la solución del problema planteado.
Tomando en la última igualdad t=60 min obtenemos que, pasada una hora, la cantidad de
sal en el tanque es igual a 10 e−3 kg ≈ 0.5 kg .
Prob.39: Una habitación de 200m3 de volumen contiene aire con un 0,15% de CO 2. Un
ventilador introduce una corriente de aire de 20m3/min que contiene un 0,04% de CO2
¿En cuánto tiempo la cantidad de CO2 en al aire de la habitación disminuirá 3 veces?
Solución: Sea Q(t) la cantidad de CO2 en la habitación en el instante t después de que
comenzó a trabajar el ventilador. Entonces, Q(t )/200 es la concentración de CO2 en la
habitación en el instante t . Por consiguiente, los 20m3 de aire que salen de la habitación en
un minuto contienen 0 , 1Q (t)m3 de CO2. Por tanto, en el transcurso de dt minutos de la
habitación saldrán0 , 1Q (t)m3 de CO2. En ese mismo intervalo de tiempo, el ventilador
entregará (0,04%/100%)20dt m3=0,008dt m3 de CO2. De esta manera, el incremento de la
cantidad de CO2 en el tiempo dt es igual a (0,008-0,1Q(t))dt , de donde resulta la ecuación
diferencial
dQ=( 0,008−0 ,1 Q ) dt
Después de integrar, obtenemos
Q ( t ) =( 0 , 08−C e−0 ,1 t ) m3
Puesto que en el instante t=0 se tiene Q=0,3m3 (es decir el 0,15 % de 200m3), entonces
C=−¿0,22. De este modo, Q=0 ,08+ 0 , 22e−0 , 1t . El instante T en que la cantidad de C02 será
0,1m3 se halla a partir de la igualdad 0 , 1=0 , 08+ 0 ,22 e−0 ,1 T . Resulta que
T =10 ln 11≈ 24 min .
Prob.40: La velocidad de enfriamiento (o calentamiento) de un cuero es proporcional a la
diferencia de temperaturas del cuerpo y del medio ambiente. La temperatura del aire que
rodea el cuerpo se mantiene constante e igual a T =20 ° C . ¿En cuánto tiempo el cuerpo se
enfriará hasta 25 ° C si en 10min su temperatura disminuyó desde 100 ° C hasta 60 ° C ?
Solución: De acuerdo con la ley señalada podemos escribir la ecuación
dT
=k ( T −T 0) …(1)
dt
donde T es la temperatura del cuerpo, T 0 es la temperatura del aire circundante, y k, el
coeficiente de proporcionalidad. En nuestro caso T 0=20 °C . Integrando la ecuación (1),
obtenemos
kt
T =20+C e
A partir de la condición T ( t ) ¿t=0 =100° C hallamos que C = 80. La otra condición
T ( t ) ¿t=0 =60 ° C permite calcular k. Así pues,

k =−0 , 1 ln 2
−0.1 tln2 −0.1 t
T ( t )=20+ 80 e =20+80 ⋅2
Tomando T =25 ° C , encontramos que el tiempo buscado es t 1=40 min .
Prob.41: En un recipiente que contiene 1kg de agua a temperatura 20 ° C se introduce un
cuerpo de aluminio de masa 0,5kg y temperatura 75 ° C . Durante un minuto, el agua se
calienta 2 °C . Teniendo en cuenta que la capacidad calorífica especifica del aluminio es
c al =0 ,2 c H , determinar cuanto tiempo deberá transcurrir para que las temperaturas del
2O

agua y el cuerpo se diferencien en 1 °C . Despreciar el calor que se pierde al calentar el


recipiente.
Solución: Por analogía con el ejemplo anterior,
dTc
=k c (T c −T a)
at
dTa
=k a (T a−T c )
dt
donde T c y T a son las temperaturas del cuerpo y del agua, respectivamente; k c y k a son
coeficientes constantes. Restando miembro a miembro la segunda fórmula de la primera y
denotando R=T c −T a ¿, podemos escribir
dR
=kR
dt
k =k c +k a

Integrando esta ecuación, hallamos que R=C e kt Puesto que en el instante inicial t=0 la
diferencia de temperatura será igual a R=55° , entonces C=55. Por tanto, R=55 ekt .
Utilicemos la ecuación de balance calorífico para determinar el coeficiente k. Según la
fórmula general tenemos que Q=cm(T k −T c ), q), donde c es la capacidad calorífica
específica del cuerpo, m su masa y T k −T c la diferencia de temperaturas. De este modo,
Q1=2 c H 2O

Q2=0 ,2 ⋅0 , 5(75−T )c H 2O

donde Q1 es la cantidad de calor que absorbe el agua, y Q2, la cantidad de calor que
desprende el cuerpo cuando se enfría hasta la temperatura T . Conforme a las condiciones
del problema Q1=Q2. De esta manera, hallamos que T =55 ° C y al cabo de un minuto
R=55° C−22 ° C−33 ° C . Entonces 33 ° C=ek 55 °C , de donde k =ln 0 , 6. Por lo tanto, la ley
de equilibrio de las temperaturas del agua y el cuerpo es R=55⋅(0 , 6)t . A partir de la
1

igualdad 1=55 ⋅(0 , 6), hallamos


ln 55
t 1= ≈ 8 min
ln 5−ln 3
Éste es el tiempo al cabo del cual la temperatura del cuerpo será 1 °C mayor que la del
agua.
Prob.42: Un pedazo de metal a una temperatura a se introduce en un horno cuya
temperatura en el transcurso de una hora crece uniformemente desde a hasta b . La
velocidad de calentamiento del metal en el horno es proporcional a la diferencia T de
temperaturas del horno y el metal siendo k el coeficiente de proporcionalidad. Hallar la
temperatura del cuerpo al cabo de una hora.
Solución: De acuerdo con las condiciones de partida
dTm
=k (T h −T m ) …(1)
dt
dTm
Donde T h y T m son las temperaturas del horno y el metal respectivamente, es la
dt
velocidad de crecimiento de la temperatura del metal. Además, a causa del crecimiento
t
uniforme de la temperatura se tiene T h=a+(b−a) , donde el tiempo t se mide en
60
minutos. Por consiguiente, la ecuación diferencial (1) se puede escribir en la forma
dTm t
=k (n+ ( n−a )−T m) …(2)
dt 60
t
Introduzcamos una nueva variable z=a+(b−a) −T m . Entonces la ecuación (2) se
60
transforma en una ecuación de variables separables
dz
=−dt
kz−(b−a)/60
Integrándola, hallamos
1
k (
ln kz −
b−a
60
1
)
=−t + lnC , k ≠ 0
k
O bien
T m=a+ t− ( 1k ) b−a
60
+C e −tk

b−a
Como T m (t ) ¿t =0=a, entonces C= . De esta manera, finalmente obtenemos que
60 k

T m=a+
b−a
60
1
(
t− (1−e−tk )
k )
Evidentemente, la temperatura del metal al cabo de una hora será
b−a −60 k
T m (60 )=a+ (1−e )
60
Prob.43: Una barca disminuye su velocidad bajo la acción de la resistencia del agua, la cual
es proporcional a la velocidad de la barca. La velocidad inicial de la barca es 1,5m/s y al
cabo de 4s es 1m/s. ¿En cuánto tiempo es la barca disminuye su velocidad hasta 1cm/s?
¿Qué distancia recorrerá la barca hasta detenerse?
Solución: Sea v (t) la velocidad de la barca en el instante t después de comenzar su
movimiento. Entonces dv /dt es su aceleración. De acuerdo con la segunda ley de Newton
dv
𝑚𝑑𝑣𝑑𝑡m =F … (1)
dt
donde F es la fuerza de resistencia del agua. Según las condiciones de partida F=kv ; por
tanto, la ecuación (1) toma la forma
dv k
= v=k 0 v ( k 0=cte . )
dt m
Integrando esta ecuación, se obtiene

𝑣𝑡=𝐶𝑒𝑘0𝑡 v ( t )=C ek t 0

Utilizando la condición complementaria v (0)=1 ,5 , encontramos que C=1 ,5 . De este


modo, la ecuación (2) toma la forma
k0 t
v ( t )=1 , 5 e
donde t se mide en segundos. Puesto que v (4)=lm/ s, a partir de. la igualdad l=1 , 5 e 4 k 0

resulta que k 0=0 , 25 ln(2 /3). Por tanto, la velocidad del movimiento de la barca varía
según la ley

()
t / 4−l
2
v ( t )= m/ s …(3)
3

Sustituyendo en (3) v=1 cm/s=0 , 01 m/s , encontramos el instante t 1:

t 1=4 1+ ( ln 0 , 01
ln 2/3 )≈ 50 s

Puesto que v (t)=ds(t )/dt , donde s(t ) es el recorrido, a partir de (3) obtenemos
4 t / 4−l
s ( t )= (2/3) + s0
ln 2/3

()
−1
−4 2
donde s0 es la constante de integración. Si s(0)=0, entonces s0 = , y la ley de
ln2 /3 3
movimiento de la barca tiene la forma

(( ) )
t
6 2 4
s ( t )= −1
ln 2/3 3

A partir de la ecuación (3) vemos que tlim v ( t)=0 ; por tanto, de la ley de movimiento de la
→+ ∞
barca hallamos
6
s1= lim s( t)= ≈ 15 m
t →+∞ ln 3/2
donde s1 es la distancia recorrida por la barca hasta detenerse.
Prob.44: En el transcurso de30 días se desintegra el 50% de la cantidad inicial de cierta
sustancia radioactiva. ¿Al cabo de cuánto tiempo quedara el 1% de la cantidad inicial?
Solución: Empleemos la ley de desintegración radiactiva: la cantidad q (t) de sustancia
radiactiva que se desintegra en la unidad de tiempo es proporcional a la cantidad Q(t ) de
dicha sustancia en el momento dado, es decir, q ( t )=kQ(t), siendo k el coeficiente de
proporcionalidad. Por consiguiente, en el intervalo de tiempo desde t hasta t+ ∆ t se
desintegra una cantidad de sustancia igual a kQ (t 1 ) ∆ t , donde t 1 ∈(t ,t +∆ t) y Q(t 1) es cierto
valor intermedio de la cantidad de sustancia entre Q(t ) y Q(t +∆ t ). Por otro lado, esta
misma cantidad es igual a Q ( t+ ∆ t ) −Q(t); por tanto, tenemos:
Q ( t+ ∆ t ) −Q ( t )=kQ(t 1)∆ t
o bien
Q ( t +∆ t )−Q ( t )
=kQ(t 1)
∆t
Considerando que la función Q es diferenciable y pasando al límite en la última igualdad
cuando ∆ t →0 , obtenemos la ecuación diferencial
dQ (t)
=kQ(t) …(1)
dt
cuya solución es Q ( t ) =C e kt . Evidentemente, la constante C representa la cantidad inicial
−1
de la sustancia, A partir de la condición 0 , 5 C=C e 30 k hallamos k = ln 2, y de la
30
condición 0 , 01 C=C e−(t / 30)ln 2 obtenemos el tiempo
ln 100
t 1= ⋅30 ≈ 200(días)
ln 2
al cabo del cual quedará solamente el 1 % de la cantidad inicial de sustancia.
La fórmula general para la cantidad de sustancia remanente es
−t /30
Q ( t ) =Q ( 0 ) 2
donde t es el tiempo medido en días.
Prob.45: Los experimentos demuestran que en el transcurso de un año de cada gramo de
radio se desintegra 0,44mg. ¿Al cabo de cuánto tiempo se desintegrará la mitad de la
cantidad inicial de radio?
Solución: Sea Q(t ) la cantidad de radio. Ésta satisface la ecuación (1) del problema
anterior y, por consiguiente, la función Q ( t ) =Q ( 0 ) ekt es la ley de desintegración del radio,
donde t es el tiempo medido en años. Para determinar el coeficiente k , utilicemos la
condición de que al cabo de t=1 año queda Q=999 ,56 mg, donde Q ( 0 )=1 g. De aquí
obtenemos que e k =0,99956 . De esta manera, Q ( t ) =Q ( 0 ) ( 999 ,56 )t . Tomando
Q ( t 1) =0 , 5Q (0), obtenemos

ln 2
t 1= ≈ 1600 años
ln 0,99956
Prob.46: Se investiga un pedazo de roca que contiene 100mg de uranio y 14mg de plomo
(el plomo proviene de la desintegrando de uranio). Se conoce que el tiempo de vida media
del uranio es igual a 4 , 5⋅10 9 años, y que de la desintegrando total de 238g de uranio se
obtienen 206g de plomo. Determinar la edad de la roca, considerando que en el momento
de su formación ella no contenía ploma y despreciando la presencia de productos
radiactivos intermedios entre el uranio y el plomo (pues estos se desintegran mucho más
rápido que el uranio)
Solución: Ante todo, determinemos la cantidad inicial de uranio en el pedazo de roca. Sea y
la cantidad de uranio totalmente desintegrado en la roca. Entonces, tomando en cuenta las
condiciones de partida, podemos escribir la proporción
y 238
=
14 206
de la cual encontramos que y=14 ⋅238/206 ≈ 16 , 2mg. Por consiguiente, la cantidad inicial
de uranio en la roca es 116,2 mg.
Calculemos ahora k . A partir de la fórmula general de la desintegración radiactiva
kt
Q ( t ) =116 , 2 e , donde Q(t ) es la cantidad de uranio que no se ha desintegrado, hallamos
que k =−ln 2/(4 ,5 ⋅109 ). Además, teniendo en cuenta que al cabo del tiempo T desde el
inicio del proceso de desintegración en la roca quedaban 100 mg de uranio, a partir de la
condición 100=116 , 2 ekT obtenemos
−1
T= ln 1,162=¿
k
9
4 , 5 ⋅10 6
¿ ln 1,162 ≈ 970 ⋅10 años
ln 2
Es decir, la edad de la roca es aproximadamente igual a 970 ⋅106 años.
Prob.47: Es bien sabido que la cantidad de luz absorbida por una capa de agua de poco
espesor es proporcional a la cantidad de luz incidente y al espesor de la capa. Suponiendo
que una capa de agua de 35cm de espesor absorbe la mitad de la luz que incide sobre ella,
determinar que parte de la luz observará una capa de 2m.
Solución: Sea I (z) la cantidad de luz que atraviesa
una capa de agua de espesor z (fig.11). De
acuerdo con las condiciones de partida, la
cantidad I ( z + ∆ z ) −I (z ) de la luz absorbida es
igual kI ( z ) ∆ z(k=const .). Por eso, la cantidad de
luz I (z) que atraviesa la capa satisface Ja
licuación diferencial.
dI
=kI
dz
cuya solución es
kz
I ( z )=I ( 0 ) e

()
z /35
1 −1 1
La condición I ( 35 ) = I (0) nos proporciona k = ln 2. Por tanto, I ( z ) =I (0) , donde
2 35 2
z se mide en centímetros. Tomando en la última igualdad z=2 m=200 cm , obtenemos

I (0) ()
I (200) 1 40/ 7
=
2
. Entonces

I ( 0 ) −I (200)
()
40
1
=1− 7
≈ 0 , 98
I (0) 2

Así pues, podemos concluir que se absorbe el 98% de la luz que incide sobre la superficie.
Prob.48: Un paracaidista salto desde una altura de 1,5km y abrió el paracaídas a una
altura de 0,5km. ¿Cuánto tiempo duró el vuelo hasta la apertura del paracaídas? Se sabe
que la velocidad límite de caída del hombre en el aire de densidad normal es de 50m/s.
Despreciar la variación de la densidad del aire respecto a la altura. La resistencia del aire
es proporcional al cuadrado de la velocidad de caída.
Solución: Conforme a la segunda ley de Newton
dv 2
m −mg−k v (k >0)
dt
Donde m es la masa del paracaidista, g ≈ 10 m/s 2 es la aceleración de la caída libre,
k =const . Separando las variables v, t e integrando, obtenemos
dv
2
=dt
g−k 0 v

1
2 √k 0 g
ln
a−v| |
a+ v
=t−lnC

O bien
a=
√ g k =k
,
k0 0 m

|a−v
a+ v
|=C e , α =2 √ kgm …(1)
αt
Como la velocidad inicial es v ( 0 )=0, entonces C=1 . De la condición lim
a+ v (t)
t →+ ∞ a−v (t ) |
=+ ∞ se
|
deduce que v (t)→ a cuando t →+ ∞. Pero las condiciones de partida exigen que

lim
| |
a+ v (t)
t →+ ∞ a−v (t )
=+ ∞ se deduce que v (t)→ a cuando t →+ ∞. Pero las de partida exigen que

lim ( vt )=50 m/ s, luego a=50 o, lo que es lo mismo,


t →+ ∞

α =0 , 4 . Ya que 0 ≤ v < a, a partir de (1) hallamos:


√ mg
k
=50 . Por consiguiente, k =
mg
2500
,

ds(t )
v ( t )=50 th ( 0 , 2 t )=
dt
Integrando esta ecuación, tenemos:
s ( t )=250 ln ch ( 0 , 2 t ) + s0 …(2)

Utilizando la condición inicial s ( 0 )=0, encontramos que s0 =0. Sustituyendo luego en (2)
s=1000, llegando a la ecuación 1000=250lnch(0 ,2 t 1 ), de donde determinamos el tiempo t 1
de caída del paracaidista hasta el instante en que el paracaídas se abre:

t 1=5 ln(e4 + √ e8−1)≈ 5 ( ln 2+ 4 ) ≈ 23 s .

Prob.49: Un balón de fútbol de peso 0,4kp es lanzada hacia arriba con una velocidad inicial
de 20m/s. La resistencia del aire es proporcional al cuadrado de la velocidad; (en el caso
considerado es igual a 0,48p para una velocidad de 1m/s). Calcular la altura máxima a la
que asciende el balón y el tiempo que demora hasta alcanzar dicha altura. ¿Cómo varían lo
resultados si despreciamos la resistencia del aire?
Solución: Según la segunda ley de Newton tenemos que:
dv 2
m =−mg−k v …(1)
dt
En el caso considerado
P 0,4
m= =
8 10
2
kp ⋅ s
k =0,00048 2
m
Por tanto, la ecuación (1) toma la forma
dv 2
=−10−0,012 v
dt
Separando las variables e integrando, obtenemos
arctg √ 0,0012 v=√ 0 ,12 ( C−t )
10
v= tg ( ( C−t ) √ 0 , 12 ) …(2)
√0 , 12
1
Dado que v ( 0 )=20, de la ecuación (2) se deduce que C= arctg ( 2 √ 0 ,12 ) . Además, de
√ 0 ,12
la ecuación (2) se deduce que v=0para t=C ≈ 1 , 75 s (lo que corresponde, evidentemente,
ds (t)
a la altura máxima). Teniendo en cuenta que v ( t )= , después de sustituir v (t) en (2) e
dt
integrar la ecuación obtenida, tenemos que
250
s ( t )=s 0 + ln|cos ( ( C−t ) √ 0 , 12 )|, s0 =const . …(3)
3
La fórmula (3) expresa la ley de movimiento del balón. Tomando en (3) s(0)=0, hallamos
−250
que s0 = ln |cos ( C √ 0 ,12 )|. Si en la ecuación (3) hacemos t=C , obtenemos la altura
3
máxima a la que asciende el balón:
125
smáx =s0 ≈ ln 1 , 48 ≈ 16 ,3
3
Dejamos a cargo del lector el caso k =0 (cuando despreciamos la resistencia del aire).
Prob.50: Supongamos que un líquido fluye de un recipiente a través de un orificio con una
velocidad igual a 0 , 6 √ 2 gh, donde g=10 m/ s2 , hes la altura del nivel del líquido sobre el
orificio.
¿En cuánto tiempo saldrá todo el líquido de un recipiente cilíndrico de diámetro 2 H=1 , 8
m y altura H=2, 45 m a través de un orificio (en el fondo del recipiente) de diámetro 2 r=6
cm? El eje del cilindro es vertical.
Solución: Seah(t) la altura del nivel del líquido en el recipiente en un momento t >0. Al
cabo del tiempo ∆ t el nivel del líquido disminuye hasta h(t+ ∆ t). Por consiguiente, del
recipiente se escapa una cantidad de líquido igual a ¿. Por otro lado, a través del orificio se
escapa una cantidad de líquido igual a π r 2 v (t 1)∆ t , donde t 1 ϵ (t , t+ ∆t ) y v (t 1 ) es cierto
valor intermedio de la velocidad con que fluye el líquido en el intervalo de tiempo
(t ,t +∆ t) . De acuerdo con la ley de conservación de la masa tenemos la igualdad.

( )
2
( ) ( ) r
h t+ ∆ t −h t =− v (t¿ ¿1) ∆ t ¿
R

Dividiendo los dos miembros de esta ecuación por ∆ t y suponiendo que la función h es
diferenciable y que la función v es continua, hacemos tender ∆ t a 0, Entonces obtenemos
la ecuación diferencial
dh 2
=−k v (t)
dt
r
k= , v=0 ,6 √2 gh
R
Cuya solución es
2
h ( t )=(C−0 ,3 √ 2 g k 2 t ) , C=const .

Por cuanto h(0)=H , entonces C=√ H . Evidentemente, h(t)=0 para

10 √ H R
2
t= ≈ 1050 s=17 ,5 min
3 √ 2 g r2
Prob.51: Resolver el problema anterior suponiendo que le eje del cilindro está colocado
horizontalmente y el orificio se encuentra en la parte inferior.
Solución: Como vemos en la figura 12, cuando el nivel del
líquido disminuye en ∆ h durante el tiempo ∆ t , a través
del orificio E fluye una cantidad de líquido.
2 H √ h (2 R−h)∆ h+o (∆ h)
Por consiguiente,

−2 H √h ( 2 R−h ) ∆ h+o ( ∆ h )=π r v (t 1)∆ t


2

de donde, como en el ejemplo anterior, obtenemos la


ecuación diferencial

−2 H √ h ( 2 R−h ) dh=π r 0 ,6 √ 2 gh dt , h ≢0
2

La solución de esta ecuación es

( ) ,C=const … (1)
2 2
3 3 π r2 √ g 3
h ( t )=2 R−( 0 , 3 t+ C )
√2 H
Puesto que h(0)=2 R , entonces C=0. Tomando h=0 en la ecuación (1), hallamos el tiempo
al cabo del cual sale todo el líquido:
3 /2
40 R H
t 1= ≈ 1040 s
9 π r 2√g

Prob.52: Un embudo tiene la forma de un cono circular de radio


R=6 cm y altura H=10 cm, con el vértice dirigido hacia abajo
¿Al cabo de cuánto tiempo toda el agua del embudo se escapará
a través de un orificio circular de diámetro 0 , 5 cm situado en su
vértice?
Solución: A partir de la figura 13 vemos que la cantidad de agua
∆ V contenida en la capa sombreada es igual a π r 2 ∆ h, con una
precisión o (∆ h). Por otra parte, a través del orificio O fluye una
cantidad de agua igual a π r 21 v (t 1)∆ t . De esta manera, tenemos
la igualdad:
2 2
−π r ∆ h=π r 1 v ( t 1 ) ∆ t +o (∆ h),
Donde
r 1=0 ,25 cm , t 1 ϵ (t ,t +∆ t)

Pasando ahora al límite cuando ∆ t →O , obtenemos la ecuación diferencial


2 2
r dh+r 1 v ( t 1 ) dt=0

v ( t )=0 , 6 √ 2 gh
hR
De la semejanza de los triángulos AMO y C O1 O se deduce la relación r = . Por tanto, la
H
ecuación (1) se puede escribir en la forma
3 /2 2
h dh+ k dt=0
2
H 2
k =0 , 6 2 r 1 √ 2 g
2

R
Integrando, obtenemos
2 5 /2 2 ~ ~
h + k t=C , C=const .
5
~ 2 5/ 2
Puesto que h(0)=H , de aquí se deduce que C= H . Así pues, la solución del problema
5
planteado tiene la forma
5 /2 5/ 2 −5 2
h −H = k t
2
Tomando en esta igualdad h=0 hallamos que

( )
5 /2
2 H
t=
5 k2

Si el tiempo al cabo del cual todo el líquido se escapa del embudo. Los cálculos dan t ≈ 27 s.
Prob.53: En un tanque rectangular de dimensiones 60 cm ×73 cm y altura 80 cm entran 1 , 8
litros de agua por segundo, En el fondo del tanque hay un orificio de área 2 , 5 cm2 .¿En
cuánto tiempo se llena el tanque?
Solución: Sea h(t) la altura del nivel del agua en el recipiente, entonces
∆ V 1 =( h ( t+ ∆ t ) −h(t ) ) 60 ⋅75 es el incremento del volumen del agua en el intervalo de
tiempo desde t hasta ∆ t . Este el incremento (o decremento) de volumen compren de el
volumen ∆ V 2 de agua que entra y el volumen ∆ V 3 que fluye a través del orificio. De esta
forma, tenemos la ecuación ∆ V 1 =∆ V 2−∆ V 3. Por cuento ∆ V 2 =1800 ∆ t ,
∆ V 3 =2 ,5 ⋅0 , 6 √ 2 gh ( t 1 ) ∆ t , t 1 ∈(t , t + ∆ t ), entonces la última ecuación se puede escribir en
la forma

4500 ( h (t + ∆t )−h(t ) ) =1800 ∆ t −2 ,5 ⋅0 , 6 √ 2 gh ( t 1 ) ∆ t … (1)


3 2
g=10 cm/s
Dividiendo la ecuación (1) por ∆ t y pasando al límite cuando ∆ t →0 , obtenemos la
ecuación diferencial
dh
30 =(12−0 , 01 √2 gh)
dt
Integrando, hallamos
t+ C=
−6000
√2 g ( √h+ 1200
√2g )
ln ( 12−0 , 01 √ 2 gh ) …(2)

Sea h(0)=0. Entonces de (2) se deduce que C=−3600 ln 12. Sustituyendo h=80en (2),
hallamos el tiempo al cabo del cual el tanque se llena:

( 3
t 1=1200 3 ln −1 ≈ 260 s
2 )
Prob.54: Una cuerda elástica de 1m de longitud se alarga kf metros bajo la acción de una
fuerza f 1 kp . Hallar el alargamiento de esta cuerda ajo la acción de su propio peso
considerando que se cuelga por uno de los extremos. La longitud de la cuerda es l y su
peso es P.
Solución: SeaU ( x )el alargamiento de una cuerda de
longitud x y U (x + ∆ x )el alargamiento de una cuerda de
longitud x +∆ x . Entonces el alargamiento del elemento de
longitud Ax es igual a la diferencia U ( x+ ∆ x )−U (x )
(fig.14). Sobre el elemento de cuerda de longitud ∆ x
actúa la fuerza de tensión f , igual al peso de la cuerda de
P P
longitud l−x−θ ∆ x, es decir, f = (l−x−θ ∆ x ), donde
l l
es el peso específico de la cuerda y θ(0<θ<1)es una
magnitud introducida con el fin de considerar la
influencia de la fuerza que actúa sobre el elemento ∆ x y
que está determinada por el peso del propio elemento. De
acuerdo con las condiciones de partida, el elemento
indicado debe alargarse kf ∆ x metros. De esta manera,
obtenemos la
ecuación
P
U ( x+ ∆ x )−U ( x )=k ( l−x−θ ∆ x ) ∆ x …(1)
l
A partir de la ecuación (1) resulta la ecuación diferencial
dU kP
= (l−x )
dx l
integrando la cual obtenemos
kP
U ( x )=C+ (2l−x )x
2l
kP
Como U (0)=0, entonces C=0. Por tanto, U ( x )= (2 l−x) x . De la última fórmula
2l
obtenemos el alargamiento de la cuerda de longitud l :
kP
U (l)=
2l
Prob.55: Hallar la presión atmosférica a una altura h sin el nivel del mar la presión es
2 2
1 kp /cm y la densidad del aire es 0,0012 g/cm
Solución: Sea P(z ) la presión del aire a una altura z . Entonces la diferencia de las
presiones P(z )−P(z +∆ z ) es igual al peso de una columna de aire de base 1 cm2 y altura
∆ z , es decir, es igual a ρ(z +θ ∆ z )g ⋅ ∆ z , donde ρ es la densidad media del aire, 0<θ <1. De
este modo, tenemos
P ( z )−P ( z + ∆ z )=ρ( z +θ ∆ z )g ⋅∆ z

y pasando al límite cuando ∆ z →0 , obtenemos la ecuación diferencial


dP
=−gρ ( z ) …(1)
dz
Según la ley de Boyle—Mariotte, la densidad del aire a temperatura constante es
proporcional a la presión, es decir,
ρ ( z )=kP (z)

Utilizando esta igualdad, a partir de (1) hallamos


dP
=−kg
P
−kgz 2
P=P 0 e kp /cm

Puesto que P=1 kp /cm2 para z=0 , entonces P=e−kgz, y como


2 2 2
0,0012 g/cm = p/cm ⋅k =1000 g /cm gk
donde g=9 , 8 m/ s2 =1 p/ g es la aceleración de la caída libre, hallamos que
−5 −1
kg=0 , 12 ⋅10 cm . De este modo, a una altura de h km la presión del aire es igual a
−0 ,12 h
P=e (kp / cm) .
Prob.56: Para detener los barcos, de éstos se lanzan unas cuerdas gruesas (amarras) que
se enrollan alrededor de unos postes enclavados en el muelle (amarraderos). Determinar
la fuerza de frenado que actúa sobre el barco si la soga da 3 vueltas alrededor del
amarradero, el coeficiente de fricción entre la soga y el poste es igual a 1/3 y el obrero del
muelle tira del extremo libre de la soga con una fuerza igual a 10kp.
Solución: Sea T (φ) la fuerza de tensión (de la soga) correspondiente al ángulo de
arrollamiento φ , y sea ∆ P la fuerza de reacción normal (del poste) que actúa sobre un
elemento de la soga de longitud ∆ S=R ∆ φ (fig.15). De la condición de equilibrio de las
fuerzas T (φ), ∆ P y T (φ+ ∆ φ), se deducen las igualdades siguientes (con una precisión o (
∆ φ ¿:
T ( φ+∆ φ ) ∆ φ+ o(∆ φ)=∆ P
T ( φ+∆ φ )−T ( φ ) +o ( ∆ φ )=∆ F … (1)

donde ∆ F es la fuerza de fricción que actúa sobre el


elemento dado. Según las condiciones de partida
∆ F=k ∆ P; por tanto, a partir de (1) obtenemos la
ecuación:
T ( φ+∆ φ )−T ( φ )=KT ( φ+∆ φ ) ∆ φ+ o(∆ φ)

de la cual, pasando al límite, fácilmente obtenemos


la ecuación diferencial
dT
=kT

Su solución es T =T 0 ekφ . Según la condición complementaria, para φ=0 se tiene T 0=10 kp .


Si φ=6 π (lo que corresponde a tres vueltas), obtenemos T =10 e 2 π ≈ 5355 kp.
Prob.57: Sobre un disco que gira dentro de un líquido actúa una fuerza de fricción
proporcional a la velocidad anular de rotación. Dicha fuerza, evidentemente, frena el
movimiento del disco. Hallar la dependencia de la velocidad angular respecto al tiempo si
la velocidad angular inicial es igual a 100 revoluciones por minuto, y al cabo de un minuto
ya es igual a 60 revoluciones por minuto.
Solución: Sea ω (t) la velocidad angular del disco. De acuerdo con la ley de conservación
del momento angular, tenemos:

I =M …(1)
dt
donde I es el momento de inercia del disco, M es el momento de las fuerzas que actúan
sobre el disco. Según las condiciones de partida M =k 0 ω(k 0=const ); entonces la ecuación
(1) toma la forma:
dω k0
=k ω , k=
dt I

Su solución es ω=ω 0 ekt , donde ω se mide en revoluciones por minuto y el tiempo, en


minutos. Teniendo en cuenta la condición inicial ω 0=100, se tiene que 60=100 e k, de donde
k
e =0 ,6 . De este modo, la dependencia requerida es
t
ω=100(0 ,6) rev /min

Prob.58: En una habitación cerrada de volumen V −m0 se encuentra un recipiente abierto


que contiene agua. Se sabe que la velocidad con que se evapora el agua es proporcional a
la diferencia entre la cantidad q de vapor de agua que satura 1 m3de aire a una temperatura
dada y la cantidad de vapor de agua q en 1m3 de aire en el instante analizado. Se supone
que las temperaturas del agua y del aire no varían, y el área de la superficie de
evaporación permanece invariable. Inicialmente, en el recipiente habían m0 gramos de
agua y en un metro cúbico de aire, q 0 gramos de vapor ¿Qué cantidad de agua quedará en
el recipiente al cabo del tiempo t ?
Solución: Sea Q(t )la cantidad de agua (medida en gramos) evaporada durante el tiempo t .
Entonces, de acuerdo con las condiciones de partida la velocidad de evaporación es
dQ
=k ( q 1−q ) …(1)
dt
dondek es el coeficiente de proporcionalidad. Hallemos el valor de q . Es evidente que
−6
Q1=(V −m0 ⋅10 )q0 es la cantidad de gramos de vapor que había en la habitación y
Q2=( V −m0 ⋅10−6 ) q 0 +Q es la cantidad de gramos de vapor contenida en la habitación en el
instante t . Es obvio que la cantidad Q2 está uniformemente distribuida en el volumen
−6 −6
V 2=V −m0 ⋅10 +Q ⋅10 T6 (la presencia del factor 10−6 se debe a que un gramo de agua
ocupa aproximadamente un volumen igual a un centímetro cúbico). Por tanto,

Q 2 ( V −m0 10 ) q 0 +Q
−6

q= = … (2)
V 2 V −10−6 (m 0 +Q)

Si V ≫ 10−6 (m0 +Q), entonces de (2) se obtiene una fórmula más simple para q:
V q0 +Q Q
q= =q 0 + …(3)
V V
Sustituyendo (3) en (1), obtenemos la ecuación diferencial
dQ
dt (
=k q1−q0 −
Q
V )
cuya solución es
−kt
V
Q ( t ) =V ( q1−q0 ) +C e

Puesto que Q ( 0 )=0, entonces C=−V (q1−q 0 ). De este modo, la cantidad de agua que
queda en el recipiente en el instante t está dada por la fórmula
−kt / v
m1 ( t )=m0 −V (q 1−q 0)(1−e )
Prob.59: La masa de un coherente con las reservas llenas de combustibles es igual a M ,
mientras que sin combustible la masa es m . La velocidad con que se expulsan los
productos de la combustión es igual a c . La velocidad inicial del cohete es igual a cero.
Hallar su velocidad después de consumirse todo el combustible. Despreciar la fuerza de
gravedad y la resistencia del aire.
Solución. Sean M (t) y v (t) la masa y la velocidad el instante t , respectivamente. Entonces
la cantidad de movimiento del cohete. Según la ley conservación de la cantidad de
movimiento tenemos que:
k ( t |∆ t )=k ( t ) +∆ p …(1)

donde p es el impulso de las fuerzas externas que actúan, en el transcurso del intervalo de
tiempo ∆ t desde t hasta t+ ∆ t . En el caso dado, este impulso es generado por la expulsión
de los productos de la combustión, cuya velocidad absoluta es c−v (t) (siendo c la
velocidad relativa de esos productos respecto al cohete). Puesto que se quema una masa
M ( t )−M (t+ ∆ t), entonces

∆ p=(c−v ( t )) ( M ( t )−M ( t+ ∆ t ) ) … (2)

Sustituyendo (2) en (1), obtenemos


M ( t+ ∆t ) v ( t+ ∆ t )−M ( t ) v ( t ) =(c −v ( t ) ) ( M ( t )−M ( t+ ∆ t ) )
O bien
M ∆ v + c ∆ M +o ( ∆ v )=0

Dividiendo por ∆ v y haciendo tender ∆ v → 0 , llegamos a la ecuación diferencial


dM
M +c =0
dv
cuya solución es
1
v ( t )=cln +c 0
M (t )

La constante c 0 se determina a partir de la condición inicial v ( t ) ¿ t=0=0, es decir, de que en


el momento inicial (o, lo que es lo mismo, para M ( t )=M ) la velocidad v es igual a cero. Por
1
tanto, de la ecuación (3) resulta que c 0=−cln . Por consiguiente,
M
M M
v=cln =cln
M (t) M −x

donde x es la masa de combustible que se quemó. Si se quema todo el combustible, es


decir, x=M −m, entonces el cohete tendrá la velocidad
M
v=cln
m
Ésta es la fórmula de Tsiolkovski.

3) Ecuaciones homogéneas y ecuaciones reducibles a homogéneas


Prob.60: ( y + √ x 2− y 2)dx−xdy =0

Solución: Aquí las funciones M ( x , y )= y √ x2 − y 2 (|x|≥| y|) y M ( x , y )=−x son homogéneas de


grado m=1, pues
M ( tx ,ty )=ty+ √t 2 (x 2− y 2)=t ¿

para t ≥ 0 (es decir, a=0, b=+∞ ) y


N ( tx , ty )=−tx=tN (x , y )

para todo t . Por consiguiente, la función dada es homogénea. Realizando el cambio de


variable y=xu , obtenemos dy =xdu+udx , y la ecuación toma la forma

( xu +|x|√1−u2 ) dx−x ( xdu+udx )=0 … (1)


Es evidente que x=0 satisface la ecuación original. Por tanto, asumiendo x ≠ 0 en (1),
obtenemos
sgnx √ 1−u dx−xdu=0
2

Separando las variables e integrando, hallamos


dx du
=
|x| √1−u 2
y
sgnxln|x|=arcsen +C
x
Teniendo en cuenta que u=±1 (esto es, y=± x ) también son soluciones de la ecuación
original, escribamos definitivamente todas las soluciones:
y
sgnxln|x|=arcsen +C , y =± x , x=0
x
Nota: Hablando con rigor, la solución x=0 se obtiene suponiendo que √ − y 2 ⋅ 0=0 para
todo valor de y

Prob.61: xy ' = y (1+lny−lnx)

Solución: Si escribimos esta ecuación en la forma

(
xdy− y 1+ ln
y
x )
dx=0

(x >0 , y >0)

y
descubrimos que las funciones M (x , y )=− y 1+ ln
x ( )
y N ( x , y )=x son homogéneas y del
mismo grado m=1. Por tanto, empleando el cambio de variable
y=xu
dy =xdu+udx
Obtenemos
xdu−ulnudx=0
Separando las variables e integrando, hallamos
lnx−ln |lnu|=lnC (u ≠ 1)

O bien
Cz
y=x e
Observemos que la solución y=x correspondiente al valor u=1, pertenece a la familia de
curvas integrales para C=0 .

Prob.62: xydx +( y 2−x ¿ ¿ 2) dy=0 M (1 , 1) ¿

Solución: En el ejemplo dado se pide hallar la curva que satisface la ecuación diferencial y
que pasa por el punto M .
Hallemos todas las soluciones de la ecuación. Haciendo el cambio de variable y=xu ,
obtenemos
u3 dx + x ( u 2−1 ) du=0
De aquí, para u ≠ 0, separando las variables e integrando, hallamos
2
∫ u u−1
3
dx
du+∫ =lnC
x
1
ln |ux|+ 2
=lnC
2u
O bien
| y|
x 2+ 2 y 2 ln =0 …(1)
C
A la familia de curvas integrales obtenidas añadimos, además, la curva y=0. Tomando en
(1) x=1, y=1 llegamos a que C=e 1/ 2. Por consiguiente, la ecuación de la curva buscada es
x 2+ y 2 ( ln y 2−1 )=0

( )
x
Prob.63: xy e y + y 2 dx−x 2 e x/ y dy =0

Solución: En este caso es cómodo hacer el cambio de variable x=uy . Entonces


dx=udy + ydu, udy + y (u e u +1)du=0. Separando las variables, obtenemos la ecuación
u
dy (u e +1)du
+ =0
y u
cuya solución es ln |x|+ e x/ y =0

Prob.64: ( 2 √ xy− y ) dx−xdy=0

Solución: Es evidente que xy ≥0 . Para eludir las soluciones triviales x=0 e y=0, asumimos
que xy >0 . Realizando el cambio de variable y=ux (u> 0), tenemos
dy =udx+ xdu

2 ( √ u sgnx−u ) dx−xdu=0 …(1)


Si x >0, después de separar las variables x y u (u ≠ 1) , hallamos
dx 1 du
∫ = ∫
x 2 √ u−u

lnx=−ln |1− √ u|+ lnC


O bien

| √|
x 1−
y
x
=C
Si x <O , haciendo x=−x 1 en (1), obtenemos

2 ( √ u+u ) d x 1 + x 1 du=0

Separando las variables e integrando, tenemos

x 1 1+
( √ )
−y
x1
=C

O bien

−x 1+ ( √ xy )=C
Reuniendo las dos respuestas (para x >0 y para x <0) en una, finalmente podemos escribir
la solución de la ecuación dada en la forma
x−√ xy=C …(2)
Señalemos que la solución y=x (u=1) satisface la ecuación (2) para x >0(C=0), mientras
que la solución y=0 no la satisface. Así pues, todas las soluciones de la ecuación dada se
expresan por la fórmula (2) con la adición de y=0. La solución x=0 se puede incluir en (2)
sólo formalmente, ya que la diferencial d (x− √ xy ) no existe para x=0 .
' y
Prob.65: x y = ycosln
x
Solución: Mediante el cambio de variable y=xu (u>0 , x ≠ 0) la ecuación dad se lleva a la
forma
'
x u =u(coslnu−1)

Asumiendo que lnu ≠ 2 kπ (u≠ e 2 ka , k ϵ Z) , obtenemos


dx du
=
d u ( coslnu−1 )

d ( lnu )
ln (|x|C )=∫ =ctg(ln √ u)
coslnu−1
O bien
ln ( Cx )=ctgln
√ y
x
Asimismo, es fácil comprobar que para todo k ϵ Z la función y=x e 2 ka, también es solución
de la ecuación dada.
Prob.66: ( 6 x + y−1 ) dx + ( 4 x+ y−2 ) dy=0

Solución: Esta ecuación es de la forma (2), p. 3.2. En el caso dado a 1=6, b 1=1, a 2=4 , b 2=1,
a 1 b 2−a2 b1=2 ≠ 0. Por consiguiente, para reducirla a una ecuación homogénea utilizamos el
cambio de variable x=u+ α , y=v + β . De este modo, obtenemos
dx=du
dy =dv
( 6 u+ v +6 α + β−1 ) du+ ( 4 u+v + 4 α + β−2 )=0
−1
La ecuación diferencial obtenida se reduce a una ecuación homogénea haciendo α = ,
2
β=4.
1
De esta manera, mediante los cambios de variable x=u− , y=v +4 la ecuación original se
2
reduce a una homogénea:
( 6 u+ v ) du+(4 u+ v)dv=0
Utilicemos ahora el cambio ya conocido: u=vℇ , du=v +dℇ + ℇdv . Tenemos
( 6 ℇ +1 ) vdℇ + ( 6 ℇ2 +5 ℇ +1 ) dv=0. Separando las variables e integrando, obtenemos
1 1
(ξ ≠− , ξ ≠− )
2 3
6 ξ+1 dv
dξ+ =0
2
6 ξ +5 ξ +1 v

1 6 ξ+ 1
∫ dξ+ ln |v|=lnC
6
( )( )
6 ξ+
1
2
ξ+
1
3
o bien
2 3 ℇ +1
(2 ℇ +1) =C
v
Regresando a las variables x e y , finalmente hallamos

(
( 2 x+ y−3 )2=C 3 x + y−
5
2 )
Las soluciones 2 x+ y−3=0 ξ= (−1
2 ) 5
y 3 x+ y− =0 ξ=
2
−1
3 ( )
) se pueden incluir en la
familia de curvas integrales obtenidas para C=0 y C=∞ , respectivamente.
Prob.67: ( 5 x−7 y +1 ) dy + ( x+ y−1 ) dx=0
Solución: Por cuanto las rectas 5 x−7 y +1=0 y x + y−1=0no son parálelas
(a ¿ ¿ 1b 2−a 2 b1 ≠ 0)¿ empleamos el cambio de variable.

x=u+ α
y=v + β
donde las constantes α y β satisfacen el sistema de ecuaciones
5 α −7 β +1=0
α + β−1=0
1
Resolviendo este sistema, obtenemos α =β= . Efectuando los cambios de variable
2
1 1
x=u+ , y=v + , llegamos a la ecuación diferencial
2 2
( 5 u−7 v ) dv+ ( u+v ) du=0
Al igual que en el ejemplo anterior, realicemos el cambio de variable u=ℇv Entonces,
du=vdξ +ξdv , ( ξ 2+ 6 ξ−7 ) dv + v ( ξ+1 ) dξ=0
Separando las variables e integrando, obtenemos
( ξ+1 )
∫ ( ξ−1 ) ( ξ +7 ) dℇ + ln (|v|C )=0
1 3
ln |ξ−1|+ ln |ξ +7|ln |vC|=0
4 4
De donde
( ξ−1 ) ( ξ−1 )3 v 4 =C
O bien
( x− y ) ( x+ 7 y−4 )3=C
Prob.68: ( 2 x+ y+ 1 ) dx −( 4 x +2 y−3 ) dy =0
Solución: Puesto que las rectas 2 x+ y+1=0 y 4 x+ 2 y−3=0 son paralelas
(a ¿ ¿ 1b 2−a 2 b1 ≠ 0)¿ , no podemos recurrir al cambio de variable empleado en los ejemplos
anteriores. Sin embargo, como los coeficientes de las variables x e y son proporcionales,
entonces podemos tomar z=2 x+ y , obteniendo de este modo dy =dx−2 dxy
5 ( z−1 ) dx−( 2 x−3 ) dz=0

Integrando esta ecuación y regresando a la variable y, obtenemos


2 y− z
2 x+ y−1=C e

( )
3
' y +2
Prob.69: y =2
x + y−1

Solución: Mediante el cambio de variable x=u+3 , y=v−2reducimos esta ecuación a la


ecuación homogénea
2
dv v
=2 2
du (u+ v)

Empleando otro cambio de variable v=uξ (u) obtenemos la ecuación diferencial de


variables separables
2
dξ 2ξ
u = −ξ
du (1+ξ )2
Integrando, hallamos
2
ξ +2 ξ +1
∫ du
u
+∫ 2
ξ(1+ ξ )
dξ=C

ln |uξ|+2 arctgξ=C
Después de regresar a las variables x e y obtenemos
y +2
−2 arctg
x−3
y +2=C e
y+ x y + x
Prob.70: ( y ¿¿ '+ 1)ln = ¿
x +3 x+3
Solución: Mediante el cambio de variable x=u−3, y=v +3, la ecuación dada se reduce a la
ecuación homogénea
( dvdu +1) ln u+u v = u+u v
la cual, a su vez, se reduce con ayuda del cambio de variable v=uξ (u) a una ecuación de
variables separables:

( u ξ ' + ξ+1 )=1+ξ


u ξ ln ( 1+ξ )=(1+ξ ) ( 1−ln ( 1+ ξ ) )
'

Separando las variables en la última ecuación e integrando, hallamos


du ln ( 1+ ξ ) dξ
∫ u
+∫
(1+ξ ) ( 1−ln ( 1+ξ ) )
=lnC

de donde, regresando a las variables x e y , obtenemos


x+ y C
ln =1+
x+3 x+ y
' y−2 y−2
Prob.71: y = +tg
x+1 x+1
Solución: Esta ecuación se reduce a una ecuación homogénea con ayuda del cambio de
variable x=u−1, y=v−2
dv v v −2u
= +tg
du u u
Al igual que en los ejemplos anteriores, el cambio de variable v=uξ (u) proporciona
'
u ξ =tg (ξ−2)
Separando las variables e integrando, obtenemos
dξ du
=
tg(ξ−2) u

ln |u|−ln|sen ( ξ−2 )|=lnC

Regresando a las variables iniciales, finalmente obtenemos

x +1=Csen ( y−2
x +1 )
x

y−2 x
Señalemos que las soluciones x +l=0(u=0)y =kr k ∈ Z , ( tg ( ξ−2 )=0 ), pertenecen a la
x+1
familia obtenida de curvas integrales para C=0 y C=∞ , respectivamente.

Prob.72: x 3 ( y ' + x ) = y 2

Solución: Esta ecuación no es homogénea; sin embargo, utilizando el cambio de variable


α
y ( z ( x ) ) , vemos que la ecuación

α z α−1 x 3 dz=( z 2 α + x 4 ) dx …(1)


se reduce a una ecuación homogénea si tomamos α =2. Haciendo z=xu en (1), obtenemos
2
( u2 −1 ) dx=2uxdu
Separando las variables e integrando, hallamos
1
ln |x|+ 2
=C
u −1
o bien
2
x
ln |x|+ 2
=C
y−u
Prob.73: 2 x 2 y ' = y 3 + xy
α
Solución: Hagamos el cambio de variable y=z ; entonces
2 α x z z =z + x z ,2 α x z dz−( z + x z ) dx=0 . Vemos que las funciones 2 α x 2 z α −1 z ' y
2 α −1 ' 3α α 2 α−1 3α α

3α α
z + x z , son homogéneas si, y sólo si, se cumple la condición
α +1=3 α=α +1
1
es decir, si α = . De esta manera, asumiendo que y ≥ 0 utilizamos el cambio de variable
2
y= √ z después de lo cual la ecuación original adquiere la forma
2
x dz−( z−x ) zdx =0
Con ayuda del cambio de variable z=xu (x), esta ecuación se transforma en la ecuación de
variables separables
2
xdu−u dx=0
Integrando, obtenemos
1
+ ln| x|=C
u
O bien
x
2
+ ln|x|=C
y
La solución y=0 pertenece a la familia de soluciones obtenida para C=∞ .
Prob.74: ydx + x (2 xy +1 ) dy=0
Solución: Empleando el cambio de variable y=z α , obtenemos la ecuación

z α dx+ x ( 2 x z α + 1 ) α z α −1 dz=0

Puesto que las funciones z α y αx z α −1 ( 2 x z α +1 ) son homogéneas y del mismo grado para
α =1, entonces la ecuación dada se reduce a una ecuación homogénea con ayuda del 1
1
cambio de variable y= . Asimismo, la ecuación examinada se puede reducir directamente
z
u(x )
a una ecuación de variables separables si se utiliza el cambio de variable y= (en el
x
caso general hay que emplear el cambio de variable y=x α u( x )¿.

De esta manera, llegamos a la ecuación


2
−2 u dx+ x ( 2 u+1 ) du=0
de donde
dx du du
= +
x u 2 u'
ln |ux|+ 21u =lnC
es decir,
1
2 xy= … (1)
lnC| y|
o bien
2 −1/ xy
y e =C …(2)
Nótese que la fórmula (2) no contiene las soluciones x=0 e y=0, mientras que la fórmula
(1) incluye la solución x=0 (para C=∞ ¿ . En cuanto a la solución y=0, no es difícil
demostrar que no existe ninguna transformación de las fórmulas (1) o (2) que permita
incluir dicha solución en la familia de soluciones obtenida.
Prob.75: 2 y ' + x =4 √ y
Solución: En este ejemplo se analiza una ecuación homogénea generalizada de grado dos,
es decir, α =2. Por consiguiente, con ayuda del cambio de variable y=x 2 u(x) (u ≥ 0) la
ecuación se reduce directamente a la ecuación de variables separables

2 xdu+ ( 1+ 4 u−4 √ u ) dx=0 , ( x ≥ 0 ) … (1)


Integrando esta ecuación, resulta

x ( 1−2 √u )=Cexp
{√ }
1
2 u−1
, u=
1
4

o bien

{ }
2
( x−2 √ y ) exp x x
=C , y =
x−2 √ y 4

Si x <0, en lugar de la ecuación (1) obtenemos la ecuación

2 xdu+ ( 1+ 4 u−4 √ u ) dx=0


la cual, después de ser integrada, nos conduce al mismo resultado obtenido para la
ecuación (1).

Prob.76: 2 ( √ x 4 y 2+1−x 2 y ) dx−x 3 dy=0

Solución: Haciendo y=z α , obtenemos la ecuación

2 ( √ x y +1−x z ) dx−α x z
4 2 2 α 3 α −1
dx=0

Vemos que los factores de las diferenciales son funciones homogéneas de un mismo grado
u(x )
sólo para α =−2. De esta manera, el cambio de variable y= 2 - reduce la ecuación
x
inicial a la ecuación de variables separables

2 √ u 2+1 dx −xdu=0
de donde hallamos
x =C ( u+ √ u + 1 )
2 2

o bien
x =C ( x y + √ 1+ x y )
2 2 4 2
Prob.77: Hallar la curva para las cual las distancias del punto de intersección de cualquiera
de sus tangentes con el eje de intersección de
cualquiera de sus tangentes con el eje de abscisas hasta
el punto de tangencia y hasta el origen de coordenadas
son iguales.
Solución: De aurd con las condiciones de partida,
|OK |=|KM| (fig.16). Por tanto,
α =2 β
Y
2tgβ
tgα= 2
1−tg β
Como
y
tgβ =
x
Entonces
2 xy
tgα= 2 2
x −y
Teniendo en cuenta el significado geométrico de la derivada, finalmente resulta
' 2 xy
y= 2 2
x −y
La ecuación diferencial obtenida es homogénea. Para resolverla utilicemos el cambio de
variable y=xu (x), el cual proporciona
' 2u
x u +u= 2
1−u
2
dx 1−u
= du
x u+ u3

Integrando, hallamos
x 2
( u +1 )=C
u
Regresando a las variables x e y ,obtenemos
2 2
x + y =Cy
Prob.78: Hallar la curva para la cual la distancia desde el origen de coordenadas hasta
cualquier tangente es igual a la abscisa del punto de tangencia (fig.17)
Solución. En la figura 17 vemos que los triángulos
ONM y OML son congruentes, pues según las
condiciones de partida |ON |=|OL| y, además, la
hipotenusa OM es común a ambos triángulos. Por
consiguiente, los ángulos ^MOL y ^MON son iguales.
Como ^
π
KON = −α , entonces ^
2
MOL=
1 π
2 2 ( )
+α . A

partir del triángulo OLM, obtenemos


y
x (
π α
=tg + = )
1+ tgα /2
4 2 1−tgα /2

α y−x 2 tg (α /2)
de donde tg = . Finalmente, utilizando la fórmula trigonométrica tgα= y
2 y+x 1−tg 2 ( α /2)
el hecho de que tgα= y ' , obtenemos la ecuación diferencial
2 2
' y −x
y=
2 xy
Mediante el cambio de variable y=xu (x), esta ecuación se reduce a la ecuación de
variables separables
dx 2 udu
+ =0
x u2 +1
cuya solución es
ln |x|+ ln ( u2+ 1 )=lnC
o bien
2 2
x + y =Cx
Prob.79: Hallar la curva para las cuales |NL|=|OM|−|ON| (fig.18), donde |OM| es el
módulo del vector de posición del punto de
tangencia y N es su abscisa
Solución: A partir de la figura 18, hallamos
|NL|= ytgα
|OM|=√ x 2+ y 2
|ON |=x
Entonces
y y ' =√ x2 + y 2−x
Esta ecuación diferencial es homogénea. Para reducirla a una ecuación de variables
separables hacemos el cambio de variable y=xu (x). Si x >0, tenemos
2
2 dx du
= 2
x √ u +1−1−u 2
2
dx du
2∫ =∫ 2 +C
x √ u + 1−1−u 2
Integrando, resulta
x (1−z )=C
Donde
z=√ u2 +1
De este modo, la familia de curvas integrales buscada tiene la forma

x−√ x 2+ y 2=C
Dejamos el caso x <0 a cargo del lector.
Prob.80: Hallar la ecuación de la curva si se sabe que el triángulo formado por los ejes
coordenados y la normal a la curva en cualquiera de sus puntos (sea M, fig.19) tiene un
área igual al del triángulo formado por el eje Ox, un segmento de la normal y la tangente
en el punto M.
Solución: Según las condiciones de partida, los triángulos
KMN y OLN tienen áreas iguales (fig. 19); por tanto, los
triángulos KOP y PML también tienen áreas iguales.
Puesto que los últimos son, además, semejantes, entonces
ellos son congruentes. Por consiguiente, |OP|=| PM| A
partir de la ecuación de la tangente Y − y= y '( X−x) ,
donde X e Y son las coordenadas de un punto de la
tangente, se deduce que
|PO|= y− y ' x
La longitud del segmento PM es
|PM |= √ x 2+ y '2 x 2
De este modo, obtenemos la ecuación diferencial homogénea
' 2 2 '2 2
( y− y x) =x + y x
o bien
2 ' 2
y −2 y y −x =0
Mediante el cambio de variable y=xu (x), esta ecuación se reduce a la ecuación de
variables separables
2 '
u −2 u u x −1=0
de donde, integrando, hallamos

x ( 1+ u2 )=2 C
O bien
( x−C )2+ y 2=C 2
Ésta es la ecuación de la curva buscada.
Prob.81: Hallar la relación entre α y β para la cual la ecuación y ' =a x α +b y β se reduce a una
ecuación homogénea mediante el cambio de variable y=2π
Solución: Efectuando el cambio de variable indicado, obtenemos
m −1 ' α β
mz z =a x + b y , ( m ≠0 ) …(1)
de donde se deduce que si a y b no son nulos, entonces para que la ecuación (1) sea
homogénea es necesario y suficiente que se cumplan las igualdades
m−1=α =mβ
De la segunda igualdad obtenemos
α
m= ,(β ≠ 0; α ≠ 0)
β
Sustituyendo m en la primera igualdad, llegamos a la relación buscada:
1 1
− =1
β α
Prob.82: Demostrar que las curvas integrales de la ecuación
2 2
( ax +by −c 1 ) dx+ ( ay−bx + c2 ) dy =0 ,(a + b ≠ 0) son espirales logarítmicas
Solución: Efectuando un desplazamiento paralelo
x=u+ α
y=v + β
del sistema de coordenadas, reducimos la ecuación inicial a la forma
( au+ bv ) du+ ( av−bu ) dv=0
Es conveniente pasar a las coordenadas polares tomando u=ρcosφ, v=ρsenφ , donde ρ , φ
son coordenadas polares con origen en el (α , β ) respecto al sistema de coordenadas Oxy.
Tenemos:
aρ ' =bρ
de donde hallamos
(b /a )φ
ρ=C e
De esta manera, hemos obtenido la familia de espirales logarítmicas buscada.
Prob.83: Hallar la forma del espejo que refleja, paralelamente a una dirección dada, todos
los rayos provenientes de un unto dado.
^ ^
Solución: En la figura 20, vemos que OMN= NMS , es decir, el ángulo de incidencia del rayo
OM es igual al ángulo de reflexión del rayo MS. Supongamos que la dirección de los rayos
MS reflejados es paralela al eje Ox. Entonces, a partir de las condiciones señaladas
obtenemos
π ^
α + β= , K MO =α
2
Por consiguiente, el triángulo KMO es isósceles, esto es, |KO|=|OM|. Evidentemente,
|OM|=√ x 2+ y 2. Hallamos la longitud del segmento KO calculando la abscisa del punto de
intersección de la tangente con el eje Ox. A partir de la ecuación de la tangente a la curva
buscada
Y = y+ y ' (X −x)
donde X, Y son las coordenadas de un punto de la tangente, hallamos la abscisa requerida:
y
X =x− '
y
Entonces |KO|=− X . De este modo, tenemos la ecuación diferencial
y
−x= √ x + y
2 2
'
y
Esta ecuación es homogénea. Empleando el cambio de variable y=xu (x), tras una serie de
cálculos sencillos obtenemos la solución
2 2
y −2 Cx−C =0
que describe una familia de parábolas (C ≠ 0).
Prob.84: Sea k 0 una raíz de la ecuación f (k )=k . Demostrar que

'
1. Si f ' (k )<1, encontramos ninguna solución de la ecuación y =f ( xy ), salvo la
solución y=k 0 x, es tangente a la recta y=k 0 x en el origen de coordenadas
2. Si f ' ( k )> 1, entonces existe una cantidad infinita de soluciones tangentes a la recta
'
y =k 0 x

Solución: Emplearemos el método de reducción al absurdo. Supongamos que


' '
()
f ( k 0 )=k 0 , f ( k 0 ) <1 , y que la ecuación y =f
y
x
tiene una solución que es tangente a la recta
y ¿ k 0 x en el origen de coordenadas. Entonces, en un entorno del origen de coordenadas
dicha solución tiene la forma
y ( x ) ¿ k 0 x + xδ ( x ) …(1)

donde δ ( x ) → 0 cuando x → 0. Sustituyendo (1) en la ecuación inicial llegamos a la


ecuación
k 0 x + x δ ' +δ=f ( k 0 +δ ( x ) ) =f ( k 0 ) + f ( ) δ ( x ) +o (δ)
' k
0

o bien
'
x δ = Aδ +o (δ )
de donde
xδ ' o(δ)
= A∨
δ δ
o(δ)
Despreciando (pues tiende a cero cuando x → 0), podemos escribir
δ
xδ '
≈A
δ
A
de donde hallamos δ ≈C |x| De la última fórmula se infiere que si C ≠ 0 , entonces δ no
tiende a cero cuando x → 0 (puesto que A< 0). Así pues, hemos llegado a una
contradicción. De este modo, la afirmación 1) es cierta.
2) En este caso no aparece contradicción alguna, ya que A> 0 y δ → 0 cuando x → 0 para
todo C.
Prob.85: Se destilan dos líquidos x e y. Se sabe que en todo instante el cociente de las
cantidades evaporadas es proporcional al cociente de las cantidades que se encuentran en
estado líquido. Determinar la dependencia entre x e y .
Solución: Sean x (t) e y (t ) las cantidades no evaporadas de los líquidos en el instante t .
Entonces x (t +∆ t ) e y ( t+ ∆ t ) son las cantidades de los líquidos no convertidas en vapor en
el instante t+ ∆ t . Por consiguiente, durante el tiempo ∆ t se evaporan las siguientes
cantidades de cada líquido:
x (t )−x ( t+ ∆ t ) e y ( t )− y (t + ∆t )

Según las condiciones de partida,


y ( t +∆ t )− y (t) y (t 1 )
=k …(1)
x ( t +∆ t )−x (t ) x (t 1)

donde k es el coeficiente de proporcionalidad, t 1 ϵ ( t , t+ ∆ t ). Asumiendo que las funciones x


e y son diferenciabas, a partir de (1) y pasando al limité cuando ∆ t →0 , obtenemos la
ecuación diferencial
dy y
=k
dx x
Integrando esta ecuación hallamos la dependencia buscada y=C x k

4) Ecuaciones lineales y ecuaciones reducibles a lineales


Prob.86: y ' + ytgx=secx
Solución: En primer lugar, hallemos la solución general de la ecuación homogénea
asociada
'
y + ytgx=0
Separando las variables e integrando, hallamos
y=Ccosx
La fórmula (1) es la solución general de la ecuación homogénea, donde C es una constante
arbitraria. Para obtener todas las soluciones de la ecuación lineal no homogénea,
consideramos C como una función de x, es decir, C=C ( x) y exigimos que la función
y=C (x ) satisfaga la ecuación original dada:
'
C cosx −Csenx+Ccosxtgx=secx
' 1
o bien C = 2 . De aquí hallamos que C (x)=tgx+C 0, donde
C 0 es una constante
cos x
arbitraria nueva. Sustituyendo la expresión de C(x) en (1), finalmente obtenemos
y=senx +C 0 cosx

Prob.87: ( 2 x+1 ) y ' =4 x+ 2 y

Solución: Resolvamos primero la ecuación homogénea asociada


( 2 x+1 ) y ' =2 y
Su solución general es y=C (2 x+ 1). Empleando el método de variación de la constante,
obtenemos

( C ' ( 2 x+ 1 )+ 2C ) ( 2 x +1 ) =4 x +2 C(2 x +1)


o bien ( 2 x+1 )2 C ' =4 x , de donde hallamos
xdx 1
C ( x )=4 ∫ +C 0=ln |2 x +1|+ +C 0
( 2 x +1 ) 2
2 x+1

Así pues, la solución de la ecuación planteada es


y= ( 2 x +1 ) ( ln |2 x+1|+ C ) +1

Prob.88: ( xy +e z ) dx−xdy=0

Solución: Asumiendo que dx ≠ 0 (para eludir la solución trivial x=0 ), escribimos la


ecuación en la forma
' z
x y −xy =e
La ecuación homogénea asociada xy '−xy =0 tiene la solución general y=C e x . Aplicando
el método de variación de la constante, obtenemos

x ( C+C ' ) e z−xC e z=e z


' 1
de donde C = , luego C=ln |x|+C 0. De este modo, todas las soluciones de la ecuación
x
lineal no homogénea están dadas por las fórmulas

y=e ( ln |x|+C ) ; x=0


z

Prob.89: ( x + y 2 ) dy= ydx , M (1 ,1)

Solución: La ecuación examinada no es lineal respecto a la variable y , pero sí lo es


respecto a x . Por eso, es apropiado considerar x como una función de y .
Asumiendo que dy ≠ 0 ( y=0 es una solución trivial), tenemos
2 dx
x+ y =y
dy
dx
La ecuación homogénea asociada y tiene la solución y general x=Cy . Aplicando el
dy
método de variación de la constante, obtenemos sucesivamente

Cy+ y 2= y ( C' y+ C ) , C ' =1, C= y +C 0

Por consiguiente, todas las soluciones de la ecuación están dadas por las fórmulas
2
x=Cy + y , y=0

Prob.90: ( 2 e y −x ) y ' =1

Solución: La ecuación planteada es lineal respecto a x y, puesto que


dy 1
=
dx dx /dy
se puede escribir en la forma
y '
2 e −x=x … (1)
La solución de la ecuación lineal homogénea x ' + x=0 es
−y
x=C e … (2)

Considerando C=C ( y ) y sustituyendo (2) en (1), obtenemos sucesivamente


−y

2 e y −C e− y =C' e −C e− y
' y
C =2 e
2y
C ( y ) =e +C 0
Finalmente,
−y y
x=C e + e

Prob.91: ( sen2 y + xctgy ) y ' =1

Solución: Esta ecuación es lineal respecto a la variable x, por eso la representamos en la


forma
' 2
x −xctgy=sen y
Aplicando el método de variación de la constante, obtenemos
x ( y )=C ( y ) seny
Donde
C ( y ) =−cosy + const
1 dy −2 x x /22

Prob.92: + ( 2−x ) lny=x (e +e )


y dx
Solución: La ecuación propuesta es de la forma (3), p.4.3. Por eso, el cambio de variable
lny=z (x) produce
' 1 dy
z=
y dx
La ecuación resultante es lineal respecto a z. Utilizando método de variación de la
constante, hallamos
2
x
−2 x
z ( x ) =C ( x ) e 2

Donde

C ( x )=∫ x ( e ) dx +C
2
x
−2 x
2x
2
+e 0

De este modo, la solución general de la ecuación inicial es

(( ) )
2 2
x −x
−2 x 1 1
lny=e 2
x− e 2 x −e 2
+C
2 2

y dy
+ √ y +1=x +1
2 2
Prob.93:
√ y −1
2 dx

Solución: Esta ecuación es de la forma (3) examinada en el p.4.3. Por consiguiente, el


cambio de variable z (x)=√ y 2+1 conduce a la ecuación:

' y' ' 2


z= , z + z=x +1
√ y −1 '
2

Aplicando el método de variación de la constante, hallamos


−z
z=C ( x ) e
Donde
C ( x )=e z ( x2 −2 x +3 ) +C 0
Así
√ y 2+ 1=x 2−2 x+3+ C e− z
es la solución requerida.
1 dy 1 y
Prob.94: − x =e
e
x
dx e
Solución: Multiplicando los dos miembros de la ecuación por e x , obtenemos una ecuación
de la forma (4), p.4.3:
dy x y
−1=e e
dx
Vemos ahora que es apropiado hacer el cambio de variable z ( x)=e− y , a partir del cual
obtenemos sucesivamente
' −y
z ( x )=−e y ' ( x )
−y ' x y
−e z −1=e e
−z 1 x
La ecuación obtenida es lineal respecto a z . Su solución es z=C e − e . De este modo, la
2
solución general de la ecuación inicial es
−y −x 1 x
−e =C e − e
2

Prob.95: 3 dy + ( 1+e x+3 y ) dx=0

Solución: Transformado esta ecuación en


dy x 3y
3 +1=−e ⋅ e
dx
Obtenemos una ecuación de la forma (4), p. 4.3, lo que sugiere el cambio de variable
3y
z ( x ) =−e . De este modo, obtenemos

z ' ( x )=−3 e−3 y y '


' x
−z −e
+1=
z z
' x
z −z=e
cuya solución es
x x
z ( x ) =C e + x e
La solución general de la ecuación diferencial examinada es

Prob.96:
√lny dy + 2
√ (lny)3=1
y dx 3(x +1)
Solución: La ecuación dada es de la forma (3), p.4.3, puesto que
d 2
(
dx 3
√ (lny)3 =√ lny )
dy
dx
2
Por tanto, haciendo el cambio de variable z ( x ) =
3
√ 3
(lny) obtenemos la ecuación lineal

' z
z+ =1
x +1
Su solución general es
1 C
z= ( x +1 ) +
2 x +1
La solución general de la ecuación inicial es
2 1 C
3
√ ( lny )3= ( x +1 ) +
2 x+1

Prob.97: ( x +1 ) ( y ' + y 2 ) =− y

Solución. Considerando que x ≠ 1, dividimos los dos miembros de la ecuación por x +1, lo
que conduce a la ecuación de Bernoulli
' y 2
y+ =− y
x +1
Dividamos ambos miembros de esta ecuación por y 2 y hagamos luego el cambio de
variable y −1=z (x ). Obtenemos
' −2
z ( x )=− y y '
' z
z− =1
x +1
La ecuación lineal obtenida se resuelve por el método de variación de la constante.
Tenemos:
z=C ( x ) ( x+ 1 ) , donde C ( x )=ln |x +1|+C 0

Concluyendo, la solución de la ecuación no lineal dada es


1
y=
( x+1)(ln|x +1|+C )

Prob.98: xydx + ( x 2 + y 2+ 1 ) dy=0

Solución: Utilizando el cambio de variable x 2=u( y), llegamos a la ecuación lineal de


primer orden
1 dy 2
y +u=−( y +1)
2 dy
C( y)
Su solución general tiene la forma u= , donde
y2

( )
2
2 y
C ( y ) =− y +1 + const .
2

Así pues, las soluciones de la ecuación inicial son


4 2 2 2
y + 2 x y +2 y =C

Prob.99: ( x 3 y−3 x 2 y + y 3 ) dx +2 x3 dy=0

Solución: Dividiendo los dos miembros de la ecuación por dx ≠ 0 (la solución x=0 es
trivial), obtenemos la ecuación de Bernoulli
dy
+ ( x −3 x ) y=− y
3 3 2 3
2x
dx
Considerando que y ≠ 0 (evidentemente, la solución y=0 es trivial), dividimos los dos
1
miembros de la ecuación por − y 3 ¿ y hacemos el cambio de variable 2 =z ( x ). Entonces
y
obtenemos
'
−2 y
3
=z ' (x )
y

x 3 z ' −( x 3−3 x 2 ) z=1


Resolviendo esta ecuación lineal, hallamos
−3 x
z=C ( x ) x e
Donde
−x
C ( x )=−e +C 0

Ahora podemos escribir todas las soluciones de la ecuación diferencial examinada:


2 x 2 3
C y e − y −x =0 , x=0 , y=0
' x xy
Prob.100: 2 y − = 2
y x −1

Solución: Multiplicando los dos miembros de esta ecuación n por y y haciendo el cambio
de variable y 2=u ( x), obtenemos la ecuación lineal
' x
u− 2
u=x … (1)
x −1
Buscamos la solución en la forma u=f (x )w(x ). Sustituyendo u y u' en (1), hallamos

( f−
' xf
2
x −1
'
)
w+w f =x

Las funciones f y w se obtienen a partir de las ecuaciones


' xf '
f− 2
=0 , w f =x
x −1
De la primera ecuación resulta que

f =C √|x −1|
2

y de la segunda ecuación se deduce que


1
w=
C
√|x 2−1|sgn ( x 2−1 ) +C 0
Por consiguiente,

u=|x −1|sgn ( x −1 ) +C √| x −1|=x −1+C √ x −1


2 2 2 2 2

de donde

y =x −1+C √|x −1|


2 2 2

3
Prob.101: y ' x seny =x y ' −2 y

Solución: Dividiendo los dos miembros de la ecuación por y ≠ 0 ( y=0 es una solución
evidente) y considerando x como una función de y, obtenemos la ecuación de Bernoulli
dx 3
2y −x=−x seny
dy

Mediante el cambio de variable x−2=z ( y ) reducimos esta ecuación a la ecuación lineal


'
y z + z=seny
C cosy
cuya solución general es z= − . De este modo, el conjunto de soluciones de la
y y
ecuación diferencial dada es
1 C cosy
y=0 , = −
x
2
y y

O bien
2 2
y + x cosy−C x =0

Prob.102: ( x 2 + y 2+ 2 x−2 y ) dx +2 ( y−1 ) dy=0

Solución: Transformemos la ecuación dada en la forma

( ( x +1 )2 + ( y−1 )2−2 ) d ( x +1 ) +dy ( y−1 )2=0


Haciendo los cambios de variable x +1=u , ( y−1 )2=v , llegamos a la ecuación lineal

Su solución general es v=C e−u−u 2+ 2u . Así pues, las soluciones de la ecuación inicial están
dadas mediante la familia siguiente:
2 2 −x
x + y −2 y=C e

Prob.103: ( e y − y' ) =2

Solución: La ecuación dada se transforma en una ecuación de Bernoulli mediante el


cambio de variable e y =z ( x ):
' 2 2
z + z=z
x
1
Su solución general tiene la forma z ( x ) = . Consiguientemente, la solución general
x (1+Cx)
de la ecuación inicial es
2
y=−ln ⁡( x +C x )
z

Prob.104: y ( x )=∫ y (t)dt + x+1


0

Solución: Derivando los dos miembros de la ecuación, obtenemos la ecuación lineal


'
y = y +1
cuya solución general es y=C e x −1. A partir de la condición inicial evidente y (0)=1,
hallamos C=2, Por consiguiente, y=2e x −1 es la solución de la ecuación planteada
z z

Prob.105: ∫ ( x−t ) y (t ) dt=2 x+∫ y ( t ) dt


0 0

Solución: Derivando dos veces ambos miembros de la ecuación dada, obtenemos


x

∫ y ( t ) dt=2+ y ( x ) , y ( x )= y ' ( x)
0

de donde hallamos y (0)=−2 e y (x )=C e x. De ia condición inicial se deduce que C=2. Así,
la función y (x )=−2 e x es la solución de la ecuación planteada.

Prob.106: ydx + xdy + y 2 ( xdy− ydx )=0


Solución: Ésta es una ecuación de Minding—Darboux. En este caso, las funciones
M (x , y )= y y N (x , y)= x son homogéneas y de grado 1. La función R(x , y)= y 2 también es
homogénea y su grado es 2. Por ello, el cambio de variable y=ux (u) nos conduce a

uxdx + x ( udx + xdu )+u x ( x ( udx+ xdu )−uxdx ) =0


2 2

O bien
2 udx+ x ( 1+ x 2 u2 ) du=0 , x=0 …(1)
Dividamos ambos miembros de la ecuación obtenida por du . Entonces resulta la ecuación
de Bernoulli
dx 2 3
2 u + x=−u x
du
Mediante el cambio de variable x=z , reducimos la ecuación de Bernoulli a la ecuación
lineal
' 2
u z =u
No es difícil comprobar que su solución general es z=u2 +Cu . Regresando a las variables x
e y , finalmente obtenemos
2
y +Cxy−1=0
Las soluciones x=0 e y=0son casos particulares de la fórmula general para C=∞ .

Prob.107: ( x 2 y+ y 3−xy ) dx + x 2 dy=0

Solución: Escribiendo esta ecuación en la forma

0 ⋅ y 3 dy + ( x 2 y + y 3 ) dy + x ( xdy− ydx )=0


vemos que es una ecuación de Minding—Darboux. Por tanto, hagamos el mismo cambio
que utilizamos antes y=ux . Tenemos:

x 3 ( u+ u3 ) dx + x 3 du=0
O bien
du
x=0 , 2
+dx =0 ,u=0
u(1+u )
u −x
De aquí se infiere que =C e . Escribiendo la solución en términos de las variables
√1+u 2

iniciales, resulta

y 2=C 2 e−2 x ( x 2 + y 2 )
x=0

Prob.108: y 2 ( x+ a ) dx+ x ( x +a )

Solución: Ésta es una ecuación de Minding—Darboux, pues puede ser representada en la


forma
2 3
y xdx + x dy +ay ( ydx−xdy ) =0
Haciendo el cambio de variable y=ux (u), obtenemos la ecuación diferencial lineal
dx x a
+ =
du u ( u+1 ) u+1

cuya solución se obtiene mediante el método de variación de la constante. La ecuación


homogénea asociada tiene la solución
u+1
x= C …(1)
u
Considerando C=C (u) como una función incógnita de u, obtenemos la ecuación
diferencial
dC
du
=a
1
(

1
u+1 ( u+1 )2 )
que tiene por solución
a
C ( u )=aln ( C0 ( u+1 ) ) +
u+1
C 0=const .

y
Sustituyendo la expresión de C (u) en (1) y teniendo en cuenta que u= , escribimos la
x
solución general en la forma
x ( a− y )
x a ( x+ y )
=C e
x+ y

Prob.109: ( 2 xy−x 2 y− y 3 ) dx=( x 2 + y 2−x 3−x y 2 ) dy

Solución: La ecuación dada es una ecuación de Minding— Darboux. Para cerciorarnos de


ello, escribámosla de la manera siguiente:

( 2 xy−x 2 y− y 3 ) dx−( x 2 + y 2−x 3−x y 2 ) dy=2 xydx−( x 2 + y 2 ) dy−( x 2 + y 2 ) ( ydx−xdy )


Haciendo y=ux (u) , llegamos a la ecuación de variables separables siguiente:

( u−u3 ) dx + x ( 1+u 2) ( x−1 ) du=0


de la cual se deduce que
x−1 u
=C
x 1−u2

Por consiguiente, en las variables "antiguas" tenemos


2 2
y ( x−1 )=C ( x − y )
Prob.110: Hallar la curva que tiene la propiedad de que el segmento del eje Ox
comprendido entre el origen de coordenadas y el punto de intersección con la tangente a
la curva en cualquier punto, es proporcional a la ordenada de ese punto.
Solución: De la ecuación
Y − y= y '( X−x)
de la tangente a la curva buscada en el punto M(x,y), donde X e Y son las coordenadas de
un punto de la tangente, se infiere que la abscisa del punto de intersección y de la tangente
y
con el eje Ox es igual a x= . De este modo, según las condiciones de partida tenemos la
y'
ecuación.
y
x− '
=ky
y
O bien
dx
y −x =−ky
dy
Todas las soluciones de la ecuación tienen la forma
x= y (C−kln| y|)

Prob.111: Hallar las curvas para las cuales el área del trapecio formado por los ejes
coordenados, la tangente a la curva en un punto dado y la ordenada del punto de
tangencia es una magnitud contante e igual a 3 a2 .
Solución: De la ecuación de la tangente (v. ej. 91) hallamos la longitud del segmento
OK :Ok= y−xy ’ (fig. 21). Sea S el área del trapecio KONM. Tenemos:
|KO|+|MN |
S= |ON|
2
De acuerdo con las condiciones de partida, S=3 a2; luego
1
3 a = ( y−x y + y ) x
2 '
2
La ecuación que obtuvimos es lineal respecto a y :
2
' −6 a
x y −2 y =
x
2
2a 2
Su solución general es y= +C x
x
Prob.112: Hallar una curva tal que en todo punto M suyo la longitud del segmento NL
(fig.22) sea igual a la media aritmética de los cuadrados de las coordenadas de dicho
punto.
Solución: Según las condiciones de partida (fig. 22), tenemos
1
|NL|+ (|ON | +|MN | )
2 2

2
Veamos el triángulo MNL y hallemos la longitud del cateto NL. Tenemos ¿ NX∨¿ yy ' . De
este modo, la ecuación diferencial de las curvas buscadas es
' 2 2
2 y y =x + y
Haciendo el cambio de variable y 2=u , llegamos a Ja ecuación lineal u '−u=x 2, cuya
integración proporciona u=C e x −x 2−2 x−2. Finalmente, obtenemos y 2=C e x −x 2−2 x−2.
Prob.113: Un tanque contiene 100l de una disolución de sal. La cantidad de sal en la
disolución es igual a 10kg. En un tanque entran 5 litros de agua por minuto y con esa
misma velocidad la mezcla pasa a otro tanque de 100l de capacidad, el cual inicialmente
estaba lleno de agua. El excedente de líquido se derrama del segundo tanque ¿Cuándo la
cantidad de sal en el segundo tanque alcanzará su valor máximo? ¿A qué igual esa
cantidad?

Solución: Sean Q1 ( t ) y Q 2(t), respectivamente, las cantidades de sal (en kilogramos) en el


primer y segundo tanques en cierto instantet desde el inicio del transvase. Entonces
5Q 1 ( t 11 ) ∆ t
es la cantidad de sal que pasa del primer tanque al segundo en el intervalo de
100
5Q 2 ( t 12) ∆t
tiempo desde t hasta t+ ∆ t , mientras que es la cantidad de sal que se derrama
100
del segundo tanque en ese mismo intervalo de tiempo, donde t 11 ∈ ( t , t+ ∆ t ) , t 12 ∈ ( t , t+ ∆t ).
Por consiguiente
5 Q1 ( t 11 ) 5 Q2 ( t 11 )
Q2 ( t+ ∆ t )=Q 2 ( t )+ ∆ t− ∆ t … (1)
100 100
es la cantidad de sal en el segundo tanque en el instante t+ ∆ t , y
5Q1 ( t 11 )
Q1 ( t+ ∆ t )=Q 1 ( t )− ∆ t …(2)
100
es la cantidad de la sal en el primer tanque en ese mismo instante. Pasando al límite en (2)
cuando ∆ t →0 , obtenemos la ecuación diferencial
d Q1
=−0 ,05 Q 1
dt

de donde después de integrar hallamos Q1=C e−0 ,05 t , asumiendo que el tiempo t se mide
en minutos. Puesto que Q1 ( 0 )=10 , entonces C=10. Así tenemos que
−0 ,05 t
Q1=C e …(3)

Pasando al límite en (1) cuando 0 y teniendo en cuenta (3), se obtiene


d Q2 −0 ,05 t
=−0 ,05 Q 2 +0 , 5 e
dt
Resolviendo esta ecuación lineal, obtenemos
−0 , 05t
Q2 ( t ) =(0 , 5 t+C )e

Puesto que en el momento inicial Q2 (0)=0, entonces C=0 . Finalmente, hallamos


−0 , 05 t
Q2 ( t ) =0 ,5 e

Haciendo un análisis de la función Q2, vemos que ésta alcanza su valor máximo para
t=20 min, donde
10
Q 2 ( 20 )= ≈ 3 , 68 kg
e
Prob.114: Al cabo de un intervalo de tiempo ∆ t (donde ∆ t →0 ) de cada gramo de radio se
desintegra 0,00044 ∆ t gramos y se forma 0,00043 ∆ t gramos de radón. A su vez, de cado
gramo de radón al cabo del tiempo ∆ t se desintegran 70 ∆ t gramos. A l comenzar el
experimento se disponía de cierta cantidad x 0 de radio. Determinar cuándo alcanzará su
máximo la cantidad de radón formado y aún no desintegrado. Asumir que el tiempo se
mide en años.
Solución: Sean P(t ) y Q(t ) las cantidades respectivas de radio y radón no desintegrados en
el momento t desde el inicio de la desintegración. Entonces, P ( t ) −P(t+ ∆ t) es la cantidad
de radio desintegrado en el intervalo de tiempo desde t hasta t+ ∆ t , mientras que
Q ( t+ ∆ t ) −Q ( t ) es la cantidad de radón formado durante ese mismo intervalo de tiempo.
Según las condiciones de partida, tenemos las ecuaciones.
P ( t ) −P ( t +∆ t )=P ( t 11 ) ⋅0,00044 ∆ t …(1)

Q ( t+ ∆ t ) −Q ( t )=P ( t 11 ) ⋅0,00043 ∆ t−Q ( t 12) 70 ∆ t …(2)

donde t 11 ∈ ( t , t+ ∆ t ) , t 12 ∈ ( t , t+ ∆t ) . Pasando al límite cuando ∆ t →0 (dividiendo


previamente por ∆ t los dos miembros de las ecuaciones (1) y (2)), obtenemos las
ecuaciones diferenciales
dP
=−0,00044 P ( t ) …(3)
dt
dQ
=0,00044 P ( t )−70 Q (t ) …( 4)
dt
La solución de la ecuación (3) es
−0,00044 t
P ( t ) =x 0 e …(5)

Sustituyendo (5) en (4) e integrando, hallamos

−70 t 0,00043 x 0 −0,00044t


Q ( t ) =C e + e
69,99956
Teniendo en cuenta la condición inicial Q(0)=0, determinamos la constante
−0,00043 x 0
C= . Finalmente tenemos
69,99956
0,00043 x 0 −0,00044 t −70 t
Q(t )= (e −e )
69,99956
Examinando la función f ( t )=e−0,00044t −e−70 t, hallamos que f ( t ) alcanza su máximo para
1 70
t= ln ≈ 0 , 17 años ≈ 62 días
69,99956 0,00044
Prob.115: Dadas dos soluciones diferentes y 1 e y 2 de una ecuación lineal de primer orden
expresar mediante ellas la solución general de la ecuación.
Solución: Para toda ecuación diferencial lineal de primer orden
'
y + P ( x ) y=Q( x )
su solución general tiene la forma
y= ( C+ α ( x )) β ( x ) … (1)

Donde α ( x ) =∫ Q ( x ) β ( x ) dx , β ( x ) =e ∮ P (x ) dx
−1

Según las condiciones de partida, haciendo uso de (1) tenemos


y 1 ( x ) = ( C 1+ α ( x ) ) β ( x )

y 2 ( x ) =( C 2+ α ( x ) ) β ( x ) …(2)
donde C 1 , C 2 son las constantes correspondientes a las soluciones y 1 , y 2 respectivamente.
Ahora, utilizando las igualdades (2), expresamos las funciones α , β mediante y 1 , y 2
obteniendo de este modo
y 1 ( x )− y 2 (x)
α (x )= (C 1 ≠ C2 )
C1−C 2

C 1 y 1 ( x )−C 2 y 2 (x)
β (x )= ( y 1 ( x )− y 2(x ))
y 1− y 2

Finalmente, sustituyendo las expresiones de α ( x ) y β ( x ) en (1), hallamos


1 ~
y=
C 1−C2
( ( C 1−C ) y 1 ( x )+ ( C−C2 ) y 2 (x))= y 2 ( x ) + C ( y 2 ( x )− y 1(x ))

~ C2−C
Siendo C= una constante arbitraria.
C 1−C 2

Prob.116: Dada la ecuación


'
y sen 2 x=2 ( y+ cosx )
π
Hallar aquella solución que permanece acotada cuando x →
2
Solución: La solución general de esta ecuación se expresa mediante Ja fórmula
1
y=Ctgx −
cosx
de la cual se deduce que

lim y (x)= lim Ctgx−


x→ π / 2 x →π/2
( 1
cosx )
es finito (igual a cero) sólo para C=1. Por consiguiente, y=tgx−secx es la solución
requerida

Prob.117: Demostrar que la ecuación x y' +ay =f ( x) (a=const >0), donde f es una función
continua tal que f → b cuando x → 0, tiene una solución única acotada cuando x → 0 . Hallar
dicha cota
Solución: La solución general de la ecuación dada se expresa mediante la fórmula
x
C 1
y= a + a ∫ f (t)|t| d (|t |)
a−1

|x| | x| 0
¿, o bien
x
C
b 1
y= a + + a ∫ ε (t )|t | d (|t|) …(1)
a −1

|x| a |x| 0
donde ε (t)→b cuando t → 0 de acuerdo con las condiciones del ejercicio. Según la
estimación

|∫ |
x
1 a −1 1
ε (t )|t | d (|t|) ≤ 0 ≤ t ≤ z|ε ( t )|→ 0 , x →0
|x|
a
0 a¿
a partir de (1) se deduce que lim y existe, está acotado sólo para C=0 y es igual a b . De
x →0 a
este modo, la solución requerida tiene la forma
x
1
∫ f (t )|t|
a −1
y= a
d (|t|)
|x| 0

Prob.118: Supongamos que en la ecuación diferencial del problema anterior a=const <0 y
f (x)→ b cuando x → 0 . Demostrar que todas las soluciones de esta ecuación tienen un
mismo límite finito cuando x → 0, la solución general de la ecuación analizada se puede
escribir en la forma

y=| x|
−a
( C +∫ f (x)|x|a−1 d (|x|))= ba +|x|−a ( C+∫ ε ( x )|x|a−1 d (|x|) )
Si la integral ∫ ε (x)|x|a−1 d (|x|) está acotada, entonces paro b todo C se tiene,
b
evidentemente, que lim y ( x )= . Si la integral no está acotada cuando x → 0 , entonces
x →0 a
aplicamos la regla de L’Ho^ pital, obteniendo

lim
∫ a−1
ε (x)|x| d (| x|)
=lim
( x )| x|
a−1

=0
a a−1
x →0 |x| x→ 0 a|x|
b
De esta manera, para todos los valores de C se tiene que lim y ( x )=
x →0 a
dx
Prob.119: Demostrar que la ecuación + x=f (t ), donde f es una función continua y
dt
|f (t)|≤ M para −∞ <t <+∞ , tiene una única solución acotada para −∞ <t <+∞ . Hallar dicha
solución y demostrar que ésta es periódica si la función f es periódica
Solución: La solución general de la ecuación dada es
t
−t
x (t )=C e +e
−t
∫ f ( τ ) eτ dτ … (1)
−∞

Tal representación es posible debido a que la integral impropia


t

∫ f ( τ ) e τ dτ
−∞
converge en virtud de la estimación

|∫ |
t
τ τ
f ( τ ) e dτ ≤ M e …(2)
−∞

De la desigualdad (2) se deduce, además, que la función


t
e
τ
∫ f ( τ ) e τ dτ
−∞

está acotada por el número M para todo t ϵ(−∞ ,+ ∞). De esta manera, la condición
necesaria (y suficiente) para que la función x sea acotada es C=0. La solución requerida
tiene la forma
t
x (t )=e
−τ
∫ f ( τ ) e τ dτ …(3)
−∞

Ahora bien, supongamos que ∀ t ϵ (−∞ ,+∞ ) t ( τ +T )=f ( τ ), donde T > 0. Entonces, a partir de
(3) hallamos
t t +T
−τ
x (t )=e ∫ f ( τ ) e dτ=x ( t ) =e
τ −(t +T )
∫ f ( τ 1 ) eτ d τ 1=x (t +T )
1

−∞ −∞

Donde τ 1 =τ +T . Por consiguiente, x es una función periódica

Nota: La condición de que f sea continua no es necesaria. Dejamos a cargo del lector
determinar la clase de funciones f para las cuales se cumple la afirmación anterior

Prob.120: Demostrar que solamente una solución de la ecuación x y' − ( 2 x 2 +1 ) y=x 2


converge a un límite finito cuando x →+ ∞. Hallar ese límite. Escribir dicha solución en
forma integral
Solución: Partimos de la solución general de la ecuación dada:

y=x e ( C+∫ e dx ) …(1)


2 −2
x x

∫ e−t dt , la suma entre paréntesis de


2

En virtud de la convergencia de la integral impropia


−∞
(1) se puede representar en la forma
x
C+∫ e dx=C 1+ ∫ e dt …(2)
−2 2
x −t

−∞

donde C 1 es cierta constante. (La igualdad (2) se comprueba fácilmente mediante


diferenciación.) De este modo, todas las soluciones de la ecuación examinada se expresan
mediante la fórmula

( )
x
x2
∫ e−t dt
2

y=x e …( 3)
−∞

Sea x →+ ∞. Entonces, a partir de (3) se deduce que y está acotada cuando x →+ ∞ sólo si
se cumple la condición
+∞
C 1=−∫ e dt=−√ π
2
−t

−∞

La última condición resulta suficiente para que lim


x →0
y también sea finito. Efectivamente,
según la regla de L'Ho^ pital tenemos
x

∫ e−t dt−√ π
2

−∞ e− x −1
lim =lim =
x→ 0 −( x +2 ) e 2
−1 −x 2 −2 −x 2
x →0 x e

Para concluir, escribamos la solución requerida:


(∫ )
x x
e dt−√ π =x ∫ e
x2 −t 2 x 2−t 2
y=x e dt
−∞ −∞

Prob.121: Hallar la solución periódica de la ecuación


' 2
y =2 y cos x−senx
Solución: En vista de la convergencia de la integral impropia
+∞

∫ exp {−t−sentcost } sentdt


−∞

la solución general de la ecuación dada se escribe en la forma


+∞
y=Cexp { x + senxcosx }−∫ exp {−t−sentcost } sentdt …(1)
x

Puesto que, evidentemente, la función exp { x +senxcosx } no es periódica y la función


+∞ +∞
y 1 ( x ) =∫ exp {−t + x−sentcost+ senxcosx } sentdt=∫ exp {−s−senscos ( s+2 x ) } sen (x+ s ¿)ds ¿
x 0

es periódica con período 2 π (esto se deduce de la identidad)


+∞
y1 ( x ) = ∫ exp {−t 1 +( x+ 2 π )−sen t 1 cos t 1+ sen ( x+ 2 π ) cos(x +2 π ) } sen t 1 d t 1 ≡ y 1 (x +2 π )
x+2 π

Donde t 1=t+ 2 π , entonces la función y es periódica tan sólo para C=0. Tomando C=0 en
(1), obtenemos la solución periódica y= y 1(x ).
Prob.122: Consideremos la ecuación
dx
+ a ( t ) x=f (t)
dt
Donde a (t)≥0 , f (t)→0 cuando t → ∞ , y las funciones a , f son continuas para t >t 0 .
Demostrar que toda solución de esta ecuación tiende a cero cuando t →+ ∞.
Solución: Haciendo uso de la solución general

x (t )=( C +∫ f ( t ) exp {∫ a ( t ) dt } dt ) exp {−∫ a ( t ) dt }

de la ecuación dada, de las condiciones de partida, y aplicando la regla de L'H o^ pital,


obtenemos

lim C+∫ f ( t ) exp {∫ a ( t ) dt } dt f ( t ) exp {∫ a ( t ) dt }


lim x ( t )=¿ t →+∞ = lim =0 ¿
t →+ ∞ exp {∫ a ( t ) dt } t →+∞ a ( t ) exp
{∫ a ( t ) dt }
Nótese que la continuidad de las funciones a y f garantiza la diferenciabilidad de las
integrales correspondientes. La condición a ( t ) ≥ C >0 se utiliza dos veces:

e∫
a ( t ) dt Ct
> e →+∞
Y

| ||f (t) f ( t)|


a( t)

C
→0
Cuando t →+ ∞.
Prob.123: Supongamos que en la ecuación del problema anterior a ( t ) ≥ C >0 y que x 0 (t) es
su solución para la condición inicial x 0 ( 0 )=b . Demostrar que ∀ ε > 0∃ δ> 0 tal que variamos
la función f y el número b en un valor menor que δ (es decir, f y b se cambian por una
función f 1 y un número b 1 tales que |f 1 ( t )−f (t )|< δ ), entonces la solución x 0 (t) varía para
t ≥ 0 en un valor menor que ε . Dicha propiedad de la solución se denomina estabilidad
respecto a las perturbaciones de acción constante.
Solución: Sea t 0=0. Entonces la solución general de la ecuación planteada se puede
escribir en la forma siguiente:

( { }) { }
t τ t

x (t )= C+∫ f ( τ ) exp ∫ a ( ξ ) dξ dτ × exp −∫ a ( τ ) dτ … (1)


0 0 0

Teniendo en cuenta las condiciones de partida, de (1) hallamos

( {∫ } ) { ∫ }
t τ t

x 0 (t)= b +∫ f ( τ ) exp a ( ξ ) dξ dτ × exp − a ( τ ) dτ …(2)


0 0 0

( {∫ } ) { ∫ }
t τ t

x 1 (t)= b1 +∫ f 1 ( τ ) exp a ( ξ ) dξ dτ ×exp − a ( τ ) dτ …(3)


0 0 0

Restando miembro a miembro de (2) la igualdad (3), obtenemos la estimación

{ } { }
t t t

|x 0 ( t ) −x1 ( t )|≤|b−b 1|exp −∫ a ( τ ) dτ +∫|f ( τ ) −f 1 ( τ )|exp −∫ a ( ξ ) dξ …(4)


0 0 τ

Dado que

|b−b 1|< δ
f ( τ )−f 1 ( τ ) <δ

{ }
t
exp −∫ a ( τ ) dτ ≤ e
−Ct

{ }
t
exp −∫ a ( ξ ) dξ ≤ e
−C(t−τ )

a partir de (4) resulta la estimación

|x 0 ( t ) −x1 ( t )|≤ δ (e −Ct


+
1
C
( 1−e−Ct ) < δ 1+ 1
C) ( )
εC
De este modo, si para un ε > 0 dado elegimos δ= , entonces para t ≥ 0 se cumple la
1+C
desigualdad

|x 0 ( t ) −x1 ( t )|<ε
la cual significa, por definición, que la solución x 0 (t) es estable respecto a las
perturbaciones de acción constante.
5) Ecuaciones diferenciales exactas. Factor integrante
Hallar las integrales generales de las ecuaciones siguientes:

Prob.124: ( xlny−x 2+ cosy ) dy + ( x 3 + ylny− y−2 xy ) dx=0

Solución: Puesto que en el conjunto

D= {( x ; y ) ∈ R2 :−∞< x< +∞ , y 0 }

Se cumple la igualdad
∂ 3
( x + ylny− y−2 xy )= ∂ ( xlny−x 2+ cosy ) =lny−2 x
∂y ∂y
entonces el primer miembro de la ecuación examinada es la diferencial total de cierta
función Φ Utilizando la fórmula (3), p.5.1, obtenemos que
x y 4
x
Φ ( x , y )=∫ ( t + ylny− y−2 ty ) +∫ costdt = + xy ( lny−1 )− y x + seny−sen y 0
3 2

0 y0 x

De este modo, la integral general de la ecuación dada escribe en la forma


4 2
x + 4 xy (lny−1 )−4 x y+ 4 seny=C

(
Prob.125: x +
1
√ y +x
2 2 ) (
dx + y−
1
y √ y 2−x 2
dy=0
)
Solución: Dado que para | y|>|x| se cumple la identidad


∂y ( 1
x+ 2 2 =

√ y −x ∂ y
y−
1
) (
y √ y −x
2 2
2 −3/ 2
=− y ( y −x )
2
)
entonces la ecuación dada es una ecuación diferencial exacta. Aplicando la fórmula (3), p.
5.1, obtenemos
x y
Φ ( x , y )=∫
0
( 1
) 1 x 1
t+ 2 2 dt +∫ tdt= ( y 2 + x 2 ) +arcsen − y 02
√ y −t y 2 0
y 2

la integral general de la ecuación analizada es


2 2 x
x + y + 2 arcsen =0
y
Resolver las ecuaciones diferenciales siguientes con ayuda de factores integrantes de la
forma μ=f (x ), o bien μ=f ( y ):

( xy ) dx+( 1x + 2xy ) dy =0
Prob.126: 1+ 2 2

y 1 2y
Solución: Tomando ω=x , M =1+ 2 y
N= + 2 en la fórmula (6), p.5.3, obtenemos
x x x

( 1 2 y dμ
+ 2
x x dx )
=μ 2 + 3 =
x x
+
(
2 4 y 2μ 1 2 y
x x x2 ) ( )
de donde hallamos que μ=x 2 (C=1). Vemos pues que la elección de la función CJ resultó
apropiada. Multiplicando ambos miembros de la ecuación dada por x 2, obtenemos la
ecuación diferencial exacta

( x 2 + y ) dx + ( x+ 2 y ) dy=0
Aplicando la fórmula (3), p.5.2, hallamos la integral general
3 2
x +3 xy + 3 y =C
Prob.127: y 2 ( x−3 y ) dx +(1+3 x y 2 )dy =0

Solución: Hagamos ω=x , M = y 2 (x−3 y ) N=1−3 x y 2 en la fórmula (6), p.5.3. Tenemos:

( 1−3 x y 2 ) dμ =2 y (x −3 y) μ
dx
1 dμ 2 y ( x−3 y)
=
μ dx 1−3 x y 2
de donde vemos que μ, no puede depender sólo de x , pues el primer miembro de la última
igualdad es una función sólo de x mientras que el segundo miembro es una función de x e
y (las variables x e y son independientes).
Probemos entonces el factor ω= y . Tenemos:
2 dμ
− y ( x−3 y ) =2 y (x−3 y)μ
dx

−y =2 μ
dx

Integrando la última ecuación hallamos μ= y−2 (C=1) .

Multiplicando los dos miembros de la ecuación inicial por y −2, llegamos a la ecuación
diferencial exacta siguiente:

( x−3 y ) dx + ( y−2−3 x ) dy=0


Haciendo uso de la fórmula (3), p.5.2, obtenemos su integral general:
2 2
x y −6 x y −2=Cy ,( y ≠ 0)

Obsérvese que la división de la ecuación inicial por y 2 condujo a la pérdida de la solución


y=0; por eso, la integral general de la ecuación examinada es
2 2
x y −6 x y −2 ≠ Cy ,( y=0 para C=∞)
Integrar las siguientes ecuaciones con ayuda de un factor integrante de la forma
μ=μ (x+ y) o bien μ=μ (x− y ):

Prob.128: ( 2 x 3+ 3 x 2 y+ y 2− y 3 ) dx + ( 2 y 3+ 3 x y 2 + x 2−x3 ) dy =0

Solución: Sea μ=μ (x− y ). Entonces tomando en la fórmula (6), p.5.3, w=x− y , obtenemos

( 2 y 3 +3 x y 2 + x 2−x 3 + ( 2 x3 +3 x 2 y+ y 2− y 3 ) ) dμ =( 6 x 2−6 y2 +2 y−2 x ) μ



Reuniendo términos semejantes

( y 3 −x3 + x 2 + y 2 +3 xy ( x + y ) ) =2 ( 3 x −3 y + y−x ) μ
2 2

Es evidente que la expresión
2 ( 3 x 2−3 y 2+ y−x )
y 3−x 3+ x2 + y 2 +3 xy (x + y )

no es una función de (x− y ); por tanto, buscaremos la función μ en la forma μ=μ (x+ y).
Análogamente obtenemos

(2 y 3 +3 x y 2 + x 2−x 3 + ( 2 x3 +3 x 2 y+ y 2− y 3 ) ) dμ =( 6 x 2−6 y2 +2 y−2 x ) μ



Factorizamos:

( 3 ( x+ y)2−( x− y ) ) =( 2−6 (x+ y) ) μ

O bien
( 3 ω 2−ω ) μ' =2(1−3 ω) μ
'
ω μ +2 μ=0
Resolviendo esta última ecuación obtenemos
1 1
μ= 2
= 2
(C=1)
ω (x+ y)
1
Multiplicando la ecuación inicial por el factor integrante , llegamos a una ecuación
( x+ y )2
diferencial exacta. Su integral general tiene la forma
x 3 2 2 3 y

∫ 2 t +3( y t +)2y − y dt +∫ 2 tdt=const ( y 0 ≠0)


0 t+ y y 0

De donde
3 3
x + y + xy=C ( x+ y)
Nótese que la solución y=− xes un caso particular de la integral general para C=∞ .

(
Prob.129: y −
ay
x )
+ x dx −ady=0

Solución: Busquemos el factor integrante en la forma μ=μ (x+ y). Tomando en la fórmula
(6), p. 5.3, w=x− y , obtenemos

( a− y + ayx −x) dωdμ =(1− ax ) μ


O bien
( ( a−x ) x − y ( x−a ) μ ' )=( x−a ) μ
( x + y ) μ ' =0
'
ω μ + μ=0
Al integrar la última ecuación resulta que
−1 −1
μ=ω =( x+ y ) ,(C=1)
1
De esta manera, multiplicando por el factor integrante reducimos la ecuación dada a
x+ y
la forma

(1− x (xay+ y ) ) dx+ x +a y dy =0


Integrándola, obtenemos
x y
dt
∫ dt +∫ x+ t
=const (x 0 ≠ 0)
x0 0

y de aquí,

| |
a
x y
e 1+ =C
x

Resolver las siguientes ecuaciones considerando que el factor integrante tiene la forma
μ=μ (xy) , μ=μ ( x 2+ y 2), o bien μ=μ ( x 2− y2 ):

Prob.130: ( x 2 + y ) dy + x ( 1− y ) dx=0

Solución: Probemos el factor μ=μ (xy) , tomando en la fórmula (6), p.5.3, w=xy . Tenemos:

(( x2 + y ) y− (1− y ) x2 ) dμ =(−x−2 x) μ

O bien
( 2 y x 2 + y 2−x 2 ) dμ +3 xμ=0

3x
Es evidente que 2 22 no se puede expresar en términos de xy , luego el factor
2 y x + y −x
integral en la forma elegida no existe.
Hagamos w=x 2 + y 2. Entonces obtenemos

( ( x2 + y ) 2 x−( 1− y ) 2 xy ) dμ =−3 x μ

De donde
2 ( x 2+ y 2) μ+3 μ=0
'
2 ω μ +3 μ=0
−3
Integrando la última ecuación hallamos que μ=ω 2 (C=1). Así pues, la multiplicación por
−3 /2
el factor( x 2 + y 2 ) transforma la ecuación inicial en la ecuación diferencial exacta

x +y
2
x (1− y )
dy + dx=0
2 3/2 3 /2
(x + y )
2
( x2 + y2 )
Integrándola, obtenemos
d (| y|)
x y
t ( 1− y ) dt
∫ 3/2 ∫
+ 2
=const ,( y 0 ≠ 0)
0 ( t 2+ y 2 ) y0 | y|
O bien

( 1− y )
( 1 1
− 2 2 − =C
1
| y| √ x + y | y| )
y 1− y
+ =C
| y| √ x2 + y2

También podría gustarte