Método da Bisseção em Scilab

Versão do Scilab: 5.5.2

clear
clc
//METODO DA BISSEÇÃO
function [resultado] = equacao(x)
    resultado = x^2+log(x)             //aqui vai a equacao de entrada
endfunction

//INTERVALO DAS RAIZES
a = 0.5;    //intervalo inicial
b = 1 ;        //intervalo final

//ERRO TOLERAVEL
e=10^(-4);            //erro
d=2*e;                   //diferença - condicao de parada do laço

x0 = 0;          // x inicial   
contIteracao = 0;

while d > e
    fa = equacao(a) 
    fb = equacao(b)  
    x = (a*fb - b*fa)/(fb - fa);
    fx = equacao(x)
    if fa*fx
        b = x
    else 
        a = x
    end
    d = abs(x-x0);
     x0 = x;            
     contIteracao = contIteracao +1;
end

disp("Raiz aproximada da funcao: ")
disp(x0)
disp("Numero de iterações: ")
disp(contIteracao)