Jerk is the rate of acceleration change in physics; that is, the time derivative of acceleration, and as such the second velocity derivative, or the third time position derivative. The jerk is important in several mechanics and acoustics applications. The Jerk vector is here resolved into tangential-normal and radial-transverse components for planar motion, and the normal component is expressed as an affine differential invariant recognized as the aberrancy. Several geometric properties of the Jerk vector are established for plane motion using known aberrancy properties of curves  .
Nonlinear third-order differential equations, known as nonlinear Jerk equations, involving the third temporal displacement derivative, are of great interest in analyzing some structures which exhibit rotating and translating movements, such as robots or machine tools, where excessive Jerk leads to accelerated wear of transmissions and bearing elements, noisy operations and large contouring errors in discontinuities (such as corners) in the machining path  .
Many authors have studied the numerical solutions of the Jerk equation, harmonic balance approach to periodic solutions is used in  , in  they have written the high-order ordinary differential equation in terms of its differential invariants. New algorithm for the numerical solutions of nonlinear third-order differential equations was used jacobi-gauss collocation method in  , He’s variational iteration method was used in  for nonlinear Jerk equations. Modified harmonic balance method was used for nonlinear Jerk equations in  . In this paper, we consider a Jerk equation of the form
With initial conditions
. The most general function of Jerk, which is invariant in time-reversal and space-reversal and only nonlinear cubic, can be written as
where the parameters are
is constants. The current work is motivated by optimizing local truncation errors in order to find a hybrid point in a two-step block method to have the most accurate solution for Jerk Equation (1). We organize this paper as follows: The next section illustrates the method derivation, Section 3 presents the analysis of the method involving order four, Section 4 presents the numerical examples showing the productivity of the new technique when it is contrasted with the different strategies proposed in the scientific writing.
To derive two steps hybrid block method with one off-step point we have a polynomial of degree 6 as follows:
with third derivative given by
Substituting Equation (4) into Equation (3) gives
By interpolating Equation (3) at
and collocating Equation (5) at
we get a system of equations written in the matrix form
Solving the above system (by using Gaussian elimination method for solving the system of linear equations  ) gives us the coefficients of the polynomial
By making the substitution
, the polynomial in Equation (3) may be written in the form:
Now, by evaluating the solution approximation at the point
Assess the approximation at the point of first derivative
Evaluate the second derivative approximation at the point
We choose to optimize the local truncation errors in the Equation (8), Equation (9) and Equation (10) to determine appropriate values for r Equation (8), Equation (9) and Equation (10). This choice at the end of the block
, which result respectively in
Determine the r values equating to zero the
coefficients in the local truncation error formulas in Equation (11), and we obtain the system
By solving Equation (12) and substituting r’s in Equation (11) we can choose the value r that gives the least truncation errors
Substituting r in the local truncation errors gives:
To obtain a two-step hybrid block solution method for solving Equation (1), we evaluate
at the points
, we get the following block of six equations
3. Characteristics of the Method
This section is presented the basic properties of the main method and analyzed it to establish their validity.
3.1. Order of the Method
We can rewrite the hybrid block method in the form
are matrices of coefficients of dimensions
is a sufficiently differentiable function, the linear difference operator
associated with the implicit two-step block hybrid method is considered in Equation (8), Equation (9), Equation (10), Equation (13), that is given
are respectively the vector columns of the matrices
in Taylor series about
Note that the proposed method has order
at least  .
3.2. Zero Stability
We can write the method as a vector As
in Equation (14), we can write the method as a vector form.
The first characteristic polynomial
, the roots of polynomials are
. Hence the block method is zero -stable  .
The block method has order
, in case
, this is a sufficient condition to be consistent with the associated block method  .
We can establish the convergence of the two-step with three points hybrid block method if and only if it is consistent and zero stable  .
3.5. Region of Absolute Stability
As we mentioned earlier, zero-stability is a concept of the numerical method behavior for
. To decide whether a numerical method will produce good results with a given value of
, we need a concept of stability that is different from zero-stability. In most numerical methods intended to solve problems of third order, the stability properties are usually analyzed by considering the linear equation given by the Dalquist test  .
This problem has bounded solutions for
that tend to zero for
. We will define the region where the numerical method reproduces the manner of the exact solutions. Let us explain the procedure for obtaining such a region. Our method has nine equations in which there are four different terms of first derivatives:
, and second derivative:
and one intermediate values
. Let us depict the procedure to gain such a region  , We eliminated these terms from the equations system by using mathematica, and get a recurrence equation in the terms
. This recurrence equation reads
Through its characteristic equation, we study the extent boundedness of their solutions to define the region of stability. The root causes of the characteristic equation
The roots of the characteristic equation must be less than 1, for the method to be stable. The stability region for the method has shown in Figure 1.
4. Numerical Examples
Example 1: Consider the Jerk equation in the following form:
With initial conditions:
, and exact solution:
Tables 1-3 show the absolute Errors at
respectively. Example 1 was solved in  , the comparison between our method and result in  is proposed in Table 4 and Table 5. It is seen in Figures 2-4 the method gives a good approximation in the interval
Example 2: Consider the Jerk equation in the following form:
With initial conditions:
, and exact solution:
Figure 2. The solution of Example 1 at
Figure 3. The solution of Example 1 at
Figure 4. The solution of Example 1 at
Table 1. Absolute error of numerical solutions at
Table 2. Absolute error of numerical solutions at
Table 3. Absolute error of numerical solutions at
Table 4. Comparison between the numerical solution in our method and the method in  ,
Table 5. Comparsion between the numerical solution in our method and the method in  ,
Tables 6-8 show the absolute Errors at
respectively. Example 2 was solved in  , the comparison between our method and result in  is proposed in Table 9 and Table 10. It is seen in Figures 5-7 the method gives a good approximation in the interval
Figure 5. The solution of Example 2 at
Figure 6. The solution of Example 2 at
Figure 7. The solution of Example 2 at
Table 6. Absolute error of numerical solutions at
Table 7. Absolute error of numerical solutions at
Table 8. Absolute error of numerical solutions at
Table 9. Comparsion between the numerical solution in our method and the method in  ,
Table 10. Comparison between the numerical solution in our method and the method in  ,
A two-step with one hybrid point was proposed and proceeded as a self-starting method for solving nonlinear Jerk equations. We considered one hybrid point and specified for approximation after optimizing local truncation errors related to the main formula. Therefore, our method’s good convergent and stability properties make it attractive for the numerical solution of nonlinear problems. The method presented is zero stable, consistent and convergence of four-algebraic order. The numerical results and figures show their efficiency and precision compared to other methods in the literature.