Пример реализации модифицированного метода

Касательных в пакете MATLAB для решения алгебраического уравнения.

function res=p_modif_kasat(f,A,B,eps);

% polyval(f,A) – вычисление значения полинома f в точке A

% polyder(f) – вычисление первой производной для полинома

% выбор точки проведения касательной

if polyval(f,A)*polyval(polyder(polyder(f)),A)>0

x=A

else

x=B;

end;

if x==A

prev_x=B

else

prev_x=A;

end;

fix=x;

% цикл уточнения корня

while abs(prev_x-x)>eps

prev_x=x;

x=x-polyval(f,x)/polyval(polyder(f),fix);

end

res=x;

Вызов функции p_modif_kasat в командном режиме. Полином необходимо задать набором соответствующих коэффициентов

>> res=p_modif_kasat([1,5,-3],0,2,0.001)

res =

0.5416

Варианты лабораторных работ