Taller Metodos Numericos
Taller Metodos Numericos
Taller Metodos Numericos
1.
% Metodo de Regula Falsi
i=1;
a(i)=a;
b(i)=b;
if (funcion(a)*funcion(b)>0)
error('la funcion no tiene raices en el intervalo')
end
fa(i)=funcion(a(i));
fb(i)=funcion(b(i));
while min(abs(fb(i)),abs(fa(i)))>tol
if (fa(i)*fb(i)==0)
if(fa(i)==0)
c(i)=a(i);
break
else
c(i)=b(i);
break
end
else
c(i)=a(i)-(fa(i)*(b(i)-a(i)))/(fb(i)-fa(i));
fc(i)=funcion(c(i));
if (fa(i)*fc(i)>0)
a(i+1)=c(i);
b(i+1)=b(i);
fa(i+1)=fc(i);
fb(i+1)=fb(i);
i=i+1;
else
a(i+1)=a(i);
b(i+1)=c(i);
fa(i+1)=fa(i);
fb(i+1)=fc(i);
i=i+1;
end
end
end
end
end
disp("xf: "+x0);
disp("porcentaje error: "+calcE)
4.
El método secante: también es un método recursivo para encontrar la raíz de los polinomios
mediante aproximaciones sucesivas.
Es similar al método Regular-falsi pero aquí no necesitamos verificar f (x 1 ) f (x 2 ) <0 una y
otra vez después de cada aproximación.
En este método,
las raíces de los vecindarios se aproximan mediante una línea secante o una cuerda a la
función f (x) .
También es una ventaja de este método que no necesitamos diferenciar la función dada f
(x) ,
como hacemos en el método de Newton-raphson .
% Metodo de la Secante
% Este metodo permite resolver la ecuacion f(x)=0, para una funcion
% continua en un intervalo [xa, xb].
clear
clc
format long
% Construccion de la funcion f(x):
f=inline('x.^2-2','x'); %#ok<DINLN>
if f(xa)*f(xb)>0
error('La función no tiene raices en el intervalo')
end
return
else
xa=xb;
ya=yb;
xb=xn;
yb=yn;
end
end
fprintf('Se alcanzo el numero maximo de iteraciones\n')
5.
clear
clc
r=3;
A=(pi*(2*r)^2)/4;
h=30/A;
i=0;
while(1>0)
i=i+1;
disp("Iteracion "+i)
fh=(pi*h^3)-9*(pi*h^2)+90;
dfh=(3*pi*h^2)-(18*pi*h);
hf=h-(fh/dfh);
err=abs(((hf-h)/hf))*100;
disp("error: "+err);
disp("h supuesta: "+h)
disp("h nueva : "+hf)
disp(" ");
if(err>0.0001)
h=hf;
else
break;
end
end
disp("error: "+err);
disp("h finnal : "+hf)