1
PRÁCTICA DE MÉTODOS NUMÉRICOS
INFORME DE LABORATORIO N°12
‘’ ECUACIONES DIFERENCIALES y MÉTODO DE RUNGE
KUTTA’’
1. OBJETIVOS:
- Representar los factores, que intervienen en la creación de Las Ecuaciones
Diferenciales.
- Método de Resolución de Ruge Kutta.
- Uso del software MATLAB
2. MARCO TEÓRICO:
2.1. ECUACIÓN DIFERENCIAL:
Una ecuación diferencial (ED) es una ecuación que relaciona de manera no
trivial a una función desconocida y una o más derivadas de esta función
desconocida con respecto a una o más variables independientes. Si la función
desconocida depende de una sola variable la ecuación diferencial se
llama ordinaria , por el contrario, si depende de más de una variable, se
llama parcial
Un ejemplo de ecuación diferencial ordinaria es:
2.1.1. ORDEN DE UNA ECUACIÓN DIFERENCIAL
El orden de una ecuación diferencial está dado por el orden mayor de su
derivada.
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
2
PRÁCTICA DE MÉTODOS NUMÉRICOS
2.1.2. GRADO DE UNA ECUACIÓN DIFERENCIAL:
El grado de una ecuación diferencial está dado por el exponente del
mayor orden de su derivada.
Ejemplos:
Determinar el orden y grado de las siguientes ecuaciones diferenciales
ordinarias.
2.1.3. SOLUCION DE UNA ECUACIÓN DIFERENCIAL
Una función que cuando se remplaza en la ecuación diferencial da una
igualdad, se llama una solución de la ecuación diferencial, por lo tanto,
resolver una ecuación diferencial es encontrar una función desconocida
que al ser sustituida en la ecuación diferencial se obtiene una igualdad.
- FUNCIÓN PRIMITIVA DE UNA ECUACIÓN DIFERENCIAL
Es una expresión equivalente a la ecuación diferencial que carece de
derivadas.
Ejemplo:
Resolver la ecuación diferencial
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
3
PRÁCTICA DE MÉTODOS NUMÉRICOS
- PROBLEMA DE VALOR INICIAL:
Un problema de valor inicial es un problema que busca determinar una
solución a una ecuación diferencia sujeta a condiciones sobre la función
desconocida y sus derivadas especificadas en un valor de la variable
independiente. Tales condiciones se llaman condiciones iniciales.
Un problema de valor de frontera es un problema que busca determinar
una solución a una ecuación diferencia sujeta a condiciones sobre la
función desconocida especificadas en dos o más valores de la variable
independiente. Tales condiciones se llaman condiciones de frontera.
Ejemplo ilustrativo
Una curva tiene la propiedad de que su pendiente en cualquier punto
(x,y) de ella es igual a 2x. Hallar la ecuación de la curva si ésta pasa por
el punto (2,5)
Solución:
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
4
PRÁCTICA DE MÉTODOS NUMÉRICOS
2.2. MÉTODO DE RUNGE KUTTA:
Es uno de los métodos más utilizados para resolver numéricamente problemas
de ecuaciones diferenciales ordinarias con condiciones iniciales es el método
de Runge-Kutta de cuarto orden, el cual proporciona un pequeño margen de
error con respecto a la solución real del problema y es fácilmente programable
en un software para realizar las iteraciones necesarias.
El método de Runge-Kutta se utiliza para resolver ecuaciones diferenciales de
la forma explícita:
O en su forma Explícita:
Y es sumamente útil para casos en los que la solución no puede hallarse por
los métodos convencionales (como separación de variables). Hay variaciones
en el método de Runge-Kutta de cuarto orden pero el más utilizado es el
método en el cual se elige un tamaño de paso h y un número máximo de
iteraciones n.
El método Runge- Kutta para este problema está dado por la siguiente
ecuación:
Para i=0,…, n-1. La solución se da a lo largo del intervalo (xo,xo+hn) , Donde:
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
5
PRÁCTICA DE MÉTODOS NUMÉRICOS
Así, siguiente valor (yi+1) es determinado por el presente valor (yi) más el
producto del tamaño del intervalo (h) por una pendiente estimada. La pendiente
un promedio ponderado de pendientes:
- k1 es la pendiente al principio del intervalo.
- k2 es la pendiente en el punto medio del intervalo, usando k1 para determinar
el valor de y en el punto xi + h/2.
- k3 es otra vez la pendiente del punto medio, pero ahora usando k2 para
determinar el valor de y
- k4 es la pendiente al final del intervalo, con el valor de y determinado por k3
Promediando las cuatro pendientes, se le asigna mayor peso a las pendientes
en el punto medio:
3. EQUIPOS, MATERIALES Y REACTIVOS:
3.1. EQUIPOS:
- Computadora personal
- Programa Microsoft Excel 2019
3.2. MATERIALES:
- Información de métodos Numéricos
- Plantillas en Excel para los problemas
3.3. REACTIVOS:
- No se usaron reactivos
4. METODOLOGÍA PROCEDIMIENTO:
4.1. PRÁCTICA DE LABORATORIO 13-03-2020:
1. En un tanque perfectamente agitado se tienen 400 L de una salmuera de
la cual están disueltas 25 kg de sal común (NaCl), en cierto momento se
hace llegar al tanque un gasto de 80 L/min de una salmuera que
contiene 0.5 kg. de sal común por litro, si se tiene un gasto de salida de
80 L/min determinar:
a) ¿Qué cantidad de sal hay en el tanque transcurridos 10 minutos?
b) ¿Qué cantidad de sal hay transcurrido un tiempo muy grande?
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
6
PRÁCTICA DE MÉTODOS NUMÉRICOS
SOLUCIÓN:
X: la cantidad de sal en Kg., en el tanque después de t minutos.
dX
La acumulación de sal en el tanque esta dado por dt y por la relación:
𝒅𝑿
= 𝒎𝒂𝒔𝒂 𝒅𝒆 𝒔𝒂𝒍 𝒒𝒖𝒆 𝒆𝒏𝒕𝒓𝒂 − 𝒎𝒂𝒔𝒂 𝒅𝒆 𝒔𝒂𝒍 𝒒𝒖𝒆 𝒒𝒖𝒆𝒅𝒂
𝒅𝒕
𝒅𝑿 𝒍𝒊𝒃 𝒌𝒈 𝒍𝒊𝒃 𝑿 𝒌𝒈
= 𝟖𝟎 ∗ 𝟎. 𝟓 − 𝟖𝟎 ∗ ∗
𝒅𝒕 𝒎𝒊𝒏 𝑳 𝒎𝒊𝒏 𝟒𝟎𝟎 𝑳
𝒅𝑿
= 𝟒𝟎 − 𝟎. 𝟐𝑿
𝒅𝒕
con las condiciones iniciales:
𝒅𝑿
= 𝟒𝟎 − 𝟎. 𝟐𝑿
𝑷. 𝑽. 𝑰 = { 𝒅𝒕
𝒙(𝟎) = 𝟐𝟓
𝒙(𝟏𝟎) =?
𝟏
𝒅𝑿 = 𝒅𝒕
𝟒𝟎 − 𝟎. 𝟐𝑿
𝟏
∫ 𝒅𝑿 = ∫ 𝒅𝒕
𝟒𝟎 − 𝟎. 𝟐𝑿
−𝟓 𝐥𝐧(𝟒𝟎 − 𝟎. 𝟐𝑿) + 𝑪𝟏 = 𝒕 + 𝑪𝟐
−𝟓 𝐥𝐧(𝟒𝟎 − 𝟎. 𝟐𝑿) = 𝒕 + 𝑪𝟐 − 𝑪𝟏
−𝟓 𝐥𝐧(𝟒𝟎 − 𝟎. 𝟐𝑿) = 𝒕 + 𝑪
𝒕 𝑪
𝐥𝐧(𝟒𝟎 − 𝟎. 𝟐𝑿) = +
−𝟓 −𝟓
−𝒕
𝐥𝐧(𝟒𝟎 − 𝟎. 𝟐𝑿) = +𝑪
𝟓
−𝒕
𝟒𝟎 − 𝟎. 𝟐𝐗 = 𝒆 𝟓 +𝑪
−𝒕
𝟎. 𝟐𝐗 = 𝟒𝟎 − 𝒆 𝟓 +𝑪
𝒕
𝟒𝟎 𝒆−𝟓+𝑪
𝐗= −
𝟎. 𝟐 𝟎. 𝟐
𝒕
𝒆−𝟓+𝑪
𝐗 = 𝟐𝟎𝟎 −
𝟎. 𝟐
𝒕
𝒆−𝟓 ∗ 𝒆𝑪
𝐗 = 𝟐𝟎𝟎 −
𝟎. 𝟐
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
7
PRÁCTICA DE MÉTODOS NUMÉRICOS
𝒕
𝒆−𝟓 ∗ 𝒆𝑪
𝐗 = 𝟐𝟎𝟎 − ; 𝑿(𝟎) = 𝟐𝟓
𝟎. 𝟐
𝟎
𝒆−𝟓 ∗ 𝒆𝑪
𝟐𝟓 = 𝟐𝟎𝟎 −
𝟎. 𝟐
𝟏 ∗ 𝒆𝑪
𝟐𝟓 = 𝟐𝟎𝟎 −
𝟎. 𝟐
𝟏 ∗ 𝒆𝑪
= 𝟐𝟎𝟎 − 𝟐𝟓
𝟎. 𝟐
𝟏 ∗ 𝒆𝑪
= 𝟏𝟕𝟓
𝟎. 𝟐
𝟏 ∗ 𝒆𝑪 = 𝟏𝟕𝟓 ∗ 𝟎. 𝟐
𝟏 ∗ 𝒆𝑪 = 𝟏𝟕𝟓 ∗ 𝟎. 𝟐
𝒆𝑪 = 𝟑𝟓
𝒕
𝒆−𝟓 ∗ 𝟑𝟓
∴ 𝐗 = 𝟐𝟎𝟎 −
𝟎. 𝟐
a) t=10min
(𝟏𝟎)
𝒆− 𝟓 ∗ 𝟑𝟓
𝐗 = 𝟐𝟎𝟎 −
𝟎. 𝟐
𝒆−𝟐 ∗ 𝟑𝟓
𝐗 = 𝟐𝟎𝟎 −
𝟎. 𝟐
𝐗 = 𝟐𝟎𝟎 − 𝟐𝟑. 𝟔𝟖𝟑𝟔𝟕𝟒𝟓𝟕
𝐗 = 𝟏𝟕𝟔, 𝟑𝟏𝟔𝟑𝟐𝟓𝟒 𝒌𝒈 𝒅𝒆 𝒔𝒂𝒍
b) t=infinito, t=100
𝒆−𝟏𝟎𝟎 ∗ 𝟑𝟓
𝐗 = 𝟐𝟎𝟎 −
𝟎. 𝟐
𝐗 = 𝟐𝟎𝟎 − 𝟔, 𝟓𝟏𝟎𝟏𝟑𝟐𝟗𝟓𝟖 ∗ 𝟏𝟎−𝟒𝟐
𝐗 = 𝟐𝟎𝟎 𝒂𝒑𝒓𝒐𝒙 𝒌𝒈 𝒅𝒆 𝒔𝒂𝒍
𝐒𝐞 𝐨𝐛𝐭𝐢𝐞𝐧𝐞 𝟐𝟎𝟎 𝐤𝐠 𝐝𝐞 𝐬𝐚𝐥
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
8
PRÁCTICA DE MÉTODOS NUMÉRICOS
EJECUTANDO EL PROBLEMA EN EXCEL
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
9
PRÁCTICA DE MÉTODOS NUMÉRICOS
4.2. PRÁCTICA DE LABORATORIO 13-03-2020:
𝒅𝒚 𝒚
=𝟐−
𝒅𝒙 𝒙
𝒅𝒚 𝒚
+ =𝟐
𝒅𝒙 𝒙
𝒅𝒚 𝟏
+ ∗𝒚=𝟐
𝒅𝒙 𝒙
EDO LINEAL DE PRIMER ORDEN:
𝒚′ (𝒙) + 𝒑(𝒙) ∗ 𝒚 = 𝒒(𝒙)
Donde: p(x)= 1/x; q(x)=2
𝑪
𝒚=𝒙+
𝒙
Si y(1)=0
𝑪
𝟎=𝟏+
𝟏
𝑪 = −𝟏
Por lo tanto:
𝟏
𝒚=𝒙−
𝒙
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
10
PRÁCTICA DE MÉTODOS NUMÉRICOS
EJECUTANDO EL PROBLEMA EN EXCEL
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
11
PRÁCTICA DE MÉTODOS NUMÉRICOS
4.3.- PRÁCTICA DE TEORÍA, ECUACIONES DIFERENCIALES CON VALOR
INICIAL:
4.3.1 Dada la condición inicial Y(o) = 0, encontrar la solución particular para
la ecuación
𝐘 ′ = 𝒆𝟐𝒙+𝟑𝒚
Aplicar también el Método de Runge Kutta 4to Orden
Hasta: Y(2), paso h= 0.01
SOLUCIÓN ANALÍTICA:
𝒚′ = 𝒆𝟐𝒙+𝟑𝒚
𝒅𝒚
= 𝒆𝟐𝒙 ∗ 𝒆𝟑𝒚
𝒅𝒙
𝟏
𝒅𝒚 = 𝒆𝟐𝒙 𝒅𝒙
𝒆𝟑𝒚
𝒆−𝟑𝒚 𝒅𝒚 = 𝒆𝟐𝒙 𝒅𝒙
∫ 𝒆−𝟑𝒚 𝒅𝒚 = ∫ 𝒆𝟐𝒙 𝒅𝒙
𝒆−𝟑𝒚 𝒆𝟐𝒙
+ 𝑪𝟏 = + 𝑪𝟐
−𝟑 𝟐
𝒆−𝟑𝒚 𝒆𝟐𝒙
= + 𝑪𝟐 − 𝑪𝟏
−𝟑 𝟐
𝒆−𝟑𝒚 𝒆𝟐𝒙
= + 𝑪𝟐 − 𝑪𝟏
−𝟑 𝟐
−𝟑𝒚
−𝟑𝒆𝟐𝒙
𝒆 = − 𝟑(𝑪𝟐 − 𝑪𝟏)
𝟐
−𝟑𝒆𝟐𝒙
𝒆−𝟑𝒚 = +𝑪
𝟐
−𝟑𝒆𝟐𝒙
−𝟑𝒚 = 𝒍𝒏 ( + 𝑪)
𝟐
−𝟑𝒆𝟐𝒙
−𝒍𝒏 ( 𝟐 + 𝑪)
𝒚= … 𝑺𝒐𝒍𝒖𝒄𝒊ó𝒏 𝒈𝒆𝒏𝒆𝒓𝒂𝒍
𝟑
Si y(0)=0
−𝟑𝒆𝟐(𝟎)
−𝒍𝒏 ( + 𝑪)
𝟐
𝟎=
𝟑
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
12
PRÁCTICA DE MÉTODOS NUMÉRICOS
−𝟑(𝟏)
𝟎 = −𝒍𝒏 ( + 𝑪)
𝟐
−𝟑
𝒆𝟎 = +𝑪
𝟐
−𝟑
𝟏= +𝑪
𝟐
𝟑
𝑪=𝟏+
𝟐
𝟓
𝑪=
𝟐
Por lo tanto:
−𝟑𝒆𝟐𝒙 𝟓
−𝒍𝒏 ( 𝟐 + 𝟐)
𝒚= … 𝑺𝒐𝒍𝒖𝒄𝒊ó𝒏 𝑷𝒂𝒓𝒕𝒊𝒄𝒖𝒍𝒂𝒓
𝟑
SOLUCIÓN POR RUNGE KUTTA 4TO ORDEN:
Nota: Mientras el método busca más iteracciones los cálculos se vuelven
más difíciles de hallar, por ende no se puede llegar a ver el xf=2
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
13
PRÁCTICA DE MÉTODOS NUMÉRICOS
4.3.2. Dada la condición inicial Y(o) = 1, encontrar la solución particular para
la ecuación
𝒅𝒚
= 𝒙𝒚 + 𝒙𝒚𝟐
𝒅𝒙
Aplicar también el Método de Runge Kutta 4to Orden
Hasta: Y(1), paso h= 0.1
𝒅𝒚
= 𝒙𝒚 + 𝒙𝒚𝟐
𝒅𝒙
𝒅𝒚
= 𝒙 ( 𝒚 + 𝒚𝟐 )
𝒅𝒙
𝟏
𝒅𝒚 = 𝒙𝒅𝒙
𝒚 + 𝒚𝟐
𝟏
∫ 𝒅𝒚 = ∫ 𝒙 𝒅𝒙
𝒚 + 𝒚𝟐
𝒙𝟐
𝒍𝒏(𝒚) − 𝒍𝒏(𝟏 + 𝒚) + 𝑪𝟏 = + 𝑪𝟐
𝟐
𝒚 𝒙𝟐
𝒍𝒏 ( )= + 𝑪𝟐 − 𝑪𝟏
𝟏+𝒚 𝟐
𝒚 𝒙𝟐
𝒍𝒏 ( )= +𝑪
𝟏+𝒚 𝟐
𝒚 𝒙𝟐
= 𝒆 𝟐 +𝑪
𝟏+𝒚
𝒙𝟐
𝒚 = (𝟏 + 𝒚)𝒆 𝟐 +𝑪
𝒙𝟐 𝒙𝟐
𝒚 = 𝒆 𝟐 +𝑪 +𝒚𝒆 𝟐 +𝑪
𝒙𝟐 𝒙𝟐
𝒚 − 𝒚𝒆 𝟐 +𝑪 = 𝒆 𝟐 +𝑪
𝒙𝟐 𝒙𝟐
𝒚 (𝟏 − 𝒆 𝟐 +𝑪 ) = 𝒆 𝟐 +𝑪
𝒙𝟐
𝒆 𝟐 +𝑪
𝒚= … 𝑺𝒐𝒍𝒖𝒄𝒊ó𝒏 𝑮𝒆𝒏𝒆𝒓𝒂𝒍
𝒙𝟐
(𝟏 − 𝒆 𝟐 +𝑪 )
Sí y(0)=1; entonces:
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
14
PRÁCTICA DE MÉTODOS NUMÉRICOS
(𝟎)𝟐
+𝑪
𝒆 𝟐
𝟏= (𝟎)𝟐
(𝟏 − 𝒆 𝟐 +𝑪 )
(𝟏 − 𝒆𝑪 ) = 𝒆𝑪
𝟏 = 𝟐𝒆𝑪
𝟏
= 𝒆𝑪
𝟐
Por lo tanto:
𝒙𝟐 +𝟏
𝒆 𝟐
𝒚= … 𝑺𝒐𝒍𝒖𝒄𝒊ó𝒏 𝑷𝒂𝒓𝒕𝒊𝒄𝒖𝒍𝒂𝒓
𝒙𝟐 +𝟏
(𝟏 − 𝒆 𝟐 )
SOLUCIÓN POR RUNGE KUTTA 4TO ORDEN:
NOTA: Podemos observación que en xf=1 nos da como resultado:
yf=4.690776751
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
15
PRÁCTICA DE MÉTODOS NUMÉRICOS
4.3.3. Dada la condición inicial Y(0.5) = 1, encontrar la solución particular
para la ecuación
Y’− 𝟐𝒙𝒚 = 𝒙
Aplicar también el Método de Runge Kutta 4to Orden
Hasta: Y(1), paso h= 0.1
𝒅𝒚
= 𝒙 + 𝟐𝒙𝒚
𝒅𝒙
𝒅𝒚
= 𝒙(𝟏 + 𝟐𝒚)
𝒅𝒙
𝟏
𝒅𝒚 = 𝒙𝒅𝒙
𝟏 + 𝟐𝒚
𝟏
∫ 𝒅𝒚 = ∫ 𝒙 𝒅𝒙
𝟏 + 𝟐𝒚
𝟏 𝒙𝟐
𝒍𝒏(𝟏 + 𝟐𝒚) + 𝑪𝟏 = + 𝑪𝟐
𝟐 𝟐
𝟏 𝒙𝟐
𝒍𝒏(𝟏 + 𝟐𝒚) = + 𝑪𝟐 − 𝑪𝟏
𝟐 𝟐
𝒍𝒏(𝟏 + 𝟐𝒚) = 𝒙𝟐 + 𝟐(𝑪𝟐 − 𝑪𝟏)
𝒍𝒏(𝟏 + 𝟐𝒚) = 𝒙𝟐 + 𝑪
𝟐 +𝑪
𝟏 + 𝟐𝒚 = 𝒆𝒙
𝟐 +𝑪
𝟐𝒚 = 𝒆𝒙 −𝟏
𝟐 +𝑪
𝒆𝒙 −𝟏
𝒚= … 𝑺𝒐𝒍𝒖𝒄𝒊ó𝒏 𝑮𝒆𝒏𝒆𝒓𝒂𝒍
𝟐
Si y(0.5)=1; entonces:
𝟐 +𝑪
𝒆𝟎.𝟓 −𝟏
𝟏=
𝟐
𝟐 +𝑪
𝟐 = 𝒆𝟎.𝟓 −𝟏
𝟐 +𝑪
𝟑 = 𝒆𝟎.𝟓
𝒍𝒏(𝟑) = 𝟎. 𝟓𝟐 + 𝑪
𝟏
𝒍𝒏(𝟑) = +𝑪
𝟒
𝟏
𝑪 = 𝒍𝒏(𝟑) −
𝟒
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
16
PRÁCTICA DE MÉTODOS NUMÉRICOS
Por lo tanto:
𝟐 +𝒍𝒏(𝟑)−𝟏
𝒆𝒙 𝟒 −𝟏
𝒚= … 𝑺𝒐𝒍𝒖𝒄𝒊ó𝒏 𝑷𝒂𝒓𝒕𝒊𝒄𝒖𝒍𝒂𝒓
𝟐
SOLUCIÓN POR RUNGE KUTTA 4TO ORDEN:
NOTA: Observamos que cuando xf=1 nos da el valor 2.67548681
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
17
PRÁCTICA DE MÉTODOS NUMÉRICOS
4.3.4 Dada la condición inicial Y(2) = 4, encontrar la solución particular para
la ecuación
𝒅𝒚
= 𝟎. 𝟏√𝒚 + 𝟎. 𝟒𝒙𝟐
𝒅𝒙
Aplicar también el Método de Runge Kutta 4to Orden.
Hasta: Y(2.5), 10 iteraciones
Observación: La ecuación diferencial es muy difícil para hallarla
analíticamente, por ello, se aplicara de frente del método de Runge kutta
4to orden porque para eso fue creado ese método numérico, para
ahorrarnos el trabajo de calcular ecuaciones diferenciales difíciles
NOTA: Observamos que cuando xf=1 nos da como valor 1.168406979
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
18
PRÁCTICA DE MÉTODOS NUMÉRICOS
4.4.- PRÁCTICA DE TEORÍA,PROBLEMA DE TANQUE CILÍNDRICO:
4.4.1. Se tiene un tanque cilíndrico abierto, de 2 metros de diámetro con
agua y tiene un agujero circular al costado de la base con un diámetro de
abertura de 0.05 metros, la altura del agua en el tanque es de 4 metros.
Determinar la altura que tendrá el líquido a los 200segundos
𝑸𝒔𝒂𝒍𝒊𝒅𝒂 = 𝑽𝑨𝑻
𝝅𝒅𝟐
𝑨𝑻 =
𝟒
Reemplazando:
𝝅𝒅𝟐
𝑸𝒔𝒂𝒍𝒊𝒅𝒂 = 𝑽 ( ) … … … … … (𝟏)
𝟒
Cálculo de la velocidad de salida:
𝑬𝒄𝒊𝒏é𝒕𝒊𝒄𝒂 = 𝑬𝒑𝒐𝒕𝒆𝒏𝒄𝒊𝒂𝒍
𝒎𝑽𝟐
= 𝒎𝒈𝒉
𝟐
𝑽 = √𝟐𝒈𝒉
Reemplazando en ecuación (1):
𝝅𝒅𝟐
𝑸𝒔𝒂𝒍𝒊𝒅𝒂 = (√𝟐𝒈𝒉) ( )
𝟒
𝒅𝑽
𝑸=
𝒅𝒕
Entonces:
𝒅𝑽 𝝅𝒅𝟐
= (√𝟐𝒈𝒉) ( )
𝒅𝒕 𝟒
𝝅𝒅𝟐
𝒅𝑽 = (√𝟐𝒈𝒉) ( ) 𝒅𝒕 … … … … . (𝟐)
𝟒
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
19
PRÁCTICA DE MÉTODOS NUMÉRICOS
Volumen del tanque:
𝑽 = 𝑨𝒃𝒂𝒔𝒆 𝒉
𝝅𝑫𝟐
𝑽=( ) 𝒉 … … … … . . (𝟑)
𝟒
Si V2= 6.28m^3 y V1=12.56m^3; entonces:
𝒅𝑽 = 𝑽𝟐 − 𝑽𝟏 = 𝟔. 𝟐𝟖𝒎𝟑 − 𝟏𝟐. 𝟓𝟔𝒎𝟑 = −𝟔. 𝟐𝟖𝒎𝟑
Es negativo porque se esta vaciando, entonces derivando ecuación (3) se
tiene:
𝝅𝑫𝟐
𝒅𝑽 = − ( ) 𝒅𝒉 … … … … (𝟒)
𝟒
Igualando las ecuaciones (2) y (4)
𝝅𝑫𝟐 𝝅𝒅𝟐
−( ) 𝒅𝒉 = (√𝟐𝒈𝒉) ( ) 𝒅𝒕
𝟒 𝟒
−𝟒 𝝅𝒅𝟐
𝒅𝒉 = ( ) (√𝟐𝒈𝒉) ( ) 𝒅𝒕
𝝅𝑫𝟐 𝟒
𝒅𝒉 √𝟐𝒈𝒉 ∗ 𝒅𝟐
= (− )
𝒅𝒕 𝑫𝟐
Si g=9.81m^2/2, d=0.05m y D=2m; entonces
𝒅𝒉 √𝟐(𝟗. 𝟖𝟏)𝒉 ∗ (𝟎. 𝟎𝟓)𝟐
= (− )
𝒅𝒕 (𝟐)𝟐
𝒅𝒉
= −𝟐. 𝟕𝟔𝟖𝟒 ∗ 𝟏𝟎−𝟑 √𝒉
𝒅𝒕
SOLUCIÓN ANALÍTICA:
𝒅𝒉
= −𝟐. 𝟕𝟔𝟖𝟒 ∗ 𝟏𝟎−𝟑 √𝒉
𝒅𝒕
𝟏
𝒅𝒉 = −𝟐. 𝟕𝟔𝟖𝟒 ∗ 𝟏𝟎−𝟑 𝒅𝒕
√𝒉
𝟏
∫ 𝒅𝒉 = − ∫ 𝟐. 𝟕𝟔𝟖𝟒 ∗ 𝟏𝟎−𝟑 𝒅𝒕
√𝒉
𝟐√𝒉 + 𝑪𝟏 = −𝟐. 𝟕𝟔𝟖𝟒 ∗ 𝟏𝟎−𝟑 𝒕 + 𝑪𝟐
𝟐√𝒉 = −𝟐. 𝟕𝟔𝟖𝟒 ∗ 𝟏𝟎−𝟑 𝒕 + 𝑪𝟐 − 𝑪𝟏
𝟏𝟎−𝟑 𝒕
√𝒉 = −𝟐. 𝟕𝟔𝟖𝟒 ∗ +𝑪
𝟐
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
20
PRÁCTICA DE MÉTODOS NUMÉRICOS
𝟐
𝟏𝟎−𝟑 𝒕
𝒉 = (−𝟐. 𝟕𝟔𝟖𝟒 ∗ + 𝑪)
𝟐
𝟐
𝟏𝟎−𝟑 𝒕
𝒉 = (−𝟐. 𝟕𝟔𝟖𝟒 ∗ + 𝑪) … . . 𝑺𝒐𝒍𝒖𝒄𝒊ó𝒏 𝑮𝒆𝒏𝒆𝒓𝒂𝒍
𝟐
Si t=0, h=4m porque el tanque al segundo 0 está completamente lleno;
entonces
𝟐
𝟏𝟎−𝟑 (𝟎)
𝟒 = (−𝟐. 𝟕𝟔𝟖𝟒 ∗ + 𝑪)
𝟐
𝟒 = (𝟎 + 𝑪)𝟐
𝑪=𝟐
Por lo tanto:
𝟐
𝟏𝟎−𝟑 𝒕
𝒉 = (−𝟐. 𝟕𝟔𝟖𝟒 ∗ + 𝟐) … . 𝑺𝒐𝒍𝒖𝒄𝒊ó𝒏 𝑷𝒂𝒓𝒕𝒊𝒄𝒖𝒍𝒂𝒓
𝟐
HALLANDO EL PROBLEMA ANALITICAMENTE:
t=200 segundos
𝟐
𝟏𝟎−𝟑 (𝟐𝟎𝟎)
𝒉 = (−𝟐. 𝟕𝟔𝟖𝟒 ∗ + 𝟐)
𝟐
𝒉 = (−𝟎. 𝟐𝟕𝟔𝟖𝟒 + 𝟐)𝟐
𝒉 = (−𝟎. 𝟐𝟕𝟔𝟖𝟒 + 𝟐)𝟐
𝒉 = (𝟏. 𝟕𝟐𝟑𝟏𝟔)𝟐
𝒉 = 𝟐. 𝟗𝟔𝟗𝟐𝟖𝟎𝟑𝟖𝟔
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
21
PRÁCTICA DE MÉTODOS NUMÉRICOS
HALLANDO EL PROBLEMA POR RUNGE KUTTA:
NOTA: Podemos observar que cuando t=200 con el método Runge Kutta
de 4to orden nos da como resultado 2.969280386
5. CONCLUSIONES:
- Se pudo identificar los valores de las ecuaciones diferenciales en las hojas
de calculo de Excel por el método de runge kutta de 4to orden
- Se pudo comprobar con éxito el calculo del método
- El uso de del software MATLAB no se pudo realizar porque no todos los
alumnos matriculados en el curso de métodos numéricos no cuentan con
dicho programa
6. RECOMEDACIONES:
- Se debe tener conocimientos básicos en el uso y aplicaciones de hojas
de cálculo, esto permite entender mejor la resolución de problemas e
idear propias formas de resolución de problemas.
- Identificar que ecuaciones diferenciales son factibles para ser usar con
el método de runge kutta de 4to roden
7. BIBLIOGRAFÍA:
- https://tarwi.lamolina.edu.pe/~duenas/rk4_2009II_pp2.pdf
- https://www.mathstools.com/section/main/Metodos_de_Runge_Kutta?lan
g=es#.XxUToChKiUk
- https://www.ugr.es/~lorente/APUNTESMNQ/cap23.pdf
ALUMNO: WILLIAMS ALEXANDER CHOQUE ORMACHEA
