Omitir e ir al contenidoIr a la página de accesibilidadMenú de atajos de teclado
Logo de OpenStax
Cálculo volumen 2

4.2 Campos de direcciones y métodos numéricos

Cálculo volumen 24.2 Campos de direcciones y métodos numéricos

Menú
Índice
  1. Prefacio
  2. 1 Integración
    1. Introducción
    2. 1.1 Aproximación de áreas
    3. 1.2 La integral definida
    4. 1.3 El teorema fundamental del cálculo
    5. 1.4 Fórmulas de integración y el teorema del cambio neto
    6. 1.5 Sustitución
    7. 1.6 Integrales con funciones exponenciales y logarítmicas
    8. 1.7 Integrales que resultan en funciones trigonométricas inversas
    9. Revisión del capítulo
      1. Términos clave
      2. Ecuaciones clave
      3. Conceptos clave
      4. Ejercicios de repaso
  3. 2 Aplicaciones de la integración
    1. Introducción
    2. 2.1 Áreas entre curvas
    3. 2.2 Determinar los volúmenes mediante el corte
    4. 2.3 Volúmenes de revolución: capas cilíndricas
    5. 2.4 Longitud del arco de una curva y superficie
    6. 2.5 Aplicaciones físicas
    7. 2.6 Momentos y centros de masa
    8. 2.7 Integrales, funciones exponenciales y logaritmos
    9. 2.8 Crecimiento y decaimiento exponencial
    10. 2.9 Cálculo de las funciones hiperbólicas
    11. Revisión del capítulo
      1. Términos clave
      2. Ecuaciones clave
      3. Conceptos clave
      4. Ejercicios de repaso
  4. 3 Técnicas de integración
    1. Introducción
    2. 3.1 Integración por partes
    3. 3.2 Integrales trigonométricas
    4. 3.3 Sustitución trigonométrica
    5. 3.4 Fracciones parciales
    6. 3.5 Otras estrategias de integración
    7. 3.6 Integración numérica
    8. 3.7 Integrales impropias
    9. Revisión del capítulo
      1. Términos clave
      2. Ecuaciones clave
      3. Conceptos clave
      4. Ejercicios de repaso
  5. 4 Introducción a las ecuaciones diferenciales
    1. Introducción
    2. 4.1 Fundamentos de las ecuaciones diferenciales
    3. 4.2 Campos de direcciones y métodos numéricos
    4. 4.3 Ecuaciones separables
    5. 4.4 La ecuación logística
    6. 4.5 Ecuaciones lineales de primer orden
    7. Revisión del capítulo
      1. Términos clave
      2. Ecuaciones clave
      3. Conceptos clave
      4. Ejercicios de repaso
  6. 5 Secuencias y series
    1. Introducción
    2. 5.1 Secuencias
    3. 5.2 Serie infinita
    4. 5.3 Las pruebas de divergencia e integral
    5. 5.4 Pruebas de comparación
    6. 5.5 Series alternadas
    7. 5.6 Criterios del cociente y la raíz
    8. Revisión del capítulo
      1. Términos clave
      2. Ecuaciones clave
      3. Conceptos clave
      4. Ejercicios de repaso
  7. 6 Serie de potencias
    1. Introducción
    2. 6.1 Series y funciones de potencia
    3. 6.2 Propiedades de las series de potencia
    4. 6.3 Series de Taylor y Maclaurin
    5. 6.4 Trabajar con la serie de Taylor
    6. Revisión del capítulo
      1. Términos clave
      2. Ecuaciones clave
      3. Conceptos clave
      4. Ejercicios de repaso
  8. 7 Ecuaciones paramétricas y coordenadas polares
    1. Introducción
    2. 7.1 Ecuaciones paramétricas
    3. 7.2 Cálculo de curvas paramétricas
    4. 7.3 Coordenadas polares
    5. 7.4 Área y longitud de arco en coordenadas polares
    6. 7.5 Secciones cónicas
    7. Revisión del capítulo
      1. Términos clave
      2. Ecuaciones clave
      3. Conceptos clave
      4. Ejercicios de repaso
  9. A Tabla de integrales
  10. B Tabla de derivadas
  11. C Repaso de Precálculo
  12. Clave de respuestas
    1. Capítulo 1
    2. Capítulo 2
    3. Capítulo 3
    4. Capítulo 4
    5. Capítulo 5
    6. Capítulo 6
    7. Capítulo 7
  13. Índice

Objetivos de aprendizaje

  • 4.2.1 Dibujar el campo de direcciones para una ecuación diferencial de primer orden dada.
  • 4.2.2 Utilizar un campo de direcciones para dibujar una curva de solución de una ecuación diferencial de primer orden.
  • 4.2.3 Utilizar el método de Euler para aproximar la solución de una ecuación diferencial de primer orden.

En el resto de este capítulo nos centraremos en diversos métodos para resolver ecuaciones diferenciales y analizar el comportamiento de las soluciones. En algunos casos es posible predecir las propiedades de una solución de una ecuación diferencial sin conocer la solución real. También estudiaremos métodos numéricos para la resolución de ecuaciones diferenciales, que pueden programarse utilizando diversos lenguajes informáticos o incluso utilizando un programa de hoja de cálculo, como Microsoft Excel.

Creación de campos de direcciones

