Добавил:
deadpigeon
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:БДЗ2 / newton
.m clc,clear
syms x
fun = x.^5 - 2*sin(x);
a = 0;
b = 1;
e = 0.001;
der = diff(fun);
der2 = diff(der);
xmin = (a+b)/2;
fmin = subs(fun,xmin);
A = [];
B = [];
while abs(subs(der,xmin))>=e
xmin = xmin - subs(der,xmin)/subs(der2,xmin);
A(end+1) = xmin;
B(end+1) = subs(fun,xmin);
end
A((end-2):end)
B((end-2):end)
xmin
fmin = subs(fun,xmin)
syms x
fun = x.^5 - 2*sin(x);
a = 0;
b = 1;
e = 0.001;
der = diff(fun);
der2 = diff(der);
xmin = (a+b)/2;
fmin = subs(fun,xmin);
A = [];
B = [];
while abs(subs(der,xmin))>=e
xmin = xmin - subs(der,xmin)/subs(der2,xmin);
A(end+1) = xmin;
B(end+1) = subs(fun,xmin);
end
A((end-2):end)
B((end-2):end)
xmin
fmin = subs(fun,xmin)