Пример вычислений по методу хорд в пакете MATHCAD.Функция hord уточняет корень уравнения Входные параметры: значения концов отрезка изоляции корня, точность решения. Выходные параметры функции: два последних приближения для корня, количество итераций для выполнения условия точности. В программном блоке используются: 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 Варианты лабораторных работ
|
на интервале [a,b] с точностью до 
.
