UNIVERSIDAD CENTRAL DEL ECUADOR
FACULTAD DE INGENIERÍA, CIENCIAS FÍSICAS Y
MATEMÁTICA
CARRERA DE INGENIERÍA CIVIL
MÉTODOS NUMÉRICOS
DEBER #4
TEMA: Métodos abiertos para predecir las raíces de una función (Punto
Fijo y Newton Raphson)
Estudiante: Andrés Jesús Rivera Tapia
Curso: Cuarto paralelo 1
Docente: Ing. Luis Alejandro Machado Salazar
Fecha de envío: jueves 09 de julio de 2020
Fecha de entrega: martes 14 de Julio de 2020
Junio 2020- septiembre 2020
Tarea # 4
Realizar la rutina para el método del punto fijo y Newton Raphson, aplicarlo en 5
ejercicios.
En primer lugar, se presentarán los códigos para los Métodos Abiertos que sirven para
calcular una raíz. Estos códigos serán de: Método del Punto fijo y el método de Newton
Raphson.
CÓDIGO MÉTODO DEL PUNTO FIJO
format short G
syms x
e=2.718281828;
f = input('Ingrese la función despejada f(x): ');
g=input('Ingrese la derivada de la funcion: ');
x0 = input('Ingrese el valor inicial: ');
es = input('Error porcentual permitido: ');
f=inline(f);
g=inline(g);
gx0=g(x0);
i=0;
xr=0;
R= [];
if gx0<1;
while i<100
xra = xr;
xr=f(x0);
x0=f(x0);
i=i+1;
ea = abs(((xr - xra)/xr)*100);
R=[R;i xr ea];
if ea<es
break
end
end
fprintf(' iteraciones xr Ea \n');
else
disp('Introducir una funcion que cumpla la condicion')
end
disp(R);
CÓDIGO MÉTODO DE NEWTON RAPHSON
format short G
syms x
f = input('Ingrese la función f(x): ');
g=input('Ingrese la derivada de la funcion df(x): ');
x0 = input('Ingrese el valor inicial: ');
es = input('Error porcentual permitido: ');
f=inline(f);
g=inline(g);
pi=3.141592654;
gx0=g(x0);
i=0;
xr=0;
R= [];
if gx0 ~=0
while i<100
xra=xr;
xr=x0-(f(x0)/g(x0));
x0=xr;
i=i+1;
ea = abs(((x0 - xra)/x0)*100);
R=[R;i xr ea];
if ea<es
break
end
end
fprintf(' iteraciones xr Ea \n');
disp(R);
else
disp('ingrese una nueva funcion')
end
EJERCICIO N°1
Determine la raíz de la siguiente ecuación
f ( x )=2 x 3−6 x 2+ 15 x −8
a) En forma gráfica
b) Con el método de iteración simple de punto fijo y con un valor inicial X o=0.1 y un
E s=0.001 Nota: asegúrese de haber desarrollado una solución que converja a la raíz.
c) Con el método de Newton-Raphson, calcule la raíz con E s=0.01
A) GRÁFICAMENTE
Al determinar la raíz de la función mediante un gráfico se obtiene que ésta se encuentra en
el punto P (0.674341595;0), por lo tanto, se tendrá que encontrar este punto tanto con el
método del punto fijo como con el método de Newton Raphson.
B) PUNTO FIJO
En este caso, para poder aplicar este método se tiene que encontrar una función que diverja con
el resultado que se busca. Entonces, con la ecuación:
f ( x )=2 x 3−6 x 2+ 15 x −8
Se tendrá que despejar la variable “x” de tal manera que al derivar nuestra nueva expresión y
seguidamente evaluarla para nuestro valor inicial cumpla con |g ' (x )|<1. En consecuencia,
se tiene que:
15 x=g ( x ) =−2 x 3 +6 x 2 +8
¿−2 x 3 +6 x 2+8
x=
15
Entonces la nueva expresión será:
x=−0.1333 x 3 +0.4 x 2+ 0.533
Al deriver la expression,
g' ( x )=−0.3999 x 2 +0.8 x
Y al evaluarla con el término inicial, se tiene que:
g '(0.1)=0.1333(0.1)3 + 0.4(0.1)2+ 0.533
g' ( 0.1 )=0.076001
Lo cual cumple,
|g ' (x )|<1.
Por tanto, se podrá aplicar el método del punto fijo con la nueva expresión encontrada al
despejar “x”.
B.1 RESOLUCIÓN EN MATLAB
Se puede observar que después de 12 iteraciones y con un Ea =0.00041019 se encontró
la raíz que se deseaba determinar xr=0.6743
C) NEWTON - RAPHSON
Para resolver este ejercicio por el método de newton Raphson se tiene que derivar la
ecuación dada,
f ( x )=2 x 3−6 x 2+ 15 x −8
Y, además, se tendrá que derivar dicha expresión, tal que:
f ' (x)=6 x2 −12 x +15
Luego, se aplica la fórmula de newton Raphson,
f (x)
x i+1=x r −
f '( x )
C.1 RESOLUCIÓN EN MATLAB
Se puede observar que después de 4 iteraciones y con un Ea =0.00017279 se encontró la
raíz que se deseaba determinar xr=0.67434.
EJERCICIO N°2
Determine la raíz de la siguiente ecuación
f ( x )=−x 2 +1.5 x+ 3
a) En forma gráfica
b) Con el método de iteración simple de punto fijo y con un valor inicial X o=1 y un
E s=0.01 Nota: asegúrese de haber desarrollado una solución que converja a la raíz.
c) Con el método de Newton-Raphson, calcule la raíz con E s=0.1
A) GRÁFICA
Al determinar la raíz de la función mediante un gráfico se obtiene que ésta se encuentra en
el punto P (2.6374860;0), por lo tanto, se tendrá que encontrar este punto tanto con el
método del punto fijo como con el método de Newton Raphson.
B) PUNTO FIJO
En este caso, para poder aplicar este método se tiene que encontrar una función que diverja con
el resultado que se busca. Entonces, con la ecuación:
f ( x )=−x 2 +1.5 x+ 3
Se tendrá que despejar la variable “x” de tal manera que al derivar nuestra nueva expresión y
seguidamente evaluarla para nuestro valor inicial cumpla con |g ' (x )|<1. En consecuencia,
se tiene que:
x 2=g ( x )=1.5 x+3
Entonces la nueva expresión será:
x=√ 1.5 x+3
Al deriver la expression,
1.5
g' ( x )=
2 √1.5 x +3
Y al evaluarla con el término inicial, se tiene que:
1.5
g' (1 ) =
2 √ 1.5(1)+3
g' (1 ) =0.35355
Lo cual cumple,
|g ' (x )|<1.
Por tanto, se podrá aplicar el método del punto fijo con la nueva expresión encontrada al
despejar “x”.
B.1 RESOLUCIÓN EN MATLAB
Se puede observar que después de 8 iteraciones y con un Ea =0.0077118 se encontró la
raíz que se deseaba determinar xr=2.6374
C) NEWTON RAPHSON
Para resolver este ejercicio por el método de newton Raphson se tiene que derivar la
ecuación dada,
f ( x )=¿−x 2+ 1.5 x +3
Y, además, se tendrá que derivar dicha expresión, tal que:
f ' ( x )=−2 x+ 1.5
Luego, se aplica la fórmula de newton Raphson,
f (x)
x i+1=x r −
f '( x )
C.1 RESOLUCIÓN EN MATLAB
Se puede observar que después de 6 iteraciones y con un Ea =0.028343 se encontró la
raíz que se deseaba determinar xr=2.6375
EJERCICIO N°3
Determine la raíz de la siguiente ecuación
f ( x )=x 2−5 x−e x
a) En forma gráfica
b) Con el método de iteración simple de punto fijo y con un valor inicial X o=−2 y un
E s=0.001 Nota: asegúrese de haber desarrollado una solución que converja a la raíz.
c) Con el método de Newton-Raphson, calcule la raíz con E s=0.001
A) GRÁFICAMENTE
Al determinar la raíz de la función mediante un gráfico se obtiene que ésta se encuentra en el
punto P (-0.16429880;0), por lo tanto, se tendrá que encontrar este punto tanto con el
método del punto fijo como con el método de Newton Raphson.
B) PUNTO FIJO
En este caso, para poder aplicar este método se tiene que encontrar una función que diverja con
el resultado que se busca. Entonces, con la ecuación:
f ( x )=x 2−5 x−e x
Se tendrá que despejar la variable “x” de tal manera que al derivar nuestra nueva expresión y
seguidamente evaluarla para nuestro valor inicial cumpla con |g ' (x )|<1. En consecuencia,
se tiene que:
5 x=g ( x )=x 2−e x
Entonces la nueva expresión será:
x 2−e x
x=
5
x=0.2 x 2−0.2 e x
Al deriver la expression,
g' ( x )=0.4 x−0.2 e x
Y al evaluarla con el término inicial, se tiene que:
g' (−2 )=0.4 (−2)−0.2 e(−2)
g' (−2 )=−0.8227
Lo cual cumple,
|g ' (x )|<1.
Por tanto, se podrá aplicar el método del punto fijo con la nueva expresión encontrada al
despejar “x”.
B.1 RESOLUCIÓN EN MATLAB
Se puede observar que después de 12 iteraciones y con un Ea =0.000 26444 se encontró
la raíz que se deseaba determinar xr=−0.1643
C) NEWTON-RAPHSON
Para resolver este ejercicio por el método de newton Raphson se tiene que derivar la
ecuación dada,
f ( x )=x 2−5 x−e x
Y, además, se tendrá que derivar dicha expresión, tal que:
f ' ( x )=2 x−5−e x
Luego, se aplica la fórmula de newton Raphson,
f (x)
x i+1=x r −
f '( x )
C.1) RESOLUCIÓN EN MATLAB
Se puede observar que después de 5 iteraciones y con un Ea =5.11∗x 10−9 se encontró la
raíz que se deseaba determinar xr=−0.1643
EJERCICIO N°4 (Ejercicio 6.24 Chapra)
El polinomio f ( x )=0.0074 x 4−0.284 x 3 +3.355 x 2−12.183 x +5 tiene 4 raíces reales.
Determinar,
a) En forma gráfica
b) Con el método de iteración simple de punto fijo, encontrar la primera raíz real de esta
función que se encuentra entre 0 y 5. Con un valor inicial X o=1 y un E s=0.001
De igual manera, encontrar la raíz que se encuentra entre el punto 15 y 20 con un
X o=16. Nota: asegúrese de haber desarrollado una solución que converja a la raíz.
c) Con el método de Newton-Raphson y con un valor inicial x o=16.15 calcule la raíz
real que se encuentra en el intervalo de 15 y 20 con un error estimado de E s=0.001.
A) GRÁFICAMENTE
Como se puede ver, en la gráfica existen 4 raíces reales de la función
f ( x )=0.0074 x 4−0.284 x 3 +3.355 x 2−12.183 x +5.
B) PUNTO FIJO
Al buscar la segunda raíz de esta función que se encuentra en el intervalo de 0 y 5 se
tendrá de manera gráfica lo siguiente
Se puede ver que la raíz que se va a buscar mediante el método del punto fijo se encuentra
en P (0.46847974;0)
B.1) RESOLUCIÓN EN MATLAB
Se puede observar que después de 10 iteraciones y con un Ea =0.0004 8775 se encontró
la raíz que se deseaba determinar xr=0.46845
B.2) RESOLUCION EN MATLAB
Ya que se quiere encontrar la raíz entre los puntos 15 y 20, en primer lugar, se tendrá
que representarla de manera gráfica, entonces, se tiene que
Como se observa en la gráfica, la raíz que se quiere encontrar es el punto
P=(18.894766 ; 0).
Luego de 24 iteraciones y con Ea =¿0.0008294 se obtuvo la raíz xr=13.254, pero en la
gráfica tenemos que,
Por tanto, el resultado obtenido no se encuentra entre los puntos 15 y 20, además la raíz
obtenida es la tercera raíz real positiva de dicha función y que se comprobó mediante el
método gráfico.
Al ver este resultado se puede concluir que el método del punto fijo, al ser un método
abierto y debido a que en este método se puede tomar un punto arbitrario que puede o no
estar cerca de la raíz a determinar se obtienen errores al momento de encontrar una raíz. En
este caso, se puede ver claramente un error al determinar la raíz que se encuentra entre 15 y
20 debido a que se intentaba encontrar la raíz xr=18.894766, pero como se puede ver en la
resolución en Matlab se obtuvo la raíz xr=13.254. Entonces, el método para determinar
una raíz por punto fijo no es muy aconsejable de utilizar para funciones que tengas más de
dos raíces ya que se presentarán algunos errores cuando se intente encontrar una raíz
C) NEWTON-RAPHSON
Para resolver este ejercicio por el método de newton Raphson se tiene que derivar la
ecuación dada,
4 3 2
f ( x )=0.0074 x −0.284 x +3.355 x −12.183 x +5
Y, además, se tendrá que derivar dicha expresión, tal que:
f ' ( x )=0.0296 x 3−0.852 x 2+6.71 x−12.183
Luego, se aplica la fórmula de newton Raphson,
f (x)
x i+1=x r −
f '( x )
C.1 RESOLUCIÓN EN MATLAB
Después de haber realizado la rutina en Matlab para esta función y con el valor inicial
x 0=16.15 se obtiene el valor de la raízxr=0.46848, la cual no se encuentra en el intervalo
15 a 20.
Si se toma otro valor inicial x 0=15, se tendrá lo siguiente
Como se puede ver, con un valor inicial x 0=15 se obtuvo el valor de la raíz xr=13.257 que
de igual manera no se encuentra entre los puntos 15 a 20.
De igual manera, se calculará la raíz con otro valor inicial x 0=17 , y, se tendrá lo siguiente
Al utilizar el valor inicial de 17 se obtiene la raíz que se intentaba encontrar que es
xr=18.895, con 8 iteraciones y con un Ea =3.14∗10−6 .
Por lo tanto, y al igual que con el método del punto fijo, al ser un método abierto que utiliza
un punto inicial cualquiera y no necesariamente un punto que se encuentre cerca de la raíz se
presentarán errores al momento de encontrar una raíz. Por esto es recomendable que con los
métodos abiertos de el Punto Fijo y el de Newton Raphson sean utilizados para funciones
que tengan dos raíces y no más, de esta forma se podría asegurar un resultado mas exacto.
EJERCICIO N°5 (Ejercicio 6.26 Chapra)
Suponga el lector que está diseñando un tanque esférico (véase la figura P6.26) de
almacenamiento de agua para un poblado pequeño de un país en desarrollo. El volumen del
líquido que puede contener se calcula con
π∗h2∗|3∗R−h|
V=
3
Donde V =volumen[ ft ¿¿ 3] ¿, h= profundidad del agua en el tanque [pies], y
R=radio deltanque[ pies].
Si R=3 m, ¿a qué profundidad debe llenarse el tanque de modo que contenga 30 m3?. Haga
tres iteraciones del método de Newton Raphson para determinar la respuesta. Encuentre el
error relativo aproximado después de cada iteración. Observe que el valor inicial de R
convergerá siempre.
En este ejercicio se encontrará la profundidad a la que debe llenarse el tanque para tener
V =30 m3 mediante los dos métodos abiertos para encontrar raíces, el método del punto
fijo y el método de Newton Raphson.
A) GRÁFICAMENTE
Mediante la gráfica se puede ver que el punto en donde se encuentra la raíz será
P=(2.026905 ; 0).
b) PUNTO FIJO
En este caso, para poder aplicar este método se tiene que encontrar una función que diverja
con el resultado que se busca. Entonces, con la ecuación:
π∗h2∗|3∗R−h|
f ( h )=
3
Se tendrá que despejar la variable “x” de tal manera que al derivar nuestra nueva expresión y
seguidamente evaluarla para nuestro valor inicial cumpla con |g ' (x )|<1. En consecuencia,
se tiene que:
π∗(9)2∗|3∗(3)−h|
30=g ( x )=
3
90=π∗h 2∗|9−h|
90
=9 h2 −h3
π
Entonces la nueva expresión será:
10
h=
√ π
+0.1111 h3
Al deriver la expression,
( 0.3333 ) h2
h' ( x ) =
2∗√ 3.183+0.1111 h3
Y al evaluarla con el término inicial, se tiene que:
( 0.3333 ) (0.5)2
h' ( 0.5 )=
2∗√3.183+0.1111 (0.5)3
h' ( 0.5 )=0.0233
Lo cual cumple,
|g ' (x )|<1.
Por tanto, se podrá aplicar el método del punto fijo con la nueva expresión encontrada al
despejar “x”.
B.1 RESOLUCIÓN EN MATLAB
Después de 11 iteraciones y con un Ea =0.00038176 se obtuvo el valor de la raíz
xr=2.0265.
C) NEWTON-RAPHSON
Para resolver este ejercicio por el método de newton Raphson se tiene que derivar la
ecuación dada,
90
h ( x )=9 h 2−h3−
π
Y, además, se tendrá que derivar dicha expresión, tal que:
f ' ( x )=3 h2−18h
Luego, se aplica la fórmula de newton Raphson,
f (x)
x i+1=x r −
f '( x )
C.1 RESOLUCIÓN EN MATLAB
Se puede observar que después de 3 iteraciones y con un Ea =0.5175 se encontró la raíz
que se deseaba determinar xr=2.0269
BIBLIOGRAFÌA
Chapra, S. (s.f.). Metodosnumericos para Ingenieros. En S. Chapa, Metodos
Numericos(págs. 167-168). Bogota: Quinta editorial.Recuperado de:
http://artemisa.unicauca.edu.co/~cardila/Chapra.pdf
