Владимир Дьяконов - Maple 9.5/10 в математике, физике и образовании

> evala(AFactor(%));

> evala(AFactors(х^2-2*у^2));

Нетрудно заметить, что разложение полинома на множители позволяет оценить наличие у него корней. Однако для этого удобнее воспользоваться специальными функциями, рассмотренными ниже.

Для вычисления действительных и комплексных корней полиномов служит уже известная нам функции solve(p, x), возвращающая список корней полинома p одной переменной. Кроме того, имеются следующие функции для вычисления корней полиномов:

roots(р)

roots(р, K)

roots(р, х)

roots(р, x, K)

Эти функции вычисляют точные корни в рациональной или алгебраической области чисел. Корни возвращаются в виде [[r1,m1], [rn, mn]], где mi — это корень полинома, a mi — порядковый номер полинома. С действиями этих функций можно разобраться с помощью приведенных ниже примеров:

> р:=х^4 1-9*х^3+31*х^2+59*х+60;

> solve(р,х);

> roots(р,х);

> roots(х^2-4,х);

> expend((х-1)*(х-2)*(х-3)*(х-4));

> roots(%,х);

С полиномами могут выполняться различные операции. Прежде всего, отметим некоторые функции, которые относятся к одному полиному:

psqrt(p) — возвращает квадрат полинома;

proot(p,n) — возвращает n-ю степень полинома;

realroot(p) — возвращает интервал, в котором находятся действительные корни полинома;

randpoly(vars, eqns) — возвращает случайный полином по переменным vars (список) с максимальной степенью eqns;

discrim(p, var) — вычисление дискриминанта полинома по переменной var;

Primitive(a) mod p — проверка полинома на примитивность (возвращает true, если полином примитивен).

Действие этих функций достаточно очевидно, поэтому ограничимся приведением примеров их использования (файл polop):

> psqrt(х^2+2*х*у+у^2);

> proot(х^3+3*х^2+3*х+1, 3);

> psqrt(x+y);

> proot(x+y, 2);

> р:=х^3-3*х^2+5*х-10;

> discrim(p,x);

> readlib(realroot):

> realroot(p);

> randpoly([x],degree=10);

> randpoly([x],degree=10);

> randpoly([x],degree=10);

> Primitive(х^4+х+1) mod 2;

Обратите внимание на то, что для использования некоторых из приведенных функций необходим вызов их из стандартной библиотеки. Для функции randpoly приведенные результаты случайны, так что, скорее всего, их повторение невозможно.

С полиномами можно выполнять обычные операции, используя для этого соответствующие операторы:

> readlib(psqrt):

> readlib(proot):

> Primitive(х^4+х+1) mod 2;

> p1:=a1*x^3+b1*x^2+c1*x+d1: p2:=а2*х^2+b2*х+с2:

> p1+p2;

> p1*p2;

> collect(%,х);

> p1/p2;

> expand(%,х);

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

> diff(p1, х);

> diff(p1, x$2);

> Int(p1,x)=int(p1,x);

> Int(p1,х=0..1)=int(p1,х=0..1);

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

> pp:=а*х^(-2)+b*x^(-1)+c*x+d+e*x^2+f*х^3;

Нетрудно показать, что с ним можно выполнять различные операции:

> рр+рр;

> рр-рр;

> pp^2

> simplify(%);

> Diff(pp, x)=diff(pp, x)

> Int(pp,x);

> int(рр,х);

Хотя Maple и не накладывает ограничений на применение степенных многочленов (полиномов) с отрицательными степенями свойства таких полиномов заметно отличаются от свойств полиномов с положительными степенями, поэтому при применении первых надо проявлять известную осторожность.

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