Los campos de direcciones (también llamados campos de pendiente) son útiles para investigar las ecuaciones diferenciales de primer orden. En particular, consideramos una ecuación diferencial de primer orden de la forma

y=f(x,y).y=f(x,y).

Un ejemplo aplicado de este tipo de ecuación diferencial aparece en la ley de enfriamiento de Newton, que resolveremos explícitamente más adelante en este capítulo. Primero, sin embargo, vamos a crear un campo de direcciones para la ecuación diferencial

T(t)=−0,4(T72).T(t)=−0,4(T72).

Aquí T(t)T(t) representa la temperatura (en grados Fahrenheit) de un objeto en el tiempo t,t, y la temperatura ambiente es 72 °F.72 °F. La Figura 4.6 muestra el campo de direcciones para esta ecuación.

Un gráfico de un campo de direcciones para la ecuación diferencial dada en los cuadrantes uno y dos. Las flechas apuntan directamente a la derecha en y = 72. Por debajo de esa línea, las flechas tienen una pendiente cada vez más positiva a medida que y disminuye. Por encima de esa línea, las flechas tienen una pendiente cada vez más negativa a medida que y aumenta. Las flechas señalan la convergencia en y = 72. Se dibujan dos soluciones: una para la temperatura inicial inferior a 72, y otra para temperaturas iniciales superiores a 72. La solución superior es una curva cóncava decreciente hacia arriba, que se acerca a y = 72 a medida que t llega al infinito. La solución inferior es una curva cóncava creciente hacia abajo, que se acerca a y = 72 a medida que t llega al infinito.
Figura 4.6 Campo de direcciones para la ecuación diferencial T ( t ) = −0,4 ( T 72 ) . T ( t ) = −0,4 ( T 72 ) . Se representan dos soluciones: una con temperatura inicial inferior a 72  ° F 72  ° F y la otra con temperatura inicial superior a 72  ° F . 72  ° F .

La idea de un campo de direcciones es el hecho de que la derivada de una función evaluada en un punto determinado es la pendiente de la línea tangente al gráfico de esa función en el mismo punto. Otros ejemplos de ecuaciones diferenciales para las que podemos crear un campo de direcciones son

y=3x+2 y4y=x2 y2 y=2 x+4y2 .y=3x+2 y4y=x2 y2 y=2 x+4y2 .

Para crear un campo de dirección, empezamos con la primera ecuación: y=3x+2 y4.y=3x+2 y4. Suponemos que (x0,y0)(x0,y0) es cualquier par ordenado, y sustituimos estos números en el lado derecho de la ecuación diferencial. Por ejemplo, si elegimos x=1yy=2 ,x=1yy=2 , sustituyendo en el lado derecho de la ecuación diferencial se obtiene

y=3x+2 y4=3(1)+2 (2 )4=3.y=3x+2 y4=3(1)+2 (2 )4=3.

Esto nos dice que si una solución de la ecuación diferencial y=3x+2 y4y=3x+2 y4 pasa por el punto (1,2 ),(1,2 ), entonces la pendiente de la solución en ese punto debe ser igual a 3.3. Para empezar a crear el campo de direcciones, ponemos un segmento de línea corto en el punto (1,2 )(1,2 ) con pendiente 3.3. Podemos hacer esto para cualquier punto del dominio de la función f(x,y)=3x+2 y4,f(x,y)=3x+2 y4, que consiste en todos los pares ordenados (x,y)(x,y) en 2 .2 . Por lo tanto, cualquier punto del plano cartesiano tiene asociada una pendiente, suponiendo que por ese punto pasa una solución de la ecuación diferencial. El campo de direcciones para la ecuación diferencial y=3x+2 y4y=3x+2 y4 se muestra en la Figura 4.7.

Un gráfico del campo de direcciones para la ecuación diferencial y' = 3 x + 2 y - 4 en los cuatro cuadrantes. En los cuadrantes dos y tres, las flechas apuntan hacia abajo y ligeramente hacia la derecha. En una línea diagonal, aproximadamente y = -x + 2, las flechas apuntan cada vez más a la derecha, se curvan y luego apuntan hacia arriba por encima de esa línea.
Figura 4.7 Campo de direcciones para la ecuación diferencial y = 3 x + 2 y 4 . y = 3 x + 2 y 4 .

Podemos generar un campo de direcciones de este tipo para cualquier ecuación diferencial de la forma y=f(x,y).y=f(x,y).

Definición

Un campo de direcciones (campo de pendiente) es un objeto matemático utilizado para representar gráficamente las soluciones de una ecuación diferencial de primer orden. En cada punto de un campo de direcciones aparece un segmento de línea cuya pendiente es igual a la de la solución de la ecuación diferencial que pasa por ese punto.

Uso de los campos de direcciones

Podemos utilizar un campo de direcciones para predecir el comportamiento de las soluciones de una ecuación diferencial sin conocer la solución real. Por ejemplo, el campo de direcciones en la Figura 4.7 sirve de guía para el comportamiento de las soluciones de la ecuación diferencial y=3x+2 y4.y=3x+2 y4.

