TAREA 5 METODOS MATEMATICOS CON PROGRAMACION
Dr. Juan Carlos García Castrejón
Sahiril Fernanda Rodriguez fuentes
Método de Newton Raphson
Ejercicio 1
Funciones: 𝑓(𝑥) = 𝑒 −𝑥 − 𝑥 𝑓 ′ (𝑥) = −𝑒 −𝑥 − 1
Iteraciones, porcentaje de error, raíz.
COMANDOS
clear clf;
clear all;
xg=[-6:1:6]
fxg=(exp (-xg))-xg
fxg2=-exp(-xg)-1
hold on
plot(xg,fxg)
plot (xg,fxg2)
grid on
x0=0
x(1)=x0
tolerancia=0.000001;
error=100;
for i=1:1:100
y(i)=((exp(-x(i))-x(i)));
i
x(i+1)=x(i)-((exp(-x(i)))-x(i))/(-exp(-x(i))-1);
error=abs(x(i)-x(i+1))*100/x(i)
if (tolerancia>error);
break
end
end
y(i+1)=((exp(-x(i+1))-x(i+1)));
plot(x,y,'*');
Ejercicio 2
Funciones: 𝑓(𝑥) = 2𝑥 3 − 11.7𝑥 2 + 17.7𝑥 − 15
𝑓 ′ (𝑥) = 6𝑥 2 − 23.4𝑥 + 17.7
Iteraciones, porcentaje de error, raíz.
COMANDOS
clear clf;
clear all;
xg=[-6:1:6]
fxg=(2*xg.^3)-(11.7*xg.^2)+(17.7*xg)+5
fxg2=(6*xg.^2)-(23.4*xg)+(17.7)
hold on
plot(xg,fxg)
plot (xg,fxg2)
grid on
x0=3
x(1)=x0
tolerancia=0.001;
error=100;
for i=1:1:inf
y(i)=((2*x(i).^3)-(11.7*x(i).^2)+(17.7*x(i))+5);
i
x(i+1)=x(i)-((2*x(i).^3)-(11.7*x(i).^2)+(17.7*x(i))+5)/((6*x(i).^2)-
(23.4*x(i))+(17.7));
error=(abs((x(i+1)-x(i))/x(i+1)))*100
if tolerancia>error
break
end
end
y(i+1)=((2*x(i+1).^3)-(11.7*x(i+1).^2)+(17.7*x(i+1))+5);
plot(x,y,'*');
Ejercicio 3
Funciones: 𝑓(𝑥) = −4.3𝑥 3 + 18𝑥 2 − 21𝑥 − 12
𝑓 ′ (𝑥) = −12𝑥 2 + 36𝑥 − 21
Iteraciones, porcentaje de error, raíz.
COMANDOS
clear clf;
clear all;
xg=[-30:1:20]
fxg=-(4.3*xg.^3)+(18*xg.^2)-(21*xg)-12
fxg2=-(12.9*xg.^2)+(36*xg)-21
hold on
plot(xg,fxg)
plot (xg,fxg2)
grid on
x0=3
x(1)=x0
tolerancia=0.05;
error=100;
for i=1:1:100
y(i)=-(4.3*x(i).^3)+(18*x(i).^2)-(21*x(i))-12;
i
x(i+1)=x(i)-(-(4.3*x(i).^3)+(18*x(i).^2)-(21*x(i))-12)/(-
(12.9*x(i).^2)+(36*x(i))-21);
error=(abs((x(i+1)-x(i))/x(i+1)))*100
if tolerancia>error
break
end
end
y(i+1)=(-(4.3*x(i).^3)+(18*x(i).^2)-(21*x(i))-12);
plot(x,y,'*');
Ejercicio 4
Funciones: 𝑓(𝑥) = sin 𝑥 + cos(1 + 𝑥 2 ) − 1
𝑓 ′ (𝑥) = −2𝑥𝑠𝑖𝑛(𝑥 2 + 1) + cos(𝑥)
Iteraciones, porcentaje de error, raíz.
COMANDOS
clear clf;
clear all;
xg=[-12.85:.01:12.85]
fxg=sin(xg)+ cos(1+(xg.^2))-1
fxg2=cos(xg)-(2*xg).*sin((xg.^2)+1)
hold on
plot(xg,fxg)
plot (xg,fxg2)
grid on
x0=0
x(1)=x0
tolerancia=0.1;
error=100;
for i=1:1:100
y(i)=sin(x(i))+ cos(1+x(i).^2)-1;
i
x(i+1)=x(i)-((sin(x(i))+ cos(1+x(i).^2)-1)/(cos(x(i))-
(2.*x(i)).*sin((x(i).^2)+1)));
error=(abs((x(i+1)-x(i))/x(i+1)))*100
if tolerancia>error
break
end
end
y(i+1)=sin(x(i))+ cos(1+x(i).^2)-1;
plot(x,y,'*');
TABLA DE COMPRACIÓN DE RESULTADOS
EJERCICIO PROGRAMA RAIZ % ERROR ITERACIONES Xo TOLERANCIA
En clase 0.5671 5.08e-10 5 0 0.000001
1
De fernanda 0.567 5.08e-10 5 0 0.000001
En clase -0.2421 1.12e-8 7 3 0.001
2
De fernanda -0.242 1.12e-8 7 3 0.001
En clase -0.4118 2.94e-5 13 3 0.05
3
De fernanda -0.4118 2.94e-5 13 3 0.05
En clase 4.27 5.59 100 0 0.1
4
De fernanda ---------- ----------- ---------- 0 0.1