Scilab - Método da Falsa posiçã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]=falsaposic(a,b,Ep,it)
 exec func.sci;
 n=2;
 x(1)=a;
 x(2)=b;
 while abs(x(n)-x(n-1)) > Ep &n<it
 n=n+1;
 x(n)=(a*func(b)-b*func(a))/(func(b)-func(a))
 if func(a)/abs(func(a))==func(x(n))/abs(func(x(n))) then
 a=x(n);
 else
 b=x(n);
 end
 end
endfunction