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

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

Int(expr,x) — инертная форма вычисления неопределенного интеграла;

Doubleint(expr,x,у,Domain) — вычисление двойного интеграла по переменным х и у по области Domain;

Tripleint(expr,x,y,z) — вычисление тройного интеграла;

intparts(f,u) — интегрирование по частям.

Ниже дан пример применения функции Tripleint пакета student:

> Tripleint(f(х,у,z),х,у,z);

> Tripleint(х*у*z^2,x=0..2,y=0..3,z=0..5);

> evalf(%);

> int(int(int(x*y*z^2,x=0..2),y=0..3),z=0..5);

Пакет student имеет три графические функции для иллюстрации интегрирования методом прямоугольников:

leftbox(f(x), x=a..b, о) или leftbox(f(x), x=a..b, n, 'shading'=, о);

rightbox(f(x), x=a..b, о) или rightbox(f(x), x=a..b, n, o);

middlebox(f(x), x=a..b, о) или middlebox(f(x), x=a..b, n, o);

Здесь f(x) — функция переменной x, x — переменная интегрирования, a — левая граница области интегрирования, b — правая граница области интегрирования, n — число показанных прямоугольников, color — цвет прямоугольников, о — параметры (см. ?plot,options).

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

showtangent(f(х), х=а)

Рисунок 4.36 показывает все эти возможности пакета student. Четыре отмеченных вида графиков здесь построены в отдельных окнах.

Рис. 4.36. Примеры иллюстративной графики пакета student

Графические средства пакета student ограничены. Но они предоставляют как раз те возможности, которые отсутствуют в основных средствах построения графиков. В Maple 9/9.5 функции пакета резко расширены и мы вернемся к их рассмотрению в Главе 9.

Пакет Student обеспечивает визуализацию ряда методов численного интегрирования: методов прямоугольников с различным расположением их, метода трапеций и метода парабол (Симпсона). Это возможно в символьном виде, например (файл intvis):

> with(student): middlesum(x*exp(-x), x=a..b);

Ниже представлено несколько примеров такой визуализации (для метода прямоугольников со средним расположением их, метода трапеций и метода Симпсона):

> with(student): middlesum(х*exp(-x), x=0..4);

> trapezoid(х*ехр(-х), x=0..4);

> simpson(x*sin(-х), х=1..4);

> evalf(%)

В последнем примере показано вычисление по представлению методом Симпсона.

Для работы с алгебраическими кривыми служит пакет расширения algcurves. Он загружается командами:

> restart;with(algcurves);

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

Weierstrassform(f,x,y,x0,y0,opt) — вычисление нормальной формы для эллиптических или гиперболических алгебраических кривых;

differentials(f, x, у, opt) — голоморфные дифференциалы алгебраических кривых;

genus(f,x,y,opt) — подлинность алгебраической кривой;

homogeneous(f,x,y,z) — создание полинома двух переменных, гомогенного в трех переменных;

homology(f, х, у) — находит канонический гомологический базис по алгоритму Треткоффа;

integral_basis(f, х, у, S) — интегральный базис алгебраического поля функции;

is_hyperelliptic(f, х, у) — тестирует кривую на ее принадлежность к гиперболической;

j_nvariant(f,x,y) — j-инвариант алгебраической кривой;

monodromy(f, х, у, opt) — вычисляет монодромию алгебраической кривой;

parametrization(f,x,y,t) — находит параметризацию для кривой с родом (даваемым функцией genuc), равным 0;