Implementação em scilab falsa posição

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

function y=f(x)
 y = x^2 + log(x);
endfunction
function [raiz,it]=falsaposicao(f,a,b,e,N)
 it=0; // começa com interações 0
 erro=1; // pode ser inicializado em 1 pois o e é na base de 10^-...
 x=a; //pode ser a ou b 
 inter = 0;
 while erro>e & it<N
 xa= x // Xanterior(guardar o ultimo x que foi calculado) recebe o x
 x=a - f(a)*(b-a)/(f(b)-f(a));
 if f(x)*f(a)<0 
 b=x;
 else
 a=x;
 end
 it=it+1; // implementando as interações 
 erro=abs((x-xa)/x);
 inter= it;
 end
 it = inter;
 raiz = x;
endfunction