In an old post, vv reported a bug in simpl/max, which has been "fixed" in Maple 2018. However, it seem that such repairs are not complete enough.
For example, suppose it is required to find the (squared) distance between the origin and a point on x^{3} - x + y^{2} = â…“ which is closest to the origin. In other words, one needs to minimize x²ï¼‹y² among the points on this curve, i.e.,

extrema(x^2 + y^2, {x^3 + y^2 - x = 1/3}, {x, y}, 's'); # in exact form

Unfortunately, an identical error message appears again:

This is about functionality introduced in Maple 2022, which however is still not well known: Integral Vector Calculus and parametrization using symbolic (algebraic) vector notation. Four new commands were added to the Physics:-Vectors package, implementing the parametrization of curves, surfaces and volumes, as well as the computation of path, surface and volume vector integrals. Those are integrals where the integrand is a scalar or vector function. The computation is done from any description (algebraic, parametric, vectorial) of the region of integration - a path, surface or volume.

The search query in the new Maple Application center is broken. There is no advanced search options and a search for mapleflow or maple flow brings up 0 results. There should at least be one found, for example the search should have at least brought up The Liquid Volume in a Partially-Filled Horizontal Tank".

Notation is one of the most important things to communicate with others in science. It is remarkable how many people use or do not use a computer algebra package just because of its notation. For those reasons, in the context of the Physics package, strong emphasis is put on using textbook notation as much as possible regarding input and output, including, for that purpose, as people here know, significant developments in Maple typesetting.

Still, for historical reasons, when using the Physics package, the labels used to refer to a coordinate system had been a single Capital Letter, as in X, Y, ...It was not possible to use, e.g. X', or x.

That has changed. Starting with the Maplesoft Physics Updates v.1308, any symbol can be used as a coordinate system label. The lines below demo this change.

Following the previous post on The Electromagnetic Field of Moving Charges, this is another non-trivial exercise, the derivation of 4-dimensional relativistic Lorentz transformations, a problem of a 3rd-year undergraduate course on Special Relativity whose solution requires "tensor & matrix" manipulation techniques. At the end, there is a link to the Maple document, so that the computation below can be reproduced, and a link to a corresponding PDF file with all the sections open.

Deriving 4D relativistic Lorentz transformations

Freddy Baudine^{(1)}, Edgardo S. Cheb-Terrab^{(2)}

(1) Retired, passionate about Mathematics and Physics

(2) Physics, Differential Equations and Mathematical Functions, Maplesoft

Lorentz transformations are a six-parameter family of linear transformations that relate the values of the coordinates of an event in one inertial reference system to the coordinates of the same event in another inertial system that moves at a constant velocity relative to the former. An explicit form of can be derived from physics principles, or in a purely algebraic mathematical manner. A derivation from physics principles is done in an upcoming post about relativistic dynamics, while in this post we derive the form of mathematically, as rotations in a (pseudo) Euclidean 4 dimensional space. Most of the presentation below follows the one found in Jackson's book on Classical Electrodynamics [1].

The problem is to find a group of linear transformations,

that represent rotations in a 4D (pseudo) Euclidean spacetime, and so they leave invariant the norm of the 4D position vector ; that is,

For the purpose of deriving the form of , a relevant property for it can be inferred by rewriting the invariance of the norm in terms of . In steps, from the above,

from where,

or in matrix (4 x 4) form, ,

where is the transpose of . Taking the determinant of both sides of this equation, and recalling that , we get

The determination of is analogous to the determination of the matrix (3D tensor ) representing rotations in the 3D space, where the same line of reasoning leads to . To exclude reflection transformations, that have and cannot be obtained through any sequence of rotations, because they do not preserve the relative orientation of the axes, the sign that represents our problem is +. To explicitly construct the transformation matrix , Jackson proposes the ansatz

Summarizing: the determination of consists of determining entering such that followed by computing the exponential of the matrix .

Determination of

In order to compare results with Jackson's book, we use the same signature he uses, , and lowercase Latin letters to represent space tensor indices, while spacetime indices are represented using Greek letters, which is already Physics' default.

(1)

Start by defining the tensor whose components are to be determined. For practical purposes, define a macro to represent the tensor and use to represent its components

(2)

(3)

