# Integrable and Superintegrable Systems with Higher Order Integrals of Motion: Master Function Formalism

Authors: Z. Alizadeh; H. Panahi
Journal: Advances in High Energy Physics
 (2018)

---

## Abstract

We construct two-dimensional integrable and superintegrable systems in terms of the master function formalism and relate them toMielnik’s andMarquette’s construction in supersymmetric quantum mechanics. For two different cases of the master functions, we obtain two different two-dimensional superintegrable systems with higher order integrals of motion.


---

## Body

## 1. Introduction

It is known from classical and quantum mechanics that a system withN degrees of freedom is called completely integrable if it allowsN functionally independent constants of the motion [1]. From the mathematical and physical point of view, these systems play a fundamental role in description of physical systems due to their many interesting properties. A system is superintegrable if one could obtain more thanN constants of the motion and if there exist 2
N
-
1 constants of the motion, the system is maximally superintegrable or just superintegrable [2–5]. Recently the study of superintegrable systems has been considered for different potentials and many researches have been studied for calculating the spectrum of these systems by different methods. In [6, 7], the spectrum of these systems has been calculated by an algebraic method using the realization of some Lie groups.For a two-dimensional quantum integrable system with HamiltonianH, there is always one operator like A
1 which commutes with Hamiltonian of the system, that is, [
H
,
A
1
]
=
0. For a quantum superintegrable system, one should define another operator such as A
2 which commutates with the Hamiltonian of system, that is, [
H
,
A
2
]
=
0, but [
A
1
,
A
2
]
≠
0. In other words, for a two-dimensional superintegrable system, there are two integrals of the motion (
A
1
,
A
2
) in addition to the Hamiltonian. The superintegrability with the second- and third-order integrals was the object of a series of articles [8–11]. The systems studied have second- and third-order integrals. Although superintegrability and supersymmetric quantum mechanics (SUSYQM) are two separated fields, many quantum systems, such as the harmonic oscillator, the Hydrogen atom, and the Smorodinsky-Winternitz potential, have both supersymmetry and superintegrable conditions [12–16]. These articles show that superintegrability is accurately connected with supersymmetry. For example, in [17], Marquette used the results obtained by Mielnik [18] and generated new superintegrable systems. Mielnik has shown that the factorization of second-order operators is not essentially unique. He has considered the Hamiltonian of the harmonic oscillator in one dimension as the simplest case:(1)
H
=
-
1
2
d
2
d
x
2
+
1
2
x
2
,where it can be factorized by two types of the first-order operators of creation and annihilation as follows:(2)
a
±
=
1
2
∓
d
d
x
+
x
,
b
±
=
1
2
∓
d
d
x
+
β
x
.For two superpartner Hamiltonians H
1 and H
2 where a
+
a
-
=
H
-
1
/
2
=
H
1 and a
-
a
+
=
H
+
1
/
2
=
H
2, he has demanded that H
2
=
b
-
b
+ and obtained the inverted product b
+
b
- as a certain new Hamiltonian:(3)
H
′
=
b
+
b
-
=
-
1
2
d
2
d
x
2
+
x
2
2
-
φ
′
x
,where φ
(
x
) is a function obtained from the general solution of Riccati equation considering β
=
x
+
φ
(
x
). The creation and annihilation operators of the third order for H
′ are described by expressions s
+
=
b
+
a
+
b
-, s
-
=
b
+
a
-
b
-, where a
+ and a
- are the creation and annihilation operators for H
2. Marquette [17] has taken the Hamiltonian H
2 in the x-axis and its superpartner H
′ given by (3) in the y-axis. Hence, he has obtained a two-dimensional superintegrable system as H
s
=
H
x
+
H
y, which can be separated in Cartesian coordinates with creation and annihilation operators a
+
(
x
), a
-
(
x
), s
+
(
y
), and s
-
(
y
). Also, he has shown that the Hamiltonian H
s possesses the following integrals of motion:(4)
K
=
H
x
-
H
y
,
A
1
=
a
+
x
s
-
y
-
a
-
x
s
+
y
,
A
2
=
a
-
x
s
+
y
+
a
+
x
s
-
y
,where these integrals are of order 2, 3, and 4 for shape invariant potentials [17].On the other hand, in [19, 20], the authors have shown that the second-order differential equations and their associated differential equations in mathematical physics have the shape invariant property of supersymmetry quantum mechanics. They have shown that by using a polynomial of a degree not exceeding two, called the master function, the associated differential equations can be factorized into the product of rising and lowering operators. The master function formalism has been used in relativistic quantum mechanics for solving the Dirac equation [21, 22].AsMielnik’s andMarquette’s method for generating superintegrable systems can be applied to other systems obtained in the context of supersymmetric quantum mechanics, hence, in this paper, we show that the supersymmetry method for obtaining the integrable and superintegrable systems can be related to master function formalism. In fact, we use the master function approach for 1-dimensional shape invariant potentials and generate 2-dimensional integrable systems. Also for a particular class of shape invariant systems, we generate 2-dimensional superintegrable systems. This class contains the harmonic oscillator, the singular harmonic oscillator, and their supersymmetric isospectral deformations.The paper is presented as follows: in Section2, we review how one can generate integrals of motion for two-dimensional superintegrable system from the creation and annihilation operators. In Section 3, we consider a particular quantum system for applying theMielnik’s andMarquette’s method and obtain a superintegrable potential separable in Cartesian coordinates. In Section 4, we briefly review the master function formalism and then in Section 5, we use this approach to obtain integrable systems and particular cases of the superintegrable systems that satisfy the oscillator-like (Heisenberg) algebra with higher order integrals of motion in terms of the master function and weight function. In Section 6, we give two examples to show how this method works in constructing oscillator-like two-dimensional superintegrable systems. Paper ends with a brief conclusion in Section 7.


