Scilab - Método da Bissecção

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

//'a' e 'b' são os extremos do intervalo onde a raiz está localizada
//Ep = precisão (0<Ep<1)
//it - número máximo de iterações
function x=bissec(a,b,Ep,it)
 exec func.sci;
 n=0;
 while abs(b-a) > Ep &n<it
 n=n+1;
 x(n)=(a+b)/2;
 if func(a)/abs(func(a))==func(x(n))/abs(func(x(n))) then
 a=x(n);
 else
 b=x(n);
 end
 end
endfunction