Package Details - Maple Help

All Products Maple MapleSim

Home : Support : Online Help : Mathematics : Algebra : Polynomials : Groebner : Package Details

Details of the Groebner Package

	Calling Sequence
	Groebner:-command(arguments) command(arguments)

Description

•	This page contains complete information about the Groebner package, with the exception of computations of Groebner bases for modules, for which users are referred to Groebner Bases for Modules. For basic information on the Groebner package, see the Groebner help page.

•	The Groebner package is a collection of commands for doing Groebner basis calculations in skew algebras like Weyl and Ore algebras and in corresponding modules like D-modules. It can also be used in the case of usual commutative polynomials.

•	Whereas in the commutative case the ideals (respectively modules) are two-sided, in the skew case we focus on left-sided ideal (respectively modules). This corresponds to skew polynomials acting on functions on the left.

•	Each command in the Groebner package can be accessed by using either the long form or the short form of the command name in the command calling sequence.

List of Groebner Package Commands

•	The following is a list of available commands.

Basis	FGLM	HilbertDimension	HilbertPolynomial
HilbertSeries	Homogenize	InitialForm	InterReduce
IsBasis	IsProper	IsZeroDimensional	LeadingCoefficient
LeadingMonomial	LeadingTerm	MatrixOrder	MaximalIndependentSet
MonomialOrder	MultiplicationMatrix	MultivariateCyclicVector	NormalForm
NormalSet	RationalUnivariateRepresentation	Reduce	RememberBasis
Solve	SPolynomial	SuggestVariableOrder	Support
TestOrder	ToricIdealBasis	TrailingTerm	UnivariatePolynomial
Walk	WeightedDegree

•	To display the help page for a particular Groebner command, see Getting Help with a Command in a Package.

•

These commands make Groebner basis calculations and related facilities available in commutative algebras of polynomials and in skew algebras. For details, see Ore_algebra:-poly_algebra, Ore_algebra:-Weyl_algebra, and Ore_algebra. Such calculations are also available in the case of algebraic and modular coefficients. For the case of modules, see Computing Groebner Bases for Modules.

•

The MonomialOrder command declares a monomial order on an algebra and returns a table that describes it. This table is needed by other commands in the package. Alternatively, a short monomial order description can be used in simple cases (see type(...,ShortMonomialOrder)). The TestOrder command checks if two monomials are in a given monomial order. Many predefined monomial orderings are available, including total degree, lexicographic, elimination, and matrix-defined orderings.

•	The commands LeadingCoefficient, LeadingMonomial, and LeadingTerm compute the leading coefficient, leading monomial, and leading term of a polynomial respectively.

•

The NormalForm command computes the normal form of a polynomial modulo an ideal. The command Reduce returns the pseudo-remainder of a polynomial in the division by a list of polynomials. The InterReduce command inter-reduces a list of polynomials. Both the NormalForm and Reduce commands are usually used after a (reduced) Groebner basis has been computed.

•

The Basis command computes a (reduced) Groebner bases either in an algebra of skew polynomials (an Ore algebra) or in an algebra of usual polynomials. It also computes Groebner bases for modules over (usual or skew) polynomial rings. The RememberBasis command makes a Groebner basis known to the system without performing any computation.

•	In the commutative case, the Solve command preprocesses a set of generators in view of solving.

•	The SPolynomial command computes the S-polynomial of two polynomials.

•	The UnivariatePolynomial command finds univariate polynomials of least degree in a polynomial ideal.

•	The MultivariateCyclicVector command is a general-purpose and parametrizable iteration algorithm based on idea behind the FGLM algorithm.

•	The commands IsZeroDimensional and IsProper decide whether an ideal has finitely many solutions, or at least one solution, respectively.

•	The commands HilbertDimension, HilbertPolynomial, and HilbertSeries compute the Hilbert dimension, Hilbert polynomial, and Hilbert series of an ideal, respectively.

•	A typical Groebner basis computation is one of the following.

1. Define the appropriate algebra with one of: Ore_algebra:-poly_algebra, Ore_algebra:-diff_algebra, Ore_algebra:-shift_algebra, Ore_algebra:-qshift_algebra, or Ore_algebra:-skew_algebra.

