Untitled
unknown
plain_text
2 years ago
511 B
10
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