octave

>> % Metodo de Newton para f(x) = x^3 - 3x + 1

>> f = @(x) x^3 - 3*x + 1;

>> f_prime = @(x) 3*x^2 - 3;

>> % Valor inicial

>> x = 2;

>> % Numero de iteraciones deseadas

>> num_iter = 5;

>> % Encabezado

>> disp('Iteracin    x               f(x)');

Iteracin    x               f(x)

>> for i = 1:num_iter

  fx = f(x);

 fx_prime = f_prime(x);

% Imprimir los valores de cada iteracin

 fprintf('%5d     %12.8f     %12.8f\n', i, x, fx);

  % Aplicar la frmula de Newton

  x = x - fx / fx_prime;

end

    1       2.00000000       3.00000000

    2       1.66666667       0.62962963

    3       1.54861111       0.06804022

    4       1.53239016       0.00121814

    5       1.53208899       0.00000042

>>