Аналитическое и численное интегрирование.

Неопределенный интеграл вычисляется с помощью 2-х команд:

1) прямого исполнения – int(f, x), где f – подынтегральная функция, x – переменная интегрирования;

2) отложенного исполнения – Int(f, x) – где параметры команды такие же, как и в команде прямого исполнения int. Команда Int выдает на экран интеграл в аналитическом виде математической формулы.

Для вычисления определенного интеграла в командах int и Int добавляются пределы интегрирования, например,

> Int((1+cos(x))^2, x=0..Pi)=

int((1+cos(x))^2, x=0..Pi);

Если в команде интегрирования добавить опцию continuous: int(f, x, continuous), то Maple будет игнорировать любые возможные разрывы подынтегральной функции в диапазоне интегрирования. Это позволяет вычислять несобственные интегралы от неограниченных функций. Несобственные интегралы с бесконечными пределами интегрирования вычисляются, если в параметрах команды int указывать, например, x=0..+infinity.

Численное интегрирование выполняется командой evalf(int(f, x=x1..x2), e), где e – точность вычислений (число знаков после запятой).

Интегралы, зависящие от параметра. Ограничения для параметров.

Если требуется вычислить интеграл, зависящий от параметра, то его значение может зависеть от знака этого параметра или каких-либо других ограничений. Рассмотрим в качестве примера интеграл , который, как известно из математического анализа, сходится при а >0 и расходится при а <0. Если вычислить его сразу, то получится:

> Int(exp(-a*x),x=0..+infinity)=

int(exp(-a*x),x=0..+infinity);

Definite integration: Can't determine if the integral is convergent.

Need to know the sign of --> a

Will now try indefinite integration and then take limits.

.

Таким способом интеграл с параметром не вычислить. Для получения явного аналитического результата вычислений следует сделать какие-либо предположения о значении параметров, то есть наложить на них ограничения. Это можно сделать при помощи команды assume(expr1), где expr1 – неравенство. Дополнительные ограничения вводятся с помощью команды additionally(expr2), где expr2 – другое неравенство, ограничивающее значение параметра с другой стороны.

После наложения ограничений на параметр Maple добавляет к его имени символ (~), например параметр a, на который были наложены некоторые ограничения, в сроке вывода будет иметь вид: a ~.

Описание наложенных ограничений параметра a можно вызвать командой about(a). Пример: наложить ограничения на параметр a такие, что a > - 1, a £3:

> assume(a>-1); additionally(a<=3);

> about(a);

Originally a, renamed a~:

is assumed to be: RealRange(Open(-1),3)

Вернемся к вычислению интеграла с параметром , которое следует производить в таком порядке:

> assume(a>0);

> Int(exp(-a*x),x=0..+infinity)=

int(exp(-a*x),x=0..+infinity);