Scilab - Método de Newton

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

//x1 é o "chute" inicial
//Ep = precisão (0<Ep<1)
//it - número máximo de iterações
function x=newton(x1,Ep,it)
 exec func2.sci
 h=10^(-8)
 x(1)=x1
 dx=(func2(x1+h)-func2(x1))/h
 x(2)=x1-(func2(x1)/dx)
 n=2
 while abs(x(n)-x(n-1))>Ep &n<(it+2)
 n=n+1
 x(n)=x(n-1)-(func2(x(n-1))/dx)
 dx=(func2(x(n)+h)-func2(x(n)))/h
 end
endfunction