Método de Newton Raphson em Scilab

Versão do Scilab: 5.5.2


clear
clc
//Método Newton-Raphson
function [resultado]=equacao(x)
    resultado=(sin(x)-x^2+4)   //colocar a equação de entrada aqui
endfunction

x0= 2.5;   //x inicial 
numeroIteracoes=0;

// Erro tolerável:
e=10^-5    //erro tolerável
d=2*e;     //diferença - condição de parada

df = numderivative(equacao, 2)
pol= df
n=size(pol,2)
y=[]

while d > e
    j=1
    for i=n:-1:1
        y(j)=[x0^(i-1)]
        j=j+1
    end

    derivada=pol*y;
    x=x0-equacao(x0)/derivada
    d=abs(x-x0)
    x0=x;
    numeroIteracoes=numeroIteracoes+1;
end

disp("Raiz da Função : ")
disp(x)
disp("Número de Iterações : ")
disp(numeroIteracoes)