Как реализовать и модифицировать метод Ньютона в MATLAB / OCTAVE?

Question

Я пытаюсь реализовать метод Ньютона.Итеративный метод порядка 2. и адаптировать его для вычисления ненулевого решения уравнения (x = 1-e ^ (- 2x) с 4 точными десятичными знаками.

Я не могу понять, как должен выглядеть результатЯ пытался реализовать метод для нахождения лучших приближений к нулю с начальной апроксимацией

function [x , res , xvec , resvec] = newton( f, df, x0, maxiter, tol )

function x(k) = newton( f, df, x0, maxiter, tol )
k(1:maxiter)
f = @x(k)
df = @(x(k))'
x(k)= x0-(f(x0)/df(x0))
x(k+1) = x(k) - (f(x(k))/df(x(k)))
end
while (iter < maxiter)
  while (abs(f(n)-(f(n-1)) > tol)
    x(1)= x0-(f(x0)/df(x0))
  elseif (iter > maxiter) 
    error('maximum amount of iterations exceeded')
  end
end
fprintf('step(k)      x(k)       |f(xk)|     |x(k+1)-x(k)|')
fprintf('------    -----------  ---------    -------------\n')
fprintf('%d %d %.4f %.4f\n', k, x(k), mod(f(xk)), mod(x(k+1)−x(k)) )

Я ожидал показать мне результат с qpproximation до 4 децимел

.