UNIVERSIDAD NACIONAL DE

INGENIERÍA
FACULTAD DE INGENIERÍA ElÉCTRICA Y
ELECTRÓNICA

PROBLEMAS RESUELTOS PARA PRIMERA PRÁCTICA

Curso:Métodos Numéricos

Código de Curso: MA195

Sección: O

Profesor: Diaz Acha Israel Manuel

Estudiantes:Aguilar Taipe Luis Fernando 20170442I

Escobar Amanca Gerson Renato 20170328A
Loayza Acosta Denis Franco 20170056A

Ciclo: 2020 - I

LIMA-PERÚ
2020

1
Facultad de Ingeniería Eléctrica y Electrónica

Método de la Bisección y el Método de la Regla Falsa

Dada la ecuación siguiente:

x2
exp(x) = cos(x) + + 1 , x ∈ [0, 1]
2

a) Muestre la existencia y unicidad de la solución de dicha ecuación en el intervalo dado.

Existencia Para ∈ [0, 1]

02
f(0) = cos(0) + - e0 + 1 f(0) = 1 > 0
2

12
f(1) = cos(1) + - e1 + 1 f(1) = -0.67797 < 0
2

Entonces para que f (x) sea continua debe cumplir el teorema de Bolzano, sea f(x) una función
cotinua en el intervalo cerrado [a, b] tal que f (a).f (b) < 0, entonces nos permite afirmar que
existe al menos un punto c ∈ [a, b] tal que f (c) = 0
Si cumple con el teorema de Bolzano, entonces cumple la existencia de la función.

Unicidad
f 0 (x) = x - sen(x) - ex
Comprobamos si es una función inyectiva, es creciente
cuando x, y ∈ A, x<y ⇒ F(x) ≤ F(y), decreciente
cuando x,y ∈ A, x<y ⇒ F(x) ≤ F(y), es decir cuando
-F es creciente monotona cuando es creciente o decreciente.
f 0 (x) = x - sen(x) - ex vemos que para cada valor que damos a x f 0 (x) < 0 (Decreciente)
entonces es inyectiva, entonces cumple la unicidad f. Inyectiva unicidad

*) Calcule además el número de iteraciones necesarias por el Método de la Bisección para que
el error sea menor a 10−6 .
Sea n = Número de iteraciones necesarias por el método de la Bisección para que el error sea
menor a 10−6
x ∈ [0, 1] , a = 0 , b = 1

2
Facultad de Ingeniería Eléctrica y Electrónica

b−a
< 10−6
2n

1−0 1
n
= n < 10−6
2 2

2n > 106

log(2n ) > log(106 )

6
n>
log(2)

n > 19.931

Entonces a partir de 20 iteraciones saldrá la solución de la función.

*) Implemente un algoritmo en Matlab que permita visualizar la gráfica de la ecuación en el

intervalo dado.

Algoritmo en Matlab
% Algoritmo para graficar intersectando las 2 funciones
x=0:0.01:1
y1 = exp(x);
y2= cos(x) + 0.5*(x.∧ 2) + 1;

plot(x,y1,’r’)% Graficando la funcion y1

hold on;
plot(x,y2,’b’)% Graficando la funcion y2
legend(’y1’,’y2’);% Intersectando en el mismo plano y1 e y2

xlabel(’Eje x’);
ylabel(’Eje y’);

grid on;
hold off;

3
Facultad de Ingeniería Eléctrica y Electrónica

4
Facultad de Ingeniería Eléctrica y Electrónica

b) Implemente un algoritmo en Matlab que permita aproximar por el método de la Bisección,

la solución de la ecuación con 6 cifras decimales exactas.
Algoritmo en Matlab
a = input(’Ingrese intervalo inicial:’);
b = input(’Ingrese intervalo final:’);
y = input(’Ingrese la funcion f(x)=’,’s’);

%Ingresamos: exp(x) - cos(x) - 0.5*(x.∧ 2) − 1;

n= 20; % Número de iteraciones

f=inline(y);
for k=1:n

