
The deterministic Degasperis-Procesi equation admits weak multi-shockpeakon solutions of the form
u(x,t)=n∑i=1mi(t)e−|x−xi(t)|−n∑i=1si(t)sgn(x−xi(t))e−|x−xi(t)|,
where sgn(x) denotes the signum function with sgn(0)=0, if and only if the time-dependent parameters xi(t) (positions), mi(t) (momenta) and si(t) (shock strengths) satisfy a system of 3n ordinary differential equations. We prove that a stochastic perturbation of the Degasperis-Procesi equation also has weak multi-shockpeakon solutions if and only if the positions, momenta and shock strengths obey a system of 3n stochastic differential equations.
Citation: Lynnyngs K. Arruda. Multi-shockpeakons for the stochastic Degasperis-Procesi equation[J]. Electronic Research Archive, 2022, 30(6): 2303-2320. doi: 10.3934/era.2022117
[1] | Qian He, Wenxin Du, Feng Shi, Jiaping Yu . A fast method for solving time-dependent nonlinear convection diffusion problems. Electronic Research Archive, 2022, 30(6): 2165-2182. doi: 10.3934/era.2022109 |
[2] | Dongdong Ruan, Xiaofeng Wang . A high-order Chebyshev-type method for solving nonlinear equations: local convergence and applications. Electronic Research Archive, 2025, 33(3): 1398-1413. doi: 10.3934/era.2025065 |
[3] | Junseok Kim . Maximum principle preserving the unconditionally stable method for the Allen–Cahn equation with a high-order potential. Electronic Research Archive, 2025, 33(1): 433-446. doi: 10.3934/era.2025021 |
[4] | Chang Hou, Hu Chen . Stability and pointwise-in-time convergence analysis of a finite difference scheme for a 2D nonlinear multi-term subdiffusion equation. Electronic Research Archive, 2025, 33(3): 1476-1489. doi: 10.3934/era.2025069 |
[5] | Shuaikang Wang, Yunzhi Jiang, Yongbin Ge . High-order compact difference methods for solving two-dimensional nonlinear wave equations. Electronic Research Archive, 2023, 31(6): 3145-3168. doi: 10.3934/era.2023159 |
[6] | Lili Li, Boya Zhou, Huiqin Wei, Fengyan Wu . Analysis of a fourth-order compact θ-method for delay parabolic equations. Electronic Research Archive, 2024, 32(4): 2805-2823. doi: 10.3934/era.2024127 |
[7] | Peng Zhou, Teng Wang . The parameter-Newton iteration for the second-order cone linear complementarity problem. Electronic Research Archive, 2022, 30(4): 1454-1462. doi: 10.3934/era.2022076 |
[8] | Yang Song, Beiyan Yang, Jimin Wang . Stability analysis and security control of nonlinear singular semi-Markov jump systems. Electronic Research Archive, 2025, 33(1): 1-25. doi: 10.3934/era.2025001 |
[9] | Jun Pan, Yuelong Tang . Two-grid H1-Galerkin mixed finite elements combined with L1 scheme for nonlinear time fractional parabolic equations. Electronic Research Archive, 2023, 31(12): 7207-7223. doi: 10.3934/era.2023365 |
[10] | Haiyan Song, Fei Sun . A numerical method for parabolic complementarity problem. Electronic Research Archive, 2023, 31(2): 1048-1064. doi: 10.3934/era.2023052 |
The deterministic Degasperis-Procesi equation admits weak multi-shockpeakon solutions of the form
u(x,t)=n∑i=1mi(t)e−|x−xi(t)|−n∑i=1si(t)sgn(x−xi(t))e−|x−xi(t)|,
where sgn(x) denotes the signum function with sgn(0)=0, if and only if the time-dependent parameters xi(t) (positions), mi(t) (momenta) and si(t) (shock strengths) satisfy a system of 3n ordinary differential equations. We prove that a stochastic perturbation of the Degasperis-Procesi equation also has weak multi-shockpeakon solutions if and only if the positions, momenta and shock strengths obey a system of 3n stochastic differential equations.
Nonlinear problems are a hot topic in modern mathematics, attracting the attention of many researchers. Problems in fields such as economics, physics, chemistry, and fluid dynamics can be transformed into nonlinear equation systems. The solution of nonlinear equation systems is difficult to discover, and this problem can be solved using some numerical methods [1,2,3]. Chen et al. proposed the alternating direction implicit (ADI) compact differential scheme, which reduces the CPU time of two-dimensional problems [4]. Yang et al. proposed an efficient compact finite difference method [5,6,7] for solving nonlinear equations and an orthogonal Gauss collocation method (OGCM) [8] for solving systems of partial differential equations.
Therefore, in-depth research on numerical methods for solving nonlinear systems has important practical significance [9]. In solving nonlinear systems, we can obtain approximate solutions with the help of various iterative methods. Multi-point iterative methods [10,11,12,13] are an effective class of methods for solving nonlinear equation systems.
Wang et al. achieved research results in the study of multi-point iterative methods, such as local convergence analysis [14,15,16] and fractal theory-based investigations of iterative methods[17,18]. Cordero et al. studied the three-step iterative method [19,20] for solving nonlinear equations.
Newton's iterative method [21] is one of the most famous iterative methods
x(n+1)=x(n)−S′(x(n))−1S(x(n)),n∈0,1,2,⋯. | (1.1) |
where S:D⊆Rn⟶Rn and the inverse of the Jacobian matrix S′(xn) is S′(xn)−1.
On the basis of Newton's method, Potra-Pták [22] proposed a third-order iterative method
{yn=xn−s(xn)s′(xn),xn+1=xn−s(xn)+s(yn)s′(xn). | (1.2) |
These methods (1.2) have the advantage of not needing to compute second-order derivatives.
Cordero et al.[23] proposed an iterative method using division difference instead of the traditional Jacobi matrix
Sn+1=Sn−[m∑h=1AhJS(ηh(Sn))]−1F(Sn), | (1.3) |
where ηh(x)=x−τhJ−1S(x)S(x) and the parameters Ah satisfy ∑mh=1Ah=1. Considering τ1=0, the following iterative method is obtained:
Sn+1=Sn−[14JS(xn)+34JS(xn−23JS(xn)−1S(xn))]−1S(xn). | (1.4) |
Kou et al. [24] proposed an iterative method for solving nonlinear equations with convergence of at least third order
{yn=xn−s(xn)s′(xn),xn+1=xn−γs(xn)+s(yn)s′(xn)−(1−γ)s(xn)2s′(xn)(s(xn)−s(yn)), | (1.5) |
where γ∈R. Let S[xn,yn]=S(yn)−S(xn)yn−xn and S(xn)2S′(xn)(s(xn)−S(yn)) can be represented asS(xn)S[xn,yn]. We obtain an iterative method (1.6) for solving nonlinear systems. The iterative method (1.5) can be represented as
x(n+1)=x(n)−γS′(x(n))−1(S(x(n))+S(y(n)))−(1−γ)S[x(n),y(n)]−1S(x(n)), | (1.6) |
where
y(n)=x(n)−S′(x(n))−1S(x(n)), |
γ∈R is a real parameter.
The contents of this paper are as follows: Section 2 analyzes the order of convergence of the extended iterative method (1.6). In Section 3, we adopt the Möbius conjugacy map to convert the iterative method (1.6) to a rational operator and examine the fixed points and critical points for this rational operator in depth. The parameter space is constructed and some specific parameters are chosen in the parameter space. Section 4 compares the iterative method (1.6) with other iterative methods through numerical experiments. Section 5 shows some fractal diagrams. We verify that the expanded method has optimal stability and convergence for the parameter γ=0.
We prove the convergence order for the iterative method (1.6) with the following theorem.
Theorem 1. Let S:D⊆Rn→Rn be a sufficiently differentiable function from an open convex set D and let x∗∈D denote the solution of S(x)=0 that makes S′ continuous and nonsingular at x∗. The error equation of the iterative method (1.6) is
e(n+1)=c22(1+γ)(e(n))3+(−3c32(1+2γ)+c2c3(3+4γ))(e(n))4+O((e(n))5). | (2.1) |
When the parameter γ=−1, the iterative method is fourth-order convergent and the error expression is
e(n+1)=(3c32−c2c3)(e(n))4+O((e(n))5). | (2.2) |
Proof. Let C(n)=1n!S′(x∗)−1S(n)(x∗),i≥2 and e(n)=x(n)−x∗. Expanding S in terms of the Taylor's series on x∗, we get
S(x(n))=S′(x∗)(e(n)+C2(e(n))2+C3(e(n))3+C4(e(n))4+C5(e(n))5)+O((e(n))6), | (2.3) |
S′(x(n))=S′(x∗)(I+2C2e(n)+3C3(e(n))2+4C4(e(n))3+5C5(e(n))4)+O((e(n))5). | (2.4) |
According to the equation S′(x(n))−1S′(x(n))=I, we can calculate the inverse of S′(x(n)) :
S′(x(n))−1=S′(x(∗))−1(I+A1e(n)+A2(e(n))2+A3(e(n))3+A4(e(n))4)+O((e(n))5)), | (2.5) |
where
A1=−2C2,A2=4C22−3C3,A3=−8C32+6C2C3+6C3C2−4C4,A4=16C42+9C23+8C2C4+8C4C2−12C22C3−12C2C3C2−12C3C22−5C5. |
Now, using(2.3)–(2.5), we get
y(n)−x∗=x(n)−x∗−S′(x(n))−1S(x(n))=e(n)+C2(e(n))2+(2C3−2C22)(e(n))3+(3C4−4C2C3−3C2C3+4C32)(e(n))4+O((e(n))5). | (2.6) |
Taylor expansion of S(y(n)) at x∗ gives
S(y(n))=S(x∗)+S′(x∗)(y(n)−x∗)+S″(x∗)2!(y(n)−x∗)2+⋯+S(n)(x∗)n!(y(n)−x∗)n, | (2.7) |
and we have
S(y(n))=S′(x(n))(C2(e(n))2−2(C22−C3)(e(n))3+(5C32−7C2C3+3C4)(e(n))4)+O((e(n))5). | (2.8) |
In the iterative method (1.6), a first-order difference quotient operator is introduced, which can be regarded as a mapping relationship S[x,y]:D×D⊂Rn×Rn→L(Rn),
S[x,y]=∫10S′(x+ρd)dρ,∀(x,d)∈Rn×Rn. |
By performing Taylor expansion on S′(x+ρd) at point x, we can obtain the following results:
∫10S′(x+ρd)dρ=S′(x)+12S′′(x)d+16S′′′(x)d2+O(d3). | (2.9) |
According to
x+d=y,d=y−x=−S′(x(n))−1S(x(n)), | (2.10) |
the following conclusion can be reached:
S[x(n),y(n)]=S′(x∗)(I+U1e(n)+U2(e(n))2+U3(e(n))3+U4(e(n))4+O((e(n))5), | (2.11) |
where
U1=C2,U2=C22+C3,U3=−2C23+3C2C3+C4,U4=4C24−8C22C3+2C32+4C2C4+C5. | (2.12) |
Based on the equation S[x(n),y(n)]−1S[x(n),y(n)]=I, the inverse of S[x(n),y(n)]−1 is
S[x(n),y(n)]−1=S′(x∗)−1(I−U1e(n)+(U21−U2)(e(n))2+(−U31+2U1U2−U3)(e(n))3+O((e(n))4)). | (2.13) |
Therefore, one has
S[x(n),y(n)]−1=S′(x∗)−1(I−C2e(n)+(C22−C3)(e(n))2+(−2C23+2C2C3−C4)(e(n))3+O((e(n))4)). | (2.14) |
The error equations for the iterative method (1.6) give us
e(n+1)=C22(1+γ)(e(n))3+(−3C32(1+2γ)+C2C3(3+4γ))(e(n))4+O((e(n))5). | (2.15) |
The equation above shows that the iterative method converges at least to the third order. In particular, when the parameter γ=−1, the error equation is as follows:
e(n+1)=(3C32−C2C3)(e(n))4+O((e(n))5). | (2.16) |
When the parameter γ=1, the method (1.6) is a fourth-order convergent method with the following iteration format:
x(n+1)=x(n)+S′(x(n))−1(S(x(n))+S(y(n)))−2S[x(n),y(n)]−1S(x(n)), | (2.17) |
where
y(n)=x(n)−S′(x(n))−1S(x(n)). | (2.18) |
Fractal theory [25,26] can be used to study the stability and convergence of iterative methods. To visualize the convergence of the proposed iterative approach, fractal maps are constructed for a variety of nonlinear functions. Such methods have been employed in numerous research studies for practical applications. Lee et al.[27] studied a dynamic view of a class of single-parameter optimal eighth-order multi-root finders under Riemannian ball Möbius conjugate mappings using fractal knowledge. The complex dynamic behavior can be used for exploring the relevant properties of rational operators connected with iterative methods, which are the unique characteristics of rational operators closely connected with iterative methods, as demonstrated by the Möbius conjugate mappings on Riemann spheres. A thorough study of the behavior of the complex dynamics at fixed points allows us to perform a careful analysis of the stability of rational operators. A parameter space that was constructed by the critical point provides us with a reference to have a better idea of the stability in the iterative method, thus enabling us to pick more appropriate parameter values to make iterative methods more stable.
An in-depth study of the method's (1.6) complex dynamic behavior is presented in this section. First, we need to construct rational operators related to the iterative method (1.6).
Now, we will analyze the iterative method (1.5) in the context of quadratic polynomial dynamics. The rational operator can be constructed on an arbitrary nonlinear function by a principle based on Riemann ball dynamics [28] and the scalar theorem [29]. Therefore, a rational operator for quadratic polynomials is designed.
Theorem 2. Let E(x)=(x−m)(x−a) be some arbitrary polynomial of quadratic form, where m and a are its roots. The corresponding rational operator De(x;γ) of the family in (1.5) applied to E(x) is
De(x;γ)=x3(γ+(1+x)3)1+3x+3x2+(1+γ)x3. | (3.1) |
Proof. Applying the iterative method (1.5) in E(x), we obtain a rational function Be associated only with m, a, and γ. Next, we construct a rational function Be by applying the Möbius transform in Be with the
k(x)=x−mx−a. | (3.2) |
The properties of the conjugate mappings are as follows (see[30]):
k(∞)=1,k(a)=∞,k(m)=0. | (3.3) |
We can get
Be=x+(−1+γ)(m+a−2x)(m−x)(−a+x)m2+a2+m(a−3x)−3ax+3x2+γ(m−x)(a−x)(m2+3ma+a2−5mx−5ax+5x2)(m+a−2x)3. | (3.4) |
and
De(x;γ)=(k∘Be∘k−1)(x)=x3(γ+(1+x)3)1+3x+3x2+(1+γ)x3. | (3.5) |
According to the theorem above, we find that the expressive form of the rational operator De(x;γ) rests on the choice of the parameter value γ. Furthermore, it is known that Be(x) and De(x;γ) are conjugate to each other, and thus the expression for De(x;γ) contains only the parameter γ. The form of De(x,γ) is affected by the parameter γ, so it is only necessary to explore its relation to γ. When γ is −1 or 0, the expressions of De(x;γ) will be further simplified by factorizing De(x;γ) as follows:
De(x;−1)=x3((1+x)3−1)1+3x+3x2, | (3.6) |
and
De(x;0)=x3(1+x)31+3x+3x2+x3. | (3.7) |
In this section, the fixed points and the stability of the rational operator De(x;γ) are analyzed. The particular value of the parameter γ will determine the number of fixed points and their stability.
De(x)−x=x(x−1)λ(x)ϑ(x), | (3.8) |
where
λ(x)=1+4x−(−6+γ)x2+4x3+x4, | (3.9) |
and
ϑ(x)=1+3x+3x2+(1+γ)x3. | (3.10) |
From solving the equation De(x)−x=0, we know that the fixed points include x=0,x=∞,x=1, and the roots of the polynomial λ(x)=1+4x−(−6+γ)x2+4x3+x4.
Here, x=0 and x=∞ are the free points and fixed points of the parameter γ, respectively. Varying values of the parameter γ lead to different expressions for the iterative method (1.6). So, for the number of strange fixed points for different values of the parameter, we have Theorem 3.
Theorem 3. For λ(x) and ϑ(x), we have
● If γ=0, (1+x)3 is the common factor of the polynomial λ(x) and ϑ(x). In this case, the operator De(x) exhibits one unique strong fixed point, namely x=1.
● If γ=9, the operator De(x) has four strange fixed points, namely x=−4.79129,x=−0.208712, the complex number x=−0.5+0.866025i and x=0.5−0.866025i .
● If γ=1, in this case, the operator De(x) exhibits two strange fixed points namely x=−2.61803,x=−0.381966.
● If γ=16, in this case, the operator De(x) exhibits four strange fixed points namely x=−5.82843,x=1 (with multiplicity 2) and x=−0.171573.
● If γ satisfies (γ−1)(γ+1)(γ−(−1)13)(γ−(−1)23)≠0, the operator De(x) has nine strange fixed points: x=1 and the eight roots of the polynomial λ(x)=0. Then the operator De(x) will have five strange fixed points: in addition to x=1, there are four other roots of the polynomial λ(x)=0.
Proof. ● With λ(x)=0 and ϑ(x)=0, we can get
(1+x)(1+x+x2)=0. | (3.11) |
Suppose that x∈C is some root such that λ(x)=0 and ϑ(x)=0. Linking these two multinomial equations, we are in a position to remove the parameter γ, giving us the equation (1+x)(1+x+x2)=0. From (1+x)=0 or (1+x+x2)=0, we get x=−1 or (−1)13 or x=(−1)23. The common factorization of λ(x)=0 and ϑ(x)=0 is (1+x). As we substitute x=−1 for λ(x;−1) and ϑ(x;−1), we get the parameter γ=0. So for the operator De(x;0), there is now one anomaly in the fixed points.
● The remainder of the proof is consistent with the approach outlined earlier.
Theorem 3 allows us to conclude that the maximum number of fixed points is 5, while the minimum is 1. The stability of strange fixed points is discussed below. We first need to calculate the first-order derivatives of the iterative method's rational operator De(x;γ) with regard to the parameter γ. The first-order derivatives of De(x;γ) are computed as follows:
D′e(x;γ)=−3x2(1+x)2δ(x)ϱ(x)2, | (3.12) |
where
δ(x)=1+γ+4x+6x2+4x3+x4+γx4, | (3.13) |
ϱ(x)=1+3x+3x2+(1+γ)x3. | (3.14) |
The stability characteristics of the fixed points described above are discussed in detail next. With x=1, for the stability of any parameter value γ : ∀δ∈C∖{−8}, the following conclusions are obtained:
1) The point x is repulsive when |248+γ|>1;
2) The point x is attractive when |248+γ|<1;
3) The point x is parabolic when |248+γ|=1;
4) The point x is a superattracting point when |248+γ|≠0.
Proof. This follows from Eq (3.1):
D′e(x;γ)=−3x2(1+x)2δ(x)ϱ(x)2, | (3.15) |
where
δ(x)=1+γ+4x+6x2+4x3+x4+γx4, | (3.16) |
and
ϱ(x)=1+3x+3x2+(1+γ)x3. | (3.17) |
Substituting x=1 into (3.15), we acquire D′e(1;γ)=248+γ. It is easy to obtain |248+γ|=1⇔|24|=|8+γ|. Let γ=x+iy and (γ∈C). Then the following equation holds: (x+8)2+y2=242. Therefore, |D′e(1;γ)|>1⇔|γ+8|<24.
The strange fixed point is the root of x=1 and the polynomial λ(x)= 1+4x−(−6+γ)x2+4x3+x4. The roots of λ(x) are denoted as λi(x), for any i ranging from 1 to 4. For γ satisfying x(x−1)≠0, we have the following:
● If x=1, taking parameter values in the region [−35,20]×[−25,25], it is a point of attraction.
● λ1(x) is exclusive and independent of the value of the parameter γ.
● λ2(x), λ3(x), and λ4(x) act as attractors for the value of γ within the small region of the complex plane [−10,12.5]×[0.96,0.98], [−2.58,2.56]×[−10.96,12.98], [−1.58,1.56]×[−2.06,1.98].
Figure 1 displays the stability regions when x = 1. The repulsive region is depicted as a gray area, while the attractive region is shown as a gold area. When the selected parameter value is inside the disk, a behaves as a repelling point, while a becomes an attractor when the γ value is outside the disk. We are generally more attentive to values inside the disk. Figures 1 and 2 illustrate the stability characteristics of the strange fixed points for the operator De(x;γ).
On the basis of the definition of critical point, the critical point of the operator De(x;γ) is found through solving the equation D′e(x;γ)=0. Clearly, the critical points of the operator are 0 and 1, and these are closely related to the roots of the quadratic polynomial E(x)=(x−m)(x−a). Apart from this, the other critical points are considered to be defined as free critical points. The following theorem shows what happens when number of free critical points changes with various values of the parameter γ.
Theorem 4. If D′e(x;γ)=0, then x=−1 and δ(x)=0 correspond to a free critical point.
● The parameter γ takes the value of 0 when x=−1,δ(x)=(1+x)4, then the rational operator De(x;γ) has four free critical points, which are x=−1 (with multiplicity 4).
● The value of parameter γ is -8 when x=1, ς(x)=(1+x)4−8(1+x4), then the rational operator De(x;γ) has two free critical points, which are x=−0.714286−0.699854i and x=−0.714286+0.699854i.
● The value of the parameter γ is 1 when x=(−1)23, ς(x)=(1+x)4+(1+x4), at which point, the rational operator De(x;γ) has four free critical points, which are x=−1 (with multiplicity 2) and x=0 (with multiplicity 2).
● For any γ(γ+8)(γ−9)(γ−1) that is not equal to 0, we can learn that there are four free critical points.
cr1=−11+γ−τ−12√Ω−ξ−Δ−Γ;
cr2=−11+γ−τ+12√Ω−ξ−Δ−Γ;
cr3=−11+γ+τ−12√Ω−ξ−Δ+Γ;
cr4=−11+γ+τ+12√Ω−ξ−Δ+Γ;
where τ=√−γ+γ2√2√1+2γ+γ2;ξ=61+γ;Δ=2+2γ1+γ;Γ=√1+2γ+γ2(−64(1+γ)3+96(1+γ)2−321+γ)4√2√−γ+γ2;Ω=8(1+γ)2. In spite of the large number of free critical points, however, they are not completely independent of each other, but they are interconnected and together influence the dynamical behavior of the iterative method, as cr1=1cr2and cr3=1cr4.
Proof. Under the assumption that x∈C indicates certain values of δ(x)=0 if we ϱ(x)=0, and combine these two equations, by eliminating the parameter γ, we get (1+x)3(−1+x3)x3=0, which can be reduced to (1+x)3(−1+x3)=0. It follows that (1+x)3 and (−1+x3), which can be a common factorization of δ(x)=0 and ϱ(x)=0. As a result, substituting x=−1 for δ(x)=0 and ϱ(x)=0 gives the parameter γ=0, at which point δ(x)ϱ(x)=−2, and there are four critical points. The equation (−1+x3)=0 is solved to solve for x=1,x=−(−1)13 and x=(−1)23. The equation (−1+x3)=0 is then solved to solve for γ=1 and γ=−8. There are four critical points in D′e(x)=6x2(1+x)2(1+2x)2 when γ=1. When γ=−8, D′e(x)=−3x2(1+x)2(7+10x+7x2)(1+4x+7x2)2 has two critical points.
The dynamic plane H and the parameter plane Q are plotted and used to visualize the dynamical behavior associated with the iterative method (1.6). In the process of plotting the parameter plane, each γ-value is located in the same connected region of the parameter space, resulting in a series of subsets of methods (1.6) with similar dynamic behavior. Thus, more stable regions can be identified in the parameter plane, allowing us to extract more members with stability from the family of iterative methods (1.6).
Q={γ∈C : the free critical point cr whose orbit converges to κq∈˜C under the influence of the operator De(x;γ)}.
H={x∈C : for a given γ∈Q, the orbital of x(γ) tends towards a certain number κh∈˜C under the action of De(x;γ)}.
By Theorem 3, there are, at most, three free independent critical points that can be observed. Furthermore, we recognize that x=−1 is actually the mapping point of the fixed point x=1, and this critical point does not show a distinctive feature on the parameter plane. Thus, we can distinguish two unique parameter planes that each carry complementary information.
When analyzing the free critical point cr1 (or cr2) as the starting point for a series of iterative methods (1.6) associated with a particular complex value, we color the complex plane according to the convergence properties of the different methods. The corresponding point is labeled in red in the complex plane if the iterative method (1.6) converges to zero. If the convergence is to a non-zero value, then it will be marked in blue. For all other cases, points will be marked green. Figure 3 depicts the parameter plane for the iterative method with the initial point set to cr1. The figure is generated within the range of γ values [−40,20]×[−30,30] and [−5,5]×[−5,5], using a grid of 1000×1000 points and computed with 50 iterations for each point. Figure 3b provides a more detailed presentation of the situation in Figure 3(a). In Figures 3 and 4, it is observed that the red region represents the iterative method converging to 0, the blue region represents convergence to infinity, and the green region indicates convergence to 0. In this final parameter selection process, we should try to avoid using green parameter values and prioritize red or blue parameter values to ensure the stability of the iterative process.
The iterative method (1.6) has a dynamic plane influenced by the parameter γ, and the regions are marked with different colors. The region where the initial point converges to 0 we color orange, and the region that converges to infinity is marked blue. The region is green if the initial point converges to the fixed point x=1. Black regions indicate no convergence to any root if the initial point converges to black. Additionally, the iteration track is set to red. The dynamic planes are built on a 400×400 point grid, with each point undergoing an iterative process of up to 20 iterations. In Figure 5, the parameter planes that correspond to some specific parameter values are shown. When γ=0, the blue and orange regions are the only ones, indicating that the initial points converge only to 0 and infinity. This parameter value is therefore considered to be a desirable value to take. For γ=−8 and γ=1, some black regions appear in Figure 5, and the initial values of these black regions do not converge to any root, suggesting that these parameter values are not ideal choices.
Figure 6: Dynamic planes for special parameter values γ. We have selected the parameter values that satisfy the condition |8+γ|>24 and divided these values into two subsets; a group of smaller values, specifically γ=20,30, and a larger group of values including γ=1300,1340,−1500. A related graphical presentation can be found in Figure 6. We plot γ=−7000,10,000; refer to Figure 7.
The dynamic plane shows only orange and blue as stable parameter regions. The larger the absolute value of the parameter, the more complex the structure on the dynamic plane tends to be. Putting these observations into context, it is possible to conclude that the corresponding iterative method shows better stability when γ=0.
In this section, we conduct a series of numerical experiments aimed at verifying the theoretical convergence and stability predictions made in previous sections for the method (2.17)(XM). The iterative methods (2.17) (XM) are compared with Homeier's method (HM) [31] and Chun's method (CM)[32] for solving some nonlinear equations. Then two numerical experiments are carried out to solve the nonlinear systems.
Homeier's method (HM) [31]:
x(n+1)=x(n)−S(x(n))2(S′(x(n))−1+S′(y(n))−1), | (4.1) |
where y(n)=x(n)−S′(x(n))−1S(x(n)).
Chun's method [32]:
x(n+1)=x(n)−32S′(x(n))−1S(x(n))+12(S′(x(n))2)−1(S(x(n))S′(y(n))), | (4.2) |
where y(n)=x(n)−S′(x(n))−1S(x(n)).
Within Tables 1 and 2, the time column demonstrates the computing time, the x0 column labels the starting point in the iterative process, ACOC [33] denotes the approximate convergence order.
si(x) | Method | x0 | |x1−x0| | |x2−x1| | |x3−x2| | |x4−x3| | ACOC | Time |
s1(x) | XM | 1.4 | 3.476e-2 | 9.725e-6 | 2.2109e-16 | 2.5973e-48 | 2.9999996 | 0.390530 |
HM | 1.4 | 3.4769e-2 | 1.333e-6 | 7.8428e-20 | 1.4607e-59 | 3.000001 | 0.398433 | |
CM | 1.4 | 3.475e-2 | 2.0184e-5 | 4.2017e-15 | 3.7902e-44 | 2.9999984 | 0.409339 | |
XM | 1.3 | 6.5158e-2 | 2.0184e-5 | 4.2017e-15 | 3.7902e-44 | 2.9999984 | 0.409339 | |
HM | 1.3 | 6.5223e-2 | 6.91e-6 | 9.9906e-18 | 3.0194e-53 | 2.9999994 | 0.406029 | |
CM | 1.3 | 6.507e-2 | 1.6018e-4 | 2.1006e-12 | 4.7361e-36 | 3.000016 | 0.418919 | |
s2(z) | XM | 1.4 | 1.2322e-1 | 2.772e-3 | 2.3441e-8 | 1.4102e-23 | 3.0004456 | 0.378724 |
HM | 1.4 | 1.2576e-1 | 2.3771e-4 | 2.9752e-12 | 5.8369e-36 | 2.9999689 | 0.633551 | |
CM | 1.4 | 1.1826e-1 | 7.7388e-3 | 1.1455e-6 | 3.6246e-18 | 3.0030309 | 0.473145 | |
XM | 1.6 | 7.3619e-2 | 3.689e-4 | 6.3353e-11 | 2.7837e-31 | 2.9999538 | 0.364632 | |
HM | 1.6 | 7.3904e-2 | 1.0183e-4 | 2.341e-13 | 2.8435e-39 | 3.0000122 | 0.487294 | |
CM | 1.6 | 7.3235e-2 | 7.7125e-4 | 1.1033e-9 | 3.2379e-27 | 2.9998068 | 0.482667 | |
s3(z) | XM | 1.3 | 1.0363e-1 | 8.6367e-4 | 3.9673e-10 | 3.8196e-29 | 3.0001193 | 0.385904 |
HM | 1.3 | 1.0449e-1 | 6.1387e-6 | 1.0142e-17 | 4.5733e-53 | 2.9999987 | 0.889935 | |
CM | 1.3 | 1.0241e-1 | 2.0802e-3 | 1.1527e-8 | 1.9478e-24 | 3.0005625 | 0.780309 | |
XM | 1.5 | 9.5067e-2 | 4.414e-4 | 5.2757e-11 | 9.0156e-32 | 2.9999442 | 0.364020 | |
HM | 1.5 | 9.5456e-2 | 5.2515e-5 | 6.3513e-15 | 1.1233e-44 | 3.0000132 | 0.815432 | |
CM | 1.5 | 9.4686e-2 | 8.2275e-4 | 7.0641e-10 | 4.4833e-28 | 2.9998082 | 0.777321 | |
s4(z) | XM | 1.8 | 2.2924e-1 | 5.4852e-3 | 4.8052e-8 | 3.2011e-23 | 3.0007836 | 0.770896 |
HM | 1.8 | 2.3536e-1 | 6.3333e-4 | 1.0604e-12 | 5.0524e-39 | 2.9992499 | 0.919149 | |
CM | 1.8 | 2.2038e-1 | 1.4339e-2 | 1.7532e-6 | 3.087e-18 | 3.0041621 | 0.869809 | |
XM | 1.9 | 1.3384e-1 | 8.8166e-4 | 1.9802e-10 | 2.2404e-30 | 3.0000956 | 0.764242 | |
HM | 1.9 | 1.3478e-1 | 5.852e-5 | 8.4809e-16 | 2.5851e-48 | 2.9999442 | 0.899065 | |
CM | 1.9 | 13271e-1 | 2.0161e-3 | 4.7188e-9 | 6.0185e-26 | 3.0004016 | 0.958593 | |
s5(z) | XM | 0.7 | 1.0509e-1 | 1.0106e-4 | 7.3026e-14 | 2.755e-41 | 3.0000094 | 0.372409 |
HM | 0.7 | 1.0503e-1 | 1.5417e-4 | 4.4925e-13 | 1.1114e-38 | 3.0000055 | 0.387347 | |
CM | 0.7 | 1.0484e-1 | 3.4457e-4 | 1.0807e-11 | 3.3333e-34 | 3.0000218 | 0.711144 | |
XM | 0.9 | 3.0153e-1 | 3.6606e-3 | 3.4953e-9 | 3.0209e-27 | 3.0005195 | 0.363470 | |
HM | 0.9 | 3.0066e-1 | 4.5271e-3 | 1.1409e-8 | 1.8204e-25 | 3.0002463 | 0.411865 | |
CM | 0.9 | 2.9467e-1 | 1.0518e-2 | 3.1085e-7 | 7.9314e-21 | 3.0011099 | 0.726080 | |
s6(z) | XM | 0.2 | 9.8796e-2 | 2.0631e-4 | 1.9522e-12 | 1.6541e-36 | 2.999995 | 0.379134 |
HM | 0.2 | 9.8925e-2 | 7.689e-5 | 3.4635e-14 | 3.1654e-42 | 3.0000016 | 0.546499 | |
CM | 0.2 | 9.8549e-2 | 4.5356e-4 | 4.8574e-11 | 5.9694e-32 | 2.9999684 | 0.509358 | |
XM | 0.6 | 4.7715e-1 | 2.1853e-2 | 2.2984e-6 | 2.6993e-18 | 2.9989692 | 0.369832 | |
HM | 0.6 | 4.8842e-1 | 1.0583e-2 | 9.0743e-8 | 5.6928e-23 | 3.0004065 | 0.509111 | |
CM | 0.6 | 4.5978e-1 | 3.9191e-2 | 3.0088e-5 | 1.4186e-14 | 2.9942634 | 0.548055 |
si(x) | γ | x0 | |x1−x0| | |x2−x1| | |x3−x2| | |x4−x3| | ACOC | Time |
s1(x) | 0 | 0.9 | 4.2417e-1 | 4.1039e-2 | 1.7406e-5 | 1.2675e-15 | 3.0060079 | 0.395263 |
1 | 0.9 | 3.2788e-1 | 1.3573e-1 | 1.6133e-3 | 2.0244e-9 | 3.0657447 | 0.411651 | |
9 | 0.9 | 4.4244e-1 | 54.806 | 5.2807 | 4.7567 | nc | 0.463676 | |
s2(x) | 0 | 1.5 | 2.5974e-2 | 2.0187e-5 | 9.006e-15 | 7.9969e-43 | 3.0000018 | 0.365875 |
1 | 1.5 | 2.5952e-2 | 4.1889e-5 | 1.6095e-13 | 9.1294e-39 | 3.000007 | 0.406378 | |
9 | 1.5 | 2.5778e-2 | 2.1550e-4 | 1.0967e-10 | 1.4441e-29 | 3.0000642 | 0.524439 | |
s3(x) | 0 | 1.4 | 4.4916e-3 | 5.6144e-8 | 1.0866e-22 | 7.8771e-67 | 3.0 | 0.384539 |
1 | 1.4 | 4.4915-3 | 1.1286e-7 | 1.7652e-21 | 6.7536e-63 | 3.0 | 0.413546 | |
9 | 1.4 | 4.4911e-3 | 5.6657e-7 | 1.1166e-18 | 8.5487e-54 | 3.0000001 | 0.703668 | |
s4(x) | 0 | 2.0 | 3.4712e-2 | 1.2797e-5 | 6.0467e-16 | 6.3787e-47 | 3.0000009 | 0.796057 |
1 | 2.0 | 3.4699e-2 | 2.6357e-5 | 1.0566e-14 | 6.8067e-43 | 3.000003 | 0.813411 | |
9 | 2.0 | 3.459e-2 | 1.3483e-4 | 7.0755e-12 | 1.022e-33 | 3.0000251 | 0.801792 | |
s5(x) | 0 | 0.6 | 5.189e-3 | 9.9857e-9 | 7.044e-26 | 2.4725e-77 | 3.0 | 0.794878 |
1 | 0.6 | 5.189e-3 | 1.9977e-8 | 1.1281e-24 | 2.0311e-73 | 3.0 | 1.063982 | |
9 | 0.6 | 5.1889e-3 | 9.9912e-8 | 7.0557e-22 | 1.022e-33 | 3.0000251 | 0.801792 | |
s6(x) | 0 | 1.6 | 1.1355 | 3.5448e-1 | 9.0195e-3 | 1.6249e-7 | 2.975614 | 0.407156 |
1 | 1.6 | 1.0323 | 4.3713e-1 | 2.9580e-2 | 1.119e-5 | 2.925868 | 0.542887 | |
9 | 1.6 | 0.20652 | 2.785e-1 | 3.5908e-1 | 3.8939e-1 | nc | 0.554967 |
We use six nonlinear equations mentioned in references [34,35] for the in-depth study.
s1(x)=x3+4x2−10,ε≈1.3652300134140968.s2(x)=x4−log(x)−5,ε≈1.5259939537536892.s3(x)=sin2(x)−x2+1,ε≈1.4044916482153412.s4(x)=cos(πx2+x2−π),ε≈0.59480888843087377−0.0000852783334727i.s5(x)=arcsin(x2−1)−12x+1,ε≈0.59480888843087377. |
Additionally, we explore population growth patterns[36]: population dynamics are analyzed by means of a first-order linear ordinary differential equation as follows:
P′(x)=cP(x)+s, | (4.3) |
where c stands for a constant birth rate and s stands for a constant migration rate. P(x) is used to represent the total population size at any random time point x. By solving the linear differential equation, we can get its generalized solution
P(x)=P0ecx+sc(ecx−1), | (4.4) |
By P0, we denote the initial population size. Examining and exploring the multiple possible values of the parameters and combining them with the initial conditions provided in the literature [37], we can accurately calculate the birth rate by using the nonlinear equation below:
s6(x)=1564−1000ex−435x(ex−1),ε≈0.1009979328831604.
Here, c=k represents the ideal birth rate we wish to determine.
Example 1. Take the nonlinear system
{x1+ex1−cos(x2)=0,3x1−x2−sin(x2)=0, | (4.5) |
where x1 and x2 are the unknowns in the nonlinear system (4.5).
Example 2. Take the nonlinear the system of equations for the Hammerstein equation, which plays a crucial role in nonlinear integral equations:
s(x)=1+(15)∫10M(x,t)s(t)3dt, | (4.6) |
where x∈C[0,1],x,t∈[0,1]. The kernel function M is
M(x,t)={(1−x)tt≤x,x(1−t)x≤t. | (4.7) |
A system of nonlinear equations is formed by discretization and solved (4.6) by approximating the integral in Eq (4.6) by the Gauss-Legendre product method.
∫10x(t)dt≈7∑i=1wjx(tj), | (4.8) |
Here, tj and wj are the nodes in the Gauss-Legende product method and the weights. Denote the approximation of s(tj) as xi,i=1,…,7
si−1−157∑j=1aijs3j=0,i=1,…,7 | (4.9) |
where
ai,j={wjtj(1−ti)j≤i,wjti(1−tj)i<j. | (4.10) |
According to Table 1, the iterative method (1.6) has high computational accuracy. Compared with other methods, the iterative method (1.6) takes less time for calculation. In Table 2, different values of stabilization parameters such as γ=0,γ=1 and γ=9, are chosen, and the iterative method (1.6) has better convergence. Tables 3 and 4 are numerical experiments for solving nonlinear systems.
γ | iter | ‖x(k+1)−x(k)‖ | ‖S(x(k+1))‖ | ACOC |
-1 | 6 | 7.188e-505 | 1.457e-2017 | 4.00000 |
0 | 7 | 1.126e-440 | 5.584e-1321 | 3.00000 |
1 | 7 | 8.525e-362 | 4.847e-1084 | 3.00000 |
9 | 8 | 1.170e-450 | 6.260e-1350 | 3.00000 |
γ | iter | ‖x(k+1)−x(k)‖ | ‖S(x(k+1))‖ | ACOC |
-1 | 5 | 2.815e-312 | 8.479e-1250 | 4.00017 |
0 | 6 | 4.212e-367 | 2.286e-1102 | 3.00001 |
1 | 6 | 6.901e-327 | 2.011e-981 | 3.00001 |
9 | 6 | 1.800e-239 | 1.786e-718 | 3.00001 |
This section draws fractal graphs for different iterative methods used to solve the nonlinear equation S1(x)=x2−1. The iterative methods we compare include AM (when parameter γ=0, iterative method (1.6) is defined as AM), K2 (when parameter γ=−1, iterative method (1.6) is defined as K2), Solaiman's method (MH), and Kou's method (KW).
Solaiman's method [38] (MH):
x(n+1)=x(n)−S′(x(n))−1S(x(n))+2(4(S′(x(n)))4−4S(x(n))(S′(x(n)))2S″(x(n))+(S(x(n)))2(S″(x(n)))2)−1S(x(n))S′(x(n))S″(x(n)), | (5.1) |
Kou's method [39] (KW) :
x(n+1)=x(n)−S′(x(n))−1(S(y(n))−S(x(n))), | (5.2) |
where y(n)=x(n)+S′(x(n))−1S(x(n)).
We created a 500×500 point grid to generate images in the region D=[−5,5]×[−5,5]. The maximum number is 25 iterations. Blue, red, and yellow areas respectively, converge to the roots of the polynomials, whereas areas that do not converge within the maximum number of iterations are in black. The number of iterations is indicated by the brightness of the colors. For example, regions with fewer iterations show brighter colors. From the fractal maps, method AM showed better convergence than the other methods. Therefore, the iterative method (AM) has better convergence and stability than the other methods.
The paper investigates the complex dynamic behavior exhibited by Kou's extended iterative method, which is used to solve systems for nonlinear equations. Strange fixed points and critical points of the iterative method (1.6) are intensively studied. We found the relevant parameter plane and select some parameter values for further analysis. Some unique parameter planes with values such as γ=0,−1 and 9 are observed. By observing the parameter plane and the dynamic plane, we discover the parameter γ=0 to have good stability. Numerical experiments and dynamic analyses both agree with each other, verifying the accuracy of the conclusion. We will study the iterative methods with higher convergence and their stability in the future.
The authors declare they have not used artificial intelligence (AI) tools in the creation of this article.
This research was supported by Educational Commission Foundation of Liaoning Province of China (No. LJ212410167008), the Key Project of Bohai University (No. 0522xn078), and the "14th Five-Year Plan" Project of Liaoning Province Education Science (No. JG24EB014).
The authors declare there are no conflicts of interest.
[1] |
L. K. Arruda, N. V. Chemetov, F. Cipriano, Solvability of the Stochastic Degasperis-Procesi Equation, J. Dyn. Differ. Equ., (2021), 1–20. https://doi.org/10.1007/s10884-021-10021-5 doi: 10.1007/s10884-021-10021-5
![]() |
[2] |
S. Hakkaev, I. D. Iliev, K. Kirchev, Stability of periodic travelling shallow-water waves determined by Newton' s equation, J. Phys. A: Math. Theor., 41 (2008), 085203. https://doi.org/10.1088/1751-8113/41/8/085203 doi: 10.1088/1751-8113/41/8/085203
![]() |
[3] |
R. Camassa, D. D. Holm, An integrable shallow water equation with peaked solitons, Phys. Rev. Lett., 71 (1993), 1661–1664. https://doi.org/10.1103/PhysRevLett.71.1661 doi: 10.1103/PhysRevLett.71.1661
![]() |
[4] |
A. Fokas, B. Fuchssteiner, Symplectic structures, their Backlund transformation and hereditary symmetries, Physica D, 4 (1982), 47–66. https://doi.org/10.1016/0362-546X(81)90025-0 doi: 10.1016/0362-546X(81)90025-0
![]() |
[5] |
J. Lenells, Classification of all travelling-wave solutions for some nonlinear dispersive equations, Phil. Trans. R. Soc. A, 365 (2007), 2291–2298. https://doi.org/10.1098/rsta.2007.2009 doi: 10.1098/rsta.2007.2009
![]() |
[6] |
T. M. Bendall, C. J. Cotter, D. D. Holm, Perspectives on the formation of peakons in the stochastic Camassa-Holm equation, Proc. R. Soc. A., 477 (2021), 20210224. https://doi.org/10.1098/rspa.2021.0224 doi: 10.1098/rspa.2021.0224
![]() |
[7] |
Y. Chen, H. Gao, B. Guo, Well-posedness for stochastic Camassa-Holm equation, J. Differ. Equ., 253 (2012), 2353–2379. https://doi.org/10.1016/j.jde.2012.06.023 doi: 10.1016/j.jde.2012.06.023
![]() |
[8] |
D. Crisan, D. D. Holm, Wave breaking for the stochastic Camassa-Holm equation, Physica D, 376 (2018), 138–143. https://doi.org/10.1016/j.physd.2018.02.004 doi: 10.1016/j.physd.2018.02.004
![]() |
[9] |
D. D. Holm, T. M. Tyranowski, Variational principles for stochastic soliton dynamics, Proc. Math. Phys. Eng. Sci., 472 (2016), 20150827. https://doi.org/10.1098/rspa.2015.0827 doi: 10.1098/rspa.2015.0827
![]() |
[10] |
G. Lv, J. Wei, G. Zou, The dependence on initial data of stochastic Camassa-Holm equation, Appl. Math. Lett., (2020), 106472. https://doi.org/10.1016/j.aml.2020.106472 doi: 10.1016/j.aml.2020.106472
![]() |
[11] |
X. Deng, A note on exact travelling wave solutions for the modified Camassa-Holm and Degasperis-Procesi equations, Appl. Math. Comput., 218 (2009), 2269–2276. https://doi.org/10.1016/j.amc.2011.07.044 doi: 10.1016/j.amc.2011.07.044
![]() |
[12] | A. Darós, L. K. Arruda Saraiva de Paiva, The Modified Camassa-Holm Equation: Wave breaking, Classification of Traveling Waves and Explicit Elliptic Peakons, arXiv preprint, (2018), arXiv: 1809.02529. |
[13] |
H. R. Dullin, G. A. Gottwald, D. D. Holm, An integrable shallow water equation with linear and nonlinear dispersion, Phys. Rev. Lett., 87 (2001), 194501. https://doi.org/10.1103/PhysRevLett.87.194501 doi: 10.1103/PhysRevLett.87.194501
![]() |
[14] |
C. Rohde, H. Tang, On the stochastic Dullin-Gottwald-Holm equation: Global existence and wave-breaking phenomena, Nonlinear Differ. Equ. Appl., 28 (2021), 1–34. https://doi.org/10.1007/s00030-021-00676-w doi: 10.1007/s00030-021-00676-w
![]() |
[15] | A. Degasperis, M. Procesi, Asymptotic integrability, in Symmetry and Perturbation Theory, (Rome, 1998), World Sci. Publishing, River Edge, NJ, (1999), 23–37. |
[16] |
G. M. Coclite, K. Karlsen, N. Risebro, Numerical schemes for computing discontinuous solutions of the Degasperis-Procesi equation, IMA J. Numer. Anal., 28 (2008), 80–105. https://doi.org/10.1093/imanum/drm003 doi: 10.1093/imanum/drm003
![]() |
[17] |
Z. Yin, On the Cauchy problem for an integrable equation with peakon solutions, Illinois J. Math., 47 (2003), 649–666. https://doi.org/10.1215/ijm/1258138186 doi: 10.1215/ijm/1258138186
![]() |
[18] |
H. Lundmark, Formation and dynamics of shock waves in the Degasperis-Procesi equation, J. Nonlinear Sci., 17 (2007), 169–198. https://doi.org/10.1007/s00332-006-0803-3 doi: 10.1007/s00332-006-0803-3
![]() |
[19] |
H. Lundmark, J. Szmigeilski, Multi-peakon solutions of the Degasperis-Procesi equation, Inverse Probl., 19 (2003), 1241–1245. https://doi.org/10.1088/0266-5611/19/6/001 doi: 10.1088/0266-5611/19/6/001
![]() |
[20] |
H. Lundmark, J. Szmigeilski, Degasperis-Procesi peakons and the discrete cubic string, Int. Math. Res. Papers, 2 (2005), 53–116. https://doi.org/10.1155/IMRP.2005.53 doi: 10.1155/IMRP.2005.53
![]() |
[21] |
A. Degasperis, D. D. Holm, A. N. W. Hone, A New Integrable Equation with Peakon Solutions, Theor. Math. Phys., 133 (2002), 1463–1474. https://doi.org/10.1023/A:1021186408422 doi: 10.1023/A:1021186408422
![]() |
[22] |
G. M. Coclite, K. H. Karlsen, On the well-posedness of the Degasperis-Procesi equation, J. Funct. Anal., 233 (2006), 60–91. https://doi.org/10.1016/j.jfa.2005.07.008 doi: 10.1016/j.jfa.2005.07.008
![]() |
[23] |
J. Escher, Y. Liu, Z. Yin, Shock waves and blow-up phenomena for the periodic Degasperis - Procesi equation, Indiana Univ. Math. J., 56 (2007), 87–117. https://doi.org/10.1512/iumj.2007.56.3040 doi: 10.1512/iumj.2007.56.3040
![]() |
[24] |
J. Szmigeilski, L. Zhou, Peakon-antipeakon interactions in the Degasperis-Procesi equation, in algebraic and geometric aspects of integrable systems and random matrices, Contemp. Math., 593 (2013), 83–107. https://doi.org/10.1090/conm/593/11873 doi: 10.1090/conm/593/11873
![]() |
[25] |
Y. Chen, H. Gao, Global existence for the stochastic Degasperis-Procesi equation, Discrete Contin. Dyn. Syst., 35 (2015), 5171–5184. https://doi.org/10.3934/dcds.2015.35.5171 doi: 10.3934/dcds.2015.35.5171
![]() |
[26] | B. Perthame, Kinetic formulation of conservation laws, Oxford University Press, 2002. |
[27] | H. A. Hoel, A numerical scheme using multi-shockpeakons to compute solutions of the Degasperis-Procesi equation, Electron. J. Differ. Equ., 2007 (2007), 1–22. |
si(x) | Method | x0 | |x1−x0| | |x2−x1| | |x3−x2| | |x4−x3| | ACOC | Time |
s1(x) | XM | 1.4 | 3.476e-2 | 9.725e-6 | 2.2109e-16 | 2.5973e-48 | 2.9999996 | 0.390530 |
HM | 1.4 | 3.4769e-2 | 1.333e-6 | 7.8428e-20 | 1.4607e-59 | 3.000001 | 0.398433 | |
CM | 1.4 | 3.475e-2 | 2.0184e-5 | 4.2017e-15 | 3.7902e-44 | 2.9999984 | 0.409339 | |
XM | 1.3 | 6.5158e-2 | 2.0184e-5 | 4.2017e-15 | 3.7902e-44 | 2.9999984 | 0.409339 | |
HM | 1.3 | 6.5223e-2 | 6.91e-6 | 9.9906e-18 | 3.0194e-53 | 2.9999994 | 0.406029 | |
CM | 1.3 | 6.507e-2 | 1.6018e-4 | 2.1006e-12 | 4.7361e-36 | 3.000016 | 0.418919 | |
s2(z) | XM | 1.4 | 1.2322e-1 | 2.772e-3 | 2.3441e-8 | 1.4102e-23 | 3.0004456 | 0.378724 |
HM | 1.4 | 1.2576e-1 | 2.3771e-4 | 2.9752e-12 | 5.8369e-36 | 2.9999689 | 0.633551 | |
CM | 1.4 | 1.1826e-1 | 7.7388e-3 | 1.1455e-6 | 3.6246e-18 | 3.0030309 | 0.473145 | |
XM | 1.6 | 7.3619e-2 | 3.689e-4 | 6.3353e-11 | 2.7837e-31 | 2.9999538 | 0.364632 | |
HM | 1.6 | 7.3904e-2 | 1.0183e-4 | 2.341e-13 | 2.8435e-39 | 3.0000122 | 0.487294 | |
CM | 1.6 | 7.3235e-2 | 7.7125e-4 | 1.1033e-9 | 3.2379e-27 | 2.9998068 | 0.482667 | |
s3(z) | XM | 1.3 | 1.0363e-1 | 8.6367e-4 | 3.9673e-10 | 3.8196e-29 | 3.0001193 | 0.385904 |
HM | 1.3 | 1.0449e-1 | 6.1387e-6 | 1.0142e-17 | 4.5733e-53 | 2.9999987 | 0.889935 | |
CM | 1.3 | 1.0241e-1 | 2.0802e-3 | 1.1527e-8 | 1.9478e-24 | 3.0005625 | 0.780309 | |
XM | 1.5 | 9.5067e-2 | 4.414e-4 | 5.2757e-11 | 9.0156e-32 | 2.9999442 | 0.364020 | |
HM | 1.5 | 9.5456e-2 | 5.2515e-5 | 6.3513e-15 | 1.1233e-44 | 3.0000132 | 0.815432 | |
CM | 1.5 | 9.4686e-2 | 8.2275e-4 | 7.0641e-10 | 4.4833e-28 | 2.9998082 | 0.777321 | |
s4(z) | XM | 1.8 | 2.2924e-1 | 5.4852e-3 | 4.8052e-8 | 3.2011e-23 | 3.0007836 | 0.770896 |
HM | 1.8 | 2.3536e-1 | 6.3333e-4 | 1.0604e-12 | 5.0524e-39 | 2.9992499 | 0.919149 | |
CM | 1.8 | 2.2038e-1 | 1.4339e-2 | 1.7532e-6 | 3.087e-18 | 3.0041621 | 0.869809 | |
XM | 1.9 | 1.3384e-1 | 8.8166e-4 | 1.9802e-10 | 2.2404e-30 | 3.0000956 | 0.764242 | |
HM | 1.9 | 1.3478e-1 | 5.852e-5 | 8.4809e-16 | 2.5851e-48 | 2.9999442 | 0.899065 | |
CM | 1.9 | 13271e-1 | 2.0161e-3 | 4.7188e-9 | 6.0185e-26 | 3.0004016 | 0.958593 | |
s5(z) | XM | 0.7 | 1.0509e-1 | 1.0106e-4 | 7.3026e-14 | 2.755e-41 | 3.0000094 | 0.372409 |
HM | 0.7 | 1.0503e-1 | 1.5417e-4 | 4.4925e-13 | 1.1114e-38 | 3.0000055 | 0.387347 | |
CM | 0.7 | 1.0484e-1 | 3.4457e-4 | 1.0807e-11 | 3.3333e-34 | 3.0000218 | 0.711144 | |
XM | 0.9 | 3.0153e-1 | 3.6606e-3 | 3.4953e-9 | 3.0209e-27 | 3.0005195 | 0.363470 | |
HM | 0.9 | 3.0066e-1 | 4.5271e-3 | 1.1409e-8 | 1.8204e-25 | 3.0002463 | 0.411865 | |
CM | 0.9 | 2.9467e-1 | 1.0518e-2 | 3.1085e-7 | 7.9314e-21 | 3.0011099 | 0.726080 | |
s6(z) | XM | 0.2 | 9.8796e-2 | 2.0631e-4 | 1.9522e-12 | 1.6541e-36 | 2.999995 | 0.379134 |
HM | 0.2 | 9.8925e-2 | 7.689e-5 | 3.4635e-14 | 3.1654e-42 | 3.0000016 | 0.546499 | |
CM | 0.2 | 9.8549e-2 | 4.5356e-4 | 4.8574e-11 | 5.9694e-32 | 2.9999684 | 0.509358 | |
XM | 0.6 | 4.7715e-1 | 2.1853e-2 | 2.2984e-6 | 2.6993e-18 | 2.9989692 | 0.369832 | |
HM | 0.6 | 4.8842e-1 | 1.0583e-2 | 9.0743e-8 | 5.6928e-23 | 3.0004065 | 0.509111 | |
CM | 0.6 | 4.5978e-1 | 3.9191e-2 | 3.0088e-5 | 1.4186e-14 | 2.9942634 | 0.548055 |
si(x) | γ | x0 | |x1−x0| | |x2−x1| | |x3−x2| | |x4−x3| | ACOC | Time |
s1(x) | 0 | 0.9 | 4.2417e-1 | 4.1039e-2 | 1.7406e-5 | 1.2675e-15 | 3.0060079 | 0.395263 |
1 | 0.9 | 3.2788e-1 | 1.3573e-1 | 1.6133e-3 | 2.0244e-9 | 3.0657447 | 0.411651 | |
9 | 0.9 | 4.4244e-1 | 54.806 | 5.2807 | 4.7567 | nc | 0.463676 | |
s2(x) | 0 | 1.5 | 2.5974e-2 | 2.0187e-5 | 9.006e-15 | 7.9969e-43 | 3.0000018 | 0.365875 |
1 | 1.5 | 2.5952e-2 | 4.1889e-5 | 1.6095e-13 | 9.1294e-39 | 3.000007 | 0.406378 | |
9 | 1.5 | 2.5778e-2 | 2.1550e-4 | 1.0967e-10 | 1.4441e-29 | 3.0000642 | 0.524439 | |
s3(x) | 0 | 1.4 | 4.4916e-3 | 5.6144e-8 | 1.0866e-22 | 7.8771e-67 | 3.0 | 0.384539 |
1 | 1.4 | 4.4915-3 | 1.1286e-7 | 1.7652e-21 | 6.7536e-63 | 3.0 | 0.413546 | |
9 | 1.4 | 4.4911e-3 | 5.6657e-7 | 1.1166e-18 | 8.5487e-54 | 3.0000001 | 0.703668 | |
s4(x) | 0 | 2.0 | 3.4712e-2 | 1.2797e-5 | 6.0467e-16 | 6.3787e-47 | 3.0000009 | 0.796057 |
1 | 2.0 | 3.4699e-2 | 2.6357e-5 | 1.0566e-14 | 6.8067e-43 | 3.000003 | 0.813411 | |
9 | 2.0 | 3.459e-2 | 1.3483e-4 | 7.0755e-12 | 1.022e-33 | 3.0000251 | 0.801792 | |
s5(x) | 0 | 0.6 | 5.189e-3 | 9.9857e-9 | 7.044e-26 | 2.4725e-77 | 3.0 | 0.794878 |
1 | 0.6 | 5.189e-3 | 1.9977e-8 | 1.1281e-24 | 2.0311e-73 | 3.0 | 1.063982 | |
9 | 0.6 | 5.1889e-3 | 9.9912e-8 | 7.0557e-22 | 1.022e-33 | 3.0000251 | 0.801792 | |
s6(x) | 0 | 1.6 | 1.1355 | 3.5448e-1 | 9.0195e-3 | 1.6249e-7 | 2.975614 | 0.407156 |
1 | 1.6 | 1.0323 | 4.3713e-1 | 2.9580e-2 | 1.119e-5 | 2.925868 | 0.542887 | |
9 | 1.6 | 0.20652 | 2.785e-1 | 3.5908e-1 | 3.8939e-1 | nc | 0.554967 |
γ | iter | ‖x(k+1)−x(k)‖ | ‖S(x(k+1))‖ | ACOC |
-1 | 6 | 7.188e-505 | 1.457e-2017 | 4.00000 |
0 | 7 | 1.126e-440 | 5.584e-1321 | 3.00000 |
1 | 7 | 8.525e-362 | 4.847e-1084 | 3.00000 |
9 | 8 | 1.170e-450 | 6.260e-1350 | 3.00000 |
γ | iter | ‖x(k+1)−x(k)‖ | ‖S(x(k+1))‖ | ACOC |
-1 | 5 | 2.815e-312 | 8.479e-1250 | 4.00017 |
0 | 6 | 4.212e-367 | 2.286e-1102 | 3.00001 |
1 | 6 | 6.901e-327 | 2.011e-981 | 3.00001 |
9 | 6 | 1.800e-239 | 1.786e-718 | 3.00001 |
si(x) | Method | x0 | |x1−x0| | |x2−x1| | |x3−x2| | |x4−x3| | ACOC | Time |
s1(x) | XM | 1.4 | 3.476e-2 | 9.725e-6 | 2.2109e-16 | 2.5973e-48 | 2.9999996 | 0.390530 |
HM | 1.4 | 3.4769e-2 | 1.333e-6 | 7.8428e-20 | 1.4607e-59 | 3.000001 | 0.398433 | |
CM | 1.4 | 3.475e-2 | 2.0184e-5 | 4.2017e-15 | 3.7902e-44 | 2.9999984 | 0.409339 | |
XM | 1.3 | 6.5158e-2 | 2.0184e-5 | 4.2017e-15 | 3.7902e-44 | 2.9999984 | 0.409339 | |
HM | 1.3 | 6.5223e-2 | 6.91e-6 | 9.9906e-18 | 3.0194e-53 | 2.9999994 | 0.406029 | |
CM | 1.3 | 6.507e-2 | 1.6018e-4 | 2.1006e-12 | 4.7361e-36 | 3.000016 | 0.418919 | |
s2(z) | XM | 1.4 | 1.2322e-1 | 2.772e-3 | 2.3441e-8 | 1.4102e-23 | 3.0004456 | 0.378724 |
HM | 1.4 | 1.2576e-1 | 2.3771e-4 | 2.9752e-12 | 5.8369e-36 | 2.9999689 | 0.633551 | |
CM | 1.4 | 1.1826e-1 | 7.7388e-3 | 1.1455e-6 | 3.6246e-18 | 3.0030309 | 0.473145 | |
XM | 1.6 | 7.3619e-2 | 3.689e-4 | 6.3353e-11 | 2.7837e-31 | 2.9999538 | 0.364632 | |
HM | 1.6 | 7.3904e-2 | 1.0183e-4 | 2.341e-13 | 2.8435e-39 | 3.0000122 | 0.487294 | |
CM | 1.6 | 7.3235e-2 | 7.7125e-4 | 1.1033e-9 | 3.2379e-27 | 2.9998068 | 0.482667 | |
s3(z) | XM | 1.3 | 1.0363e-1 | 8.6367e-4 | 3.9673e-10 | 3.8196e-29 | 3.0001193 | 0.385904 |
HM | 1.3 | 1.0449e-1 | 6.1387e-6 | 1.0142e-17 | 4.5733e-53 | 2.9999987 | 0.889935 | |
CM | 1.3 | 1.0241e-1 | 2.0802e-3 | 1.1527e-8 | 1.9478e-24 | 3.0005625 | 0.780309 | |
XM | 1.5 | 9.5067e-2 | 4.414e-4 | 5.2757e-11 | 9.0156e-32 | 2.9999442 | 0.364020 | |
HM | 1.5 | 9.5456e-2 | 5.2515e-5 | 6.3513e-15 | 1.1233e-44 | 3.0000132 | 0.815432 | |
CM | 1.5 | 9.4686e-2 | 8.2275e-4 | 7.0641e-10 | 4.4833e-28 | 2.9998082 | 0.777321 | |
s4(z) | XM | 1.8 | 2.2924e-1 | 5.4852e-3 | 4.8052e-8 | 3.2011e-23 | 3.0007836 | 0.770896 |
HM | 1.8 | 2.3536e-1 | 6.3333e-4 | 1.0604e-12 | 5.0524e-39 | 2.9992499 | 0.919149 | |
CM | 1.8 | 2.2038e-1 | 1.4339e-2 | 1.7532e-6 | 3.087e-18 | 3.0041621 | 0.869809 | |
XM | 1.9 | 1.3384e-1 | 8.8166e-4 | 1.9802e-10 | 2.2404e-30 | 3.0000956 | 0.764242 | |
HM | 1.9 | 1.3478e-1 | 5.852e-5 | 8.4809e-16 | 2.5851e-48 | 2.9999442 | 0.899065 | |
CM | 1.9 | 13271e-1 | 2.0161e-3 | 4.7188e-9 | 6.0185e-26 | 3.0004016 | 0.958593 | |
s5(z) | XM | 0.7 | 1.0509e-1 | 1.0106e-4 | 7.3026e-14 | 2.755e-41 | 3.0000094 | 0.372409 |
HM | 0.7 | 1.0503e-1 | 1.5417e-4 | 4.4925e-13 | 1.1114e-38 | 3.0000055 | 0.387347 | |
CM | 0.7 | 1.0484e-1 | 3.4457e-4 | 1.0807e-11 | 3.3333e-34 | 3.0000218 | 0.711144 | |
XM | 0.9 | 3.0153e-1 | 3.6606e-3 | 3.4953e-9 | 3.0209e-27 | 3.0005195 | 0.363470 | |
HM | 0.9 | 3.0066e-1 | 4.5271e-3 | 1.1409e-8 | 1.8204e-25 | 3.0002463 | 0.411865 | |
CM | 0.9 | 2.9467e-1 | 1.0518e-2 | 3.1085e-7 | 7.9314e-21 | 3.0011099 | 0.726080 | |
s6(z) | XM | 0.2 | 9.8796e-2 | 2.0631e-4 | 1.9522e-12 | 1.6541e-36 | 2.999995 | 0.379134 |
HM | 0.2 | 9.8925e-2 | 7.689e-5 | 3.4635e-14 | 3.1654e-42 | 3.0000016 | 0.546499 | |
CM | 0.2 | 9.8549e-2 | 4.5356e-4 | 4.8574e-11 | 5.9694e-32 | 2.9999684 | 0.509358 | |
XM | 0.6 | 4.7715e-1 | 2.1853e-2 | 2.2984e-6 | 2.6993e-18 | 2.9989692 | 0.369832 | |
HM | 0.6 | 4.8842e-1 | 1.0583e-2 | 9.0743e-8 | 5.6928e-23 | 3.0004065 | 0.509111 | |
CM | 0.6 | 4.5978e-1 | 3.9191e-2 | 3.0088e-5 | 1.4186e-14 | 2.9942634 | 0.548055 |
si(x) | γ | x0 | |x1−x0| | |x2−x1| | |x3−x2| | |x4−x3| | ACOC | Time |
s1(x) | 0 | 0.9 | 4.2417e-1 | 4.1039e-2 | 1.7406e-5 | 1.2675e-15 | 3.0060079 | 0.395263 |
1 | 0.9 | 3.2788e-1 | 1.3573e-1 | 1.6133e-3 | 2.0244e-9 | 3.0657447 | 0.411651 | |
9 | 0.9 | 4.4244e-1 | 54.806 | 5.2807 | 4.7567 | nc | 0.463676 | |
s2(x) | 0 | 1.5 | 2.5974e-2 | 2.0187e-5 | 9.006e-15 | 7.9969e-43 | 3.0000018 | 0.365875 |
1 | 1.5 | 2.5952e-2 | 4.1889e-5 | 1.6095e-13 | 9.1294e-39 | 3.000007 | 0.406378 | |
9 | 1.5 | 2.5778e-2 | 2.1550e-4 | 1.0967e-10 | 1.4441e-29 | 3.0000642 | 0.524439 | |
s3(x) | 0 | 1.4 | 4.4916e-3 | 5.6144e-8 | 1.0866e-22 | 7.8771e-67 | 3.0 | 0.384539 |
1 | 1.4 | 4.4915-3 | 1.1286e-7 | 1.7652e-21 | 6.7536e-63 | 3.0 | 0.413546 | |
9 | 1.4 | 4.4911e-3 | 5.6657e-7 | 1.1166e-18 | 8.5487e-54 | 3.0000001 | 0.703668 | |
s4(x) | 0 | 2.0 | 3.4712e-2 | 1.2797e-5 | 6.0467e-16 | 6.3787e-47 | 3.0000009 | 0.796057 |
1 | 2.0 | 3.4699e-2 | 2.6357e-5 | 1.0566e-14 | 6.8067e-43 | 3.000003 | 0.813411 | |
9 | 2.0 | 3.459e-2 | 1.3483e-4 | 7.0755e-12 | 1.022e-33 | 3.0000251 | 0.801792 | |
s5(x) | 0 | 0.6 | 5.189e-3 | 9.9857e-9 | 7.044e-26 | 2.4725e-77 | 3.0 | 0.794878 |
1 | 0.6 | 5.189e-3 | 1.9977e-8 | 1.1281e-24 | 2.0311e-73 | 3.0 | 1.063982 | |
9 | 0.6 | 5.1889e-3 | 9.9912e-8 | 7.0557e-22 | 1.022e-33 | 3.0000251 | 0.801792 | |
s6(x) | 0 | 1.6 | 1.1355 | 3.5448e-1 | 9.0195e-3 | 1.6249e-7 | 2.975614 | 0.407156 |
1 | 1.6 | 1.0323 | 4.3713e-1 | 2.9580e-2 | 1.119e-5 | 2.925868 | 0.542887 | |
9 | 1.6 | 0.20652 | 2.785e-1 | 3.5908e-1 | 3.8939e-1 | nc | 0.554967 |
γ | iter | ‖x(k+1)−x(k)‖ | ‖S(x(k+1))‖ | ACOC |
-1 | 6 | 7.188e-505 | 1.457e-2017 | 4.00000 |
0 | 7 | 1.126e-440 | 5.584e-1321 | 3.00000 |
1 | 7 | 8.525e-362 | 4.847e-1084 | 3.00000 |
9 | 8 | 1.170e-450 | 6.260e-1350 | 3.00000 |
γ | iter | ‖x(k+1)−x(k)‖ | ‖S(x(k+1))‖ | ACOC |
-1 | 5 | 2.815e-312 | 8.479e-1250 | 4.00017 |
0 | 6 | 4.212e-367 | 2.286e-1102 | 3.00001 |
1 | 6 | 6.901e-327 | 2.011e-981 | 3.00001 |
9 | 6 | 1.800e-239 | 1.786e-718 | 3.00001 |