## 2. Two-Dimensional Superintegrable System and Its Integrals of Motion

According to [17, 23, 24], for a two-dimensional Hamiltonian separable in Cartesian coordinates as(5)
H
x
,
y
,
p
x
,
p
y
=
H
x
x
,
p
x
+
H
y
y
,
p
y
,where the creation and annihilation operators (polynomial in momenta) A
+
(
x
), A
-
(
x
), A
+
(
y
), and A
-
(
y
) satisfy the following equations:(6)
H
x
,
A
-
x
=
-
λ
x
A
-
x
,
H
y
,
A
-
y
=
-
λ
y
A
-
y
,
H
x
,
A
+
x
=
λ
x
A
+
x
,
H
y
,
A
+
y
=
λ
y
A
+
y
,one can show that the operators f
1
=
A
+
m
(
x
)
A
-
n
(
y
) and f
2
=
A
-
m
(
x
)
A
+
n
(
y
) commute with the Hamiltonian H: that is,(7)
H
,
f
1
=
H
,
f
2
=
0
,if(8)
m
λ
x
-
n
λ
y
=
0
,
m
,
n
∈
Z
+
.Also the following sums of f
1 and f
2 commute with the Hamiltonian(9)
I
1
=
A
+
m
x
A
-
n
y
-
A
-
m
x
A
+
n
y
,
I
2
=
A
+
m
x
A
-
n
y
+
A
-
m
x
A
+
n
y
;that is, I
1 and I
2 are the integrals of motion. The order of these integrals of motion depends on the order of the creation and annihilation operators. On the other hand, the Hamiltonian H possesses a second-order integral as K
=
H
x
-
H
y, such that the integral I
2 is the commutator of I
1 and K. Thus, the Hamiltonian H is a superintegrable system and H, I
1, and K are its integrals of motion.


## 3. Mielnik-Marquette Method and Superintegrable Model Obtained from Shifted Oscillator Hamiltonian