2. Define the appropriate monomial order on the algebra with Groebner:-MonomialOrder.

3. Call Groebner:-Basis, or in simple cases call Groebner:-Basis with the short syntax.

•	A typical test of membership in an ideal is:

1. Perform a Groebner basis computation to find a Groebner basis of the ideal (usually with tdeg);

2. Call Groebner:-Reduce to test whether it returns zero.

•	A typical use of Groebner bases for elimination is:

1. Compute a Groebner basis with respect to an elimination monomial order (with the lexdeg syntax);

2. Select those polynomials in which the indeterminate to be eliminated has disappeared.

Note: Elimination of a single indeterminate between two polynomials can also be performed using Ore_algebra:-skew_elim.

•

When relevant, most commands of the Groebner package accept algebraic numbers and algebraic functions denoted by I, radicals, or RootOfs in their the inputs. Exceptions are Groebner:-MultiplicationMatrix, Groebner:-NormalSet, and Groebner:-Solve. Also when relevant, they accept PolynomialIdeal data structures in their input.

•	You can use the Groebner package in conjunction with the PolynomialIdeals package.

Note: Both packages have the following exports of the same name and almost identical functionality: HilbertDimension, IsProper, IsZeroDimensional, and UnivariatePolynomial. After loading both packages using with, these names refer to the corresponding exports of the most recently loaded package, but you can still access the exports with the same name in the other package using the long form, that is, by prepending the package name.

Examples

Short Versus General Syntax for Term Orders

>	$with (Groebner) &colon;$

Simple Groebner basis calculation using the short syntax:

>	$G ≔ map (numer, [\frac{s^{2}}{2 s + 1} x^{2} y - (s + 1) y, \frac{3 s - 1}{3} x y^{2} - \frac{1}{s - 1} x]) &colon;$

>	$Basis (G, tdeg (x, y))$

$[3 s^{2} x^{2} + (- 6 s^{4} - s^{3} + 7 s^{2} + s - 1) y^{2}, (3 s^{2} - 4 s + 1) y^{3} - 3 y, (3 s^{2} - 4 s + 1) x y^{2} - 3 x]$

(1)

More complicated one using the general syntax:

>	$with (Ore_algebra) &colon;$

>	$A ≔ poly_algebra (s, x, y, characteristic = 3, rational = s) &colon;$

>	$T ≔ MonomialOrder (A, tdeg (x, y)) &colon;$

>	$Basis (G, T)$

$[y^{2}, s^{2} x^{2} y + s^{2} y + 2 y]$

(2)

Lauricella Functions and Partial Differential Systems

The following example deals with Lauricella functions. It illustrates Groebner basis computation in D-modules with parameters in the ground field.

>	$n ≔ 2 &colon;$

>	$S ≔ add (x [i] D [i], i = 1 .. n)$

$S ≔ x_{1} D_{1} + x_{2} D_{2}$

(3)

>	$A ≔ skew_algebra (comm = \{a, b, seq (c [i], i = 1 .. n)\}, seq (diff = [D [i], x [i]], i = 1 .. n), polynom = \{seq (x [i], i = 1 .. n)\}) &colon;$

>	$P ≔ skew_product (S + a, S + b, A)$

$P ≔ D_{1}^{2} x_{1}^{2} + 2 D_{1} D_{2} x_{1} x_{2} + (a + b + 1) x_{1} D_{1} + D_{2}^{2} x_{2}^{2} + (a + b + 1) x_{2} D_{2} + a b$

(4)

The following operators introduce $n$ other parameters.

>	$for i to n do l [i] ≔ skew_product (D [i], x [i] D [i] + c [i] - 1, A) - P end do$

$l_{1} ≔ D_{1}^{2} x_{1} + D_{1} c_{1} - D_{1}^{2} x_{1}^{2} - 2 D_{1} D_{2} x_{1} x_{2} - (a + b + 1) x_{1} D_{1} - D_{2}^{2} x_{2}^{2} - (a + b + 1) x_{2} D_{2} - a b$

$l_{2} ≔ D_{2}^{2} x_{2} + D_{2} c_{2} - D_{1}^{2} x_{1}^{2} - 2 D_{1} D_{2} x_{1} x_{2} - (a + b + 1) x_{1} D_{1} - D_{2}^{2} x_{2}^{2} - (a + b + 1) x_{2} D_{2} - a b$

