Codigo Biseccion Octave
Codigo Biseccion Octave
Codigo Biseccion Octave
Este programa toma como entrada una función f, un intervalo [a,b] en el que buscar la raíz, y una
tolerancia tolerance. Devuelve una estimación de la raíz y el número de iteraciones necesarias
para encontrar la raíz.
Para utilizar el programa, simplemente define una función f en Octave y llámala junto con los
parámetros a, b y tolerance. Por ejemplo, si quisieras encontrar una raíz de la función f(x) = x^2 - 2
en el intervalo [0,2] con una tolerancia de 0.001, podrías hacer lo siguiente:
scssCopy code
f = @(x) x^2 - 2; [a, b] = [0, 2]; tolerance = 0.001; [root, num_iterations] = bisection(f, a, b,
tolerance); printf("La raíz es %f\n", root); printf("Se necesitaron %d iteraciones.\n",
num_iterations);
Copy code
% Initialize variables
fa = f(a);
fb = f(b);
num_iterations = 0;
if fc == 0
root = c;
return
elseif sign(fc) == sign(fa)
a = c;
fa = fc;
else
b = c;
fb = fc;
end
num_iterations = num_iterations + 1;
end
f = @(x) x^2 - 2;
[a, b] = [0, 2];
tolerance = 0.001;
La raíz es 1.414062
Se necesitaron 10 iteraciones.