Общие сведения

Огромное количество задач вычислительной математики связано с решением нелинейных алгебраических уравнений, а также систем таких уравнений. При этом необходимость решения нелинейных уравнений возникает зачастую на промежуточных шагах, при реализации фрагментов более сложных алгоритмов (к примеру, при расчетах дифференциальных уравнений при помощи разностных схем и т. п.).

Численное решение нелинейного уравнения

Алгоритм приближенного решения уравнения f(x)=0 состоит из двух этапов:

1. нахождения промежутка, содержащего корень уравнения (или начальных приближений для корня);

2. получения приближенного решения с заданной точностью с помощью функции root.

Если после многих итераций MathCAD не находит подходящего приближения, то появится сообщение (отсутствует сходимость).

Эта ошибка может быть вызвана следующими причинами:

· уравнение не имеет корней;

· корни уравнения расположены далеко от начального приближения;

· выражение имеет комплексный корень, но начальное приближение было вещественным.

Чтобы установить причину ошибки, исследуйте график f(x). Он поможет выяснить наличие корней уравнения f(x) = 0 и, если они есть, то определить приблизительно их значения. Чем точнее выбрано начальное приближение корня, тем быстрее будет root сходиться.

Нахождение корней полинома

Для нахождения корней выражения, имеющего вид v₀+v₁x+… v_n-1x^n-1 +v_nxⁿ, лучше использовать функцию polyroots, нежели root. В отличие от функции root, функция polyroots не требует начального приближения и возвращает сразу все корни, как вещественные, так и комплексные.

ФункцияPolyroots(v) - возвращает корни полинома степени n. Коэффициенты полинома находятся в векторе v длины n + 1. Возвращает вектор длины n, состоящий из корней полинома.

Решение систем уравнений