Shevchenko / 2sem / LR9_2

clc
clear
f = @(x) 0.1*x.^3 + 20*cos(x) + 10;
h = 0.1; %Шаг производной
df = @(x) (f(x+h) - f(x-h))/(2*h); %Вычисление производной
tol = 1e-6; %Точность
iterlimit = 200; %Предел итераций
ezplot(f); %Создание графика
grid on
x0 = [];

function[root, iter] = newton(f, df, x0, tol, iterlimit);
iter = 0;
while abs(f(x0)) > tol && iter < iterlimit
x1 = x0 - f(x0)/df(x0);
x0 = x1;
iter = iter + 1;
end
if abs(f(x0)) <= tol
root = x0;
else
root = [];
end
end

for x = -6.5:h:6.5
if f(x)*f(x+h) < 0
x0 = [x0, x];
end
end


tic
for i=1:length(x0);
[root, iter] = newton(f, df, x0(i), tol, iterlimit);
if ~isempty(root);
disp(['Корень при х0 = ', num2str(x0(i))]);
disp(['х = ', num2str(root)]);
disp(['f(x) = ', num2str(f(root))]);
disp(['Количество итераций = ', num2str(iter)]);
disp(' ');
else
disp('Превышено число итераций');
end
end
toc