bissecao

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

a=0;
b=1;
p=0.5;
for i=0:4
 p = (a+b)/2;
 disp([a,b,p,f(a),f(b),f(p)]);
 if (f(a) * f(p) < 0 ) then
 b=p;
 else
 a=p;
 end
 if (f(p) == 0) then
 break;
 end
 i=i+1;
 
end
disp(p)