Para utilizar un campo de direcciones, empezamos eligiendo cualquier punto del campo. El segmento de línea en ese punto sirve como señal que nos indica la dirección que debemos tomar a partir de ahí. Por ejemplo, si una solución de la ecuación diferencial pasa por el punto (0,1),(0,1), entonces la pendiente de la solución que pasa por ese punto viene dada por y=3(0)+2 (1)4=–2.y=3(0)+2 (1)4=–2. Ahora supongamos que xx aumenta ligeramente, digamos que a x=0,1.x=0,1. Utilizando el método de aproximaciones lineales se obtiene una fórmula para el valor aproximado de yy para x=0,1.x=0,1. En particular,

L(x)=y0+f(x0)(xx0)=12 (x0)=12 x.L(x)=y0+f(x0)(xx0)=12 (x0)=12 x.

Sustituyendo x=0,1x=0,1 en L(x)L(x) da un valor yy aproximado de 0,8.0,8.

En este punto la pendiente de la solución cambia (de nuevo según la ecuación diferencial). Podemos seguir avanzando, recalculando la pendiente de la solución a medida que damos pequeños pasos hacia la derecha, y observando el comportamiento de la solución. La Figura 4.8 muestra un gráfico de la solución que pasa por el punto (0,1).(0,1).

Un gráfico del campo de direcciones para la ecuación diferencial y' = 3 x + 2 y - 4 en los cuatro cuadrantes. En los cuadrantes dos y tres, las flechas apuntan hacia abajo y ligeramente hacia la derecha. En una línea diagonal, aproximadamente y = -x + 2, las flechas apuntan cada vez más a la derecha, se curvan y luego apuntan hacia arriba por encima de esa línea. Se muestra la solución que pasa por el punto (0, 1). Se curva hacia abajo a través de (-5, 10), (0, 2), (1, 0) y (3, -10).
Figura 4.8 Campo de direcciones para la ecuación diferencial y = 3 x + 2 y 4 y = 3 x + 2 y 4 con la solución que pasa por el punto ( 0 , 1 ) . ( 0 , 1 ) .

La curva es el gráfico de la solución del problema de valor inicial

y=3x+2 y4,y(0)=1.y=3x+2 y4,y(0)=1.

Esta curva se denomina curva solución que pasa por el punto (0,1).(0,1). La solución exacta de este problema de valor inicial es

y=32 x+5414e2 x,y=32 x+5414e2 x,

y el gráfico de esta solución es idéntico a la curva en la Figura 4.8.

Punto de control 4.7

Cree un campo de direcciones para la ecuación diferencial y=x2 y2 y=x2 y2 y dibuje una curva solución que pase por el punto (–1,2 ).(–1,2 ).

Medios

Vaya a este sitio web para ver más sobre los campos de pendiente.

Consideremos ahora el campo de direcciones para la ecuación diferencial y=(x3)(y2 4),y=(x3)(y2 4), se muestra en la Figura 4.9. Este campo de direcciones tiene varias propiedades interesantes. En primer lugar, en y=−2y=−2 y y=2 ,y=2 , aparecen guiones horizontales en todo el gráfico. Esto significa que si y=–2,y=–2, entonces y=0.y=0. Sustituyendo esta expresión en el lado derecho de la ecuación diferencial se obtiene

(x3)(y2 4)=(x3)((−2)2 4)=(x3)(0)=0=y.(x3)(y2 4)=(x3)((−2)2 4)=(x3)(0)=0=y.

Por lo tanto, y=−2y=−2 es una solución de la ecuación diferencial. De la misma manera, y=2 y=2 es una solución de la ecuación diferencial. Estas son las únicas soluciones de valor constante de la ecuación diferencial, como podemos ver en el siguiente argumento. Supongamos que y=ky=k es una solución constante de la ecuación diferencial. Entonces y=0.y=0. Sustituyendo esta expresión en la ecuación diferencial se obtiene 0=(x3)(k2 4).0=(x3)(k2 4). Esta ecuación debe ser cierta para todos los valores de x,x, por lo que el segundo factor debe ser igual a cero. Este resultado produce la ecuación k2 4=0.k2 4=0. Las soluciones de esta ecuación son k=−2k=−2 y k=2 ,k=2 , que son las soluciones constantes ya mencionadas. Estas son las llamadas soluciones de equilibrio de la ecuación diferencial.

Un campo de direcciones para la ecuación diferencial dada. Las flechas son horizontales y apuntan hacia la derecha en y = -4, y = 4 y x = 6. Cuanto más cerca están las flechas de x = 6, más horizontales se vuelven. Cuanto más lejos, más verticales son. Las flechas apuntan hacia abajo para y > 4 y y < 4, -4 < y < 4 y x > 6 y y < -4 y x < 6. En todas las demás áreas, las flechas apuntan hacia arriba.
Figura 4.9 Campo de direcciones para la ecuación diferencial y = ( x 3 ) ( y 2 4 ) y = ( x 3 ) ( y 2 4 ) mostrando dos soluciones. Estas soluciones están muy próximas, pero una de ellas está apenas por encima de la solución de equilibrio x = −2 x = −2 y la otra está apenas por debajo de la misma solución de equilibrio.

Definición

Consideremos la ecuación diferencial y=f(x,y).y=f(x,y). Una solución de equilibrio es cualquier solución de la ecuación diferencial de la forma y=c,y=c, donde cc es una constante.

Para determinar las soluciones de equilibrio de la ecuación diferencial y=f(x,y),y=f(x,y), iguale el lado derecho a cero. Una solución de equilibrio de la ecuación diferencial es cualquier función de la forma y=ky=k tal que f(x,k)=0f(x,k)=0 para todos los valores de xx en el dominio de f.f.