c=(a+b)/2;

e=abs((b-a)/2);

A(k,:)=[k a b c f(c) e];

if f(a)*f(c)<0

b=c;

else

a=c;

end

end

printf(’\n \tk \ta \tb \tc \tf(c) \terror \n’,’%f’)

5
Facultad de Ingeniería Eléctrica y Electrónica

disp(A)
printf(’Solucion:\n c=%6f\’,c)
printf(’f(c)=%6f\n’,f(c))
printf(’error=%6f\n’,e)

*) Muestre el algoritmo y los resultados de su ejecución mediante una tabla de datos detalla-
da. Verifique si el número de iteraciones en dicha aproximación de la solución coincide con el
número calculado en la parte (a).

6
Facultad de Ingeniería Eléctrica y Electrónica

Las iteraciones para esta parte son para n ≥ 20, pero en la parte de Matlab sale la solución
para n = 20.

¿Cuales son sus conclusiones con respecto a este método?

Concluimos que haciendo con Matlab es más rápido, porque hacer 20 iteraciones a mano nos
costaría más tiempo, porque el programa nos ayuda tener la solución en poco tiempo.

c) Implemente un algoritmo en Matlab que permita aproximar por el Método de la Regla Falsa,
la solución de la ecuación con 6 cifras decimales exactas.
Algoritmo en Matlab

7
Facultad de Ingeniería Eléctrica y Electrónica

function falsa(h,a,b,to1)

%% ’ exp(x) − cos(x) − 0.5 ∗ (x.∧ 2) − 1’ va en cadena

clc
f = inline(h);
c = (a*f(b)-b*f(a))/(f(b)-f(a));
n = 0;
error = abs(f(c));
fprintf(’\t n \t \t \ta \t \t c\t \t \t b \t \t error\n ’)

while(error>to1)

c= (a*f(b) - b*f(a))/(f(b)-f(a));

fprintf(’\t %d \t \t \t %f \t \t \t \t %f\t \t \t \t \t %f \t \t \t \t %f\n’, n,a,b,c,error’)

if(f(a)*f(c)<0)

b=c;

else a=c;

end

error=abs(f(c));

n=n+1;
end

*)Muestre el algoritmo y los resultados de la ejecución del algoritmo mediante una tabla deta-
llada.

8
Facultad de Ingeniería Eléctrica y Electrónica

d) Determine el numero de iteraciones que fueron necesarias en la parte (c) y compare con las
de la parte (b).
En la parte (c) el número de iteraciones para hallar la solución es de 8 iteraciones; pero en la
parte (b) el numero de iteraciones es de 20 iteraciones.

*) ¿Cual de los dos métodos fue más rápido y a que se debe eso? Describa las ventajas y
desventajas de un método con respecto al otro después de esta experiencia.

9
Facultad de Ingeniería Eléctrica y Electrónica

Método de la Bisección
Ventajas Desventajas
• Siempre resulta convergente • Converge de forma lenta.
• Brinda la posibilidad de hallar una raíz
• Permite establecer cierto límite de error.
aunque halla otras en el intervalo.
• Se necesita que f(x) sea de carácter conti- • No sirve para establecer raíces de natura-
nuo en el intervalo previsto. leza compleja.

Método de la Falsa Posición

Ventajas Desventajas
• Con respecto al método de la bi- • Aunque el método de la falsa posición parecía ser siem-
sección mejora la elección del in- pre la mejor opción de los que usan intervalos hay casos
tervalo. donde funcionan deficiente.
• Este método es tan seguro como
• Puede fallar en dar un intervalo que contenga a la raíz
el método de la bisección.

Método del Punto Fijo y el Método de Newton-Rapshon

Dada la función siguiente:
f (x) = 2x3 + x2 − 2x ,x ∈ [0, 1]
a) Muestre que la función f (x) es cóncava hacia arriba y creciente en el intervalo dado.
Solución:
f 0 (x) = 6x2 + 2x
Para hallar los puntos críticos hacemos f 0 (x) = 0
f 0 (x) = 6x2 + 2x = 0
x(6x + 2) = 0
1
x=0∨x=-
3