In this section, for reviewing the Mielnik-Marquette method, we consider shifted oscillator Hamiltonian as(10)
H
=
-
d
2
d
x
2
+
1
4
ω
2
x
-
2
b
ω
2
-
ω
2
.We introduce the following first-order operators:(11)
a
-
=
d
d
x
+
1
2
ω
x
-
b
,
a
+
=
-
d
d
x
+
1
2
ω
x
-
b
,where the supersymmetric partner Hamiltonians are calculated as(12)
H
1
=
a
-
a
+
=
H
+
ω
,
H
2
=
a
+
a
-
=
H
.It is obvious that H
1 and H
2 have the shape invariant properties. Now, according to (2), we define the new operators b
- and b
+ such that(13)
H
1
=
H
+
ω
=
b
-
b
+
.The above equation gives the following Riccati equation:(14)
β
2
+
β
′
=
1
4
ω
2
x
2
-
b
ω
x
+
b
2
+
ω
2
,where a particular solution is(15)
β
x
=
β
0
x
=
1
2
ω
x
-
b
.Now, if we consider(16)
β
x
=
β
0
x
+
φ
x
,then we can obtain the following first-order linear inhomogeneous equation:(17)
z
′
+
-
2
β
0
z
=
1
,where z
=
1
/
φ
(
x
). After solving the above equation, we get(18)
φ
x
=
1
z
x
=
e
-
ω
/
2
x
2
+
2
b
x
π
/
2
ω
e
2
b
2
/
ω
erf
⁡
ω
/
2
x
-
2
b
/
ω
+
C
,where C is the constant of integration. Using the function φ
(
x
), we obtain(19)
H
′
=
b
+
b
-
=
H
1
-
φ
′
x
,where its creation and annihilation operators are given by the following expressions:(20)
s
+
=
b
+
a
+
b
-
,
s
-
=
b
+
a
-
b
-
.According to Marquette method, we take the x-axis for Hamiltonian H
1 and the y-axis for its superpartner H
′ and we have the following two-dimensional superintegrable system:(21)
H
s
=
H
x
+
H
y
=
H
1
+
H
′
=
-
d
2
d
x
2
-
d
2
d
y
2
+
1
4
ω
2
x
-
2
b
ω
2
+
1
4
ω
2
y
-
2
b
ω
2
-
ω
-
d
φ
d
y
.This Hamiltonian possesses the integral of motion given by (4), which are of order 2, 3, and 4.


## 4. The Master Function Formalism

According to [19, 20], the general form of the differential equation in master function approach is written as(22)
A
x
Φ
n
′
′
+
A
x
w
x
′
w
x
Φ
n
′
x
-
n
A
x
w
x
′
w
x
′
+
n
n
-
1
2
A
′
′
x
Φ
n
x
=
0
,where A
(
x
) as master function is at most a second-order polynomial and w
(
x
) is the nonnegative weight function in interval (
a
,
b
). By differentiating (22) m times and then multiplying it by (
-
1
)
m
A
m
/
2
(
x
), we get the following associated second-order differential equation in terms of the master function and weight function:(23)
A
x
Φ
n
,
m
′
′
+
A
x
w
x
′
w
x
Φ
n
,
m
′
x
+
-
1
2
n
2
+
n
-
m
2
A
′
′
+
m
-
n
A
x
w
′
x
w
x
′
-
m
2
4
A
′
x
2
A
x
-
m
2
A
′
x
w
′
x
w
x
Φ
n
,
m
x
=
0
,where(24)
Φ
n
,
m
x
=
-
1
m
A
m
/
2
d
d
x
m
Φ
n
x
.Changing the variable d
x
/
d
r
=
A
(
x
) and defining the new function Ψ
n
m
(
r
)
=
A
1
/
4
(
x
)
w
1
/
2
(
x
)
Φ
n
,
m
(
x
), one can obtain the Schrodinger equation as(25)
-
d
2
d
r
2
Ψ
n
m
r
+
v
m
x
r
Ψ
n
m
r
=
E
n
,
m
Ψ
n
m
r
,
m
=
0,1
,
2
,
…
,
n
,where the most general shape invariant potential is(26)
v
m
x
r
=
-
1
2
A
x
w
′
x
w
x
′
-
2
m
-
1
4
A
′
′
x
+
1
4
A
x
A
x
w
′
x
w
x
2
+
m
2
A
′
x
w
′
x
w
x
+
4
m
2
-
1
16
A
′
2
x
A
x
,and the energy spectrum E
(
n
,
m
) is as(27)
E
n
,
m
=
-
n
-
m
+
1
A
x
w
′
x
w
x
′
+
1
2
n
+
m
A
′
′
x
.According to [19, 20] the first-order deferential operators are written as(28)
A
±
=
∓
d
d
r
+
W
m
x
r
,where the superpotential W
m
(
x
(
r
)
) is expressed in terms of the master function A
(
x
) and weigh function w
(
x
) as(29)
W
m
x
r
=
-
A
x
w
′
x
/
2
w
x
+
2
m
-
1
/
4
A
′
x
A
x
.The Hamiltonian H
1 and H
2 called the superpartner Hamiltonians are written as(30)
H
1
=
A
-
A
+
=
-
d
2
d
r
2
+
W
m
2
r
+
W
m
′
r
=
-
d
2
d
r
2
+
v
1
r
,
H
2
=
A
+
A
-
=
-
d
2
d
r
2
+
W
m
2
r
-
W
m
′
r
=
-
d
2
d
r
2
+
v
2
r
,where v
1
(
r
) and v
2
(
r
) are called the partner potentials in the concept of supersymmetry in nonrelativistic quantum mechanics. Furthermore, if the partner potentials have the same shape and differ only in parameters, then potentials v
1
(
r
) and v
2
(
r
) are called the shape invariant potentials that satisfy (31)
v
1
r
,
a
0
=
v
2
r
,
a
1
+
R
a
1
,where R
(
a
1
) is independent of any dynamical variable and a
1 is a function of a
0. Potentials which satisfy in this condition are exactly solvable, although shape invariance is not the most general integrability or superintegrability condition.