Una característica importante de las soluciones de equilibrio es si se acercan o no a la línea y=ky=k como asíntota para valores grandes de x.x.

Definición

Consideremos la ecuación diferencial y=f(x,y),y=f(x,y), y supongamos que todas las soluciones de esta ecuación diferencial están definidas para xx0.xx0. Supongamos que y=ky=k sea una solución de equilibrio de la ecuación diferencial.

  1. y=ky=k es una solución asintóticamente estable de la ecuación diferencial si existe ε>0ε>0 tal que para cualquier valor c(kε,k+ε)c(kε,k+ε) la solución del problema de valor inicial
    y=f(x,y),y(x0)=cy=f(x,y),y(x0)=c

    se acerca a kk como xx se acerca al infinito.
  2. y=ky=k es una solución asintóticamente inestable de la ecuación diferencial si existe ε>0ε>0 tal que para cualquier valor c(kε,k+ε)c(kε,k+ε) la solución del problema de valor inicial
    y=f(x,y),y(x0)=cy=f(x,y),y(x0)=c

    nunca se acerca a kk como xx se acerca al infinito.
  3. y=ky=k es una solución asintóticamente semiestable de la ecuación diferencial si no es ni asintóticamente estable ni asintóticamente inestable.

Ahora volvemos a la ecuación diferencial y=(x3)(y2 4),y=(x3)(y2 4), con la condición inicial y(0)=0,5.y(0)=0,5. El campo de direcciones para este problema de valores iniciales, junto con la solución correspondiente, se muestra en la Figura 4.10.

Un campo de direcciones para la ecuación diferencial dada. Las flechas son horizontales y apuntan hacia la derecha en y = -4, y = 4 y x = 6. Cuanto más cerca están las flechas de x = 6, más horizontales se vuelven. Cuanto más lejos, más verticales son. Las flechas apuntan hacia abajo para y > 4 y y < 4, -4 < y < 4 y x > 6 y y < -4 y x < 6. En todas las demás áreas, las flechas apuntan hacia arriba. Se grafica una solución que pasa por (0, 0,5). Comienza a lo largo de y = -4 en el cuadrante tres, aumenta de -4 a 4 entre x = -1 y 1, y termina yendo a lo largo de y = 4 en el cuadrante 1.
Figura 4.10 Campo de direcciones para el problema de valor inicial y = ( x 3 ) ( y 2 4 ) , y ( 0 ) = 0,5 . y = ( x 3 ) ( y 2 4 ) , y ( 0 ) = 0,5 .

Los valores de la solución de este problema de valor inicial se mantienen entre y=−2y=−2 y y=2 ,y=2 , que son las soluciones de equilibrio de la ecuación diferencial. Sin embargo, un problema de valor inicial que comienza con 2 <y<2 2 <y<2 nunca puede cruzar las soluciones de equilibrio y=2 y=2 y y=−2y=−2. Por lo tanto, dado que y2 4<0y2 4<0 y para x>3,y'=(x3)(y2 4)<0x>3,y'=(x3)(y2 4)<0, yy es decreciente y por eso se acerca a y=2 y=2 . Por lo tanto, y=2 y=2 es una solución asintóticamente estable de la ecuación diferencial.

¿Qué ocurre cuando el valor inicial es inferior a y=−2?y=−2? Este escenario se ilustra en la Figura 4.11, con el valor inicial y(0)=−3.y(0)=−3.

Un campo de direcciones para la ecuación diferencial dada. Las flechas son horizontales y apuntan hacia la derecha en y = -4, y = 4 y x = 6. Cuanto más cerca están las flechas de x = 6, más horizontales se vuelven. Cuanto más lejos, más verticales son. Las flechas apuntan hacia abajo para y > 4 y y < 4, -4 < y < 4 y x > 6 y y < -4 y x < 6. En todas las demás áreas, las flechas apuntan hacia arriba. Se grafica una solución que va a lo largo de y = -4 en el cuadrante 3 y curvas entre x = -1 y x = 0 para ir al infinito negativo a lo largo del eje y.
Figura 4.11 Campo de direcciones para el problema de valor inicial y = ( x 3 ) ( y 2 4 ) , y ( 0 ) = −3 . y = ( x 3 ) ( y 2 4 ) , y ( 0 ) = −3 .

Podemos ver que para los valores iniciales y<2 y<2 por x>3,y'(x3)(y2 4)>0x>3,y'(x3)(y2 4)>0 y yy es creciente y por eso se acerca a y=2 y=2 . Esto reafirma que y=2 y=2 es una solución asintóticamente estable de la ecuación diferencial.

Ejemplo 4.8

Estabilidad de una solución de equilibrio

Cree un campo de direcciones para la ecuación diferencial y=(y3)2 (y2 +y2 )y=(y3)2 (y2 +y2 ) e identifique cualquier solución de equilibrio. Clasifique cada una de las soluciones de equilibrio como estable, inestable o semiestable.

Análisis

Es posible hallar las soluciones de equilibrio de la ecuación diferencial igualando el lado derecho a cero y resolviendo para y.y. Este enfoque da las mismas soluciones de equilibrio que las que vimos en el campo de direcciones.

Punto de control 4.8