(5)

>	$M ≔ MonomialOrder (A, lexdeg ([seq (D [i], i = 1 .. n)], [seq (x [i], i = 1 .. n)])) &colon;$

The result is a third-order differential system of four equations.

>	$Basis ([seq (l [i], i = 1 .. n)], M)$

$[2 D_{1} D_{2} x_{1} x_{2} + D_{2}^{2} x_{1} x_{2} + D_{2}^{2} x_{2}^{2} + a D_{1} x_{1} + x_{2} D_{2} a + b D_{1} x_{1} + x_{2} D_{2} b - x_{1} D_{1} c_{1} - D_{2}^{2} x_{2} + c_{2} D_{2} x_{1} + a b + x_{1} D_{1} - D_{2} c_{2} + x_{2} D_{2}, D_{1}^{2} x_{1} - D_{2}^{2} x_{2} + D_{1} c_{1} - D_{2} c_{2}, 2 D_{2}^{3} x_{1}^{2} x_{2} - 4 D_{2}^{3} x_{1} x_{2}^{2} + 2 D_{2}^{3} x_{2}^{3} + 2 D_{1} D_{2} x_{1}^{2} a - D_{2}^{2} x_{1} x_{2} a + 3 D_{2}^{2} x_{2}^{2} a + 2 D_{1} D_{2} x_{1}^{2} b - D_{2}^{2} x_{1} x_{2} b + 3 D_{2}^{2} x_{2}^{2} b - 2 D_{1} D_{2} x_{1}^{2} c_{1} - 4 D_{1} D_{2} x_{1}^{2} c_{2} - 4 D_{2}^{3} x_{1} x_{2} - 4 D_{2}^{3} x_{2}^{2} + D_{2}^{2} x_{1} x_{2} c_{1} - 3 c_{1} D_{2}^{2} x_{2}^{2} + 2 D_{2}^{2} x_{1}^{2} c_{2} - 6 D_{2}^{2} x_{1} x_{2} c_{2} + x_{1} D_{1} a^{2} + x_{2} D_{2} a^{2} - 2 a x_{1} D_{1} b + 2 a b D_{2} x_{1} + 4 a x_{2} D_{2} b - 2 a D_{1} D_{2} x_{1} - 3 a D_{2}^{2} x_{2} - 3 c_{1} x_{2} D_{2} a - 3 a c_{2} D_{2} x_{1} + x_{1} D_{1} b^{2} + x_{2} D_{2} b^{2} - 2 b D_{1} D_{2} x_{1} - 3 b D_{2}^{2} x_{2} - 3 c_{1} x_{2} D_{2} b - 3 b c_{2} D_{2} x_{1} + 2 D_{1} D_{2} c_{1} x_{1} + 4 D_{1} D_{2} c_{2} x_{1} + 6 D_{1} D_{2} x_{1}^{2} - x_{1} D_{1} c_{1}^{2} + 2 D_{2}^{3} x_{2} + 3 D_{2}^{2} c_{1} x_{2} - 4 D_{2}^{2} c_{2} x_{1} - 2 D_{2}^{2} x_{2} c_{2} + 2 D_{2}^{2} x_{1}^{2} - 5 D_{2}^{2} x_{1} x_{2} + 9 D_{2}^{2} x_{2}^{2} + c_{1} c_{2} D_{2} x_{1} + a^{2} b + a b^{2} - 2 a b D_{2} - 3 c_{1} a b - a D_{2} c_{2} + 2 a D_{2} x_{1} + 6 x_{2} D_{2} a - b D_{2} c_{2} + 2 b D_{2} x_{1} + 6 x_{2} D_{2} b - 6 D_{1} D_{2} x_{1} + 2 x_{1} D_{1} c_{1} + 2 D_{2}^{2} c_{2} - 4 D_{2}^{2} x_{1} - 11 D_{2}^{2} x_{2} + 3 D_{2} c_{1} c_{2} - 3 c_{1} x_{2} D_{2} - 5 c_{2} D_{2} x_{1} + 3 a b - 2 a D_{2} - 2 b D_{2} - x_{1} D_{1} + 2 D_{2}^{2} - 3 D_{2} c_{2} + 2 D_{2} x_{1} + 5 x_{2} D_{2} - 2 D_{2}, 2 D_{1} D_{2}^{2} x_{2}^{2} - D_{2}^{3} x_{1} x_{2} + 3 D_{2}^{3} x_{2}^{2} - a D_{1} D_{2} x_{1} + 2 a D_{1} D_{2} x_{2} + a D_{2}^{2} x_{2} - b D_{1} D_{2} x_{1} + 2 b D_{1} D_{2} x_{2} + b D_{2}^{2} x_{2} - 2 D_{1} D_{2}^{2} x_{2} + D_{1} D_{2} c_{1} x_{1} - 4 D_{1} D_{2} c_{1} x_{2} + 2 D_{1} D_{2} c_{2} x_{1} + D_{2}^{3} x_{2} - 2 D_{2}^{2} c_{1} x_{2} - D_{2}^{2} c_{2} x_{1} + 4 D_{2}^{2} x_{2} c_{2} + 2 a b D_{1} - a b D_{2} - 2 a D_{1} c_{1} + 2 a D_{2} c_{2} - 2 b D_{1} c_{1} + 2 b D_{2} c_{2} - 2 D_{1} D_{2} c_{2} - 3 D_{1} D_{2} x_{1} + 6 D_{1} D_{2} x_{2} + 2 D_{1} c_{1}^{2} + D_{2}^{2} c_{2} - D_{2}^{2} x_{1} + 5 D_{2}^{2} x_{2} - 2 D_{2} c_{1} c_{2} + 2 a D_{1} - a D_{2} + 2 b D_{1} - b D_{2} - 4 D_{1} c_{1} + D_{2}^{2} + 4 D_{2} c_{2} + 2 D_{1} - D_{2}]$

