Пример вычислений по методу хорд в пакете MATHCAD.

Функция hord уточняет корень уравнения на интервале [a,b] с точностью до .

Входные параметры: значения концов отрезка изоляции корня, точность решения.

Выходные параметры функции: два последних приближения для корня, количество итераций для выполнения условия точности.

В программном блоке используются: f(x) – функция исходного уравнения ; f2(x) – функция второй производной для f(x). Обе функции должны быть заданы до программного блока как функции пользователя.

Неподвижная точка метода хорд определяется по следующему правилу: неподвижен тот конец, для которого знак функции f(x) совпадает со знаком ее второй производной f”(x) и ее значение заносится в fix.

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

Касательных в пакете 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

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

Номер варианта	Уравнение	Номер варианта	Уравнение