Cree un campo de direcciones para la ecuación diferencial y=(x+5)(y+2 )(y2 4y+4)y=(x+5)(y+2 )(y2 4y+4) e identifique cualquier solución de equilibrio. Clasifique cada una de las soluciones de equilibrio como estable, inestable o semiestable.

Método de Euler

Consideremos el problema de valor inicial

y=2 x3,y(0)=3.y=2 x3,y(0)=3.

Integrando ambos lados de la ecuación diferencial se obtiene y=x2 3x+C,y=x2 3x+C, y resolviendo para CC se obtiene la solución particular y=x2 3x+3.y=x2 3x+3. La solución de este problema de valor inicial aparece como la parábola en la Figura 4.13.

Un gráfico en el rango [-1,4] para x y y. La parábola de apertura ascendente dada se dibuja con vértice en (1,5, 0,75). Los puntos individuales se representan en (0, 3), (0,5, 1,5), (1, 0,5), (1,5, 0), (2, 0), (2,5, 0,5) y (3, 1,5) con segmentos de línea que los conectan.
Figura 4.13 Método de Euler para el problema de valor inicial y = 2 x 3 , y ( 0 ) = 3 . y = 2 x 3 , y ( 0 ) = 3 .

El gráfico rojo está compuesto por segmentos de línea que se aproximan a la solución del problema de valor inicial. El gráfico comienza con el mismo valor inicial de (0,3).(0,3). Entonces la pendiente de la solución en cualquier punto está determinada por el lado derecho de la ecuación diferencial, y la longitud del segmento de la línea se determina aumentando el valor xx por 0,50,5 cada vez (el tamaño del paso). Este enfoque es la base del método de Euler.

Antes de enunciar el método de Euler como un teorema, consideremos otro problema de valor inicial:

y=x2 y2 ,y(–1)=2 .y=x2 y2 ,y(–1)=2 .

La idea de los campos de direcciones también puede aplicarse a este problema para estudiar el comportamiento de su solución. Por ejemplo, en el punto (–1,2 ),(–1,2 ), la pendiente de la solución está dada por y=(–1)2 2 2 =−3,y=(–1)2 2 2 =−3, por lo que la pendiente de la línea tangente a la solución en ese punto también es igual a −3.−3. Ahora definimos x0=–1x0=–1 y y0=2 .y0=2 . Como la pendiente de la solución en este punto es igual a −3,−3, podemos utilizar el método de aproximación lineal para aproximar yy cerca de (–1,2 ).(–1,2 ).

L(x)=y0+f(x0)(xx0).L(x)=y0+f(x0)(xx0).

Aquí x0=−1,y0=2 ,x0=−1,y0=2 , y f(x0)=−3,f(x0)=−3, por lo que la aproximación lineal se convierte en

L(x)=2 3(x(–1))=2 3x3=−3x1L(x)=2 3(x(–1))=2 3x3=−3x1

Ahora elegimos un tamaño de paso. El tamaño del paso es un valor pequeño, normalmente 0,10,1 o menos, que sirve de incremento para x;x; se representa con la variable h.h. En nuestro ejemplo, supongamos que h=0,1.h=0,1. Incrementando x0x0 por hh da nuestro siguiente valor xx:

x1=x0+h=−1+0,1=−0,9.x1=x0+h=−1+0,1=−0,9.

Podemos sustituir x1=−0,9x1=−0,9 en la aproximación lineal para calcular y1.y1.

y1=L(x1)=−3(−0,9)1=1,7.y1=L(x1)=−3(−0,9)1=1,7.

Por lo tanto, el valor aproximado yy de la solución cuando x=−0,9x=−0,9 es y=1,7.y=1,7. Podemos entonces repetir el proceso, utilizando x1=−0,9x1=−0,9 y y1=1,7y1=1,7 para calcular x2 x2 y y2 .y2 . La nueva pendiente está dada por y=(−0,9)2 (1,7)2 =−2,08.y=(−0,9)2 (1,7)2 =−2,08. Primero, x2 =x1+h=−0,9+0,1=−0,8.x2 =x1+h=−0,9+0,1=−0,8. Utilizando la aproximación lineal se obtiene

L(x)=y1+f(x1)(xx1)=1,72,08(x(−0,9))=1,72,08x1,872=−2,08x0,172.L(x)=y1+f(x1)(xx1)=1,72,08(x(−0,9))=1,72,08x1,872=−2,08x0,172.

Por último, sustituimos x2 =−0,8x2 =−0,8 en la aproximación lineal para calcular y2 .y2 .

y2 =L(x2 )=−2,08x2 0,172=−2,08(−0,8)0,172=1,492.y2 =L(x2 )=−2,08x2 0,172=−2,08(−0,8)0,172=1,492.

Por tanto, el valor aproximado de la solución de la ecuación diferencial es y=1,492y=1,492 cuando x=−0,8.x=−0,8.

Lo que acabamos de mostrar es la idea en la que se basa el método de Euler. Repitiendo estos pasos se obtiene una lista de valores para la solución. Estos valores se muestran en la Tabla 4.2, redondeados a cuatro decimales.

nn 00 11 2 2 33 44 55
xnxn −1−1 −0,9−0,9 −0,8−0,8 −0,7−0,7 −0,6−0,6 −0,5−0,5
ynyn 2 2 1,71,7 1,4921,492 1,33341,3334 1,20461,2046 1,09551,0955
nn 66 77 88 99 1010
xnxn −0,4−0,4 −0,3−0,3 −0,2−0,2 −0,1−0,1 00
ynyn 1,00041,0004 1,91641,9164 1,84141,8414 1,77461,7746 1,71561,7156
Tabla 4.2 Usar el método de Euler para aproximar soluciones a una ecuación diferencial