(6)

Legendre Polynomials and Differential-Difference Linear Systems

The following example deals with Legendre polynomials: from a differential equation and a mixed differential-difference equation that define them, elimination of $Dx$ yields a difference equation that vanishes on them.

>	$unassign ('n','P') &colon;$

>	$DE ≔ (1 - x^{2}) {Dx}^{2} - 2 x Dx + n (n + 1)$

$DE ≔ (- x^{2} + 1) {Dx}^{2} - 2 x Dx + n (n + 1)$

(7)

>	$RDE ≔ (1 - x^{2}) Dx Sn + (n + 1) x Sn - (n + 1)$

$RDE ≔ (- x^{2} + 1) Dx Sn + (n + 1) x Sn - n - 1$

(8)

>	$with (Ore_algebra) &colon;$

>	$A ≔ skew_algebra (diff = [Dx, x], shift = [Sn, n], polynom = \{n, x\}) &colon;$

>	$with (Groebner) &colon;$

>	$M ≔ MonomialOrder (A, lexdeg ([Dx], [n, x, Sn])) &colon;$

>	$G ≔ Basis ([DE, RDE], M)$

$G ≔ [- {Sn}^{2} n^{2} + 2 Sn n^{2} x - 4 {Sn}^{2} n + 7 x Sn n - 4 {Sn}^{2} + 6 x Sn - n^{2} - 3 n - 2, - Dx Sn n + Dx n x - Dx Sn + x Dx + n^{2} + 2 n + 1, Dx {Sn}^{2} n + 2 Dx {Sn}^{2} - Dx Sn x - 2 Sn n^{2} - Dx n - 6 Sn n - Dx - 5 Sn, Dx Sn x^{2} - x Sn n - Dx Sn - x Sn + n + 1, {Dx}^{2} x^{2} - {Dx}^{2} + 2 x Dx - n^{2} - n]$

(9)

Continuing this application, we derive a proof that Legendre polynomials satisfy

>	$Sum (orthopoly [P] (n, z) u^{n}, n = 0 .. \infty) = \frac{1}{sqrt (1 - 2 z u + u^{2})}$

$\sum_{n = 0}^{\infty} orthopoly [P] (n, z) u^{n} = \frac{1}{\sqrt{u^{2} - 2 z u + 1}}$

(10)

Indeed, it follows that the summand $orthopoly [P] (n, z) u^{n}$ in the identity that we want to prove satisfies

