Вычисление корней кубического полинома

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

Векторные и матричные операторы и функции системы MatCAD позволяют решать широкий круг задач линейной алгебры. К примеру, если задана матрица Аи вектор В для системы линейных уравнений в матричной форме А*Х:=В, то вектор решения можно получить из очевидного выражения

Х:=А-1.

Поскольку решение систем линейных уравнений довольно распространенная задача, то для этого в MatCAD, начиная с шестой версии введена встроенная функция lsolve(A,B), которая возвращает вектор Хдля системы линейных уравнений А*Х=Впри заданной матрице коэффициентов А и векторе свободных членов В. Если уравнений n, размерность вектора В должна быть n, а размер матрицыАn* n.

Приведем пример решения системы линейных уравнений (обратите внимание, что решаемая система имеет матрицу с комплексными коэффициентами):

 

 

  А:=   0.24 -0.08
0.09 -0.15
0.04 -0.08

 

Матрица коэффициентов системы

 

  В:=  

 

Вектор свободных членов

 

Х:=А-1ВРешение системы

 

 

  Х=   1.909
3.195
5.0450

 

 

Результаты решения

 

 

  Х1=   1.909
3.195
5.0450

 

Решение с применением функции lsolve:

 

Х1:=lsolve(A,B)

 

 

Поиск корня нелинейного уравнения с помощью функции

Root

 

Многие уравнения, например трансцендентные, и системы из них не имеют аналитических решений. Однако они могут решаться численными методами с заданной погрешностью(не более значения, заданного системной переменной TOL). Для простейших уравнений вида F(x)=0 решение находится с помощью функции

root (выражение, Имя_переменной)

Эта функция возвращает значение переменной, при котором выражение дает 0. Функция реализует вычисления итерационным методом, причем можно задать начальное значение переменной. Это особенно полезно, если возможно несколько решений. Тогда выбор решения определяется выбором начального значения переменной. Приведем пример решения кубического полинома.

Вычисление корней кубического полинома

 

Коэффициенты полинома

 

F(x):= Задание полинома