Códigos Parcial
Códigos Parcial
Códigos Parcial
Có digo Regula
function [r,rep]=reg(a,b,n,tol)
for i= 1:n
m=(funcion(b)-funcion(a))/(b-a);
corte=funcion(a)-(m*a);
c=(-corte)/m;
if funcion(a)*funcion(c)<0
b=c;
else
a=c;
end
if abs(funcion(c))<tol
r = c;
rep=i;
break
end
end
end
function y=funcion(x)
y=sin(x);
end
Có digo Secante
function [r,rep]=secante(a,b,n,tol)
format long
for i = 1:n
m=(funcion(b)-funcion(a))/(b-a);
corte=funcion(a)-(m*a);
c=(-corte)/m;
if abs(funcion(c))<tol
r=c;
rep=i;
break
end
a=b;
b=c;
end
end
function y=funcion(x)
y=cos(x);
end
Có digo Newton
function [r,rep]=new(a,n,tol)
for i = 1:n
m=derivada(a);
corte=funcion(a)-(m*a);
c = (-corte)/m;
if abs(funcion(c))<tol
r=c;
rep=i;
break
end
a=c;
end
end
function y=funcion(x)
y=cos(x);
end
function z=derivada(x1)
z=-sin(x1);
end
Código Punto Fijo
function [r,rep]=fijo(a,n,tol)
for i = 1:n
a=funcion(a);
if abs((a-funcion(a)))<tol
r=a;
rep=i;
break
end
end
end
function y=funcion(x)
y=x^(2)-4;
end