Teorema 4.1

Método de Euler

Consideremos el problema de valor inicial

y=f(x,y),y(x0)=y0.y=f(x,y),y(x0)=y0.

Para aproximar una solución a este problema utilizando el método de Euler, defina

xn=x0+nhyn=yn1+hf(xn1,yn1).xn=x0+nhyn=yn1+hf(xn1,yn1).
(4.2)

Aquí h>0h>0 representa el tamaño de paso y nn es un número entero, empezando por 1.1. El número de pasos dados se cuenta con la variable n.n.

Normalmente hh es un valor pequeño, por ejemplo 0,10,1 o 0,05.0,05. Cuanto menor sea el valor de h,h, más cálculos se necesitan. Cuanto mayor sea el valor de h,h, menos cálculos se necesitan. Sin embargo, la compensación resulta en un menor grado de exactitud para un tamaño de paso mayor, como se ilustra en la Figura 4.14.

Dos gráficos de la misma parábola, y = x ^ 2 - 3 x + 3. El primero muestra el método de Euler para el problema de valor inicial dado con un tamaño de paso de h = 0,05, y el segundo muestra el método de Euler con un tamaño de paso de h = 0,25. El primero tiene los puntos (0, 3), (0,5, 1,5), (1, 0,5), (1,5, 0), (2, 0), (2,5, 0,5) y (3, 1,5) trazados con segmentos de línea que los conectan. El segundo tiene los puntos (0, 3), (0,25, 2,25), (0,5, 1,625), (0,75, 1,125), (1, 0,75), (1,25, 0,5), (1,5, 0,375), (2, 0,5), (2,25, 0,75), (2,5, 1,125), (2,75, 1,625) y (3, 2,25) trazados con segmentos de línea que los conectan.
Figura 4.14 Método de Euler para el problema de valor inicial y = 2 x 3 , y ( 0 ) = 3 y = 2 x 3 , y ( 0 ) = 3 con (a) un tamaño de paso de h = 0,5 ; h = 0,5 ; y (b) un tamaño de paso de h = 0,25 . h = 0,25 .

Ejemplo 4.9

Uso del método de Euler

Consideremos el problema de valor inicial

y=3x2 y2 +1,y(0)=2 .y=3x2 y2 +1,y(0)=2 .

Utilice el método de Euler con un tamaño de paso de 0,10,1 para generar una tabla de valores de la solución para valores de xx entre 00 y 1.1.

Medios

Para obtener más información sobre el método de Euler utilice esta miniaplicación.

Punto de control 4.9

Consideremos el problema de valor inicial

y=x3+y2 ,y(1)=–2.y=x3+y2 ,y(1)=–2.

Utilizando un tamaño de paso de 0,1,0,1, genere una tabla con los valores aproximados de la solución del problema de valor inicial para los valores de xx entre 11 y 2 .2 .

Medios

Visite este sitio web para ver una aplicación práctica del material de esta sección.

Sección 4.2 ejercicios

Para los siguientes problemas, utilice el siguiente campo de direcciones de la ecuación diferencial y=–2y.y=–2y. Dibuje el gráfico de la solución para las condiciones iniciales dadas.

Un campo de direcciones con flechas horizontales que apuntan hacia la derecha en 0. Las flechas sobre el eje x apuntan hacia abajo y hacia la derecha. Cuanto más lejos del eje x, más pronunciadas son las flechas, y cuanto más cerca del eje x, más planas son las flechas. Asimismo, las flechas situadas bajo el eje x apuntan hacia arriba y hacia la derecha. Cuanto más lejos del eje x, más pronunciadas son las flechas, y cuanto más cerca del eje x, más planas son las flechas.
66.

y ( 0 ) = 1 y ( 0 ) = 1

67.

y ( 0 ) = 0 y ( 0 ) = 0

68.

y ( 0 ) = –1 y ( 0 ) = –1

69.

¿Existe algún equilibrio? ¿Cuáles son sus estabilidades?

Para los siguientes problemas, utilice el siguiente campo de direcciones de la ecuación diferencial y=y2 2 y.y=y2 2 y. Dibuje el gráfico de la solución para las condiciones iniciales dadas.

Un campo de direcciones con flechas horizontales en y = 0 y y = 2. Las flechas apuntan hacia arriba para y > 2 y para y < 0. Las flechas apuntan hacia abajo para 0 < y < 2. Cuanto más cerca están las flechas de estas líneas, son más horizontales, y cuanto más lejos, son más verticales.
70.

y ( 0 ) = 3 y ( 0 ) = 3

71.

y ( 0 ) = 1 y ( 0 ) = 1

72.

y ( 0 ) = –1 y ( 0 ) = –1

73.

¿Existe algún equilibrio? ¿Cuáles son sus estabilidades?

Dibuje el campo de direcciones para las siguientes ecuaciones diferenciales, luego resuelva la ecuación diferencial. Dibuje su solución sobre el campo de direcciones. ¿Su solución sigue las flechas de su campo de direcciones?

74.

y = t 3 y = t 3

75.

