Он дополняет команду rsolve ,
которая находится в основной библиотеке.
Пакет LREtools содержит следующий
список команд:
REcontent REcreate
REplot REprimpart REreduceorder
REtoDE REtodelta
REtoproc constcoeffsol
delta
dispersion hypergeomsols
polysols ratpolysols
shift
Рассмотрим пример, в котором
рекуррентное уравнение решается как командой
rsolve так и при помощи команд пакета LREtools , чтобы
проиллюстрировать разницу в методах. Отметим,
однако, что далеко не все уравнения, решаемые
командами LREtools решаются командой rsolve .
Пусть задано уравнение
> restart;rec :=
a(n+2)-(2*n+1)*a(n+1)/n+n*a(n)/(n-1)=n*(n+1):
rsolve(rec,a(n));
Так решает рекуррентное уравнение
команда rsolve:
Теперь воспользуемся пакетом LREtools .
Вначале необходимо создать структуру данных RESol
для представления решений рекуррентного
уравнения при помощи команды REcreate
> re1 := LREtools[REcreate](rec,
a(n), {});
И только после этого вводится команда
вывести решение, например, в виде полинома
(возможно также в виде рациональной дроби или
гипергеометрической функции).
> LREtools[polysols](re1);
Так решается линейное рекуррентное
уравнение при помощи пакета LREtools .
|