Regula Falsi
Publicado por
Carmina (1 intervención) el 15/04/2015 16:46:35
]Hola. Aquí colgué el código para la Regula Falsi. Le faltan detalles, por ejemplo, no muestra el resultado de la aproximación y el error, pero eso se le puede agregar. Espero te sirva. Saludos.
% Método de la Regula Falsi.
a(1)=input('Extremo inferior del intervalo: ');
b(1)=input('Extremo superior del intervalo: ');
m=(funcion(b(1))*a(1)-funcion(a(1))*b(1))/(funcion(b(1))-funcion(a(1)));
e=abs(funcion(m));
if funcion(a(1))*funcion(m) < 0
a(2) = a(1);
b(2) = m;
else
a(2) = m;
b(2) = b(1);
end
n=2;
while e > 10^(-3) %(La tolerancia la podés ingresar al comienzo)
m=(funcion(b(n))*a(n)-funcion(a(n))*b(n))/(funcion(b(n))-funcion(a(n)));
e=abs(funcion(m));
if funcion(a(n))*funcion(m) < 0
a(n+1) = a(n);
b(n+1) = m;
else
a(n+1) = m;
b(n+1) = b(n);
end
n=n+1;
end