Si x ∈ < −∞, − 31 >

f 0 (x) = 6x2 + 2x > 0

Entonces f (x) es creciente en < −∞, − 13 >

10
Facultad de Ingeniería Eléctrica y Electrónica

Si x ∈ <- 13 , 0>

f 0 (x) = 6x2 + 2x < 0

Entonces f (x) es decreciente en <- 13 , 0 >

Si x ∈ < 0, +∞ >

f 0 (x) = 6x2 + 2x > 0

Entonces f (x) es creciente en < 0, +∞ >

*) Muestre que la ecuación f (x) = 0 tiene exactamente una solución en el intervalo dado.

Solución:
f 00 (x) = 12x + 2

Hallamos los puntos de inflexión en f 00 (x) = 0

1
f 00 (x) = 12x + 2 = 0 x = - (P.I)
6

Para x < − 16

f 00 (x) = 12x + 2 < 0

Entonces f (x) es cóncava hacia abajo en < −∞, − 21 >

Para x > - 16

f 00 (x) = 12x + 2 >0

Entonces f 0 (x) es cóncava hacia arriba en < − 61 , ∞ >

*) Implemente un algoritmo en Matlab que permita visualizar la gráfica de la función f(x) en

11
Facultad de Ingeniería Eléctrica y Electrónica

el intervalo dado. Muestre el algoritmo y la gráfica obtenida.

Solución:
Para demostrar que f (x) tiene solo una raíz en el intervalo [0, 1]

f (x) = 2x3 + 2x − 2 ; f(0) = -2 < 0

f (1) = 1 > 0

f (0).f (1) < 0

Puesto que f (x) es continua y diferenciable en todo x, entonces es continua en [0, 1] y diferen-
ciable en <0,1>.

Entonces existe c ∈ <0,1> tal que f (c) = 0.

Para demostrar que es única usamos el teorema Rolle, por contradicción supongamos que tiene
2 raíces a y b, entonces f (a) = f (b) y como es diferenciable en <a,b> por lo tanto existe un
número "c" entre a y b tal que f 0 (c) = 0, pero f 0 (x) > 0, entonces f 0 (x) no puede ser cero, lo
que da una contradicción, por lo tanto no tiene una raíz en ese intervalo.
Código en Matlab:
clear
clc
close all
x=[-100:1:100];
dim=get(0,’ScreenSize’);
figure(’position’,[dim(1) dim(2) dim(3) dim(4)]);
y=((2)*x.*x.*x)+(*x.*x)-2;
plot(x,y)
axis([-5 5 -20 180])
daspect([1 1 1])
grid on
title(’Grafica inciso a’)
xlabel(’eje x’) % etiqueta para el nombre del eje x

12
Facultad de Ingeniería Eléctrica y Electrónica

ylabel(’eje y’) % etiqueta para el nombre del eje y

13
Facultad de Ingeniería Eléctrica y Electrónica

…
x2 3
b) Muestre que g(x) = 1−
es una contracción en el intervalo [0, 1]. Muestre que g(x) es
2
cóncava hacia abajo y decreciente en dicho intervalo. Muestre que el punto fijo de g(x) es la
solución de f (x) = 0. Implemente un algoritmo en Maltab que permita visualizar la gráfica de
la ecuación g(x) = x en el intervalo dado. Muestre el algoritmo y la gráfica obtenida.

Solución: …
3 x2
Sea g(x) = 1− , x ∈ [0, 1] , se sabe que g(x) > 0 ∀ x ∈ [0, 1]
2

Luego:

|g(x) − g(y)| |g(x) − g(y)| |g 2 (x) + g(x)g(y) + g 2 (y)|

⇒
|x − y| |x − y| |g 2 (x) + g(x)g(y) + g 2 (y)|

