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