Untitled
Anonymous
plain_text
02/15/2024 1:41 PM
684 B
22
Indexable
clc;
f = @(x) x^3 + 3*x + 2;
f_prime = @(x) 3*x^2 + 3;
x0 = 0;
tolerable_error = 0.0001;
% my input
nvalue = input('my number for iterations: ');
x = x0;
for iterations = 1:nvalue
x_new = x - f(x) / f_prime(x);
disp(x_new)
if abs(x_new - x) < tolerable_error
fprintf('root found at x = %f after %d iterations\n', x_new, iterations);
break;
end
x = x_new;
end
if iterations >= nvalue
disp('Newton-Raphson not convergen number.');
end
Editor is loading...
Leave a Comment