INSTITUTO POLITÉCNICO NACIONAL
Unidad Profesional Interdisciplinaria de Biotecnología
MÉTODOS NUMÉRICOS
TAREA 2:
MÉTODO DE BISECCIÓN
GRUPO: 4BM2
PROFESORES:
Valdés Lozano Alejandra
EQUIPO 1
INTEGRANTES DEL EQUIPO:
ANGELES MONTES LIBERTAD
Juan Manuel Zamora Zúñiga
MARTÍNEZ ORTA MADAÍ
RIVERA ESQUIVEL DAVID
RODRÍGUEZ ALBARRÁN SAÚL
FECHA DE ENTREGA: 1 DE MARZO DE 2024
INTRODUCCIÓN
La determinación de las raíces de una ecuación es uno de los problemas más antiguos
en matemáticas, su importancia radica en que si podemos determinar las raíces de una
ecuación también podemos determinar máximos y mínimos, valores propios de
matrices, resolver sistemas de ecuaciones lineales y diferenciales.
Los métodos de aproximaciones sucesivas, son aquellos que a partir de una primera
aproximación a la solución de una ecuación y aplicando en forma repetida una fórmula
de recurrencia obtienen una mejor aproximación a la solución esperada.
El método de bisección es una técnica fundamental en métodos numéricos para
encontrar aproximaciones de raíces de ecuaciones no lineales. Su principio se basa en
el teorema del valor intermedio, que establece que si una función es continua en un
intervalo y cambia de signo en ese intervalo, entonces debe existir al menos una raíz
en ese intervalo.
En este caso el método de bisección se implementa mediante un proceso iterativo en el
que se divide sucesivamente un intervalo en dos mitades, evaluando la función en el
punto medio para determinar en qué mitad se encuentra la raíz. Este proceso se repite
hasta que se alcanza una aproximación suficientemente precisa de la raíz.
Directamente gracias a un sistema de programación puede implementar el método de
bisección de manera eficiente para resolver una variedad de problemas de búsqueda
de raíces en ecuaciones no lineales.
OBJETIVOS
Identificar y desarrollar los pasos del algoritmo de bisección para la solución de
ecuaciones no lineales de una
variable.
Calcular el error en el algoritmo de bisección para la solución de ecuaciones no lineales
de una variable.
CONTENIDO
APUNTE DE CLASE
Página | 2
Página | 3
Página | 4
i a r b F(a) F(r) F(b) Error
1 4.65 4.6750 4.70 + - - 0.0250
2 4.65 4.6625 4.6750 + - - 0.0125
3 4.65 4.6562 4.6625 + + - 0.0063
4 4.6562 4.6594 4.6625 + - - 0.0031
5 4.6562 4.6578 4.6594 + + - 0.0016
6 4.6578 4.6586 4.6594 + - - 0.0008
La mejor aproximación a la raíz es 4.66 con 0.001 de error máximo o de tolerancia
Página | 5
Página | 6
i a r b F(a) F(r) F(b) Error
1 37.73 37.7350 37.74 - + + 0.0050
2 37.7335 37.7325 37.74 - - + 0.0025
3 37.7325 37.7338 37.74 - - + 0.0012
4 37.7325 37.7344 37.74 - - + 0.0006
La mejor aproximación a la raíz es 37.73 m/s con 0.001 de error máximo o de
tolerancia
Página | 7
Página | 8
Página | 9
TAREA
1. Encuentre la raíz positiva más pequeña de las siguientes ecuaciones mediante el
método de bisección con una tolerancia de 0.001 y considerando la longitud del
intervalo indicado en cada caso.
a)
i a r b f(a) f(r) f(b) e
1 1.55 1.9427 1.55+pi/4 + - - 0.3927
2 1.55 1.7463 1.9427 + - - 0.1963
3 1.55 1.6482 1.7463 + - - 0.982
4 1.55 1.5991 1.6482 + - - 0.0491
5 1.55 1.5745 1.5991 + - - 0.0245
6 1.55 1.5623 1.5745 + + - 0.0123
7 1.5623 1.5684 1.5745 + + - 0.0061
8 1.5684 1.5715 1.5745 + - - 0.0031
9 1.5684 1.5699 1.5715 + + - 0.0015
10 1.5699 1.5715 1.5715 + + - 0.0008
La raíz positiva más pequeña debe ser de 1.57 con una tolerancia de 0.001
Página | 10
CÓDIGO UTILIZADO
clear
clc
f=inline('tan(x)-x+1');
f=vectorize(f);
v=-5:0.001:5;
plot(v,f(v));
grid
a=1.55;b=1.55+pi/4;tol=0.001;
r(1)=(a+b)/2
fa=f(a)
fr=f(r(1))
fb=f(b)
e(1)=(b-a)/2
b=r(1)
r(2)=(a+b)/2
fa=f(a)
fr=f(r(2))
fb=f(b)
e(2)=(b-a)/2
b=r(2)
r(3)=(a+b)/2
fa=f(a)
fr=f(r(3))
fb=f(b)
e(3)=(b-a)/2
b=r(3)
r(4)=(a+b)/2
fa=f(a)
fr=f(r(4))
fb=f(b)
e(4)=(b-a)/2
Página | 11
b=r(4)
r(5)=(a+b)/2
fa=f(a)
fr=f(r(5))
fb=f(b)
e(5)=(b-a)/2
b=r(5)
r(6)=(a+b)/2
fa=f(a)
fr=f(r(6))
fb=f(b)
e(6)=(b-a)/2
a=r(6)
r(7)=(a+b)/2
fa=f(a)
fr=f(r(7))
fb=f(b)
e(7)=(b-a)/2
a=r(7)
r(8)=(a+b)/2
fa=f(a)
fr=f(r(8))
fb=f(b)
e(8)=(b-a)/2
b=r(8)
r(9)=(a+b)/2
fa=f(a)
fr=f(r(9))
fb=f(b)
e(9)=(b-a)/2
a=r(9)
r(10)=(a+b)/2
Página | 12
fa=f(a)
fr=f(r(10))
fb=f(b)
e(10)=(b-a)/2
Página | 13
Página | 14
b)
i a r b f(a) f(r) f(b) e
1 0 0.5 1 - - + 0.5
2 0.5 0.75 1 - + + 0.25
3 0.5 0.6250 0.75 - + + 0.125
4 0.5 0.5625 0.6250 - + + 0.0625
5 0.5 0.5312 0.5625 - - + 0.0312
6 0.5312 0.5469 0.5625 - + + 0.0156
7 0.5312 0.5391 0.5469 - - + 0.0078
8 0.5391 0.5430 0.5469 - + + 0.0017
9 0.5391 0.5410 0.5430 - - + 0.0020
10 0.5410 0.5420 0.5430 - + + 0.001
La raíz positiva más pequeña debe ser de 0.54 m con una tolerancia de 0.001
Página | 15
CÓDIGO UTILIZADO
clc
clear
f=inline('sin(x)-0.3*exp(x)');
f=vectorize(f);
v=0:0.001:1;
plot(v,f(v));
grid
a=0;b=1;tol=0.001;
r(1)=(a+b)/2
fa=f(a)
fr=f(r(1))
fb=f(b)
e(1)=(b-a)/2
a=r(1)
r(2)=(a+b)/2
fa=f(a)
fr=f(r(2))
fb=f(b)
e(2)=(b-a)/2
b=r(2)
r(3)=(a+b)/2
fa=f(a)
fr=f(r(3))
fb=f(b)
e(3)=(b-a)/2
Página | 16
b=r(3)
r(4)=(a+b)/2
fa=f(a)
fr=f(r(4))
fb=f(b)
e(4)=(b-a)/2
b=r(4)
r(5)=(a+b)/2
fa=f(a)
fr=f(r(5))
fb=f(b)
e(5)=(b-a)/2
a=r(5)
r(6)=(a+b)/2
fa=f(a)
fr=f(r(6))
fb=f(b)
e(6)=(b-a)/2
b=r(6)
r(7)=(a+b)/2
fa=f(a)
fr=f(r(7))
fb=f(b)
e(7)=(b-a)/2
a=r(7)
r(8)=(a+b)/2
Página | 17
fa=f(a)
fr=f(r(8))
fb=f(b)
e(8)=(b-a)/2
b=r(8)
r(9)=(a+b)/2
fa=f(a)
fr=f(r(9))
fb=f(b)
e(9)=(b-a)/2
a=r(9)
r(10)=(a+b)/2
fa=f(a)
fr=f(r(10))
fb=f(b)
e(10)=(b-a)/2
Página | 18
c)
i a r b f(a) f(r) f(b) e
1 0.05 0.1 0.15 - - + 0.05
2 0.1 0.125 0.15 - + + 0.025
3 0.1 0.1125 0.125 - + + 0.0125
4 0.1 0.1063 0.1125 - + + 0.0062
5 0.1 0.1031 0.1063 - + + 0.0031
6 0.1 0.1016 0.1031 - - + 0.0016
7 0.1016 0.1023 0.1031 - + + 0.0008
La raíz positiva más pequeña debe ser de 0.1023 con una tolerancia de 0.001
CÓDIGO UTILIZADO
clc
clear
f=inline('16*x^5-20*x^3+x^2+5*x-0.5');
f=vectorize(f)
v=0:0.001:1;
plot(v,f(v))
grid
a=0.05;b=0.15;tol=0.001
n=ceil(log((b-a)/tol)/log(2))
r(1)=(a+b)/2
fa=f(a)
fr=f(r(1))
fb=f(b)
e=(b-a)/2
a=r(1)
r(2)=(a+b)/2
fa=f(a)
fr=f(r(2))
Página | 19
fb=f(b)
e(2)=(b-a)/2
b=r(2)
r(3)=(a+b)/2
fa=f(a)
fr=f(r(3))
fb=f(b)
e(3)=(b-a)/2
b=r(3)
r(4)=(a+b)/2
fa=f(a)
fr=f(r(4))
fb=f(b)
e(4)=(b-a)/2
b=r(4)
r(5)=(a+b)/2
fa=f(a)
fr=f(r(5))
fb=f(b)
e(5)=(b-a)/2
b=r(5)
r(6)=(a+b)/2
fa=f(a)
fr=f(r(6))
fb=f(b)
e(6)=(b-a)/2
a=r(6)
r(7)=(a+b)/2
fa=f(a)
fr=f(r(7))
fb=f(b)
e(7)=(b-a)/2
Página | 20
Página | 21
Página | 22
Página | 23
d)
i a r b f(a) f(r) f(b) e
1 0.35 0.4 0.45 - + + 0.05
2 0.35 0.375 0.4 - - + 0.025
3 0.375 0.3875 0.4 - + + 0.0125
4 0.375 0.3812 0.3875 - + + 0.0063
5 0.375 0.3781 0.3812 - - + 0.0031
6 0.3781 0.3797 0.3812 - + + 0.0016
7 0.3781 0.3789 0.3797 - + + 0.0008
La raíz positiva más pequeña debe ser de 0.37 con una tolerancia de 0.001
CÓDIGO UTILIZADO
clc
clear
f=inline('log(x)-0.2*x^2+1');
f=vectorize(f)
v=0:0.001:1;
plot(v,f(v))
grid
a=0.35;b=0.45;tol=0.001
n=ceil(log((b-a)/tol)/log(2))
r(1)=(a+b)/2
fa=f(a)
fr=f(r(1))
fb=f(b)
e=(b-a)/2
b=r(1)
r(2)=(a+b)/2
fa=f(a)
Página | 24
fr=f(r(2))
fb=f(b)
e(2)=(b-a)/2
a=r(2)
r(3)=(a+b)/2
fa=f(a)
fr=f(r(3))
fb=f(b)
e(3)=(b-a)/2
b=r(3)
r(4)=(a+b)/2
fa=f(a)
fr=f(r(4))
fb=f(b)
e(4)=(b-a)/2
b=r(4)
r(5)=(a+b)/2
fa=f(a)
fr=f(r(5))
fb=f(b)
e(5)=(b-a)/2
a=r(5)
r(6)=(a+b)/2
fa=f(a)
fr=f(r(6))
fb=f(b)
e(6)=(b-a)/2
b=r(6)
r(7)=(a+b)/2
fa=f(a)
fr=f(r(7))
fb=f(b)
Página | 25
e(7)=(b-a)/2
Página | 26
Página | 27
Página | 28
2. Se está diseñando un tanque esférico para almacenar agua para un poblado
pequeño en un país en desarrollo. El volumen del líquido que puede contener se
calcula con
Donde V= volumen (m3), h=profundidad del agua en el tanque (m), y R=radio
del tanque (m). Si R=3m, ¿a qué profundidad debe llenarse el tanque de modo
que contenga 30 m3? Maneje una tolerancia de 0.0001.
i a r b f(a) f(r) f(b) Error
1 2 2.0200 2.04 - - + 0.0200
2 2.0200 2.0300 2.04 - + + 0.0100
3 2.0200 2.0250 2.0300 - - + 0.0050
4 2.0250 2.0275 2.0300 - + + 0.0025
5 2.0250 2.0263 2.0275 - - + 0.0012
6 2.0263 2.0269 2.0275 - - + 0.0006
7 2.0269 2.0272 2.0275 - + + 0.0003
8 2.0269 2.0270 2.0272 - + + 0.0002
9 2.0269 2.0270 2.0270 - + + 0.0001
La profundidad del tanque debe ser de 2.027 m con una tolerancia de 0.0001
Página | 29
CÓDIGO UTILIZADO
clc
clear all
f=inline('(pi*h.^2)*(((3*3)-h)/3)-(30)');
f=vectorize(f);
h=-10:0.01:10;
plot(h,f(h))
grid on
a=2; b=2.04; tol=0.0001;
n=ceil(log((b-a)/tol)/log(2))
r(1)=(a+b)/2
fa=f(a)
fr=f(r(1))
fb=f(b)
e(1)=(b-a)/2
a=r(1)
r(2)=(a+b)/2
fa=f(a)
fr=f(r(2))
fb=f(b)
e(2)=(b-a)/2
b=r(2)
r(3)=(a+b)/2
fa=f(a)
fr=f(r(3))
fb=f(b)
Página | 30
e(3)=(b-a)/2
a=r(3)
r(4)=(a+b)/2
fa=f(a)
fr=f(r(4))
fb=f(b)
e(4)=(b-a)/2
b=r(4)
r(5)=(a+b)/2
fa=f(a)
fr=f(r(5))
fb=f(b)
e(5)=(b-a)/2
a=r(5)
r(6)=(a+b)/2
fa=f(a)
fr=f(r(6))
fb=f(b)
e(6)=(b-a)/2
a=r(6)
r(7)=(a+b)/2
fa=f(a)
fr=f(r(7))
fb=f(b)
e(7)=(b-a)/2
b=r(7)
Página | 31
r(8)=(a+b)/2
fa=f(a)
fr=f(r(8))
fb=f(b)
e(8)=(b-a)/2
b=r(8)
r(9)=(a+b)/2
fa=f(a)
fr=f(r(9))
fb=f(b)
e(9)=(b-a)/2
3. Por un canal trapezoidal fluye agua a una tasa de Q=20 m3/s. La profundidad
crítica y para dicho canal satisface la siguiente ecuación:
Donde g=9.8m/s2, Ac=área de la sección transversal (m2), y B=el ancho del canal
Página | 32
en la superficie (m). Para este caso, el ancho y el área de la sección transversal
se relacionan con la profundidad y por medio de
Resuelva para la profundidad crítica empleando un intervalo de longitud 1 y que
su error caiga por debajo del 0.0001.
i a r b f(a) f(r) f(b) e
1 1.5 1.55 1.6 - + + 0.05
2 1.5 1.525 1.55 - + + 0.025
3 1.5 1.5125 1.525 - - + 0.0125
4 1.5125 1.5187 1.525 - + + 0.0062
5 1.5125 1.5156 1.5187 - + + 0.0031
6 1.5125 1.5141 1.5156 - - + 0.0016
7 1.5141 1.5148 1.5156 - + + 0.0008
8 1.5141 1.5145 1.5148 - - + 0.0004
9 1.5145 1.5156 1.5148 - + + 0.0002
10 1.5145 1.5146 1.5156 - + + 0.0001
La profundidad para dicho canal será de 1.514 m con un error debajo de 0.0001
CÓDIGO
clc
clear
f=inline('(x.^6)+(18*x.^5)+(108*x.^4)+(216*x.^3)-((9600+3200*x)/9.8)');
f=vectorize(f)
v=0:0.01:2;
plot(v,f(v))
Página | 33
grid
a=1.5;b=1.6;tol=0.0001
n=ceil(log((b-a)/tol)/log(2))
r(1)=(a+b)/2
fa=f(a)
fr=f(r(1))
fb=f(b)
e=(b-a)/2
b=r(1);
r(2)=(a+b)/2
fa=f(a)
fr=f(r(2))
fb=f(b)
e(2)=(b-a)/2
b=r(2)
r(3)=(a+b)/2
fa=f(a)
fr=f(r(3))
fb=f(b)
e(3)=(b-a)/2
a=r(3)
r(4)=(a+b)/2
fa=f(a)
fr=f(r(4))
fb=f(b)
e(4)=(b-a)/2
b=r(4)
r(5)=(a+b)/2
fa=f(a)
fr=f(r(5))
fb=f(b)
e(5)=(b-a)/2
Página | 34
b=r(5)
r(6)=(a+b)/2
fa=f(a)
fr=f(r(6))
fb=f(b)
e(6)=(b-a)/2
a=r(6)
r(7)=(a+b)/2
fa=f(a)
fr=f(r(7))
fb=f(b)
e(7)=(b-a)/2
b=r(7)
r(8)=(a+b)/2
fa=f(a)
fr=f(r(8))
fb=f(b)
e(8)=(b-a)/2
a=r(8)
r(9)=(a+b)/2
fa=f(a)
fr=f(r(9))
fb=f(b)
e(9)=(b-a)/2
b=r(9)
r(10)=(a+b)/2
fa=f(a)
fr=f(r(10))
fb=f(b)
e(10)=(b-a)/2
b=r(10)
r(11)=(a+b)/2
Página | 35
fa=f(a)
fr=f(r(11))
fb=f(b)
e(11)=(b-a)/2
Página | 36
Página | 37
Página | 38
Página | 39
Página | 40
Página | 41
4. La corriente I en amperes cuando se aplica un voltaje E de una batería a través de
un varistor es
Hallar el voltaje necesario para producir una corriente de
a) 0.750 A
Trabaje con una tolerancia de 0.001 y proponga la longitud del intervalo inicial.
CÓDIGO
clc
clear
f=inline('0.0001*exp(5) - 0.1273*exp(2) + x-0.750');
f=vectorize(f);
v=-10:0.001:10;
plot(v,f(v))
grid
a=1.6;b=1.8;tol=0.001;
r(1)=(a+b)/2
fa=f(a)
fr=f(r(1))
fb=f(b)
e(1)=(b-a)/2
b=r(1)
r(2)=(a+b)/2
fa=f(a)
fr=f(r(2))
Página | 42
fb=f(b)
e(2)=(b-a)/2
a=r(2)
r(3)=(a+b)/2
fa=f(a)
fr=f(r(3))
fb=f(b)
e(3)=(b-a)/2
a=r(3)
r(4)=(a+b)/2
fa=f(a)
fr=f(r(4))
fb=f(b)
e(4)=(b-a)/2
b=r(4)
r(5)=(a+b)/2
fa=f(a)
fr=f(r(5))
fb=f(b)
e(4)=(b-a)/2
b=r(5)
r(6)=(a+b)/2
fa=f(a)
fr=f(r(6))
fb=f(b)
e(4)=(b-a)/2
Página | 43
b=r(6)
r(7)=(a+b)/2
fa=f(a)
fr=f(r(7))
fb=f(b)
e(4)=(b-a)/2
b=r(7)
r(8)=(a+b)/2
fa=f(a)
fr=f(r(8))
fb=f(b)
e(4)=(b-a)/2
Página | 44
i a r b f(a) f(r) f(b) e
1 1.6 1.7 1.8 - + + 0.1
2 1.6 1.65 1.7 - - + 0.05
3 1.65 1.6750 1.7 - - + 0.0250
4 1.6750 1.6875 1.7 - + + 0.0125
5 1.6750 1.6813 1.6875 - + + 0.0062
6 1.6750 1.6781 1.6813 - + + 0.0031
7 1.6750 1.6766 1.6781 - + + 0.0016
8 1.6750 1.6758 1.6766 - - + 0.0008
La raíz positiva más pequeña debe ser de 1.6758 con una tolerancia de 0.001
Página | 45
Página | 46
b) 1.50 A
clc
clear
f=inline('0.0001*exp(5) - 0.1273*exp(2) + x-1.5');
f=vectorize(f);
v=-10:0.001:10;
plot(v,f(v))
grid
a=2.42;b=2.43;tol=0.001;
r(1)=(a+b)/2
fa=f(a)
fr=f(r(1))
fb=f(b)
e(1)=(b-a)/2
a=r(1)
r(2)=(a+b)/2
fa=f(a)
fr=f(r(2))
fb=f(b)
e(2)=(b-a)/2
b=r(2)
r(3)=(a+b)/2
fa=f(a)
fr=f(r(3))
fb=f(b)
e(2)=(b-a)/2
b=r(3)
r(4)=(a+b)/2
fa=f(a)
fr=f(r(4))
fb=f(b)
e(4)=(b-a)/2
Página | 47
i a r b f(a) f(r) f(b) e
1 2.42 2.4250 2.43 - - + 0.0050
2 2.4250 2.4275 2.43 - + + 0.0025
3 2.4250 2.4263 2.4275 - + + 0.0013
4 2.4250 2.4256 2.4263 - - + 0.0006
La raíz positiva más pequeña debe ser de 2.4256 con una tolerancia de 0.001
Página | 48
Página | 49
CONCLUSIONES
Con los 2 métodos cerrados para la búsqueda de raíces pudimos identificar sus
desventajas y ventajas que tienen cada uno de estos, siendo el método gráfico la base
para el método de bisección ya que es importante saber identificar el intervalo donde se
encuentra nuestra raíz .
Al analizar tanto el método gráfico como el método de bisección para la búsqueda de
raíces, podemos afirmar que el método gráfico es fundamental para identificar el
intervalo inicial de búsqueda, mientras que la bisección, a pesar de su relativa lentitud,
ofrece la ventaja de asegurar convergencia a una solución dentro de ese intervalo. Esto
hace que la bisección sea una elección apropiada cuando se busca una solución
aproximada con alta certeza en la resolución de ecuaciones no lineales.
Este método se utiliza para encontrar raíces de ecuaciones no lineales al dividir
iterativamente un intervalo en dos partes y determinar en cuál de ellas cambia el signo
de la función. Aunque puede ser relativamente lento en comparación con otros
métodos, la bisección tiene la ventaja de garantizar convergencia a una solución dentro
del intervalo inicial dado, lo que lo hace adecuado cuando se necesita una solución
aproximada con un alto grado de certeza.
El método de la biseccion debe tener un intervalo de valores definido sobre el cual
realizar la búsqueda de las raices, antes de encontrar el punto medio del intervalo, se
debe comprobar que exista una raíz en el intervalo dado. Esto se obtiene al evaluar
ambos extremos del intervalo en la función y multiplicándolos. Si el resultado es menor
a cero quiere decir que la función cruza en el eje de las abscisas, por lo tanto una
existe una raiz. Confirmado eso, el método puede seguir y de lo contrario, un nuevo
intervalo debe ser seleccionado.
BIBLIOGRAFÍA
Iriarte, R., Borras, H., Duran, R. (2018) Apuntes de Métodos
Numéricos. Facultad de Ingeniería. UNAM.
Nakamura,S. (1992).Métodos Numéricos Aplicados con [Link].
[Link] Prentice Hall.
Página | 50