y = e t y = e t

76.

d y d x = x 2 cos x d y d x = x 2 cos x

77.

d y d t = t e t d y d t = t e t

78.

d x d t = cosh ( t ) d x d t = cosh ( t )

Dibuje el campo de direcciones para las siguientes ecuaciones diferenciales. ¿Qué puede decir sobre el comportamiento de la solución? ¿Existen equilibrios? ¿Qué estabilidad tienen estos equilibrios?

79.

y = y 2 1 y = y 2 1

80.

y = y x y = y x

81.

y = 1 y 2 x 2 y = 1 y 2 x 2

82.

y = t 2 sen y y = t 2 sen y

83.

y = 3 y + x y y = 3 y + x y

Haga coincidir el campo de direcciones con las ecuaciones diferenciales dadas. Explique sus selecciones.

Un campo de direcciones con flechas que apuntan hacia abajo y hacia la derecha en los cuadrantes dos y tres. Después de cruzar el eje y, las flechas cambian de dirección y apuntan hacia la derecha.


Un campo de direcciones con flechas horizontales que apuntan hacia la izquierda en los cuadrantes dos y tres. Al cruzar el eje y, las flechas cambian y apuntan hacia arriba en los cuadrantes uno y cuatro.


Un campo de direcciones con flechas horizontales que apuntan hacia la derecha en el eje x. Arriba, las flechas apuntan hacia abajo y a la derecha, y debajo, las flechas apuntan hacia arriba y a la derecha. Cuanto más lejos del eje x, más verticales se vuelven las flechas.


Un campo de direcciones con flechas horizontales en los ejes x y y. Las flechas apuntan hacia abajo y hacia la derecha en los cuadrantes uno y tres. Apuntan hacia arriba y hacia la derecha en los cuadrantes dos y cuatro.


Un campo de direcciones con flechas que apuntan hacia arriba en los cuadrantes dos y tres, hacia la derecha en el eje y, y hacia abajo en los cuadrantes uno y cuatro.
84.

y = −3 y y = −3 y

85.

y = −3 t y = −3 t

86.

y = e t y = e t

87.

y = 1 2 y + t y = 1 2 y + t

88.

y = t y y = t y

Haga coincidir el campo de direcciones con las ecuaciones diferenciales dadas. Explique sus selecciones.

Un campo de direcciones con flechas horizontales que apuntan hacia la derecha en los ejes x y y. En los cuadrantes uno y tres, las flechas apuntan hacia arriba, y en los cuadrantes dos y cuatro, hacia abajo.


Un campo de direcciones con flechas horizontales que apuntan hacia la derecha en los ejes x y y. En los cuadrantes uno y tres, las flechas apuntan hacia arriba y hacia la derecha, y en los cuadrantes dos y cuatro, las flechas apuntan hacia abajo y hacia la derecha.


Un campo de direcciones con flechas horizontales que apuntan hacia la derecha en los ejes x y y. En los cuadrantes dos y tres, las flechas apuntan hacia abajo, y en los cuadrantes uno y cuatro, hacia arriba.


Un campo de direcciones con flechas horizontales que apuntan hacia la derecha en el eje x. Las flechas apuntan hacia arriba y hacia la derecha en todos los cuadrantes. Cuanto más cerca están las flechas del eje x, más horizontales son, y cuanto más lejos están, más verticales son.


Un campo de direcciones con flechas horizontales en el eje y. Las flechas también son más horizontales cerca de y = 1,5, y = -1,5, y el eje y. Para y > 1,5 y x < 0, para y < -1,5 y x < 0, y para -1,5 < y < 1,5 y x > 0-, las flechas apuntan hacia abajo. Para y> 1,5 y x > 0, para y < -1,5, para y < -1,5 y x > 0, y para -1,5 < y < 1,5 y x < 0, las flechas apuntan hacia arriba.
89.

y = t sen y y = t sen y

90.

y = t cos y y = t cos y

91.

y = t tan y y = t tan y

92.

y = sen 2 y y = sen 2 y

93.

y = y 2 t 3 y = y 2 t 3

Estime las siguientes soluciones mediante el método de Euler con n=5n=5 pasos sobre el intervalo t=[0,1].t=[0,1]. Si es capaz de resolver el problema de valor inicial exactamente, compare su solución con la solución exacta. Si no puede resolver el problema de valor inicial, se le proporcionará la solución exacta para que la compare con el método de Euler. ¿Qué precisión tiene el método de Euler?

94.

y = −3 y , y ( 0 ) = 1 y = −3 y , y ( 0 ) = 1

95.

y = t 2 y = t 2

96.

y=3ty,y(0)=1.y=3ty,y(0)=1. La solución exacta es y=3t+4et3y=3t+4et3

97.

y=y+t2 ,y(0)=3.y=y+t2 ,y(0)=3. La solución exacta es y=5et2 t2 2 ty=5et2 t2 2 t

98.

y = 2 t , y ( 0 ) = 0 y = 2 t , y ( 0 ) = 0

99.

[T] y=e(x+y),y(0)=−1.y=e(x+y),y(0)=−1. La solución exacta es y=ln(e+1ex)y=ln(e+1ex) grandes.

100.

y=y2 ln(x+1),y(0)=1.y=y2 ln(x+1),y(0)=1. La solución exacta es y=1(x+1)(ln(x+1)1)y=1(x+1)(ln(x+1)1) grandes.

