Informe Parcial ME
Informe Parcial ME
Informe Parcial ME
2021-2
Integrantes: Yuranny Yeraldin Morillo Santacruz - Sara Hoyos Guisao - José Alberto Gómez
Echeverri
Problema [Tomado de Barrera.R & Díaz.L (2020). Métodos numéricos en Excel y Matlab Con
aplicaciones en ingeniería (1ra ed). Universidad De Antioquia.]
Como práctica académica, usted opta por la modalidad de trabajo social o servicio comunitario, y su
primera tarea es ayudar a diseñar un tanque de almacenamiento de agua (esférico) para un pequeño
poblado en las afueras de cierta cabecera municipal. Dadas las características del lugar, es imperativo
3
que el tanque esférico contenga 30 𝑚 de agua, pero las condiciones semiprecarias (limitada
tecnología y recursos disponibles) solo le permiten medir o controlar con certeza la altura del líquido
(nivel) dentro del tanque. Afortunadamente, de los primeros cursos en la universidad, usted sabe bien
que la siguiente ecuación le será muy útil:
(1)
Con V = volumen del tanque esférico (m3), h = altura o nivel del líquido en el tanque (m), r = radio
(m) y π = 3.1416. Utilice alguno de los métodos numéricos vistos en este capítulo para determinar
(redondeado a cuatro cifras decimales) cuál debe ser la altura o nivel del líquido que garantice el
3
volumen de 30𝑚 para r=2.08m.
Observación: para este ejercicio, el problema numérico presenta tres raíces reales, pero solo una se
ajusta a la física del problema que exige h ≤ 2r.
Resultados.
A continuación en la Tabla 1, se encuentran los resultados obtenidos utilizando cada método:
Tabla 1. Resultados
Discusión.
Para comenzar a evaluar todos lo métodos, se tomó la ecuación (1) dada por el problema, en esta se
2 3𝑟−ℎ
hizo un despeje teniendo asi la ecuacion πℎ ( 3
)-V=0 (2) esto por facilidad en los cálculos,
seguido de esto y con los datos entregados se hace el reemplazo en (2) y se tiene que
2 3(2.08)−ℎ
3. 1416ℎ ( 3
)-30=0 (3) pero podemos decir que esta ecuación es la función f que depende de
2 3(2.08)−ℎ
h, f(h)=3. 1416ℎ ( 3
)-30 (4).
Se procede a crear la tabla en excel y a seguir los pasos del método de Bipartición o Bisección
explicado en Solución de ecuación no lineal ☺ Método de Bisección/Bipartición ☺. Ejemplo en
Excel(2019, 1 diciembre)YouTube (Enlace: Solución de ecuación no lineal ☺ Método de
Bisección/Bipartición ☺ Ejemplo en Excel).
Para el ejercicio se toma valores iniciales 𝑎0=2 y 𝑏0 =4, al realizar 20 iteraciones, se observa que F(𝑎0
) es negativo, F(ℎ0) es positivo y F(𝑏0) es positivo, es decir, el cambio de signo (la raíz) se presenta
dentro del subintervalo [𝑎0,𝑥0] ([2, 3]) y, por lo tanto, el subintervalo [ℎ0, 𝑏0] ([3, 4]) puede
descartarse. De este modo, el nuevo intervalo de búsqueda [𝑎1, 𝑏1] se puede definir como 𝑎1 = 2 ,
comenzando una nueva etapa de cálculo para este nuevo intervalo [𝑎1, 𝑏1] = [2, 3]. Sucesivamente se
siguen realizando los cálculos hasta que el intervalo de búsqueda sea pequeño lo que permite que se
encuentre la raíz con una aproximación razonable, siendo en la iteración 20 el valor de la altura de
2,951636 metros lo que cumple con las condiciones del problema y al graficar la función con ayuda
de un recurso informático como es symbolab se puede observar que los valores son semejantes.
Figura 1. Gráfica de la función del problema 1 en función de la altura. [tomada de:Calculadora
gráfica - Symbolab]
Para hallar la raíz que soluciona el problema por el método de falsa posición se procede a crear la
tabla en excel y seguir los pasos del método explicado en Solución de ecuación no lineal ☺ Método
Falsa posición/Regla falsa(Regula falsi) ☺ Ejemplo en Excel(2019, 1 diciembre).YouTube. (Enlace:
Solución de ecuación no lineal ☺ Método Falsa posición/Regla falsa(Regula falsi) ☺ Ejemplo en
Excel ).
Para hallar la solución el problema se toma valores iniciales 𝑎0=-0,1 m y 𝑏0 =3,5 m, luego de digitar
los datos y las operaciones correspondientes se decide realizar 20 iteraciones para ver el
comportamiento. Se observa que en la sexta iteración la función evaluada en h (altura) se vuelve cero,
lo que significa que se encontró una raíz del problema es decir la solución ya que cumple con las
condiciones que presenta el problema 1. A partir de la séptima iteración los valores tienden a ser los
mismos en algunos casos y en otros cambia su valor a 0.
Para hallar la raíz que soluciona el problema se emplea la ecuación (4) y por el método de
sustituciones sucesivas o de punto fijo el cual suele ser de rápida convergencia, se procede a crear la
tabla en excel seguir los pasos del método explicado en Solución de ecuación no lineal ☺ Método
Sustituciones sucesivas/Punto fijo ☺ Ejemplo en Excel(2019, 2 diciembre).YouTube. (Enlace:
Solución de ecuación no lineal ☺ Método Sustituciones sucesivas/Punto fijo ☺ Ejemplo en Excel ).
Para hallar la solución el problema se toma como valor inicial ℎ0=-1 el cual no corresponde al
Se procede a crear la tabla en excel y a seguir los pasos del método de secante explicado en Solución
de ecuación no lineal ☺ Método Secante ☺ Ejemplo en Excel(2019, 2 diciembre)YouTube (Enlace:
Solución de ecuación no lineal ☺ Método Secante ☺ Ejemplo en Excel).
Para el ejercicio se toma valores iniciales 6 metros y 7 metros, al realizar 8 iteraciones el valor
evaluado en función de la altura es cero y su ℎ8 = 5, 16687828 𝑚𝑒𝑡𝑟𝑜𝑠 el cual es un valor que
corresponde a una raíz de la función, lo que quiere decir que es una solución para el problema pero
presenta el mismo problema con el anterior método y es que no puede ser la solución del problema
porque no cumple con la condición h ≤ 2r. Si se realizan más iteraciones aparece un error (#¡DIV/0!)
el cual quiere decir que la función al evaluarla será indeterminada.
Método Newton: para el desarrollo del método, se tomó 3 metros como valor inicial o valor semilla,
valor que se obtuvo reemplazando el radio escogido por el grupo de 2.08 m y el volumen requerido de
30 m3 para la construcción del tanque, este valor fue verificado con la gráfica obtenida en una
calculadora virtual, se encontró convergencia a las 3 iteraciones, se usó un criterio de convergencia
εf =1x10-6 y el resultado obtenido para la altura (h) fue de 2.9516 m, valor que cumple con el
parámetro h<2r definido en el problema propuesto.