Next, from (see above in Formulation of the problem) one can derive the form of . To work algebraically with representing matrices, set these symbols as noncommutative

(4)

From

(5)

it follows that

(6)

(7)

Expanding the exponential using , and taking into account that the matrix product can be rewritten as, the left-hand side of (7) can be written as

(8)

Multiplying by

(9)

Recalling that , and that for any matrix , ,

(10)

(11)

To allow for the combination of the exponentials, now that everything is in tensor notation, remove the noncommutative character of

(12)

(13)

Since every tensor component of this expression is real, taking the logarithm at both sides and simplifying tensor indices

(14)

(15)

So the components of

(16)

satisfy (15). Using TensorArray the components of that tensorial equation are

(17)

Simplifying taking these equations into account results in the form of we were looking for

(18)

This is equation (11.90) in Jackson's book [1]. By eye we see there are only six independent parameters in , or via

(19)

(20)

This number is expected: a rotation in 3D space can always be represented as the composition of three rotations, and so, characterized by 3 parameters: the rotation angles measured on each of the space planes . Likewise, a rotation in 4D space is characterized by 6 parameters: rotations on each of the three space planes, parameters and , and rotations on the spacetime planes, parameters . Define now using (18) for further computing with it in the next section

(21)

Determination of

From the components of in (18), the components of can be computed directly using the command. Then, following Jackson's book, in what follows we also derive a general formula for in terms of and shown in [1] as equation (11.98), finally showing the form of as a function of the relative velocity of the two inertial systems of references.

An explicit form of in the case of a rotation on the plane can be computed by taking equal to zero all the parameters in (19) but for and substituting in

(22)

(23)

Computing the matrix exponential,

(24)

(25)

This is formula (4.2) in Landau & Lifshitz book [2]. An explicit form of in the case of a rotation on the plane can be computed by taking equal to zero all the parameters in (19) but for

(26)

(27)

(28)

Rewriting

Following Jackson's notation, for readability, redefine the 6 parameters entering as

(29)

(Note in the above the surrounding backquotes '...' to prevent a premature evaluation of the left-hand sides; that is necessary when using the command.) With this redefinition, becomes

(30)

where each parameter is related to a rotation angle on one plane. Any Lorentz transformation (rotation in 4D pseudo-Euclidean space) can be represented as the composition of these six rotations, and to each rotation, corresponds the matrix that results from taking equal to zero all of the six parameters but one.

The set of six parameters can be split into two sets of three parameters each, one representing rotations on the planes, parameters , and the other representing rotations on the planes, parameters . With that, following [1], (30) can be rewritten in terms of four 3D tensors, two of them with the parameters as components, the other two with matrix as components, as follows:

(31)

(32)

The 3D tensors and satisfy the commutation relations

(33)

(34)

(35)

(36)

The matrix components of the 3D tensor , related to rotations on the planes, are

(37)

(38)

(39)

The matrix components of the 3D tensor , related to rotations on the 3D space planes, are

(40)

(41)

(42)

Verifying the commutation relations between and

The tensor is now expressed in terms of these objects as

(50)

where the right-hand side, without free indices, represents the matrix form of . This notation makes explicit the fact that any Lorentz transformation can always be written as the composition of six rotations

(51)

(52)

(53)

which is the same as the starting point (30)

The transformation , where , as a function of the relative velocity of two inertial systems

As seen in the previous subsection, in , the second term, , corresponds to 3D rotations embedded in the general form of 4D Lorentz transformations, and is the term that relates the coordinates of two inertial systems of reference that move with respect to each other at constant velocity . In this section, is rewritten in terms of that velocity, arriving at equation (11.98) of Jackson's book [1]. The key observation is that the 3D vector , can be rewritten in terms of , where and c is the velocity of light (for the rationale of that relation, see [2], sec 4, discussion before formula (4.3)).

Use a macro - say ub - to represent the atomic variable (this variable can be entered as `#mover(mi("β"),mo("ˆ")`. In general, to create atomic variables, see the section on Atomic Variables of the page 2DMathDetails).

(54)

(55)

With these two definitions, and excluding the rotation term we have

(56)

(57)

(58)

From this expression, the form of can be obtained as in (24) using and simplifying the result taking into account that is a unit vector