101.

y=2 x,y(0)=0,y=2 x,y(0)=0, La solución exacta es y=2 x1ln(2 )y=2 x1ln(2 ) grandes.

102.

y=y,y(0)=−1.y=y,y(0)=−1. La solución exacta es y=ex.y=ex.

103.

y=−5t,y(0)=–2.y=−5t,y(0)=–2. La solución exacta es y=52 t2 2 y=52 t2 2

Las ecuaciones diferenciales pueden utilizarse para modelar epidemias de enfermedades. En el siguiente conjunto de problemas, examinamos el cambio de tamaño de dos subpoblaciones de personas que viven en una ciudad: los individuos infectados y los individuos susceptibles a la infección. SS representa el tamaño de la población susceptible, e II representa el tamaño de la población infectada. Suponemos que si una persona susceptible interactúa con una persona infectada, existe una probabilidad cc que la persona susceptible se infecte. Cada persona infectada se recupera de la infección a una tasa rr y vuelve a ser susceptible. Consideramos el caso de la gripe, en el que suponemos que nadie muere por la enfermedad, por lo que suponemos que el tamaño total de la población de las dos subpoblaciones es un número constante, N.N. Las ecuaciones diferenciales que modelan estos tamaños de población son

S=rIcSIyI=cSIrI.S=rIcSIyI=cSIrI.

Aquí cc representa la tasa de contacto y rr es la tasa de recuperación.

104.

Demuestre que, por nuestra suposición de que el tamaño total de la población es constante (S+I=N),(S+I=N), se puede reducir el sistema a una única ecuación diferencial en I:I=c(NI)IrI.I:I=c(NI)IrI.

105.

Suponiendo que los parámetros sean c=0,5,N=5,c=0,5,N=5, y r=0,5,r=0,5, dibuje el campo de direcciones resultante.

106.

[T] Utilice un programa de computadora o una calculadora para calcular la solución del problema de valor inicial y=ty,y(0)=2 y=ty,y(0)=2 utilizando el método de Euler con el tamaño de paso dado h.h. Halle la solución en t=1.t=1. A modo de pista, aquí se presenta el "pseudocódigo" de cómo escribir un programa de computadora para realizar el Método de Euler para y=f(t,y),y(0)=2 :y=f(t,y),y(0)=2 :

Cree la función f(t,y)f(t,y)

Defina los parámetros y(1)=y0,t(0)=0,y(1)=y0,t(0)=0, tamaño de paso h,h, y número total de pasos, NN

Escriba un bucle for:

para k=1hasta Nk=1hasta N

fn=f(t(k),y(k))fn=f(t(k),y(k)) grandes.

y ( k+1 ) = y ( k ) + h*fn y ( k+1 ) = y ( k ) + h*fn

t ( k+1 ) = t ( k ) + h t ( k+1 ) = t ( k ) + h

107.

Resuelva el problema de valor inicial para la solución exacta.

108.

Dibuje el campo de direcciones

109.

h = 1 h = 1

110.

[T] h=10h=10

111.

[T] h=100h=100

112.

[T] h=1.000h=1.000

113.

[T] Evalúe la solución exacta en t=1.t=1. Haga una tabla de errores para el error relativo entre la solución del método de Euler y la solución exacta. ¿Cuánto cambia el error? ¿Puede explicarlo?

Consideremos el problema de valor inicial y=–2y,y(0)=2 .y=–2y,y(0)=2 .

114.

Demuestre que y=2 e−2xy=2 e−2x resuelve este problema de valor inicial.

115.

Dibuje el campo de direcciones de esta ecuación diferencial.

116.

[T] A mano o con calculadora o computadora, aproxime la solución mediante el método de Euler en t=10t=10 utilizando h=5.h=5.

117.

[T] Con una calculadora o una computadora, aproxime la solución mediante el método de Euler a t=10t=10 utilizando h=100.h=100.

118.

[T] Grafique la respuesta exacta y cada aproximación de Euler (para h=5h=5 y h=100)h=100) en cada hh en el campo de direcciones ¿Qué observa?

Solicitar una copia impresa

As an Amazon Associate we earn from qualifying purchases.

Cita/Atribución

¿Desea citar, compartir o modificar este libro? Este libro utiliza la Creative Commons Attribution-NonCommercial-ShareAlike License y debe atribuir a OpenStax.

Información de atribución
  • Si redistribuye todo o parte de este libro en formato impreso, debe incluir en cada página física la siguiente atribución:
    Acceso gratis en https://openstax.org/books/c%C3%A1lculo-volumen-2/pages/1-introduccion
  • Si redistribuye todo o parte de este libro en formato digital, debe incluir en cada vista de la página digital la siguiente atribución:
    Acceso gratuito en https://openstax.org/books/c%C3%A1lculo-volumen-2/pages/1-introduccion
Información sobre citas

© 2 mar. 2022 OpenStax. El contenido de los libros de texto que produce OpenStax tiene una licencia de Creative Commons Attribution-NonCommercial-ShareAlike License . El nombre de OpenStax, el logotipo de OpenStax, las portadas de libros de OpenStax, el nombre de OpenStax CNX y el logotipo de OpenStax CNX no están sujetos a la licencia de Creative Commons y no se pueden reproducir sin el previo y expreso consentimiento por escrito de Rice University.