## 5. Integrable and Superintegrable Systems Obtained from the Master Function Formalism

In this section, we try to relate the Mielnik-Marquette method to the master function approach. Hence, we define the following new operators:(32)
B
±
=
∓
d
d
r
+
ω
r
,where ω
(
r
) as the new superpotential must be related to the general form of the master function superpotential W
m
(
x
(
r
)
). Their product yields to Hamiltonians as(33)
B
-
B
+
=
-
d
2
d
r
2
+
ω
2
r
+
ω
′
r
,
B
+
B
-
=
-
d
2
d
r
2
+
ω
2
r
-
ω
′
r
.Now if we demand A
-
A
+
=
B
-
B
+ then we can obtain the following Riccati equation in terms of master function:(34)
ω
2
r
+
ω
′
r
=
W
m
2
r
+
W
m
′
r
,where a particular solution is ω
(
r
)
=
W
m
(
r
). The general solution can be obtained like(35)
ω
r
=
W
m
r
+
λ
r
,which yields(36)
λ
2
r
+
2
W
m
r
λ
r
+
λ
′
r
=
0
.We consider the transformation f
(
r
)
=
1
/
λ
(
r
) and obtain a first-order linear inhomogeneous differential equation as(37)
f
′
r
-
2
W
m
r
f
r
=
1
,where the general solution is(38)
f
r
=
exp
⁡
2
∫
W
m
r
d
r
C
+
∫
exp
⁡
2
∫
W
m
r
′
d
r
′
d
r
,where C is constant. Hence, (39)
ω
r
=
W
m
r
+
e
-
∫
2
W
m
r
d
r
C
+
∫
e
∫
2
W
m
r
′
d
r
′
d
r
.Using the function f
(
r
) given by (38), the superpartner Hamiltonian is given by(40)
H
′
=
H
2
-
λ
′
r
=
-
d
2
d
r
2
+
W
m
2
r
-
W
m
′
r
-
λ
′
r
,which is the general form of Hamiltonian in terms of master function. Now if we catch H
r
=
H
2 and H
r
′
=
H
′ (the Hamiltonian H
′ is thus given in terms of the variable r
′ vertical to r), then we obtain a new two-dimensional integrable Hamiltonian as(41)
H
s
=
H
r
+
H
r
′
=
-
d
2
d
r
2
-
d
2
d
r
′
2
+
W
m
2
r
+
W
m
2
r
′
-
W
m
′
r
-
W
m
′
r
′
-
λ
′
r
′
.Therefore, we have obtained the general form of the 2-dimensional integrable Hamiltonian in terms of master function which can be separated in radial coordinates. This separation of variable implies the existence of a second-order integral as K
=
H
r
-
H
r
′. Hence, H
s is an integrable system. Now, for generating superintegrable systems, we can obtain the creation and annihilation operators for H
′ from H
2 as(42)
S
+
=
B
+
A
+
B
-
,
S
-
=
B
+
A
-
B
-
,where A
± and B
± are given in (28) and (32). It is necessary to mention that these ladder operators satisfy (6) only for a particular class of shape invariant systems and in general, the 2-dimensional system H
s, obtained from a given master function, is not a superintegrable system. In other words, we cannot obtain 2-dimensional superintegrable system for all of the shape invariant cases given in [19, 20]. This particular class contains the harmonic oscillator, the singular harmonic oscillator, and their supersymmetric isospectral deformations.Hence, if relation (6) exists, according to (9) we can obtain the integrals of motion for Hamiltonian (41) as(43)
K
=
H
r
-
H
r
′
,
A
1
=
A
+
m
r
S
-
n
r
′
-
A
-
m
r
S
+
n
r
′
,
A
2
=
A
+
m
r
S
-
n
r
′
+
A
-
m
r
S
+
n
r
′
.In the next section, we apply this formalism for some particular cases of shape invariant potentials in terms of master function.