|g 3 (x) − g 3 (y)|
⇒ , operando
|x − y||g 2 (x) + g(x)g(y) + g 2 (y)|
2
x2
|(1 − 2
) − (1 − y2 )|
x2 2/3 2 2 y 2 1/3
, simplificando
|x − y||(1 − 2
) + (1 − y2 )2/3 + ((1 − x2 )(1 − 2
)) |

1 |x − y||x + y|
⇒ x2 2/3 y 2 2/3 x2 y 2 1/3
2 |x − y||(1 − ) + (1 − ) + ((1 − )(1 − )) |
2 2 2 2

1 |x − y| 1 |x + y|
⇒ x2 2/3 y 2 2/3 x2 y 2 1/3
= ...(1)
2 |(1 − ) + (1 − ) + ((1 − )(1 − )) | 2 h(x, y)
2 2 2 2

Evaluando:
0 ≤ |x+y| ≤ 2 ...(2)

ã2/3
x2
Å
1
√
3
≤ 1− ≤ 1 ...(3)
4 2
ã2/3
y2
Å
1
√
3
≤ 1− ≤ 1 ...(4)
4 2
ã2/3 Å ã2/3
x2 y2
Å
1
√
3
≤ 1− 1− ≤ 1 ...(5)
4 2 2

Sumando (3),(4) y (5) se obtiene

3
√3
≤ |h(x,y)| ≤ 3
4

14
Facultad de Ingeniería Eléctrica y Electrónica

√
4
1 1 4
≤ ≤ ...(6)
3 |h(x, y)| 3

De (2) y (6) en (1)

√
3
1 |x + y| 4
0≤ ≤
2 |h(x, y)| 3

Se observa que:
|g(x) − g(y)|
0≤ ≤1
|x − y|

∴ Por lo tanto g(x) es continua.

Concavidad y decreciente

… ã−2/3
x2 x2
Å
3 0 1
g(x) = 1− ⇒ g (x) = - x 1 − , igualamos g 0 (x) = 0
2 3 2

1 Ä 2−x2 ä−2/3 √ √
⇒ g 0 (x) = 0 ⇒ − x 2 = 0, Puntos críticos x = { 0,- 2, 2}
3

Como x ∈ [0,1] ⇒ tomamos como punto crítico x = 0 ...(1)

ã−2/3 Å
x2 x2 −1
Å ã
00 1 2 2
Hallamos g (x) = − 1− 1 + x (1 − )
3
|{z} 2 3 2
| {z }| {z }
<0 >0 >0

Se concluye que para g(x) ∈ [0,1]; g”(x) <0

1
∴ Por lo que g(x) es cóncava hacia abajo y g(0) = 1, g(1) = √
3
, por ende g(0) es un máximo
2
relativo y es decreciente.

Punto Fijo

… de g(x) con la función identidad

Es la intersección
3 x2 x2
g(x) = x ⇒ 1 − ⇒1- = x3 ⇒ 2x3 + x2 − 2 = 0
2 2

∴ Por lo tanto se verifica que el punto fijo de g(x) es f(x) = 0

15
Facultad de Ingeniería Eléctrica y Electrónica

Código en Matlab
x = 0 : 0.001 : 1;
g = (1 − ((x.∧ 2)/2)).∧ (1/3);
y = x;
plot(x,g, ’r’);
grid on;
hold on;
plot(x,y, ’b’);
title(’Gráfica inciso b’);
xlabel(’eje x’);
ylabel(’eje g(x)’);

16
Facultad de Ingeniería Eléctrica y Electrónica

c) Implemente un algoritmo en Matlab que permita aproximar por el método del Punto fijo, la
solución de la ecuación g(x) = x con 6 cifras decimales exactas usando como punto de partida
a x0 = 1. Muestre el algoritmo y los resultados de la ejecución del algoritmo mediante una
tabla detallada.
Solución:
Código en Matlab:
function puntofijo

fprintf(’Punto Fijo ’);

gx=input(’Ingrese la funcion: ’, ’s’);
g=inline(gx);
x=input(’Ingrese el x inicial: ’);
maxiter=100;
error=input(’Ingrese el valor tolerado: ’);
i=1;
fprintf(’i \t x \t \t error \n’);