>	$G ≔ map (numer @ normal, eval (convert (G, set), Sn = \frac{Sn}{u})) union \{u Du - n\}$

$G ≔ \{u Du - n, {Dx}^{2} x^{2} - {Dx}^{2} + 2 x Dx - n^{2} - n, Dx Sn x^{2} - x Sn n - Dx Sn - x Sn + n u + u, Dx n x u - Dx Sn n + x Dx u + n^{2} u - Dx Sn + 2 n u + u, Dx {Sn}^{2} n - Dx Sn x u - Dx n u^{2} - 2 Sn n^{2} u + 2 Dx {Sn}^{2} - Dx u^{2} - 6 Sn n u - 5 Sn u, 2 Sn n^{2} x u - {Sn}^{2} n^{2} + 7 x Sn n u - n^{2} u^{2} - 4 {Sn}^{2} n + 6 x Sn u - 3 n u^{2} - 4 {Sn}^{2} - 2 u^{2}\}$

(11)

in the new algebra $C (z, u) [n] [Dz, Du, Sn]$ . (Since the index of summation is $n$ , we consider polynomials in this variable to allow for its elimination.)

>	$with (Ore_algebra)$

$[Ore_to_DESol, Ore_to_RESol, Ore_to_diff, Ore_to_shift, annihilators, applyopr, diff_algebra, dual_algebra, dual_polynomial, poly_algebra, qshift_algebra, rand_skew_poly, reverse_algebra, reverse_polynomial, shift_algebra, skew_algebra, skew_elim, skew_gcdex, skew_pdiv, skew_power, skew_prem, skew_product]$

(12)

>	$A ≔ skew_algebra (diff = [Dx, x], diff = [Du, u], shift = [Sn, n], polynom = n) &colon;$

>	$M ≔ MonomialOrder (A, lexdeg ([n], [Dx, Du, Sn])) &colon;$

>	$GB ≔ Basis (G, M)$

$GB ≔ [- u x Du Sn + Dx Sn x^{2} + u^{2} Du - Dx Sn + u, - u^{2} {Du}^{2} + {Dx}^{2} x^{2} - 2 u Du - {Dx}^{2} + 2 x Dx, - u^{2} x {Du}^{2} - Du Dx u x^{2} + u {Du}^{2} Sn + u Dx Du - 3 u x Du - Dx x^{2} + Du Sn + Dx - x, - u Du + n]$

(13)

The previous equations are satisfied by the summand. Setting $Sn = 1$ in them yields equations satisfied by the sum. A final Groebner basis calculation returns these equations:

>	$GB ≔ Basis (eval (remove (has, GB, n), Sn = 1), M) &colon;$

>	$collect (GB, \{Du, Dx\}, distributed, factor)$

$[(u - x) u Du + (x - 1) (x + 1) Dx + u, u (u^{2} - 2 u x + 1) {Du}^{2} + (4 u^{2} - 5 u x + 1) Du + 2 u - x]$

(14)

Solving this simple system of PDEs yields the announced closed form.

Maple

Maple Add-Ons

MapleSim

MapleSim Add-Ons

Systems Engineering

Consulting Services

Maple T.A. and Möbius

Education

Industries

Automotive and Aerospace

Robotics

Machine Design & Industrial Automation

Other

Application Areas

Product Pricing

Purchasing

Institutional Student Licensing

Maplesoft Elite Maintenance (EMP)

Support

Product Training

Online Product Help

Webinars & Events

Publications

Content Hubs

Examples & Applications

Community

About Maplesoft

Media Center

User Community

Contact

Online Help

All Products Maple MapleSim

Maple

Powerful math software that is easy to use

Maple Add-Ons

MapleSim

Advanced System Level Modeling

MapleSim Add-Ons

Systems Engineering

Consulting Services

Maple T.A. and Möbius

Education

Industries

Automotive and Aerospace

Robotics

Machine Design & Industrial Automation

Other

Application Areas

Product Pricing

Purchasing

Institutional Student Licensing

Maplesoft Elite Maintenance (EMP)

Support

Product Training

Online Product Help

Webinars & Events

Publications

Content Hubs

Examples & Applications

Community

About Maplesoft

Media Center

User Community

Contact

Online Help

All Products Maple MapleSim