## 6. Examples of Two-Dimensional Superintegrable Systems as a Result of Master Function Approach

In this section, we would apply the master function formalism of the previous section for two examples and show how these results allow us to obtain 2-dimensional superintegrable systems with higher order integrals.Example 1.
LetA
(
x
)
=
1; then according to [19], w
(
x
)
=
e
-
(
β
/
2
)
x
2 that x
=
r
-
2
α
/
β, β
>
0 and the interval is (
-
∞
,
+
∞
). Using (29), we obtain the superpotential as(44)
W
m
r
=
-
β
2
r
-
2
α
β
.According to (27), the energy spectrum is as(45)
E
=
n
-
m
+
1
,and also the ladder operators given by (28) related to (44) satisfy a Heisenberg algebras (6). Now, substituting expression W
m
(
r
) in (38) yields the following relation in terms of the error function:(46)
λ
r
=
e
-
β
r
2
/
2
+
2
α
r
C
+
π
/
2
β
e
2
α
2
/
β
erf
⁡
β
/
2
r
-
2
α
/
β
,and so(47)
ω
r
=
W
m
r
+
λ
r
=
β
2
r
-
2
α
β
+
e
β
r
2
/
2
+
α
r
C
+
π
/
2
β
e
2
α
2
/
β
erf
⁡
β
/
2
r
-
2
α
/
β
.Substituting this expression in (40) and (41) yields the family of superpartner H
′ and a two-dimensional superintegrable Hamiltonian H
s, respectively, as(48)
H
′
=
H
2
-
λ
′
r
,
H
s
=
-
d
2
d
r
2
-
d
2
d
r
′
2
+
β
2
4
r
-
2
α
β
2
+
r
′
-
2
α
β
2
-
β
-
λ
′
r
′
,where(49)
λ
′
r
′
=
β
r
′
+
α
e
β
/
2
r
′
2
+
α
r
′
C
+
π
/
2
β
e
2
α
2
/
β
erf
⁡
β
/
2
r
′
-
2
α
/
β
-
e
β
/
2
r
′
2
+
α
r
′
e
2
α
2
/
β
e
-
1
/
2
β
r
′
-
2
α
/
β
2
C
+
π
/
2
β
e
2
α
2
/
β
erf
⁡
β
/
2
r
′
-
2
α
/
β
2
,
(50)
H
2
=
-
d
2
d
r
2
+
β
2
4
r
-
2
α
β
2
-
β
2
.It is seen that 2-dimensional superintegrable Hamiltonian (48) is the same as (21). We can find the general form of the operators S
+ and S
- in terms of the master function for this oscillator-like potential as follows:(51)
S
+
=
-
d
3
d
r
3
-
W
m
d
2
d
r
2
+
-
2
ω
′
-
W
m
′
+
ω
2
d
d
r
+
-
ω
′
′
-
W
m
′
ω
-
W
m
ω
′
+
ω
ω
′
+
W
m
ω
2
,
S
-
=
d
3
d
r
3
-
W
m
d
2
d
r
2
+
-
2
ω
′
-
W
m
′
-
ω
2
d
d
r
+
ω
′
′
-
W
m
′
ω
-
W
m
ω
′
+
ω
ω
′
+
W
m
ω
2
.Thus, we have obtained a 2-dimensional superintegrable system with integrals given by (43) as(52)
K
=
H
r
-
H
r
′
,
A
1
=
A
+
r
S
-
r
′
-
A
-
r
S
+
r
′
,
A
2
=
A
+
r
S
-
r
′
+
A
-
r
S
+
r
′
,where(53)
K
=
-
d
2
d
r
2
+
d
2
d
r
′
2
+
W
m
2
r
-
W
m
2
r
′
-
W
m
′
r
+
W
m
′
r
′
+
λ
′
r
′
,
A
1
=
2
W
m
r
d
3
d
r
′
3
+
2
W
m
r
′
d
3
d
r
d
r
′
2
-
2
-
2
ω
′
r
′
-
W
m
′
r
′
d
2
d
r
d
r
′
-
2
-
W
m
′
r
′
ω
r
′
-
W
m
r
′
ω
′
r
′
+
ω
r
′
ω
′
r
′
+
W
m
r
′
ω
2
r
′
d
d
r
-
2
W
m
r
ω
2
r
′
d
d
r
′
+
2
W
m
r
ω
′
′
r
′
,
A
2
=
-
2
d
4
d
r
d
r
′
3
+
2
ω
2
r
′
d
2
d
r
d
r
′
-
2
W
m
r
W
m
r
′
d
2
d
r
′
2
-
2
ω
′
′
r
′
d
d
r
+
2
W
m
r
-
2
ω
′
r
′
-
W
m
′
r
′
d
d
r
′
+
2
W
m
r
-
W
m
′
r
′
ω
r
′
-
W
m
r
′
ω
′
r
′
+
ω
r
′
ω
′
r
′
+
W
m
r
′
ω
2
r
′
.These integrals are of order 2, 3, and 4.Example 2.
According to [19] for A
(
x
)
=
x, we have w
(
x
)
=
x
α
e
-
β
x, x
=
r
2
/
4, α
>
-
1, β
>
0 and the interval is [
0
,
+
∞
). Now, using (29), the superpotential and the energy spectrum are as(54)
W
m
r
=
-
1
r
α
+
m
-
1
2
+
β
4
r
,
E
=
β
n
-
m
+
1
.This system has also the ladder operators that satisfy the form of (6); hence, substituting expression W
m
(
r
) in (38) yields the following relation in terms of Whittaker function:(55)
f
r
=
β
-
α
+
m
α
+
m
r
2
-
2
α
+
2
m
-
1
e
1
/
4
β
r
2
C
+
e
-
1
/
8
β
r
2
1
α
+
m
+
1
β
r
2
4
α
/
2
+
m
/
2
M
1
/
2
α
+
1
/
2
m
,
1
/
2
α
+
1
/
2
m
+
1
/
2
1
4
β
r
2
+
β
r
2
4
α
/
2
+
m
/
2
-
1
M
1
/
2
α
+
1
/
2
m
+
1
,
1
/
2
α
+
1
/
2
m
+
1
/
2
1
4
β
r
2
,where the Whittaker function M
μ
,
ν
(
z
) is the solution of the following differential equation:(56)
y
′
′
+
-
1
4
+
μ
z
+
1
/
4
-
ν
2
z
2
y
=
0
. It can be also defined in terms of the confluent hypergeometric function as(57)
M
μ
,
ν
z
=
e
-
1
/
2
z
z
1
/
2
+
ν
F
1
1
1
2
+
ν
-
μ
,
1
+
2
ν
,
z
.The family of superpartner Hamiltonians H
′ and the two-dimensional superintegrable Hamiltonian H
s are thus calculated by (40) and (41), respectively. The creation and annihilation operators for the Hamiltonian H
2 are as(58)
M
+
r
=
A
+
2
r
A
-
r
,
M
-
r
=
A
+
r
A
-
2
r
,where A
±
(
r
) is given in (28) and from (42), we have the creation and annihilation operators of the Hamiltonian H
′ as(59)
R
+
r
′
=
B
+
r
′
M
+
r
′
B
-
r
′
,
R
-
r
′
=
B
+
r
′
M
-
r
′
B
-
r
′
,where B
±
(
r
) is given by (32) and (39). We can also find the integrals of motion of the Hamiltonian H
s from (43) as(60)
K
=
H
r
-
H
r
′
,
A
1
=
M
+
r
R
-
r
′
-
M
-
r
R
+
r
′
,
A
2
=
M
+
r
R
-
r
′
+
M
-
r
R
+
r
′
,that are of the order 2, 7, and 8.


## 7. Conclusion

In this article, we have used the results obtained by Mielnik in the concept of SUSYQM and related it to master function formalism for constructing two-dimensional integrable and superintegrable systems with higher order integrals of motion. From this procedure, we have generated the superintegrable systems for two different cases of master functionsA
(
x
)
=
1 and A
(
x
)
=
x and have shown that the higher integrals of motion are in order 2, 3, 4 and 2, 7, 8, respectively.



---

*Source: 5647148-2018-06-21.xml*