17
Facultad de Ingeniería Eléctrica y Electrónica

fprintf(’0 \t %f \n’,x);

while(i<=maxiter)
a=x;
x=g(a);
e=abs(x-a);
fprintf(’%d \t %1.6f \t \t %f \n’, i,x,e);

if(e<=error)
i=maxiter+1;
end
i=i+1;
end
fprintf(’La aproximacion es %1.6f con un error de: %f ’, i,x,e);

end
fprintf(’i \t x \t \t error \n’);

18
Facultad de Ingeniería Eléctrica y Electrónica

f (xn )
d) Explique de forma geométrica porque la sucesión xn+1 = xn - , x0 = 1, converge a la
f 0 xn
solución de la ecuación f (x) = 0 en el intervalo [0, 1].

Solución:
Para que converja ubicamos f (x0 ) en el gráfico, luego trazamos una recta tangente la cual pasa
por ele eje x muy cerca de una raíz de f (x), luego si tomamos otro punto entre 0 y 1 y trazamos
la tangente correspondiente se aproxima más a la raíz de f (x). Por lo tanto converge en [0, 1].
f (x0 )
f 0 (x0 ) =
x0 − x1

f (x0 )
x1 = x0 -
f 0 (x0 )

19
Facultad de Ingeniería Eléctrica y Electrónica

e) Implemente un algoritmo en Matlab que permita aproximar por el método de Newton-

Rapson, la solución de la ecuación f (x) = 0 con 6 cifras decimales exactas usando como ponto
inicial x0 = 1. Muestre los resultados de la ejecución del Algoritmo mediante una tabla de
datos detallada.

Código en Matlab
function newton raphson

clear
clc
i=1;
x=1;
e=1;
tol=1e-6;

20
Facultad de Ingeniería Eléctrica y Electrónica

fprint(’La funcion es : 2*(x∧ 3) + (x∧ 2) − 2\n’);

fprintf(’i \t \t x \t \t \t xn \t \t error \n’);

while e>tol

f=2 ∗ (x∧ 3) + (x∧ 2) − 2;

df=6 ∗ (x∧ 2) + (2 ∗ x);
xn=x-f/df;
e=abs((xn-x)/xn);
fprintf(’% d \t %f \t %f \t %f \n’ ,i,x,xn,e)
x=xn;
i=i+1;

end

21
Facultad de Ingeniería Eléctrica y Electrónica

f) Compare el método del punto Fijo y el método de Newton para el cálculo de la menor raíz de
la ecuación f (x) = 0. ¿Cuál de los dos métodos fue más rápido y a que se debe eso?. Describa
las ventajas y desventajas de un método con respecto del otro después de la experiencia.

RESPECTO A ESTA FUNCIONES, EL MÉTODO DE NEWTON RAPHSON ES MAS EFI-

CIENTE QUE EL DE PUNTO FIJO YA QUE SE NECESITAN MENOS ITERACIONES
PARA APROXIMARSE A LA RAÍZ, ESTO DEBIDO A LA NATURALEZA DE LA FUN-
CIÓN Y LA EXACTITUD DE SU VALOR INICIAL.

PERO HAY CASOS DONDE NEWTON RAPHSON SE COMPORTA DE MANERA MUY

DEFICIENTE DEBIDO A SU LENTA CONVERGENCIA, OTRA DESVENTAJA ES QUE
ES NECESARIO APOYARSE DE UNA GRÁFICA Y TENER UN VALOR INICIAL CER-
CANO A LA RAÍZ.

LA VENTAJA DE PUNTO FIJO ES QUE CUANDO CONVERGE, ES DE MUCHA PRECI-

SIÓN, NO NECESITA DE UN INTERVALO PARA FUNCIONAR SINO DE ÚNICAMENTE
UN PUNTO PERTENECIENTE AL INTERVALO DONDE ESTÉ LA RAÍZ. LA DESVEN-
TAJA DEL PUNTO FIJO ES QUE NO GARANTIZA LA CONVERGENCIA.

22