UNIVERSIDAD CATÓLICA SANTO

TORIBIO DE MOGROVEJO
MÉTODOS ITERATIVOS DE APROXIMACIÓN PARA
ECUACIONES NO LINEALES

MÉTODOS NUMÉRICOS

ALUMNA: FARROÑAN SANCHEZ, JHEIMY JANELLY

DOCENTE: VILLEGAS SANTAMARIA, LUIS MIGUEL

SEPTIEMBRE 2022, CHICLAYO

1
 Metodo de Biseccion

Ejercicio: en el siguiente ejercicio tenemos una viga uniforme sujeta a una

carga distribuida uniformemente que crece en forma lineal. La ecuación para
la curva elástica resultante es la siguiente:
ω0
−X 5 + 2L2 X 3 − L4 X


y=
120EIL
Para determinar el punto de máxima flexión (es decir, el valor de x donde
dx/dy = 0. Después, sustituimos ya este valor en la ecuación a fin de determinar
en valor de la flexión de máxima. En sus cálculos, utilice los valores siguientes
para los parámetros: L = 5 cm, E = 50KN/cm2 , I = 300 cm4 y ω0 = 9 ∗
106 KN/cm
Datos: L = 5 cm
E = 50
I = 300
ω0 = (300)(120)(5)(50) = 9 ∗ 106
Solución:
9 106


−x5 + 50x3 − 625x


y=
120(5)(300)(5)
y = −x4 + 50x3 − 625x
Luego:
dy
= −5x4 + 150x2 − 625
dx
Hacer:
dy
= f (x) ⇒ f (x) = −5x4 + 150x2 − 625
dx
1. Hacer f (x) = 0

f (x) = 0
−5x + 150x2 − 625 = 0
4

3. Determinar el intervalo de confianza:

x ... 1 2(a) 3( b) 4 ...
F(x) ... −480 −105 320 495 ...

4. Los valores iniciales son a = (2) y b = (3), tales que f (2) = −105 < 0y
f (3) = 320 > 0 tiene signos opuestos, es decir:

• f (2) ∗ f (3) = (−105)(320) = −33600 < 0

5. La primera aproximación es:

2
 2+3
C0 = = 2.5 ⇒ f (2.5) = 117.187500 > 0
2
2 + 2.5
C1 = = 2.25 ⇒ f (2.25) = 6.230468 > 0
2
2 + 2.25
C2 = = 2.125 ⇒ f (2.125) = −49.610595 < 0
2
2.125 + 2.25
C3 = = 2.1875 ⇒ f (2.1875) = −21.715164 < 0
2
2.1875 + 2.25
C4 = = 2.21875 ⇒ f (2.21875) = −7.744603 < 0
2
2.21875 + 2.25
C5 = = 2.234375 ⇒ f (2.234375) = −0.757122 < 0
2
2.234375 + 2.25
C6 = = 2.2421875 ⇒ f (2.242187) = 2.736723 > 0
2
2.234375 + 2.2421875
C7 = = 2.238281 ⇒ f (2.238281) = 0.989805 > 0
2
2.234375 + 2.238281
C8 = = 2.236328 ⇒ f (2.236328) = 0.116341 > 0
2
2.234375 + 2.236328
C9 = = 2.235351 ⇒ f (2.235351) = −0.320390 > 0
2
2.235351 + 2.236328
C10 = = 2.235839 ⇒ f (2.235839) = −0.102024 < 0
2
iter a b c f(a) F(c) producto error
0 2 3 2.5 (−) (+) (−) −...
1 2 2.5 2.25 (−) (+) (−) 11.11%
2 2 2.25 2.125 (−) (−) (+) 5.88%
3 2.125 2.25 2.1875 (−) (−) (+) 2.85%
4 2.1875 2.25 2.21875 (−) (−) (+) 1.40%
5 2.21875 2.25 2.234375 (−) (−) (+) 0.83%
6 2.234375 2.25 2.242187 (−) (+) (−) 0.34%
7 2.234375 2.242187 2.238281 (−) (+) (−) 0.17%
8 2.234375 2.238281 2.236328 (−) (+) (−) 0.08%
9 2.234375 2.236328 2.235351 (−) (−) (+) 0.04%
10 2.235351 2.236328 2.235839 (−) (−) (+) 0.02

6. La raíz obtenida es: 2.235839

Método del punto fijo Ejercicio: en el siguiente ejercicio tenemos una viga uni-
forme sujeta a una carga distribuida uniformemente que crece en forma lineal.
La ecuación para la curva elástica resultante es la siguiente:
ω0
−X 5 + 2L2 X 3 − L4 X


y=
120EIL
Para determinar el punto de máxima flexión (es decir, el valor de x donde
dx/dy = 0. Después, sustituimos ya este valor en la ecuación a fin de determinar
en valor de la flexión de máxima. En sus cálculos, utilice los valores siguientes

3
para los parámetros: L = 5 cm, E = 50KN/cm2 , I = 300 cm4 y ω0 = 9 ∗
106 KN/cm
Datos:
L = 5 cm
E = 50
I = 300
ω0 = (300)(120)(5)(50) = 9 ∗ 106
Solución:
9 106


−x5 + 50x3 − 625x


y=
120(5)(300)(5)
y = −x4 + 50x3 − 625x
Luego:
dy
= −5x4 + 150x2 − 625
dx
Hacer:
dy
= g(x) ⇒ g(x) = −5x4 + 150x2 − 625
dx
1. Hacer g(x) = 0

g(x) = 0
−5x + 150x2 − 625 = 0
4

3. Determinar el intervalo de confianza:

x ... 1 2(a) 3( b) 4 ...
g(x) ... −480 −105 320 495 ...

4. Hallar los posibles despejes:

−5x4 + 150x2 − 625 = 0

r r
5x4 + 625 5x4 + 625 x3
⇒x= ⇒ g(x) = ⇒ g ′ (x) = q
150 150 4 +625
15 5x 150

5. Reemplazamos el punto medio del intervalo de confianza aquel que cumpla

con el criterio de convergencia ósea |g ′ (x0 ) ≤ 1| será apto para iterar
   
3 3



x

 ⇒ |g ′ (2.1)| =  q 2.1
 
   
|g ′ (x0 )| =  q 
 15 5x4 +625   15 5(2.1)4 +625 
  
150 150

= |0.281365| = 0.281365 ≤ 1
Apto para iterar

6. Iterar el despeje apto

4
 r
5(2.1)4 + 625
g (x1 ) ⇒ g(2.1) = = 2.194296
150
r
5(2.194296)4 + 625
g (x2 ) ⇒ g(2.194296) = = 2.222488
150
r
5(2.222488)4 + 625
g (x3 ) ⇒ g(2.222488) = = 2.231577
150
r
5(2.231577)4 + 625
g (x4 ) ⇒ g(2.231577) = = 2.234577
150
r
5(2.231577)4 + 625
g (x5 ) ⇒ g(2.231577) = = 2.235570
150
r
5(2.235570)4 + 625
g (x6 ) ⇒ g(2.235570) = = 2.235902
150
r
5(2.235902)4 + 625
g (x7 ) ⇒ g(2.235902) = = 2.236012
150
r
5(2.236012)4 + 625
g (x8 ) ⇒ g(2.236012) = = 2.236049
150
r
5(2.236049)4 + 625
g (x9 ) ⇒ g(2.236049) = = 2.236061
150
r
5(2.236061)4 + 625
g (x10 ) ⇒ g(2.236061) = = 2.236065
150
i xi g1 (xi ) = xi+1 |g1 (xi ) − g1 (xi−1 )|
1 x0 = 2.1 g1 (x1 ) = 2.194296 −...
2 x1 = 2.194296 g1 (x2 ) = 2.222488 0.048192
3 x2 = 2.222488 g1 (x3 ) = 2.231577 9.089∗10−3
4 x3 = 2.231577 g1 (x4 ) = 2.234575 2.998∗10−3
5 x4 = 2.234575 g1 (x5 ) = 2.235570 9.95∗ 10−4
6 x5 = 2.235570 g1 (x6 ) = 2.235902 3.32∗ 10−4
7 x6 = 2.235902 g1 (x7 ) = 2.236012 1.1∗ 10−4
8 x7 = 2.236012 g1 (x8 ) = 2.236049 3.7∗ 10−5
9 x8 = 2.236049 g1 (x9 ) = 2.236061 1.2∗ 10−5
10 x9 = 2.236061 g1 (x10 ) = 2.236065 4∗ 10−6

7. La raíz obtenida es: 2.236065

Método de falsa posición: Ejercicio: en el siguiente ejercicio tenemos una viga

uniforme sujeta a una carga distribuida uniformemente que crece en forma lineal.
La ecuación para la curva elástica resultante es la siguiente:
ω0
−X 5 + 2L2 X 3 − L4 X


y=
120EIL
Para determinar el punto de máxima flexión (es decir, el valor de x donde
dx/dy = 0. Después, sustituimos ya este valor en la ecuación a fin de determinar

5
en valor de la flexión de máxima. En sus cálculos, utilice los valores siguientes
para los parámetros: L = 5 cm, E = 50KN/cm2 , I = 300 cm4 y ω0 = 9 ∗
106 KN/cm
Datos: L = 5 cm
E = 50
I = 300
ω0 = (300)(120)(5)(50) = 9 ∗ 106
Solución:
9 106


−x5 + 50x3 − 625x


y=
120(5)(300)(5)
y = −x4 + 50x3 − 625x
Luego:
dy
= −5x4 + 150x2 − 625
dx
Hacer:
dy
= f (x) ⇒ f (x) = −5x4 + 150x2 − 625
dx
1. Hacer f (x) = 0

f (x) = 0
−5x + 150x2 − 625 = 0
4

3. Determinar el intervalo de confianza:

x ... 1 2(a) 3( b) 4 ...
F(x) ... −480 −105 320 495 ...

4. Los valores iniciales son a = (2) y b = (3), tales que f (2) = −105 < 0y

f (3) = 320 > 0 tiene signos opuestos, es decir:

• f (2) ∗ f (3) = (−105)(320) = −33600 < 0

5. La primera aproximación de la raíz se toma igual a:

6
 f (b)(b − a)
c=b−
f (b) − f (a)
(320)(1)
C0 = 3 − = 2.247058 ⇒ f (2.247058) = 4.915196 > 0
(320) − (−105)
(4.915196)(0.247058)
C1 = 2.247058 − = 2.236010
(4.915196) − (−105)
⇒ f (2.236010) = −0.025557 < 0
(4.915196)(0.011047)
C2 = 2.247058 − = 2.236067
(4.915196) − (−0.025557)
⇒ f (2.236067) = 2.02874 ∗ 10−7 > 0



2.02874 ∗ 10−7 5.7151 ∗ 10−5
C3 = 2.236067 − = 2.236067
(2.02874 ∗ 10−7 ) − (−0.025557)
⇒ f (2.236067) = 2.02874 ∗ 10−7 > 0



2.02874 ∗ 10−7 5.7151 ∗ 10−5
C4 = 2.236067 − = 2.236067
(2.02874 ∗ 10−7 ) − (−0.025557)
⇒ f (2.236067) = 2.02874 ∗ 10−7 > 0



2.02874 ∗ 10−7 5.7151 ∗ 10−5
C5 = 2.236067 − = 2.236067
(2.02874 ∗ 10−7 ) − (−0.025557)
⇒ f (2.236067) = 2.02874 ∗ 10−7 > 0



2.02874 ∗ 10−7 5.7151 ∗ 10−5
C6 = 2.236067 − = 2.236067
(2.02874 ∗ 10−7 ) − (−0.025557)
⇒ f (2.236067) = 2.02874 ∗ 10−7 > 0



2.02874 ∗ 10−7 5.7151 ∗ 10−5
C7 = 2.236067 − = 2.236067
(2.02874 ∗ 10−7 ) − (−0.025557)
⇒ f (2.236067) = 2.02874 ∗ 10−7 > 0



2.02874 ∗ 10−7 5.7151 ∗ 10−5
C8 = 2.236067 − = 2.236067
(2.02874 ∗ 10−7 ) − (−0.025557)
⇒ f (2.236067) = 2.02874 ∗ 10−7 > 0



2.02874 ∗ 10−7 5.7151 ∗ 10−5
C9 = 2.236067 − = 2.236067
(2.02874 ∗ 10−7 ) − (−0.025557)
⇒ f (2.236067) = 2.02874 ∗ 10−7 > 0



2.02874 ∗ 10−7 5.7151 ∗ 10−5
C9 = 2.236067 − = 2.236067
(2.02874 ∗ 10−7 ) − (−0.025557)
⇒ f (2.236067) = 2.02874 ∗ 10−7 > 0
1. Hallando el criterio de parada o el primer valor del error en %
 
 2.236010 − 2.247058 
 ∗ 100% = 0.49%
 2.236010

7
 iter a b c f (a) F (c) producto error
0 2 3 2.247058 (−) (+) (−) ···
1 2 2.247058 2.236010 (−) (−) (+) 0.49%
2 2.236010 2.236067 2.236067 (−) (+) (+) 2.54%
3 2.236010 2.236067 2.236067 (−) (+) (−) 0%
4 2.236010 2.236067 2.236067 (−) (+) (−) 0%
5 2.236010 2.236067 2.236067 (−) (+) (−) 0%
6 2.236010 2.236067 2.236067 (−) (+) (−) 0%
7 2.236010 2.236067 2.236067 (−) (+) (−) 0%
8 2.236010 2.236067 2.236067 (−) (+) (−) 0%
9 2.236010 2.236067 2.236067 (−) (+) (−) 0%
10 2.236010 2.236067 2.236067 (−) (+) (−) 0%
Método de Newton-Raphson Ejercicio: en el siguiente ejercicio tenemos una
viga uniforme sujeta a una carga distribuida uniformemente que crece en forma
lineal. La ecuación para la curva elástica resultante es la siguiente:
ω0
−X 5 + 2L2 X 3 − L4 X


y=
120EIL
Para determinar el punto de máxima flexión (es decir, el valor de x donde
dx/dy = 0. Después, sustituimos ya este valor en la ecuación a fin de determinar
en valor de la flexión de máxima. En sus cálculos, utilice los valores siguientes
para los parámetros: L = 5 cm, E = 50KN/cm2 , I = 300 cm4 y ω0 = 9 ∗
106 KN/cm
Datos:
L = 5 cm
E = 50
I = 300
ω0 = (300)(120)(5)(50) = 9 ∗ 106
Solución:
9 106


−x5 + 50x3 − 625x


y=
120(5)(300)(5)
y = −x4 + 50x3 − 625x
Luego:
dy
= −5x4 + 150x2 − 625
dx
Hacer:
dy
= f (x) ⇒ f (x) = −5x4 + 150x2 − 625
dx
1. Hacer f (x) = 0

f (x) = 0
−5x + 150x2 − 625 = 0
4

3. Determinar el intervalo de confianza:

8
 x ... 1 2(a) 3( b) 4 ...
F(x) ... −480 −105 320 495 ...
4. Los valores iniciales son a = (2) y b = (3), tales que f (2) = −105 < 0y
f (3) = 320 > 0 tiene signos opuestos, es decir:
• f (2) ∗ f (3) = (−105)(320) = −33600 < 0
3. Donde C0 = 3.6 ∈ [3, 4]
4. Hallar el criterio de convergencia donde a = (3)yb = (4)
• f (3) ∗ f (4) = (−0.959969)(1.385425) = −1.329966 < 0
• f (x) = 20x2 + 300(x) → f ′′ (x) = −60x2 + 300
′

Son continuas y no se anulan en [2; 3]

• Sea C0 = 2.1 ∈ [2; 3] ⇒ f ′ (2.1) ∗ f ′′ (2, 1) ⇒
= (718, 2)(35.4) = 760.20 > 0
Cumple con el criterio de convergencia
f (C0 )
5. Aplicamos el algoritmo Cn+1 = C0 − f ′ (C0 )

Sea:
f (C0 ) f (2.1)
C0 = 2.1 ⇒ C1 = C0 − = 2.1 − ′ = 2.236563
f ′ (C0 ) f (2.1)
f (2.236563)
C1 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)
f (2.236563)
C2 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)
f (2.236563)
C3 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)
f (2.236563)
C4 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)
f (2.236563)
C5 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)
f (2.236563)
C6 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)
f (2.236563)
C7 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)
f (2.236563)
C8 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)
f (2.236563)
C9 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)
f ′ (2.236563)
C10 = 2.236563 ⇒ C2 = − ′ = 2.236567
f (2.236563)

9
n Cn |Cn − Cn−1 |
0 C0 = 2.1 ···
1 C1 = 2.236567 1.136567
2 C2 = 2.236567 0.000000
3 C3 = 2.236567 0.000000
4 C4 = 2.236567 0.000000
5 C5 = 2.236567 0.000000
6 C6 = 2.236567 0.000000
7 C7 = 2.236567 0.000000
8 C8 = 2.236567 0.000000
La raíz obtenida es 2.236567

10
 CONCLUSIONES:
Para el ejercicio:

−5x4 + 150x2 − 625 = 0

Métodos Intervalo raíz funciona
Método de Bisección [2; 3] 2.235839 Si
Método de Punto fijo [2; 3] 2.236065 Si
Método de Newton Raphson [2; 3] 2.236067 Si
Método de falsa Posición [2; 3] 2.236067 si

11