The usual unary and binary ring operations are available for multivariate polynomials.
For polynomial rings over fields division by elements of the coefficient field are allowed (with the result in the original polynomial ring). The operator div has slightly different semantics from the univariate case: if b divides a, that is, if there exists a polynomial q∈P such that a = b.q∈P then q will be the result of a div b, but if such polynomial does not exist an error results.
The list belows contains the general ring element predicates. Also, the IsDivisibleBy function allows a divisibility test, and the IsAlgebraicallyDependent function determines if a set of ring elements is algebraically dependent. Note that not all functions are available for every coefficient ring.
Given elements a, b in a multivariate polynomial ring P, this function returns whether the polynomial a is divisible by b in P, that is, if and only if there exists q∈P such that a=q.b. If true is returned, the quotient polynomial q is also returned.
Returns true iff the set S of multivariate polynomials is algebraically dependent.
Many of the functions in this subsection come in three different forms: one in which no variable is specified, which usually returns values in the coefficient ring, and two in which a particular variable is referred, either by name or by number, and these usually return values in the polynomial ring itself.
Given a multivariate polynomial f with coefficients in R, this function returns a sequence of `base' coefficients, that is, a sequence of elements of R occurring as coefficients of the monomials in f. Note that the monomials are ordered, and that the sequence of coefficients corresponds exactly to the sequence of monomials returned by Monomials(f).
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns a sequence of coefficients with respect to a given variable v=xi, that is, the function returns a sequence of elements of P that form the coefficients of the powers of v (in ascending order) when f is regarded as a polynomial ∑j cj xij; note that the variable xi itself will not occur in the coefficients. There are two ways to indicate with respect to which variable the coefficients are to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns the coefficient of vk=xik, that is, the function returns the element of P that forms the coefficient of the k-th power of xi, when f is regarded as a polynomial ∑j cj xij; note that the variable xi itself will not occur in the coefficient. There are two ways to indicate with respect to which variable the coefficient is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a multivariate polynomial f with coefficients in R, this function returns the leading coefficient of f as an element of R; this is the coefficient of the leading monomial of f, that is, the first among the monomials occurring in f with respect to the ordering of monomials used in P.
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns the element of P that forms the coefficient of the largest power of v=xi that occurs with non-zero coefficient in f, when f is regarded as a polynomial ∑j cj xij; note that the variable xi itself will not occur in the coefficient. There are two ways to indicate with respect to which variable the leading coefficient is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a multivariate polynomial f, return the length of f, i.e., the number of terms of f.
Given a multivariate polynomial f with coefficients in R, this function returns the trailing coefficient of f as an element of R; this is the coefficient of the trailing monomial of f, that is, the last among the monomials occurring in f with respect to the ordering of monomials used in P.
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns the element of P that forms the coefficient of the least power of v=xi that occurs with non-zero coefficient in f, when f is regarded as a polynomial ∑j cj xij; note that the variable xi itself will not occur in the coefficient. There are two ways to indicate with respect to which variable the leading coefficient is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a multivariate polynomial f with coefficients in Q, a rational function field or a field of fractions of an order of a number or function field, return the LCM of the denominators of the coefficients of f (CoefficientDenominator) or the product of that LCM with f itself (CoefficientNumerator).
Given a multivariate polynomial f and a monomial m, both in P∈R[x1, ..., xn], this function returns the coefficient with which m occurs in f as an element of R.
Given a multivariate polynomial f∈P, this function returns a sequence of monomials, that is, a sequence of monomial elements of P occurring in f. Note that the monomials in P are ordered, and that the sequence of monomials corresponds exactly to the sequence of coefficients returned by Coefficients(f).
Given a multivariate polynomial f∈P, this function returns parallel sequences C and M of the coefficients and monomials, respectively, of f. Thus this function is equivalent to calling Coefficients and Monomials separately, but is more efficient (particularly for large polynomials) since only one scan of the polynomial needs to be done.
Given a multivariate polynomial f∈P this function returns the leading monomial of f, that is, the first monomial element of P that occurs in f, with respect to the ordering of monomials used in P.
Given a multivariate polynomial f∈P, this function returns the sequence of (non-zero) terms of f as elements of P. The terms are ordered according to the ordering on the monomials in P. Consequently the i-th element of this sequence of terms will be equal to the product of the i-th element of the sequence of coefficients and the i-th element of the sequence of monomials.
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns a sequence of terms with respect to a given variable v=xi, that is, the function returns a sequence of elements of P that form the terms (ascending order) of f regarded as a polynomial ∑j cj xij. There are two ways to indicate with respect to which variable the terms are to be ordered: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns the k-th term of f (with k≥0), that is, the function returns the term of f involving the k-th power of xi, when f is regarded as a polynomial ∑j cj xij. There are two ways to indicate with respect to which variable the term is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a multivariate polynomial f∈P, this function returns the leading term of f as an element of P; this is the product of the leading monomial and the leading coefficient that is, the first among the monomial terms occurring in f with respect to the ordering of monomials used in P.
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns the element of P that forms the leading term of f when f is regarded as a polynomial ∑j cj xij. Thus it is the term involving the largest power of xi that occurs with non-zero coefficient. There are two ways to indicate with respect to which variable the leading coefficient is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a multivariate polynomial f∈P, this function returns the trailing term of f as an element of P; this is the last among the monomial terms occurring in f with respect to the ordering of monomials used in P.
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns the element of P that forms the trailing term of f when f is regarded as a polynomial ∑j cj xij. Thus it is the term involving the least power of xi that occurs with non-zero coefficient. There are two ways to indicate with respect to which variable the leading coefficient is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a single term f (a polynomial having exactly one term) in a polynomial ring of rank n, return the exponents of the monomial of f, as a sequence of length n of integers. (The coefficient of f is ignored; it need not be 1.)
Given a multivariate polynomial ring P = R[x1, ..., xn], and a sequence E of non-negative integers, return the monomial (x1)E[1] ... (xn)E[n] in P. This function is a semi-inverse of Exponents.
Given a length-k sequence C of coefficients in a ring R and a length-k sequence M of monomials of a polynomial ring R, return the multivariate polynomial f∈R whose coefficients are C and monomials are M. (Thus for any f∈R, Polynomial(Coefficients(f), Monomials(f)) equals f.)
> R<x, y, z> := PolynomialAlgebra(RationalField(), 3); > f := (2*x+y)*z^3+11*x*y*z+x^2*y^2; > f; x^2*y^2 + 11*x*y*z + 2*x*z^3 + y*z^3 > Coefficients(f); [ 1, 11, 2, 1 ] > Monomials(f); [ x^2*y^2, x*y*z, x*z^3, y*z^3 ] > CoefficientsAndMonomials(f); [ 1, 11, 2, 1 ] [ x^2*y^2, x*y*z, x*z^3, y*z^3 ] > Terms(f); [ x^2*y^2, 11*x*y*z, 2*x*z^3, y*z^3 ] > Coefficients(f, y); [ 2*x*z^3, 11*x*z + z^3, x^2 ] > Terms(f, 2); [ 2*x*z^3, 11*x*y*z + y*z^3, x^2*y^2 ] > MonomialCoefficient(f, x*y*z); 11 > LeadingTerm(f); x^2*y^2 > LeadingTerm(f, z); 2*x*z^3 + y*z^3 > LeadingCoefficient(f, z); 2*x + y > Polynomial([1, 2, 3], [x*y, y, z^2]); x*y + 2*y + 3*z^2
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns the degree of f in vk=xik, that is, the function returns the degree of f when it is regarded as a polynomial ∑j cj xij. The resulting integer is thus the largest power of xi occurring in any monomial of f. There are two ways to indicate with respect to which variable the degree is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P). If f is the zero polynomial, the return value is always -1.
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns the total degree of f, which is the maximum of the total degrees of all monomials that occur in f. The total degree of a monomial m is the sum of the exponents of the indeterminates that make up m. Note that this ignores the weights on the variables if there are any (see the section on graded polynomial rings below). If f is the zero polynomial, the return value is -1.
Given a multivariate polynomial f∈P = R[x1, ..., xn], this function returns the leading total degree of f, which is the total degree of the leading monomial of f. If f is the zero polynomial, the return value is -1.
Given a multivariate polynomial f∈R[x1, ..., xn], this function returns whether f is in fact a univariate polynomial in one of its indeterminates x1, ..., xn. If true is returned, then the function also returns a univariate version u of f and (the first) i such that f is univariate in xi. Note that there will only be ambiguity about i if f is a constant polynomial. The univariate polynomial u will be an element of R[x] with the same coefficients as f.
Given a multivariate polynomial f∈R[x1, ..., xn], this function returns whether f is in fact a univariate polynomial in xi. If true is returned, then the function also returns a univariate version u of f, which will be an element of the univariate polynomial ring R[x] with the same coefficients as f. The indeterminate xi should either be specified as a (polynomial) argument v or as an integer i.
Given a multivariate polynomial f∈R[x1, ..., xn], which is known to be a univariate polynomial in xi for some i with 1≤i≤n, return a univariate version u of f, which will be an element of the univariate polynomial ring R[x] with the same coefficients as f.
> P<x,y> := PolynomialRing(GF(5), 2); > f := x^2 - y + 3; > g := y^3 - x*y + x;If we compute the resultant in either variable of the two polynomials, then we can apply UnivariatePolynomial to this to obtain a univariate version of it, from which we can compute the roots.
> ry := Resultant(f, g, y); > ry; 4*x^6 + x^4 + x^3 + 3*x^2 + 2*x + 3 > Roots(UnivariatePolynomial(ry)); [ <3, 1> ] > Evaluate(f, x, 3); 4*y + 2 > Evaluate(g, x, 3); y^3 + 2*y + 3 > GCD($1, $2); y + 3 > rx := Resultant(f, g, x); > rx; y^6 + 4*y^3 + 3*y + 3 > Roots(UnivariatePolynomial(rx)); [ <2, 1> ] > Evaluate(f, y, 2); x^2 + 1 > Evaluate(g, y, 2); 4*x + 3 > GCD($1, $2); x + 2
Given a multivariate polynomial f∈P, return the derivative of f with respect to the variable v=xi, as an element of P. There are two ways to indicate with respect to which variable the derivative is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a multivariate polynomial f∈P and an integer k > 0, return the k-th derivative of f with respect to the variable v=xi, as an element of P. There are two ways to indicate with respect to which variable the derivative is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Given a multivariate polynomial f∈P over a field of characteristic zero, return the formal integral of f with respect to v=xi as an element of P. There are two ways to indicate with respect to which variable the integral is to be taken: either one specifies i, the integer 1≤i≤n that is the number of the variable (upon creation of P, corresponding to P.i) or the variable v itself (as an element of P).
Creates the matrix with (i, j)'th entry the partial derivative of the i'th polynomial in the list with the j'th indeterminate of its parent ring.
Given an element f of a polynomial ring P=R[x1, ..., xn] and a sequence or tuple s of ring elements of length n, return the value of f at s, that is, obtained by substituting xi=s[i]. If the elements of s can be lifted into the coefficient ring R, then the result will be an element of R. If the elements of s cannot be lifted to the coefficient ring, then an attempt is made to do a generic evaluation of f at s. In this case, the result will be of the same type as the elements of s.
Given an element f of a multivariate polynomial ring P and a ring element r return the value of f when the variable v=xi is evaluated at r. If r can be coerced into the coefficient ring of P, the result will be an element in P again. Otherwise the other variables of P must be coercible into the parent of r, and the result will have the same parent as r.
Let K be a field, and P=K[x1, ..., xn] a multivariate polynomial ring over K; let v=xi be the i-th indeterminate of P. Given a sequence I of elements of K (the interpolation points) and a sequence V of elements of P (the interpolation values), both sequences of length k>0, return the unique polynomial f∈P of degree less than k in the variable xi such that f(I[j]) = V[j], for j=1, ..., k. The variable xi may not occur anywhere in the values V. There are two ways to indicate with respect to which variable to interpolate: either one specifies i, the integer 1≤i≤n that is the number of the variable or the variable v itself (as a polynomial).
> Q := RationalField(); > P<x, y, z> := PolynomialRing(Q, 3); > f := Interpolation([Q | 1, 2, 3], [y, z, y + z], 1); > f; x^2*y - 1/2*x^2*z - 4*x*y + 5/2*x*z + 4*y - 2*z > [ Evaluate(f, 1, v) : v in [1, 2, 3] ]; [ y, z, y + z ]
The quotient of the multivariate polynomial f by g in R[x1, ..., xn], provided the result lies in P again. Here R must be a domain. If a polynomial q in P exists such that f=q.g then it will be returned, but if does not exist an error results.
The reductum of a polynomial f, which is the polynomial obtained by removing the leading term of f.
The reductum of a multivariate polynomial f∈R[x1, ..., xn] obtained by removing the leading term with respect to the variable v=xi. Here either v must be specified as a polynomial, or xi must be specified by providing the integer i, with 1≤i≤n.
We provide two basic tools that deal with polynomial diagonalization.
The symmetric bilinear form (as a matrix) of a multivariate polynomial of degree 2.
The diagonal form of the multivariate polynomial of degree 2. Also returns the transformation matrix.
> Q := RationalField(); > PR<x, y, z> := PolynomialRing(Q, 3); > g := 119/44*x^2 - 93759/41440*x*y + 390935/91427*x*z > + 212/243*x - 3/17*y^2 + 52808/172227*y*z > - 287/227*y + 537/934*z^2 - 127/422*z; > SymmetricBilinearForm(g); [ 119/44 -93759/82880 390935/182854 106/243] [ -93759/82880 -3/17 26404/172227 -287/454] [390935/182854 26404/172227 537/934 -127/844] [ 106/243 -287/454 -127/844 0] > DiagonalForm(g); 119/44*x^2 - 15798558582429/4*y^2 + 34932799628335074761085292707227419544217/934*z^2 - 176588732861018934524371210556883645619275217398116147234837710457404146371/2 > > bl := SymmetricBilinearForm(g); > NBL := Matrix(PR, bl); > D, T := OrthogonalizeGram(bl); > NT := Matrix(PR, T); > C := Matrix(PR, [[x,y,z,1]]); > NC := C * NT; > NCT := Transpose(NC); > (NC * NBL * NCT)[1][1] eq DiagonalForm(g); trueThe last few statements demonstrate how the polynomial's diagonal form is obtained from its symmetric bilinear form. Note also that since the polynomial g is not homogeneous its symmetric bilinear form is given on four variables, the fourth variable being a homogenizing variable.