Porque de poquito a poquito se llega .... Metodo de Biseccion

Gabo, que es la que? Aqui esta el codigo de matlab que implementa el metodo de biseccion para hallar raices de una funcion continua bla bla bla tu sabes ..... Ojala te sirva y recuerda que debes definir f(x) para que funcione. Cuidate brother .....

function [root, iter] = bisection(x0, x1,maxiter)

iter=0;
xnext = (x1+x0)/2;

if f(x0)*f(x1)<0
while itereps
if f(x0)*f(xnext)<0
x1=xnext;
else
x0=xnext;
end
iter = iter+1;
xnext = (x1+x0)/2;
end
iter
else
'Error: el metodo no funciona para este caso.'
end

root = x0;

return
end

Bonito verda' ... bueno ahi los dejo ...

Happy Valentine's!!

Aunque se que este dia no es mas que otro mas, igual a los demas, y que realmente no hay nada que celebrar, me es inevitable la idea de aprovechar, ya que todo el mundo esta en el "mood", para mandarles un saludo a todos los amigos que pueda, y por si alguno se siente mal porque no lo llame o algo pues que sepa que al menos me paso por la mente un instante. A Abe, Santhia, Eric, Pringle, Lori, Lorena, Enily, Danny, Ricky, Melissa Scarcliff, Miraya, a LOA xoxoxo, a Aldo Perez ca..b..on, siempre te voy a querer como un hermano, a Nandy, Pedrito, Damian, los quiero gente de corazon. A todo el corillo de Colombia: Angie, Edwin, KTA !!! , Gabo, el Diablo, Sindy, Lili, Julian, Leonid, El chino, Rafa, la Mona, William, Ana, ademas a Jehu y a Trilce, ...tambien les mando un abrazo grandote. Saludos tambien a to los panitas de Mayaguez a Coco, Gabriel, Francely, JeanLuis, Karen, Ahmed y to' el corillo de Western. A mis housemates: Juan, Christian, Piwi, Tony y a mi ex-housemates for life Kevo y Carnation los quiero chorro de ca....b...nes. Un besote tambien para los vieja escuela de la UPR, a Beca .. muakata un besote para usted. Saludos a Jose Roberto, a Villa, Nilmar, al amore de la mia vita Iara kookie my love, a Vane cool, a Mario Caldari ..... pa mi tu siempre vas a ser familia, a una super amiga de esas que ya no se encuentran muchas, Maristela for you here's some love !! A to el mundo de la AECC: Osvaldo, Emilio, Mon, Poly, Pablo CUCA y todos los que no veo hace tiempo a Suleika, JuanJO(de Marbella), Edgar Rivera, Cuco, Chucho, Danny Chui, Angel Vecino, Ivy, Taira, Marina y toda la gente buena de Santa Rita, paz!!! Tambien un saludo super especial a una persona que llevo siempre muy cerquita del corazon, Txio besotes!! Saludos tambien y un abrazote gigantesco a mi madrina Laurita ... te quiero vieja... nunca te olvides. A mis viejos los amo, los amo, los amo!! e igual a mis hermanso Alanis, Joselyn, Nelson, Luigi y Rafael, felicidades a todos!!!

Gabo, tieni qui, il "Metodo de Newton", come ti ho promesso.

function [flag, root, iter] = Newton(x0, maxiter, tol)

iter=0;
xn = x0;
xnext = xn-f(xn)/fp(xn);
iter = iter+1;

NewtonErrors(iter)=abs(xnext-xn);

while((iter < maxiter & abs(xnext-xn)) >= tol)
xn = xnext;
xnext = xn-f(xn)/fp(xn);
iter = iter+1;
NewtonErrors(iter) = abs(xnext-xn);
end

if (iter) < (maxiter)
root=xnext;
flag=0;
else
flag=1;
end

iteraciones = 1:iter;
plot(iteraciones, NewtonErrors,'r+-');
grid on;
xlabel('Iteraciones');
ylabel('Error del metodo de Newton');
title('Iteraciones vs. Error obtenido');
return
end

function [f_x]= f(x)
f_x = 2*x*cos(2*x)-(x-2)^2;
end

function [fp_x] = fp(x)
fp_x = 2*(cos(2*x) - 2*x*sin(2*x)-x+2);
end


Recuerda que para usar tienes que cambiar f(x) en el codigo y lo mismo con la derivada de f(x). Este codigo tambien grafica el comportamiento del error para que asi puedas ver cuan rapido converge. Ah y una cosa importante si al correr el programita flag=1 el resultado almacenado en root es irrelevante.

Bueno, ojala te sirva de algo.... Cuidate boy e parliamo dopo.... ciao