1.1. Models in Theoretical Physics
The subject we consider belongs to the modeling approach in theoretical physics ... “whose goal is to understand specific phenomena by developing either a mathematical or computational model. You begin this by choosing phenomena to study. Then you choose an approach to representing the phenomena; can you represent it as particle? a field? or some continuous distribution of matter? Then you choose a mathematical formulation .... You then adapt your approach to the mathematical formulation, thus developing a mathematical representation of your phenomena. You then use physical, mathematical, and/or computational arguments and methods to make predictions in the form of tables, plots, and/or formulas. By studying these results in different circumstances, you can extend our understanding of the phenomena. This is the most direct method of doing theoretical physics; it is a straight application of mathematical or computational methods. It is certainly the most structured way of doing theoretical physics ... A body of models linked by physical argument, derivation methods, and/or computer simulations constitutes a physical theory ... The first step in understanding any physics is to try to simplify the situation by removing all complications and then by working out all of the consequences of the situation. The particle is this kind of simplification. For such a simple explanation, it is very rich in principles and consequences. Once we have studied many simple ideas, we need to make them more realistic by reintroducing some of the complications that we removed in the process of simplification.” see G. E. Hrabovsky  .
Feynman’s advice  : “I think equation guessing might be the best method to proceed to obtain the laws for the part of physics which is presently unknown ... The problem is not to find the best or most efficient method to proceed to a discovery, but to find any method at all.”
Thus, we put forward the method for making simplified, serial models (S-models):
1) Considering a phenomenon, we chose the fundamental model that we will simplify;
2) Modifying it by inserting an auxiliary parameter for expansion by the Taylor series,
3) The S-models we make up of the Taylor polynomials, omitting auxiliary parameter.
There are three illustrative examples of such simplified, serial modeling:
a) Given physical phenomenon described by a complicated scalar function , we insert two auxiliary parameters and , and use the Taylor series
with to formulate the following two-term S-model
This linear model is simpler than and suggests two, yet simpler one-term models:
and . (I.3)
· Note that may not be equal to its Taylor series even if the series converges!
b) Given the convolution integral
we use the Taylor series expansion of . The first m + 1 terms provide the S-model
where are constants that are defined by the moments of specified by equation (I.4).
c) We can interpret the partial differential equations as the simplified serial models obtained from the Boltzmann integro-differential transport equations, see [  , Sect.4.3.1].
· As noted by Dirac, “A great deal of my work is just playing with equations and seeing what they give …” “A good deal of my research in physics has consisted in not setting out to solve some particular problem, but simply examining mathematical equations of a kind that physicists use and trying to fit them together in an interesting way, regardless of any application that the work may have. It is simply a search for pretty mathematics. It may turn out later to have an application. Then one has good luck.” see  .
· Simplified, serial models are convenient for detailed study of phenomena. Selecting a particular element of the fundamental model for inserting the auxiliary parameter, we concentrate our attention on the specific aspect of the considered phenomenon.
· A wide variety of fundamental models implies the same S-model! Thus, one can collate a variety of simplified models without knowing the actual fundamental model.
It makes sense to consider initially various simplified models for a phenomenon in order to infer eventually its fundamental model, see  .
In 2012, we presented modeling of the long-term, asymptotic dynamics of a point-like object where and is large, (arXiv: 1205.2920 [physics. gen-ph]). There we introduced some elements and concepts that we will now clarify with simple examples. In this paper, we consider simplified, serial models of the long-term dynamics, of a particular type of rigid bodies whose motion interacts with the surrounding medium. We name them “pointlike objects” (POs) and define them:
1) A PO is a classical extended real object that is moving through a medium. The response of PO to an external force is specified solely by the trajectory of a single point, which we name “the PO-position”. Therefore, the POs are a sort of rigid bodies.
2) The external force determines the acceleration of the PO-position.
3) The PO-velocity eventually stops changing after the cessation of the external force.
4) So far, one has modeled the dynamics of a PO by using so-called “point mass”, whose acceleration one has specified by Newton’s second law of motion, dividing the force acting on the point mass by its mass, its sole kinetic constant. However, a PO may be actually a very complicated real object such as a satellite, train, ship, chain, particle ..., thus we may oversimplify the situation by using a point mass to model its dynamics. Therefore, we approach this question by considering simplified serial models of the long-term, PO dynamics (LT-dynamics for short) for a small and slowly changing external force with a variable auxiliary parameter . Expanding the LT-dynamics in powers of , we obtain with simplified, serial models: the Taylor polynomials in terms of the external force and its derivatives, which we name “the LT-formulas”.
In Section 2, we consider the damped harmonic oscillator driven by the small and slowly changing external force with auxiliary parameter in order to calculate illustrative S-models for LT-dynamics where and is large, say . As the fundamental model, we introduce a special type of Newton’s second law formula for acceleration of PO, which we name the N-formula, where . It generalizes Newton’s second law by explicitly specifying the PO-acceleration by a possibly nonlinear transform of the external force. We calculate the N-formula for the driven damped harmonic oscillator from its differential equation of motion. By the expansion of this N-formula in powers of , we obtain the corresponding LT-formulas, the Taylor polynomials in time derivatives the external force . These LT-formulas suggest that by eliminating iteratively the higher time derivatives of the trajectory from the PO differential equation of motion, we can calculate the corresponding LT-formulas up to any order of , without solving the differential equation.
In Section 3, we consider S-models for the LT-dynamics of PO whose motion interacts with the surrounding medium. We calculate iteratively exemplary LT-formulas without solving the relevant differential equations: i) the Riccati differential equation for velocity at the quadratic drag force, ii) a second-order cubic nonlinear differential equation for the driven damped oscillator, and iii) the Lorentz-Abraham-Dirac relativistic differential equation.
In Section 4, we draw coclusion and comment on:
1) Specific S-models for the LT-dynamics.
2) Dynamic properties of PO, which serial models for relativistic LT-dynamics imply.
3) Applications of LT-formulas for predicting the long-term acceleration dependence of PO on the external force and its time derivatives: this dependence determines the relative significance of kinetic constants of the PO equation of motion for the LT-dynamics.
4) Figure 1 shows the block cycle diagram that illustrates the relations between the introduced concepts and techniques.
· Since physical theories employ mathematical models to describe and predict physical phenomena, our knowledge depends on the models available to that end. To increase their scope we use the Taylor series expansion in order to introduce simplified, serial models. Using the first few terms of a Taylor series, these models provide some information about difficult and complex problems. In this paper, we use them to describe the long-term dynamics of a special type of rigid bodies at small and slowly changing external forces.
2. LT-Dynamics of the Driven Damped Harmonic Oscillator
To illustrate the mathematical framework of the proposed S-models for of the LT-dynamics, we will now consider the S-models for the LT-dynamics of the damped harmonic oscillator that is driven by an external force with its magnitude and rate of change determined by the auxiliary parameter .
2.1. N-Formula for the Driven Damped Harmonic Oscillator
We calculate the N-formula from the differential equation of motion for the particular kind of PO. This PO is based on the point mass that is moving along the x-axis under the influence of the external force , , and initially at rest at . The point mass is attached to the zero-length spring with the force constant , and slowed down by the frictional force with the non-negative viscous damping coefficient c such that c > 0 when . Thus, the PO-position satisfies the differential equation of motion for the driven damped harmonic oscillator:
where the kinetic constant m specifies the PO inertial force, whereas the kinetic constants c and k specify the interaction between the surrounding medium, and PO-velocities’ and PO-position respectively. Therefore the PO-trajectory
with the Green function
is named “the un-damped angular frequency”, and is named “the damping ratio”. The PO-acceleration as a function of the external force is given by the following linear N-formulas:
if ; (4)
if ; (5)
if ; (6)
They transform the external force into PO-acceleration by integral operators; cf. the N-formulas (16), (17), and (18).
When the external force , the PO equation of motion (1) implies:
if ; (7)
if ; (8)
if ; (9)
where the four constants and are determined by .
Thus, after the cessation of the external force the PO-velocity eventually stops changing, and the properties of the external force within every finite period of time have negligible effects on the PO-velocity as because the PO differential equation of motion (1) is linear. If , then there is self-acceleration.
· To further illustrate the relationship between linear N-formulas and linear PO equations of motions we calculate N-formula for a PO made out of two connected point masses of equal mass , which are located on the x-axis, initially resting at points and : thus and . They are connected by the zero-length spring with the force constant k/2 > 0. The point mass with the trajectory is accelerated by the external force , , and slowed down by the frictional force with the viscous damping coefficient . Whereas the point mass with the trajectory is only slowed down by the frictional force . The differential equations of motion for this system of two connected point masses are
and . (10)
Therefore, the velocity of the first point mass is given by the formula:
if ; (11)
If c > 0, we may consider as the PO-trajectory, and the time differentiation of the formula (11) or (12) defines directly the corresponding N-formula. According to the differential equations of motion (10), the differential equation of motion for the PO-trajectory is as follows:
The kinetic constants m and k specify the PO, whereas the kinetic constant c specifies the interaction between the PO-velocity and the surrounding medium.
2.2. S-models for the Linear LT-Dynamics
We now calculate approximate models for the linear LT-dynamics implied by the linear N-formulas (4) and (5), when the external force i.e. is small and . Taylor series expansion of N-formulas (4) and (5) in powers of gives:
as if ;(14)
as if ; (15)
provided for . Thus, both responses of acceleration to the external force are slow, not instantaneous.
Whereas the differential equation of motion (1) implies that and each :
if ; (16)
if ; (17)
if . (18)
Above formulas (14) and (15) show how using the S-models for the LT-dynamics specified by the PO equation of motion (1) with , we can obtain the long-term PO-acceleration expressed as a sum of time derivatives of the small and slowly changing external force , say,
provided for ; e.g. if , with . Such S-model for the LT-dynamics under small and slowly changing external force we name as the LT-formula of the order of . The real constants kn we name as the LT-constants. About the LT-formula (19) we note that
as , . (20)
If for , and , , then each LT‑formula (19) of the order implies iteratively the following novel differential equation of the order for the LT-dynamics: the LT-equation
and vice versa. When N = 4 and , the constants of the LT-equation (21) and the LT-constants of the LT-formula (19) are related iteratively:
, and ; (22)
, and . (23)
2.3. On Elements and Usage of Serial Models
1) The Equations (1), (4), (14), (19), and (21) illustrate the mathematical relations between the following elements of the proposed S-models:
a) PO differential equation of motion, b) N-formula, c) LT-dynamics, d) LT-formulas, and e) LT-equations; see Figure 1 the block cycle diagram in Section 4.
2) The LT-formulas (19) contain only one time derivative of the PO-trajectory, the LT-equations (21) contain only one time derivative of the external force, while the PO differential equations of motion are is not so constrained, cf. equation (13).
3) The Equations (14) and (15) exemplify how on expressing the LT-constants kn of the LT-formula (19) in terms of the kinetic constants m, c, and k, the S-models bring out the relative significance of these kinetic constants for the LT-dynamics.
4) The early time, , start-up dynamics of the N-formula (4) is given by
But the way in which the LT-formula (14) of the LT-dynamics depends on the kinetic constants m, c, and k is fundamentally different.
5) The Equations (4) and (14) suggest that the N-formula and the corresponding LT-formulas may differ significantly. Nevertheless, when the external force is small and slowly changing, we may use the LT-formula to calculate the approximate long-term PO-trajectories.
6) In contrast to the differential equation of motion (1) that depends continuously on the kinetic constants m, c, and k, in view of the Equations (14) and (15), the corresponding LT-formulas may not. The same applies for the LT-equations (21), which in general do not determine the PO differential equation of motion (1).
7) The Equations (15) to (18) and the estimate (20) suggest that without solving a PO equation of motion like (1) or (13), we can calculate the corresponding LT-formulas up to any order of λ by eliminating iteratively the higher time derivatives of trajectory. So assuming that the estimate (20) is correct, we calculate in that way from the differetial equation of motion (13) for the PO consisting of two connected point masses the corresponding LT-formula:
8) The formulas (7) to (9) show that for the particular PO, its initial response to an external force becomes exponentially negligible. Moreover, if the PO-acceleration is cyclic then there is no start-up, early time dynamics: in this case, the long-term dynamics and dynamics of PO are identical!
3. Exemplary LT-Formulas
In this section, where with , we calculate iteratively from the assorted differential equations some exemplary LT-formulas in terms of the time derivatives of a function of the external force . Any LT-formula or LT-equation that we do not calculate from some PO equation of motion we designate as a hypothetical one.
3.1. LT-Formulas for the Strong String
Using the strong string, we generalize the PO differential equation of motion (1):
where ; and if . Thus the three kinetic constants c, and k1 and k3, specify the reaction of the surrounding medium to the PO-velocity and PO-position respectively. To calculate the Taylor polynomial of the differential equation (26) we rewrite it by Cardano’s formula:
with , , and . We presume that the equation (26) is still a PO equation of motion if and that the estimates
are true also if , calculate the Taylor polynomial of as a function of r, and eliminate the time derivatives . This way we obtain the LT-formula:
at . (29)
This formula exemplifies a new type of S-models for the LT-dynamics with polynomials in time derivatives of the function of the external force .
3.2. The Quadratic Drag Force
Let us consider a PO with mass , which is moving along the x-axis through a fluid at relatively large velocity under the influence of the external force , and slowed down by the Lord Rayleigh type of the quadratic drag force , . The PO-mass m specifies the PO inertial force, whereas the constant specifies the reaction force of the surrounding fluid. So we presume that the PO-velocity satisfies the following Riccati differential equation:
Hypothesizing that the differential Equation (30) is an LT-equation and that the estimates
are true also for , we get iteratively a new type of serial models for LT-dynamics. They are polynomials in time derivatives of :
We can generalize the calculated hypothetical LT-formula (32) by adding the frictional force to the hypothetical LT-equation (30).
3.3. The Nonlinear Relativistic LT-Formulas
Let us consider LT-formulas of a relativistic N-formula. We base them on relativistic point mass, which is located at and moving with velocity , under the influence of the external force with the dimensionless auxiliary parameter . We define the external four-force
with , (33)
where . We will use the PO four-velocity ; and the metric with the signature , so that . We introduce an additional four-force , which specifies the properties of PO-dynamics in the case of the external four-force , and formulates the relativistic N-formula:
with , (34)
where t/γ is the proper time. As and , we may rewrite the relativistic N-formula (34) as
Generalizing the linear LT-formula (19), we model the dependence of the four-force on the external four-force by a relativistic polynomial in time derivatives of the external four-force, to get the hypothetical relativistic LT-formula:
where the real parameters are independent of the external four-force . We name them “the LT-constants” when we use them for a particular PO to specify its hypothetical relativistic LT-formula (36) up to the order of inclusive.
· If , then N-formula (35) is invariant under time reversal and the hypothetical relativistic LT-formula (36) has .
· The N-formulas (17) and (18) themselves suggest two hypothetical relativistic LT-formulas: [rhs(36)]  and [rhs(36)]  .
· By eliminating iteratively all the time derivatives except the one of from the hypothetical relativistic LT-formula (36), we get a hypothetical relativistic LT-equation, analogous to the LT-equation (21).
3.4. The Relativistic LT-Formulas for an Electrified PO
Presuming that the PO is electrified by a pointlike charge, we follow Schott  and express the additional four-force as the difference
where , , is the long-term intensity of the four-momentum emitted by the Liénard-Wiechert potentials with the cyclically moving singularity at , see (  , §6.6), and is the time derivative of an “acceleration four-momentum ”. Therefore, for an electrified PO we rewrite the relativistic N-formula (34) as
The kinetic constant d specifies the magnitude of radiation reaction force of the surrounding vacuum, which opposes the acceleration of an electrified PO. According to Dirac  , the acceleration four-momentum may be any four-function of the time derivatives ; and the N-formula (38) for an electrified PO conserves the four-momentum by the Equation (39). Furthermore, Bhabha  pointed out that if is such that the cross product
is a total differential with respect to the proper time, then such an electrified PO conserves the angular four-momentum!
Inspired by the LT-equation (21), we assume that the nth time derivative is of the order as , and then model the time derivative in the relativistic N-formula (38) by the relativistic polynomials in , subject to the Bhabha condition (40), cf. (  , Ch.9) and  . Accordingly, the hypothetical, relativistic LT-equation for an electrified PO equals the LT-equation
up to the order of inclusive, disregarding the Bhabha condition (40). Assuming the Bhabha condition (40), we gave such a hypothetical, relativistic LT-equation up to the order of inclusive, see  . On eliminating iteratively from the (41) we get
which is the LT-equation such as the LT-equation (36) up to the order of . Up to the order of inclusive, we gave such relativistic LT-equation in (  , Sect.11.4).
3.5. On Equation of Motion for a Charged Particle
1) In 1892, H. A. Lorentz started an ongoing quest to take account of the radiation reaction force (the effect of the loss of four-momentum by the electromagnetic radiation) by the classical equation of motion for a charged particle. In 1938, Dirac assumed that an electron is such a simple thing that the lowest order hypothetical LT-equation (41) ought to be the correct equation of motion with , and no additional polynomial terms are needed, cf.  . The equation (41) is named the Lorentz-Abraham-Dirac equation. Since it exhibits the self-acceleration, it has baffled the mathematical physicists ever since the Dirac invented it. Since we obtained it as the first known hypothetical relativistic LT-equation about the LT-dynamics of an electrified PO, we puzzled out its significance. An equation about long-term dynamics, where , is not necessarily an equation of motion, where , the initial dynamics is missing.
2) In 2008, Rohrlich  stated that the physically correct equation of motion for a classical charged particle is the lowest order relativistic LT-formula (42) for an electrified PO with , provided , see  for comment.
There is a century old discussion with an infinite number of proposals about the appropriate equation of motion for an electrified PO, see e.g.    and the references cited therein. Nevertheless, we do not consider here the theoretical problem about the properties of a real object and conditions under which it behaves largely like PO, and we may idealize it as PO by constructing an appropriate PO equation of motion.
4.1.1. Multi-Simplified Serial Model
To simplify an implicit relation we modify it by inserting six parameters:
where three auxiliary parameters , , and point to the variables , , and the constant c respectively, whereas the parameters , , and specify their starting values. The Taylor series to nth order in auxiliary parameters , , and provides the nth order multi-simplified, serial model.
4.1.2. Specific S-Models for LT-Dynamics
Usage of S-models requires selecting and identifying the relevant aspects of the LT-dynamics. Instead of using the modified external force , we might obtain more detailed information about the significance of kinetic constants for LT-dynamics of a given PO by using the modified external force with two dimensionless, auxiliary parameters λ and ε, to control separately the magnitude and the rate of change of the external force. Considering the same PO, we can start also with various simplified equations of motion. Thereby we can obtain Taylor polynomials in time derivatives of different functions of the external force, e.g. for the equation (26) there are three possibilities: , , or . The resulting LT-formulas provide information about the significance of individual components of PO for its LT-dynamics !
4.1.3. General Relativistic Properties of the LT-Dynamics
According to the hypothetical relativistic LT-formula (36), the LT-dynamics of a PO is specified up to the order of inclusive by the four LT-constants and . According to Einstein, the first LT-constant is . In addition, within Dirac’s classical theory of radiating electrons, the second LT-constant is . Therefore, in general, we expect that the LT-constant is determined by the intensity of the loss of the PO four momentum in response to the time derivative of the external force, it is not negative as a real PO may provide only a finite amount of the four-momentum. Therefore, the first two terms of the relativistic LT-formula (36) provide for a classical electrified PO the definitions of its mass and charge by their role in the relativistic LT-dynamics.
The physical interpretation of two third-order LT-constants and is open. Under the Bhabha condition (40), they are related as follows:
and , (44)
where is a real constant.
Using the LT-formula (36) with increasing , we can simulate how at a given precision of our observations the number of the observable LT-constants increases with the magnitude and rate of change of the external four-force : at very small and slowly changing external force, just the PO-mass can be determined. Due to the multitude of actual POs, we see no physical reason to believe that there is only a limited number of the independent relativistic LT-constants.
4.1.4. Application of the LT-Formula
1) The N-formulas (17) and (18) point out exceptional cases with POs whose LT-formula equals their N-formula.
2) As we can hardly ever obtain the exact solutions in closed form for the PO equation of motion, the corresponding LT-formulas provide a welcome source of information about LT-dynamics. Whenever each following term of the LT-formula is essentially smaller than the preceding one, we may expect that this formula will provide appropriate information about the LT-dynamics. The LT-constants provide the significance of the individual kinetic constants of the PO equation of motion for LT-dynamics.
3) A hypothetical LT-equation implies hypothetical LT-formula of the same order of λ: We may always use it to calculate approximations to the long-term PO-trajectories because LT-formula cannot exhibit self-acceleration, on the contrary to PO-equation.
4) Certain PO moving through a medium might present a too complicated system to create its exact equation of motion. But we may still be able to somewhat describe its LT-dynamics by a hypothetical LT-equation which balances the external force with the sum of the inertial force and an S-model of the interaction force between PO and the surroundings, cf. the hypothetical LT-equation (30) and the Lorentz-Abraham-Dirac Equation (41).
5) To make a hypothetical LT-formula we can use as a generic ansatz some multivariate polynomials in time derivatives of an appropriate function of the external force such as given by the Equations (19), (29), (32) or (36), and extract the values of their parameters by multiple linear regressions from data about acceleration the of long-term PO-trajectories, see  . Thus without knowing an adequate N-formula, we can make appropriate LT-formulas for predicting the LT-dynamics of a given PO.
6) There are many real systems consisting of POs, each of which is treated as the point mass with acceleration specified by Newton’s second law, e.g. in astronomy, and in classical mechanics. Using appropriate LT-formulas instead of Newton’s second law, we could take into account not only the PO-masses but also some of the additional kinetic properties of these POs. That way we might get better dynamic models of such systems.
7) Basing models of continuous mechanical medium on the laws about the interaction of the point masses and Newton’s second law (see [  , Sect. 4.4]), and using an appropriate LT-formula instead of Newton’s second law, one might get better models.
8) Remaining subjects: a) Memory: PO-acceleration after cessation of the external force, b) Passive damping: energy dissipation by parts of a composite PO. c) Composite dynamics: the LT-dynamics of a PO in terms of the LT-formulas of its parts.
We presented a particular type of simplified models. Using few terms of the Taylor series, these models provide some information about difficult and complex problems. In this paper, they describe the long-term dynamics of a special type of bodies at small and slowly changing external forces. We introduced:
1) A pointlike physical object PO, a classical extended object whose motion interacts with the surrounding medium. PO response to an external force is aptly specified solely by the trajectory of a single point, whose velocity eventually stops changing after the cessation of the external force.
2) A particular type of PO-dynamics: if the PO-acceleration is cyclic then LT-dynamics equals to the basic dynamics, there is no start-up.
3) The key mathematical element is the N-formula. It generalizes Newton’s second law formula for acceleration by explicitly specifying the PO-acceleration by the external force.
4) The simplified models of the LT-dynamics under a small and slowly changing external force, i.e. the LT-formulas, approximate the long-term dependence of PO-acceleration at a given time instant by the Taylor polynomial in time derivatives of the external force at the same time instant.
Given an ordinary PO differential equation of motion, we can calculate iteratively the corresponding LT-formula of any order of λ! Nevertheless, the LT-formula does not need to imply the original PO equation of motion.
Each LT-formula implies iteratively a novel differential equation of the same order of λ about the LT-dynamics (LT-equation), and vice versa. Such an LT-equation provides certain information about the LT-dynamics, and may exhibit self-acceleration. Different POs may have identical LT-formulas of the same order of λ.
To illustrate the mathematical framework and usage of serial models of the long-term dynamics of the pointlike objects, we considered the driven damped harmonic oscillator in Section 2. As Figure 1 the block cycle diagram shows, we illuminated the relations between various elements and techniques for modeling the PO long-term dynamics, where and , i.e. is large. Thus, the PO long-term dynamics is actually a generalization of the asymptotic dynamics, where .
4.3. Main Points
1) Simplification through Taylor series expansion. We consider simplification of mathematical models by Taylor series expansion of their elements. According to Planck  “the simpler the presentation of a particular law of Nature, the more general it is though at the same time, which formula to take as the simpler, is a problem which cannot always be confidently and finally decided.”
2) Modifications of an initial ansatz through its parameter. The tacit basis of our approach to model making is the fact that when the given ansatz model is an invertible function of a particular parameter; we can modify it through this parameter into an arbitrary ansatz.
Example: Given the three parameter ansatz
Figure 1. The simplified modeling of the long-term PO dynamics.
, , (45)
then for any function we get
3) To model small and slowly changing force we use with small variable . The results are useful with if function in itself is small and slowly changing.
4) We introduced pointlike objects and serial models as the modeling elements. To illustrate them we provided various examples of heuristic approach to modeling the pointlike objects dynamics. This might prevent future mix-up over the proper usage of the long-term equations and various differential equations for modeling dynamics of the pointlike objects. We put forward this novel approach as well we knew, yet “Proof of the pudding is in eating”.
We express our thanks to S. Ribarič and to the JAMP reviewer for many helpful comments.
We appreciate any comments, references, suggestions, and viewpoints.