f = @(x) exp(-x) - x;
df = @(x) -exp(-x) - 1;
 
% Initial guess
a = 0;
 
% Number of iterations
n = 5;
 
% Newton-Raphson Iteration
for i = 1:n
    r = a - (f(a) / df(a));
    tol(i) = abs(a - r);
 
    % Display root at each step
    disp(['Iteration ', num2str(i), ': Root estimate = ', num2str(r)]);
 
    a = r;
end
 
disp(['Final root estimate after ', num2str(n), ' iterations:']);
disp(r);
 
disp('Error (tolerance) at each iteration:');
disp(tol);
# your code goes here

ZiA9IEAoeCkgZXhwKC14KSAtIHg7CmRmID0gQCh4KSAtZXhwKC14KSAtIDE7CgolIEluaXRpYWwgZ3Vlc3MKYSA9IDA7CgolIE51bWJlciBvZiBpdGVyYXRpb25zCm4gPSA1OwoKJSBOZXd0b24tUmFwaHNvbiBJdGVyYXRpb24KZm9yIGkgPSAxOm4KICAgIHIgPSBhIC0gKGYoYSkgLyBkZihhKSk7CiAgICB0b2woaSkgPSBhYnMoYSAtIHIpOwogICAgCiAgICAlIERpc3BsYXkgcm9vdCBhdCBlYWNoIHN0ZXAKICAgIGRpc3AoWydJdGVyYXRpb24gJywgbnVtMnN0cihpKSwgJzogUm9vdCBlc3RpbWF0ZSA9ICcsIG51bTJzdHIocildKTsKICAgIAogICAgYSA9IHI7CmVuZAoKZGlzcChbJ0ZpbmFsIHJvb3QgZXN0aW1hdGUgYWZ0ZXIgJywgbnVtMnN0cihuKSwgJyBpdGVyYXRpb25zOiddKTsKZGlzcChyKTsKCmRpc3AoJ0Vycm9yICh0b2xlcmFuY2UpIGF0IGVhY2ggaXRlcmF0aW9uOicpOwpkaXNwKHRvbCk7CiMgeW91ciBjb2RlIGdvZXMgaGVyZQ==