Algoritmo de Newton Raphson de Cálculo Numérico en Scilab.

newton_raphson.sce

function [raizaproximada, iteraciones] = newton(funcion, aproximacioninicial, Tolerancia, Iteracionesmaximas)
    aprox= aproximacioninicial;
    while 1
        raizaproximada = aprox - funcion(aprox)/numderivative(funcion, aprox);
        errorabsoluto = abs((raizaproximada-aprox)/raizaproximada)*100;
        aprox = raizaproximada;
        iteraciones = iteraciones + 1;
        if(errorabsoluto < Tolerancia | iteraciones  == Iteracionesmaximas ) then
        break
        end
    end
endfunction

//Datos

function y = primera(x)
    y = %e^x + 2^-x + 2*cos(x) - 6;
endfunction


A = 1;
B = 2;
Tolerancia = 0.00001;
Iteracionesmaximas = 100;
[raizaproximada, iteraciones] = newton(primera, A-B, Tolerancia, Iteracionesmaximas)