Refine
Year of publication
- 2008 (4)
Document Type
- ZIB-Report (4)
Language
- English (4)
Has Fulltext
- yes (4)
Is part of the Bibliography
- no (4)
Keywords
- computer algebra (3)
- Crack (2)
- REDUCE (2)
- discrete equations (2)
- integrable systems (2)
- large polynomial systems (2)
- Euclid's algorithm (1)
- FORM (1)
- Reduce (1)
- conservation laws (1)
Institute
- Mathematical Optimization (4) (remove)
The paper describes a method for solution of very large overdetermined algebraic polynomial systems on an example that appears from a classification of all integrable 3-dimensional scalar discrete quasilinear equations $Q_3=0$ on an elementary cubic cell of the lattice ${\mathbb Z}^3$. The overdetermined polynomial algebraic system that has to be solved is far too large to be formulated. A probing' technique which replaces independent variables by random integers or zero allows to formulate subsets of this system. An automatic alteration of equation formulating steps and equation solving steps leads to an iteration process that solves the computational problem.
An algorithmic method using conservation law multipliers is introduced that yields necessary and sufficient conditions to find invertible mappings of a given nonlinear PDE to some linear PDE and to construct such a mapping when it exists. Previous methods yielded such conditions from admitted point or contact symmetries of the nonlinear PDE. Through examples, these two linearization approaches are contrasted.
We classify all integrable 3-dimensional scalar discrete affine linear equations $Q_3=0$ on an elementary cubic cell of the lattice ${\mathbb Z}^3$. An equation $Q_3=0$ %of such form is called integrable if it may be consistently imposed on all $3$-dimensional elementary faces of the lattice ${\mathbb Z}^4$. Under the natural requirement of invariance of the equation under the action of the complete group of symmetries of the cube we prove that the only ontrivial(non-linearizable) integrable equation from this class is the well-known dBKP-system.
The purpose of this paper is twofold. An immediate practical use of the presented algorithm is its applicability to the parametric solution of underdetermined linear ordinary differential equations (ODEs) with coefficients that are arbitrary analytic functions in the independent variable. A second conceptual aim is to present an algorithm that is in some sense dual to the fundamental Euclids algorithm, and thus an alternative to the special case of a Gr\"{o}bner basis algorithm as it is used for solving linear ODE-systems. In the paper Euclids algorithm and the new dual version' are compared and their complementary strengths are analysed on the task of solving underdetermined ODEs. An implementation of the described algorithm is interactively accessible at http://lie.math.brocku.ca/crack/uode.