Loading [MathJax]/extensions/TeX/mathchoice.js
Research article Special Issues

Design and dynamical behavior of a fourth order family of iterative methods for solving nonlinear equations

  • In this paper, a new fourth-order family of iterative schemes for solving nonlinear equations has been proposed. This class is parameter-dependent and its numerical performance depends on the value of this free parameter. For studying the stability of this class, the rational function resulting from applying the iterative expression to a low degree polynomial was analyzed. The dynamics of this rational function allowed us to better understand the performance of the iterative methods of the class. In addition, the critical points have been calculated and the parameter spaces and dynamical planes have been presented, in order to determine the regions with stable and unstable behavior. Finally, some parameter values within and outside the stability region were chosen. The performance of these methods in the numerical section have confirmed not only the theoretical order of convergence, but also their stability. Therefore, the robustness and wideness of the attraction basins have been deduced from these numerical tests, as well as comparisons with other existing methods of the same order of convergence.

    Citation: Alicia Cordero, Arleen Ledesma, Javier G. Maimó, Juan R. Torregrosa. Design and dynamical behavior of a fourth order family of iterative methods for solving nonlinear equations[J]. AIMS Mathematics, 2024, 9(4): 8564-8593. doi: 10.3934/math.2024415

    Related Papers:

    [1] Zhen Lin . On the sum of the largest Aα-eigenvalues of graphs. AIMS Mathematics, 2022, 7(8): 15064-15074. doi: 10.3934/math.2022825
    [2] Bilal A. Rather, Fawad Ali, Nasim Ullah, Al-Sharef Mohammad, Anwarud Din, Sehra . Aα matrix of commuting graphs of non-abelian groups. AIMS Mathematics, 2022, 7(8): 15436-15452. doi: 10.3934/math.2022845
    [3] Alfredo Sotelo-Pejerrey . Traces of certain integral operators related to the Riemann hypothesis. AIMS Mathematics, 2023, 8(10): 24971-24983. doi: 10.3934/math.20231274
    [4] Shuo Wang, Juan Liu, Xindong Zhang . Properties of solutions for fractional-order linear system with differential equations. AIMS Mathematics, 2022, 7(8): 15704-15713. doi: 10.3934/math.2022860
    [5] Efruz Özlem Mersin . Sturm's Theorem for Min matrices. AIMS Mathematics, 2023, 8(7): 17229-17245. doi: 10.3934/math.2023880
    [6] Samer Al-Ghour, Hanan Al-Saadi . Soft weakly connected sets and soft weakly connected components. AIMS Mathematics, 2024, 9(1): 1562-1575. doi: 10.3934/math.2024077
    [7] Meraj Ali Khan, Ali H. Alkhaldi, Mohd. Aquib . Estimation of eigenvalues for the α-Laplace operator on pseudo-slant submanifolds of generalized Sasakian space forms. AIMS Mathematics, 2022, 7(9): 16054-16066. doi: 10.3934/math.2022879
    [8] Jun He, Yan-Min Liu, Jun-Kang Tian, Ze-Rong Ren . A note on the inclusion sets for singular values. AIMS Mathematics, 2017, 2(2): 315-321. doi: 10.3934/Math.2017.2.315
    [9] Xiaofei Cao, Yuyue Huang, Xue Hua, Tingyu Zhao, Sanzhang Xu . Matrix inverses along the core parts of three matrix decompositions. AIMS Mathematics, 2023, 8(12): 30194-30208. doi: 10.3934/math.20231543
    [10] Xiaomei Zhang, Xiang Chen . Uniqueness of difference polynomials. AIMS Mathematics, 2021, 6(10): 10485-10494. doi: 10.3934/math.2021608
  • In this paper, a new fourth-order family of iterative schemes for solving nonlinear equations has been proposed. This class is parameter-dependent and its numerical performance depends on the value of this free parameter. For studying the stability of this class, the rational function resulting from applying the iterative expression to a low degree polynomial was analyzed. The dynamics of this rational function allowed us to better understand the performance of the iterative methods of the class. In addition, the critical points have been calculated and the parameter spaces and dynamical planes have been presented, in order to determine the regions with stable and unstable behavior. Finally, some parameter values within and outside the stability region were chosen. The performance of these methods in the numerical section have confirmed not only the theoretical order of convergence, but also their stability. Therefore, the robustness and wideness of the attraction basins have been deduced from these numerical tests, as well as comparisons with other existing methods of the same order of convergence.



    Let

    G=(V(G),E(G))

    be a simple undirected graph on n vertices such that V(G) and E(G) denote the vertex set and edge set of G, respectively. Let V(G)={v1,v2,,vn}, for viV(G), and dG(vi)=d(vi)=di denotes the degree of vertex vi and the set of vertices adjacent to vi, denoted by N(vi), refers to the neighborhood of vi. The maximum and the minimum degree of G are denoted by Δ=Δ(G) and δ=δ(G), respectively. λ(M), λk(M) and λmin(M) denote the largest eigenvalue, the k-th largest eigenvalue, and the smallest eigenvalue of a matrix M, respectively. The set of all eigenvalues for a matrix M together with their multiplicities is called the M-spectrum. A real symmetric matrix M is called positive semidefinite if λmin(M)0. Denote by Kn, Pn, Cn and Ks,ns the complete graph, path, cycle, and complete bipartite graph with n vertices, respectively. Let A(G) and D(G) be the adjacency matrix and the diagonal matrix of the degrees of G, respectively. The signless Laplacian Q(G) of G is defined as

    Q(G)=D(G)+A(G).

    The Laplacian L(G) of G is defined as

    L(G)=D(G)A(G).

    We denote the eigenvalues of Q(G) and L(G) by ρk(G) and λk(G), respectively, where 1kn (we drop G when it is clear from the context). In particular, the largest eigenvalue of Q(G) is called the signless Laplacian spectral radius of G, denoted by ρ(G), and the largest eigenvalue of L(G) is called the Laplacian spectral radius of G, denoted by λ(G).

    In [1], Nikiforov introduced the Aα-matrix of G, which is the convex linear combination of A(G) and D(G) defined by

    Aα(G)=αD(G)+(1α)A(G),0α1.

    Clearly,

    A0(G)=A(G),  A1(G)=D(G)and2A12(G)=Q(G).

    Thus, the matrix Aα(G) is a generalization of the adjacency matrix and the signless Laplacian matrix. The largest eigenvalue of Aα(G) is called the Aα-spectral radius of G, denoted by ρα(G). Nikiforov [1] had studied the matrix Aα(G), and he has investigated many properties on Aα(G), including bounds on the k-th largest (especially, the largest, the second largest, and the smallest) eigenvalue of Aα(G), the positive semidefiniteness of Aα(G), etc. From then on, the study of Aα-spectra and the Aα-spectral radius for the graph has attracted the attention of many researchers. In [2], Nikiforov et al. gave an upper bound of the spectral radius of Aα(TΔ), where TΔ is the tree of maximal degree Δ. Also, he obtained several bounds on the spectral radius of Aα of general graphs. The graphs with

    λk(Aα(G))=αn1,2kn

    had been characterized by Lin et al. in [3]. In [4], Guo and Zhou gave upper bounds for the Aα-spectral radius for some graphs under some conditions. Guo and Zhang [5]; obtained a sharp upper bound on the Aα-spectral radius for α[12,1), and proved that for two connected graphs G1 and G2, ρα(G1)>ρα(G2) under certain conditions. For further related studies, one may see [6,7,8].

    The research on L(G) has shown that it is a remarkable matrix in many respects. Since; L(G) is just the difference between D(G) and A(G), then to understand to what extent each of the summands A(G) and D(G) determines the properties of L(G), and motivated by the above works, we introduce the Aα-matrix defined by

    Aα(G)=αD(G)+(α1)A(G),0α1. (1.1)

    Through our study, several facts indicate that the study of the Aα(G) family is of a unique importance. First of all, obviously

    A(G)=A0(G),D(G)=A1(G)andL(G)=2A12(G).

    Since A12(G) is essentially equivalent to L(G), then the matrix Aα(G) can be regarded as a generalization of the Laplacian matrix L(G). Since

    A1(G)=A1(G)=D(G);

    unless otherwise specified, we only consider the case of 0α<1 in this paper.

    In this paper, we study some properties on Aα(G) and obtain some bounds of the largest eigenvalue of Aα(G). For the relation between the Q-spectrum and L-spectrum of a connected graph G, it is well known that G is bipartite if and only if Q(G) and L(G) has the same spectrum ([9, Proposition 1.3.10]), and we extend this relation on the Aα- and Aα-spectrum. In [10], Akbari et al. proved that, for a graph G of order n and any real number β, if 0<β1 or 2β3, then

    nk=1ρβknk=1λβk;

    and if 1β2, then

    nk=1ρβknk=1λβk.

    In this work we present the relation between the sum of powers of Aα and Aα eigenvalues, which is a generalization of the relation between the sum of powers of Q(G) and L(G) eigenvalues in [10].

    The rest of this paper is organized as follows. In Section 2, we give some preliminaries and lemmas used later. In Section 3, we derive new basic properties of Aα(G). In Section 4, we present the Aα-spectra of the complete graphs and the complete split graphs. Section 5 gives some bounds of the largest eigenvalue of Aα(G). In Section 6, we deduce the relation between the Aα- and Aα-spectral radius. Finally, in Section 7, we prove the relation between the sum of powers of Aα and Aα eigenvalues.

    Although Weyl's inequalities have been known for almost a century (see, e.g., [11]), their equality case was first established by So in [12], and for convenience we state Weyl and So's complete theorem as follows:

    Theorem 2.1. Let A and B be Hermitian matrices of order n, and let 1in and 1jn. Then,

    λi(A)+λj(B)λi+jn(A+B),ifi+jn+1, (2.1)
    λi(A)+λj(B)λi+j1(A+B),ifi+jn+1. (2.2)

    In either of these inequalities, the equality holds if and only if there exists a nonzero n-vector that is an eigenvector to each of the three eigenvalues involved.

    A simplified version of (2.1) and (2.2) gives

    λk(A)+λmin(B)λk(A+B)λk(A)+λ(B). (2.3)

    We shall need the following lemmas for our new results.

    Lemma 2.2. [3] Let G be a graph of order n. If eE(G) and α12, then,

    ρi(Aα(G))ρi(Aα(Ge)),for1in.

    Lemma 2.3. [13] Let B be a real symmetric nonnegative irreducible matrix and λ be the largest eigenvalue of B. ZRn. If ZtBZ=λ and Z=1, then BZ=λZ.

    Lemma 2.4. [5] Let G be a connected graph with n4 vertices and m edges. If α[12,1), then,

    ρα(G)max{αΔ(G),(1α)(mn12)}+2α.

    Equality holds if and only if α=n1n+1 and G=Kn.

    For a graph G of order n, suppose that λ is an eigenvalue of Aα(G) and x is an eigenvector of Aα(G) with respect to λ. We use x(v) to denote the entry of x corresponding to the vertex vV(G). It is clear that the system of eigenequations for the matrix Aα(G) is

    λx(v)=αdG(v)x(v)+(α1)uN(v)x(u). (3.1)

    If G is a graph of order n with Aα(G)=Aα, and x is a real vector, the quadratic form Aαx,x can be represented in several equivalent ways:

    Aαx,x=uvE(G)(αx(u)2+2(α1)x(u)x(v)+αx(v)2), (3.2)
    Aαx,x=(2α1)uV(G)x(u)2dG(u)+(1α)uvE(G)(x(u)x(v))2, (3.3)
    Aαx,x=αuV(G)x(u)2dG(u)+2(α1)uvE(G)x(u)x(v). (3.4)

    Each of these representations can be useful in proofs.

    Now, we give some of the spectral properties of the Aα-matrix. Let us call the largest eigenvalue of Aα(G) the Aα-spectral radius of G, and denote it as λα(G). Let us also denote the smallest eigenvalue of Aα(G) as μα(G). Since Aα(G) is a real symmetric matrix, and by using Rayleigh's principle, the following result holds:

    Proposition 3.1. If α[0,1) and G is a graph of order n, then

    λα(G)=maxx2=1Aα(G)x,xandμα(G)=minx2=1Aα(G)x,x. (3.5)

    Moreover, if x is a unit n-vector, then,

    λα(G)=Aα(G)x,x

    if and only if x is an eigenvector to λα(G), and

    μα(G)=Aα(G)x,x

    if and only if x is an eigenvector to μα(G).

    By using these relations, the following result is evident:

    Proposition 3.2. If α[0,1) and G is a graph, then

    λα(G)=max{λα(H):HisacomponentofG},
    μα(G)=min{μα(H):HisacomponentofG}.

    It is clear that if G is a d-regular graph of order n, then

    Aα(G)=αdIn+(α1)A(G),

    and so there is a linear correspondence between the spectra of Aα(G) and of A(G),

    λk(Aα(G))=αd+(α1)λnk+1(A(G)), 1kn. (3.6)

    In particular, if k=n, then

    μα(G)=(2α1)d

    for any α[0,1]. As a consequence of Weyl's inequality (2.3), the following result is immediate:

    Proposition 3.3. If α[0,1] and G is a graph with

    A(G)=AandAα(G)=Aα,

    then,

    αδ+(α1)λnk+1(A)λk(Aα)αΔ+(α1)λnk+1(A).

    An important property of the Laplacian L(G) is that it is positive semidefinite. This is certainly not true for Aα(G) if α<12 and G is a regular graph, but if α12, then Aα(G) is like L(G). We give this result in the following:

    Proposition 3.4. If α12, and G is a graph, then Aα(G) is positive semidefinite, and if α>12 and G has no isolated vertices, then Aα(G) is positive definite. Moreover, if α<12 and G is a regular graph, then Aα(G) is not positive semidefinite.

    Proof. Let x be a nonzero vector. If α12, then for any edge uvE(G), we see that

    Aα(G)x,x(2α1)(x(u)2+x(v)2)+(1α)(x(u)x(v))20. (3.7)

    Hence, Aα(G) is positive semidefinite. Now, assume that G has no isolated vertices. Choose a vertex u with x(u)0 and let uvE(G). Then, (3.7) becomes a strict inequality, and so Aα(G) is positive definite. Finally, suppose that G is a dregular graph, α<12, and let A be its adjacency matrix, then λ(A)=d. Thus

    μα(G)=(2α1)d<0.

    By the above proposition we get the following lemma which gives a relation between the Aα eigenvalues of G and the Aα eigenvalues of spanning subgraphs of G.

    Lemma 3.5. Let G be a graph of order n and let α[12,1). If G=Ge, where eE(G), then,

    λi(G)λi(G)for1in.

    Proof. Let e=uv such that u,vV(G). It is easy to see that

    Aα(G)=Aα(G)+M,

    where M is the matrix of order n indexed by the vertices of G having (u,v)th and (v,u)th entries both equal to α1, and the (u,u)th and (v,v)th entries both equal to α and all other entries equal to zero, hence M is an Aα-matrix of a graph containing only one edge. Since α[12,1), then Aα(G), Aα(G) and M are positive semidefinite and Weyl's inequalities (2.3) imply that λi(G)λi(G).

    Equality (3.6) and the fact the eigenvalues of A(Kn) are {n1,1,,1} give the spectrum of Aα(Kn) as follows:

    Proposition 4.1. The eigenvalues of Aα(Kn) are

    μα(Kn)=(2α1)(n1)andλk(Aα(Kn))=α(n2)+1for1kn1.

    If SV(G), then we use G[S] to denote the subgraph of G induced by S. Recall that G[S] is an independent set if no two vertices of S are adjacent and G[S] is a clique if it is a complete subgraph of G. The graph Kr(nr)K1 is called a complete split graph, denoted by CSr,nr. The work in the following proposition is motivated by the proof of [3, Proposition 2.4].

    Proposition 4.2. Let G be a graph with Aα(G), and α[0,1). Let SV(G) and |S|=k. Suppose that dG(u)=d for each vertex uS, and N(v){w}=N(w){v} for any two vertices v,wS. Then, we have the following statements:

    (i) If G[S] is a clique, then α(d1)+1 is an eigenvalue of Aα(G) with multiplicity at least k1.

    (ii) If G[S] is an independent set, then αd is an eigenvalue of Aα(G) with multiplicity at least k1.

    Proof. Let S={v1,v2,,vk}. Clearly, d1==dk=d. Let z1,z2,,zk1 be vectors such that

    {zi(v1)=1,zi(vi+1)=1,zi(v)=0,if vV(G){v1,vi+1},

    for i=1,,k1. Assume that G[S] is a clique. It is easy to obtain that

    Aα(G)zi=(α(d11)+1,0,,0,α(1di)1,0,,0)=(α(d1)+1)zi,

    for i=1,,k1. Hence, α(d1)+1 is an eigenvalue of Aα(G) and z1,z2,,zk1 are eigenvectors of Aα(G) corresponding to α(d1)+1. In addition, since z1,z2,,zk1 are linearly independent, the multiplicity of α(d1)+1 is at least k1. Now, supposing that G[S] is an independent set, we have

    Aα(G)zi=(αd1,0,,0,αdi,0,,0)=αdzi,

    for i=1,,k1. Since z1,z2,,zk1 are linearly independent, it follows that αd is an eigenvalue of Aα(G) with multiplicity at least k1. Thus, the proof is completed.

    Consider an n×n real symmetric matrix

    S=(S11S12S1tS21S22S2tSt1St2Stt),

    whose rows and columns are partitioned according to a partitioning P1, P2, , Pt of {1,2,,n}. The quotient matrix B of the matrix S is the t×t matrix whose entries are the average row sums of the blocks Sij of S. The partition is equitable if each block Sij of S has constant row sum.

    Lemma 4.3. [14] Let B be an equitable quotient matrix of a symmetric real matrix S. If λ is an eigenvalue of B, then λ is also an eigenvalue of S.

    Now, we can determine all Aα-eigenvalues of CSr,nr as follows:

    Proposition 4.4. The Aα-spectrum of CSr,nr contains α(n2)+1 with multiplicity r1, αr with multiplicity nr1, and the remaining two Aα-eigenvalues are

    α(n+2(r1))+1±(α(n+2(r1))+1)2+4r(12α)(nr+α(r1))2. (4.1)

    Proof. We can write Aα(CSr,nr) as

    Aα(CSr,nr)=((1α)Jr,r+(αn1)Ir(1α)Jr,nr(1α)Jnr,rαrInr,nr).

    Then, the quotient matrix of Aα(CSr,nr) is equitable and it can be written in the form

    B(Aα(CSr,nr))=((nr)α+r1(1α)(nr)(1α)rαr).

    Thus, by Lemma 4.3, the eigenvalues of B(Aα(CSr,nr)) are eigenvalues of Aα(CSr,nr), and according to Proposition 4.2, we get the result.

    In this section we give a few general bounds on λα(G).

    Proposition 5.1. Let G be a graph, with

    Δ(G)=Δ,  A(G)=AandD(G)=D.

    The following inequalities hold for λα(G):

    λα(G)αΔ+(α1)λ(A), (5.1)
    λα(G)(2α1)λ(A). (5.2)

    Proof. Inequality (5.1) follows by Weyl's inequalities (2.3) because

    αΔ+(α1)λ(A)=αλ(D)+(α1)λ(A)=λ(αD)+λmin((α1)A)λ(αD+(α1)A)=λα(G).

    To prove the inequality (5.2), let H be a component of G such that λ(A)=λ(A(H)). Let x be a positive unit vector to λ(A(H)). For every edge uvE(H), the AM-GM inequality implies that

    2x(u)x(v)=2αx(u)x(v)+2(1α)x(u)x(v)α(x(u)+x(v))22+2(1α)x(u)x(v)=12(αx(u)2+2(α1)x(u)x(v)+αx(v)2)+(32α)x(u)x(v).

    Summing this inequality over all edges uvE(H), and using (3.2), we get

    λ(A)=λ(A(H))=A(H)x,x12Aα(H)x,x+12(32α)λ(A);

    and then we get

    2λ(A)(32α)λ(A)Aα(H)x,xλα(G);

    hence

    (2α1)λ(A)λα(G).

    Having inequality (5.2) in hand, if α12, then every lower bound of λ(A) gives a lower bound on λα(G), but if α<12, then every upper bound of λ(A) gives a lower bound on λα(G). We mention just two such bounds.

    Corollary 5.2. Let G be a graph such that α12. If G is of order n and has m edges, then

    λα(G)(2α1)2mn.

    Corollary 5.3. Let G be a connected graph such that α<12. If G is of order n and has m edges, then

    λα(G)(2α1)2mn+1.

    Proposition 5.4. Let G be a graph of order n. If α(12,1], then,

    λα(G)α(n2)+1.

    Moreover,

    λα(G)=α(n2)+1

    with multiplicity k1 if G has k vertices of degree n1.

    Proof. Applying Lemma 3.5 leads to

    λα(G)λα(Kn)=α(n2)+1.

    If G has k vertices of degree n1, then it follows from Proposition 4.2 that α(n2)+1 is an eigenvalue of Aα(G) with multiplicity at least k1, and since

    λα(G)α(n2)+1,

    we get

    λα(G)=α(n2)+1

    with multiplicity k1.

    Let G be a connected graph. Merris [15] pointed out λ(L(G))ρ(Q(G)), and the equality holds if G is a bipartite graph. In the next result, we generalize this result to the Aα- and Aα-spectral radius of a connected graph.

    Theorem 6.1. Let G be a graph of order n, α(0,1), λα(G)=λα and ρα(G)=ρα. We have λαρα. Moreover, if G is connected, then the equality holds if and only if G is bipartite.

    Proof. Let

    V(G)={v1,v2,,vn},  x=(x1,x2,,xn)tRn

    be an arbitrary vector such that x=1. Let

    y=(y1,y2,,yn)tRn

    be a unit eigenvector of Aα(G) belonging to λα and

    y=(y1,y2,,yn)tRn

    be a unit eigenvector of Aα(G) belonging to ρα. Let

    |y|=(|y1|,|y2|,,|yn|)t.

    First, we prove that λαρα.

    λα=maxxtAαx=maxxt(αD+(α1)A)x=max(αni=1x2idi+2(α1)vivjE(G)xixj)=ytAαy=αni=1y2idi+2(α1)vivjE(G)yiyj

    and

    ρα=maxxtAαx=maxxt(αD+(1α)A)x=max(αni=1x2idi+2(1α)vivjE(G)xixj)=ytAαy=αni=1yi2di+2(1α)vivjE(G)yiyj.

    Thus,

    λα=αni=1y2idi+2(α1)vivjE(G)yiyjαni=1y2idi+2(1α)vivjE(G)|yiyj|=|yt|Aα|y|maxxtAαx=ρα. (6.1)

    Now, if G is bipartite, then the matrix Aα and the matrix Aα are similar by a diagonal matrix D with diagonal entries ±1, that is, Aα=DAαD1. Therefore, Aα and Aα have the same spectrum, and thus we get λα=ρα. Finally, when G is connected and λα=ρα, all inequalities (6.1) must be equalities. By Lemma 2.3 and the equality

    |yt|Aα|y|=ρα;

    we know that |y| is an eigenvector of Aα belonging to ρα. So, |y|=±y. Using the Perron-Frobenius' theorem for Aα(G), we have y>0, |y|=y, and |yi|>0, i=1,2,,n.

    Since

    αni=1y2idi+2(α1)vivjE(G)yiyj=αni=1y2idi+2(1α)vivjE(G)|yiyj|;

    we get

    vivjE(G)yiyj=vivjE(G)|yiyj|,

    hence, |yiyj|=yiyj when vivjE(G). Therefore, yiyj<0 if vivjE(G).

    Let

    U={vi:yi>0}andW={vj:yj<0}.

    For each edge e=vivj, we have yiyj<0, one of the vertices of edge e is in U, and the other is in W. So, G is a bipartite graph.

    Remark 6.2. If G is not bipartite and is not connected, and λα is the largest eigenvalue of Aα for a bipartite connected component of G, then the equality in Theorem 6.1 still holds.

    Example 6.3. Take

    G=C3C4.

    Then, G is not bipartite and not connected. Now, A14(G) has a spectrum 2[1],1.25[2],0.5[2],1[2] (where λ[i] means the eigenvalue λ is repeated i times in the spectrum). On the other hand, A14(G) has a spectrum 2[2],0.5[2],0.25[2],1[1]. Then, we have

    λ14(G)=ρ14(G)=2.

    Note that 2 is the largest eigenvalue of A14(C4) as well.

    Now, we introduce the relationship between the Aα- and Aα-spectra of bipartite graphs, which is a generalization [9, Proposition 1.3.10], and it follows from the proof of the above theorem.

    Corollary 6.4. Let G be a connected graph. Then, G is bipartite if and only if the Aα-spectrum and Aα-spectrum are equal.

    Remark 6.5. In fact, if G is bipartite and is not connected, the Aα-spectrum still equals the Aα-spectrum.

    Example 6.6. Take G=P3C4. Then G is bipartite and not connected. We have that A34(G) has a spectrum 2[1],1.64039[1],1.5[2],1[1],0.75[1],0.609612[1]. In contrast, A34(G) has the same spectrum as A34(G), although G is not connected.

    According to Corollary 6.4 and [1, Propositions 38 and 39], we get the following two results:

    Corollary 6.7. The Aα-spectrum and the Aα-spectrum of the complete bipartite graph Ka,b are equal, that is, if ab1 and α(0,1), the eigenvalues of Aα(Ka,b) are

    λα(Ka,b)=12(α(a+b)+α2(a+b)2+4ab(12α)),
    μα(Ka,b)=12(α(a+b)α2(a+b)2+4ab(12α)),
    λk(Aα(Ka,b))=αafor1<kb,
    λk(Aα(Ka,b))=αbforb<k<a+b.

    Corollary 6.8. The Aα-spectrum and the Aα-spectrum of the star K1,n1 are equal, that is, the eigenvalues of Aα(K1,n1) are

    λα(K1,n1)=12(αn+α2n2+4(n1)(12α)),
    μα(K1,n1)=12(αnα2n2+4(n1)(12α)),
    λk(Aα(K1,n1))=αfor1<k<n.

    Indeed, many practical results of Theorem 6.1 can be deduced, and here are some of them.

    Proposition 6.9. Let G be a connected graph with n4 vertices and m edges. If α[12,1), then

    λα(G)<max{αΔ(G),(1α)(mn12)}+2α. (6.2)

    Proof. Let

    λα(G)=λαandρα(G)=ρα.

    Then, Theorem 6.1 and Lemma 2.4 lead to

    λαmax{αΔ(G),(1α)(mn12)}+2α. (6.3)

    Suppose that the equality in (6.3) holds, thus ρα=λα, and so G is bipartite and

    ρα=max{αΔ(G),(1α)(mn12)}+2α.

    Therefore, by Lemma 2.4, G=Kn, and thus G is bipartite if and only if n=2, but n4, and hence the inequality is strict.

    Theorem 6.1 and [2, Theorem 2] lead directly to the next result.

    Proposition 6.10. If T is a tree of order n and α[0,1], then

    λα(T)nα+n2α2+4(n1)(12α)2.

    The equality holds if and only if T is the star K1,n1.

    By Theorem 6.1 and [1, Proposition 20], we get the next result.

    Proposition 6.11. If G is a graph with no isolated vertices, then

    λα(G)maxuV(G){αd(u)+1αd(u)uvE(G)d(v)}.

    If α(12,1) and G is connected, the equality holds if and only if G is a regular bipartite graph.

    In this part we give two explicit expressions for the sums and the sum of squares of the eigenvalues of Aα and Aα, considering [1, Propositions 34 and 35].

    Proposition 6.12. If G is a graph of order n and has m edges, then

    ni=1λi(Aα(G))=ni=1λi(Aα(G))=trAα(G)=trAα(G);

    where

    ni=1λi(Aα(G))=αuV(G)d(u)=2αm.

    A similar formula for the sum of the squares of the Aα and Aα-eigenvalues is given as follows:

    Proposition 6.13. If G is a graph of order n and has m edges, then

    ni=1λ2i(Aα(G))=ni=1λ2i(Aα(G))=trA2α(G)=trA2α(G);

    where

    ni=1λ2i(Aα(G))=α2uV(G)d2(u)+2(1α)2m.

    Proof. Let

    Aα=Aα(G),  A=A(G)andD=D(G).

    Calculating the square A2α and taking its trace, we find that

    trA2α=tr(α2D2+(1α)2A2+α(α1)DA+α(α1)AD)=α2trD2+(1α)2trA2α(1α)trDA+α(α1)trAD=α2uV(G)d2(u)+2(1α)2m.

    Pirzada et al. [16] introduced the sum of the βth powers of the Aα eigenvalues of G as

    Sαβ(G)=ni=1ρβi.

    Now, we have the notation

    Sαβ(G)=ni=1λβi

    for the sum of the βth powers of the Aα eigenvalues of G. The following theorem is a generalization [10, Theorem 2].

    Theorem 7.1. Let G be a graph of order n>1, α[12,1), and let β be a real number.

    (i) If 0<β1 or 2β3, then

    Sαβ(G)Sαβ(G).

    (ii) If 1β2, then

    Sαβ(G)Sαβ(G).

    For β(0,1)(2,3), the equality holds in (i) if and only if G is a bipartite graph. Moreover, for β(1,2), the equality holds in (ii) if and only if G is a bipartite graph.

    Proof. We recall that, for any real number r, the binomial series

    \sum\limits_{k = 0}^{\infty}\binom{r}{k}x^k

    converges to (1+x)^r if |x| < 1 . This also remains true for x = -1 if r > 0 (see, e.g., [17, p. 419]). By Lemma 3.5, we find that,

    \lambda_{\alpha^-}(G)\leq\lambda_{\alpha^-}(K_n) = \alpha(n-2)+1 \, \, \, {\text{and}} \, \, \, \rho_{\alpha}(G)\leq\rho_{\alpha}(K_n) = n-1.

    Since A_\alpha is positive semidefinite when \alpha \in [\frac{1}{2}, 1) , then if \rho_i > 0 , we have

    \lvert\frac{\rho_i}{n} -1\rvert\leq\lvert\frac{n-1}{n}-1\rvert = \lvert\frac{-1}{n}\rvert < 1

    and if \rho_i = 0 , we get

    \frac{\rho_i}{n}-1 = -1.

    Therefore,

    \begin{align*} \frac{S^{\alpha}_{\beta}(G)}{n^\beta}& = {(\frac{\rho_1}{n})}^\beta+ \cdots +{(\frac{\rho_n}{n})}^\beta\\ & = \sum\limits_{k = 0}^{\infty} \binom{\beta}{k}(\frac{\rho_1}{n}-1)^k+\cdots+ \sum\limits_{k = 0}^{\infty} \binom{\beta}{k}(\frac{\rho_n}{n}-1)^k\\ & = \sum\limits_{k = 0}^{\infty} \binom{\beta}{k}tr\left(\frac{1}{n}(\alpha D+(1-\alpha) A)-I\right)^k. \end{align*}

    Also, since A_{\alpha^-} is positive semidefinite when \alpha \in [\frac{1}{2}, 1) , then if \lambda_i > 0,

    \begin{align*} \left\lvert\frac{\lambda_i}{n}-1\right\rvert\leq\left\lvert\frac{\alpha(n-2)+1}{n}-1\right\rvert = \left\lvert(\alpha -1)+\frac{1-2\alpha}{n}\right\rvert < 1 \end{align*}

    and if \lambda_i = 0 , we have

    \frac{\lambda_i}{n}-1 = -1.

    Thus, in a similar manner as above, we obtain that

    \frac{S^{\alpha^-}_{\beta}(G)}{n^\beta} = \sum\limits_{k = 0}^{\infty} \binom{\beta}{k}tr\left(\frac{1}{n}(\alpha D+(\alpha-1) A)-I\right)^k.

    We claim that

    {\text{if}}\, \, k\, \, {\text{is even, }}\, \, \, tr(\alpha D+(1-\alpha) A-n I)^k\leq tr(\alpha D+(\alpha-1) A-n I)^k;
    {\text{if}}\, \, k\, \, {\text{is odd, }}\, \, \, tr(\alpha D+(1-\alpha) A-n I)^k\geq tr(\alpha D+(\alpha-1) A-n I)^k.

    When ((\alpha D-n I)+(1-\alpha)A)^k and ((\alpha D-n I)+(\alpha-1) A)^k are expanded in terms of the powers of \alpha D-n I and (1-\alpha) A , respectively, the terms appearing in both expansions, regardless of their signs, are the same. To prove this claim, we identify the sign of each term in both expansions. Consider the terms in the expansion of ((\alpha D-n I)+(1-\alpha) A)^k , where there are exactly j factors equal to \alpha D-n I , for some j = 0, \, 1, \, \cdots, \, k . The sign of the trace for each of these terms is (-1)^j or 0 because all entries of \alpha D-n I and (1-\alpha) A are non-positive and non-negative, respectively. On the other hand, in each term in the expansion of ((\alpha D-n I)+(\alpha-1) A)^k all factors are matrices with non-positive entries, hence the sign of the trace of each term is (-1)^k or 0 . Therefore, the claim has been proven.

    Now, note that if \, \, \, 0 < \beta < 1 or \, 2 < \beta < 3 , then the sign of \binom{\beta}{k} is (-1)^{k-1} , except that \binom{\beta}{2} > 0 , for 2 < \beta < 3 . According to this, for \, \, \, 0 < \beta < 1 and every k ,

    \binom{\beta}{k}tr(\alpha D+(1-\alpha) A-n I)^k\geq \binom{\beta}{k}tr(\alpha D+(\alpha-1) A-n I)^k.

    This inequality remains true for \, 2\leq\beta\leq 3 as

    tr(\alpha D+(1-\alpha) A-n I)^2 = tr(\alpha D+(\alpha-1) A-n I)^2,

    since tr A_\alpha^2 = tr A_{\alpha^-}^2. Thus, Part (i) is proved. For \, 1 < \beta < 2 , the sign of \binom{\beta}{k} is (-1)^{k-1} , except that \binom{\beta}{1} > 0 . Since tr A_\alpha = tr A_{\alpha^-}, we have

    tr(\alpha D+(1-\alpha) A-n I) = tr(\alpha D+(\alpha-1) A-n I),

    and so part (ii) is similarly proved.

    Now, we examine the equality case. Since A_\alpha and A_{\alpha^-} are similar if G is bipartite, it follows that the equality holds in both (i) and (ii) . Since for any positive integer i , tr A^i equals the total number of closed walks of length i in G , then if G is not bipartite, there exists an odd integer r such that tr A^r > 0 (see [18, Lemma 2.5]). Hence,

    tr(\alpha D+(1-\alpha) A-n I)^r > tr(\alpha D+(\alpha-1) A-n I)^r;

    and so the inequalities in both (i) and (ii) are strict.

    We know that the A_\alpha -spectra and A_{\alpha^-} -spectra of K_n are \{(n-1)^{[1]}, (\alpha n-1)^{[n-1]}\} and \{((2\alpha-1)(n-1))^{[1]}, (\alpha(n-2)+1)^{[n-1]}\} , respectively. Therefore,

    S^{\alpha}_{\beta}(K_n)-S^{\alpha^-}_{\beta}(K_n) = (n-1)^\beta +(n-1)(\alpha n-1)^\beta -((2\alpha-1)(n-1))^\beta - (n-1)(\alpha(n-2)+1)^\beta.

    By Theorem 7.1 and based on our numerical experiments, we propose the following conjecture:

    Conjecture 7.2. For every \alpha\in (\frac{1}{2}, 1) and each integer n\geq 3 , we have

    S^{\alpha}_{\beta}(K_n)-S^{\alpha^-}_{\beta}(K_n)\geq 0

    for any \beta \in [0, 1]\cup[2, \infty) and

    S^{\alpha}_{\beta}(K_n)-S^{\alpha^-}_{\beta}(K_n)\leq 0

    for any \beta \in (-\infty, 0)\cup(1, 2) .

    In this paper, we have introduced the A_{\alpha^-} -matrix of a graph G , which is a generalization of the Laplacian matrix L(G) , and we have studied the basic properties of A_{\alpha^-} ; and derived some bounds for its spectral radius. Furthermore, we have determined the A_{\alpha^-} -spectra for the complete graph and the complete split graph. Building upon previous results, we have extended findings related to the spectral radius of L(G) and Q(G) matrices to the A_\alpha - and A_{\alpha^-} -spectral radius. Specifically, in Theorem 6.1, we have generalized Merris' [15] observation that \lambda(L(G))\leq \rho(Q(G)) with equality holding for bipartite graphs. Additionally, we have extended the known relation that G is bipartite if and only if Q(G) and L(G) share the same spectrum, as demonstrated in Corollary 6.4. Finally, in Theorem 7.1; we have generalized a relation established by S. Akbari et al. in [10]; which relates the sum of powers of the eigenvalues of Q(G) and L(G) . In conclusion, these findings have implications for various applications involving graph analysis.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors are grateful to the referees for their valuable suggestions and comments.

    The authors declare no conflicts of interest.



    [1] F. Ahmad, F. Soleymani, F. Khaksar Haghani, S. Serra-Capizzano, Higher order derivative-free iterative methods with and without memory for systems of nonlinear equations, Appl. Math. Comput., 314 (2017), 199–211. http://dx.doi.org/10.1016/j.amc.2017.07.012 doi: 10.1016/j.amc.2017.07.012
    [2] S. Amat, S. Busquier, S. Plaza, Review of some iterative root-finding methods from a dynamical point of view, Series A: Mathematical Sciences, 10 (2004), 3–35.
    [3] I. Argyros, Á. Magreñán, Iterative methods and their dynamics with applications: a contemporary study, Boca Raton: CRC Press, 2017.
    [4] S. Artidiello, A. Cordero, J. Torregrosa, M. Vassileva, Optimal high-order methods for solving nonlinear equations, J. Appl. Math., 2014 (2014), 591638. http://dx.doi.org/10.1155/2014/591638 doi: 10.1155/2014/591638
    [5] S. Artidiello, A. Cordero, J. Torregrosa, M. Vassileva, Two weighted-order classes of iterative root-finding methods, Int. J. Comput. Math., 92 (2015), 1790–1805. http://dx.doi.org/10.1080/00207160.2014.887201 doi: 10.1080/00207160.2014.887201
    [6] R. Behl, A. Cordero, S. Motsa, J. Torregrosa, A new efficient and optimal sixteenth-order scheme for simple roots of nonlinear equations, Bull. Math. Soc. Sci. Math. Roumanie, 60 (2017), 127–140.
    [7] P. Blanchard, Complex analytic dynamics on the Riemann sphere, Bull. Am. Math. Soc., 11 (1984), 85–141.
    [8] P. Blanchard, The dynamics of Newton's method, Proceedings of Symposia in Applied Mathematics, 49 (1994), 139–154.
    [9] D. Brkić, A note on explicit approximations to Colebrook's friction factor in rough pipes under highly turbulent cases, Int. J. Heat Mass Tran., 93 (2016), 513–515. http://dx.doi.org/10.1016/j.ijheatmasstransfer.2015.08.109 doi: 10.1016/j.ijheatmasstransfer.2015.08.109
    [10] B. Campos, J. Canela, P. Vindel, Dynamics of Newton-like root finding methods, Numer. Algor., 93 (2023), 1453–1480. http://dx.doi.org/10.1007/s11075-022-01474-w doi: 10.1007/s11075-022-01474-w
    [11] P. Chand, F. Chicharro, N. Garrido, P. Jain, Design and complex dynamics of Potra–Pták-Type optimal methods for solving nonlinear equations and its applications, Mathematics, 7 (2019), 942. http://dx.doi.org/10.3390/math7100942 doi: 10.3390/math7100942
    [12] F. Chicharro, A. Cordero, J. Torregrosa, Drawing dynamical parameters planes of iterative families and methods, Sci. World J., 2013 (2013), 780153. http://dx.doi.org/10.1155/2013/780153 doi: 10.1155/2013/780153
    [13] F. Chicharro, A. Cordero, N. Garrido, J. Torregrosa, On the choice of the best members of the Kim family and the improvement of its convergence, Math. Method. Appl. Sci., 43 (2020), 8051–8066. http://dx.doi.org/10.1002/mma.6014 doi: 10.1002/mma.6014
    [14] C. Chun, Construction of Newton-like iterative methods for solving nonlinear equations, Numer. Math., 104 (2006), 297–315. http://dx.doi.org/10.1007/s00211-006-0025-2 doi: 10.1007/s00211-006-0025-2
    [15] C. Chun, B. Neta, J. Kozdon, M. Scott, Choosing weight functions in iterative methods for simple roots, Appl. Math. Comput., 227 (2014), 788–800. http://dx.doi.org/10.1016/j.amc.2013.11.084 doi: 10.1016/j.amc.2013.11.084
    [16] A. Cordero, J. García-Maimó, J. Torregrosa, M. Vassileva, P. Vindel, Chaos in King's iterative family, Appl. Math. Lett., 26 (2013), 842–848. http://dx.doi.org/10.1016/j.aml.2013.03.012 doi: 10.1016/j.aml.2013.03.012
    [17] A. Cordero, J. Torregrosa, Variants of Newton's method using fifth-order quadrature formulas, Appl. Math. Comput., 190 (2007), 686–698. http://dx.doi.org/10.1016/j.amc.2007.01.062 doi: 10.1016/j.amc.2007.01.062
    [18] A. Cordero, J. Torregrosa, On the design of optimal iterative methods for solving nonlinear equations, In: Advances in iterative methods for nonlinear equations, Cham: Springer, 2016, 79–111. http://dx.doi.org/10.1007/978-3-319-39228-8_5
    [19] R. Devaney, A first course in chaotic dynamical systems: theory and experiment, New York: Chapman and Hall/CRC, 2020. http://dx.doi.org/10.1201/9780429280665
    [20] R. Devaney, The Mandelbrot set, the Farey tree, and the Fibonacci sequence, American Mathematical Monthly, 106 (1999), 289–302. http://dx.doi.org/10.2307/2589552 doi: 10.2307/2589552
    [21] P. Fatou, Sur les équations fonctionnelles, Bull. Soc. Math. France, 47 (1919), 161–271. http://dx.doi.org/10.24033/BSMF.998 doi: 10.24033/BSMF.998
    [22] Y. Geum, Study on the dynamical analysis of a family of third-order multiple zero finders, Eur. J. Pure Appl. Math., 16 (2023) 2775–2785. http://dx.doi.org/10.29020/nybg.ejpam.v16i4.4986 doi: 10.29020/nybg.ejpam.v16i4.4986
    [23] P. Jarratt, Some fourth order multipoint iterative methods for solving equations, Math. Comput., 20 (1966), 434–437.
    [24] G. Julia, Mémoire sur l'iteration des fonctions rationnelles, J. Math. Pure. Appl., 1 (1918), 47–245.
    [25] M. Kansal, H. Sharma, Analysis of optimal iterative methods from a dynamical point of view by studying their stability properties, J. Math. Chem., 62 (2024), 198–221. http://dx.doi.org/10.1007/s10910-023-01523-2 doi: 10.1007/s10910-023-01523-2
    [26] M. Khirallah, A. Alkhomsan, Convergence and Stability of Optimal two-step fourth-order and its expanding to sixth order for solving non linear equations, Eur. J. Pure Appl. Math., 15 (2022), 971–991. http://dx.doi.org/10.29020/nybg.ejpam.v15i3.4397 doi: 10.29020/nybg.ejpam.v15i3.4397
    [27] H. Kung, J. Traub, Optimal order of one-point and multi-point iteration, J. ACM, 21 (1974), 643–651. http://dx.doi.org/10.1145/321850.321860 doi: 10.1145/321850.321860
    [28] G. Layek, An introduction to dynamical systems and chaos, New Delhi: Springer, 2015. http://dx.doi.org/10.1007/978-81-322-2556-0
    [29] A. Ostrowski, Solutions of equations and systems of equations, New York: Academic Press, 1966.
    [30] M. Petković, B. Neta, L. Petković, J. Džunić, Multipoint methods for solving nonlinear equations: a survey, Appl. Math. Comput., 226 (2014), 635–660. http://dx.doi.org/10.1016/j.amc.2013.10.072 doi: 10.1016/j.amc.2013.10.072
    [31] S. Qureshi, A. Soomro, A. Shaikh, E. Hincal, N. Gokbulut, A novel multistep iterative technique for models in medical sciences with complex dynamics, Comput. Math. Method. M., 2022 (2022), 7656451. http://dx.doi.org/10.1155/2022/7656451 doi: 10.1155/2022/7656451
    [32] K. Sayevand, R. Erfanifar, H. Esmaeili, On computational efficiency and dynamical analysis for a class of novel multi-step iterative schemes, Int. J. Appl. Comput. Math., 6 (2020), 163. http://dx.doi.org/10.1007/s40819-020-00919-x doi: 10.1007/s40819-020-00919-x
    [33] M. Scott, B. Neta, C. Chun, Basin attractors for various methods, Appl. Math. Comput., 218 (2011), 2584–2599. http://dx.doi.org/10.1016/j.amc.2011.07.076 doi: 10.1016/j.amc.2011.07.076
    [34] J. Traub, Iterative methods for the solution of equations, New York: Chelsea Publishing Company, 1982.
    [35] S. Yaseen, F. Zafar, H. Alsulami, An efficient Jarratt-type iterative method for solving nonlinear global positioning system problems, Axioms, 12 (2023), 562. http://dx.doi.org/10.3390/axioms12060562 doi: 10.3390/axioms12060562
    [36] T. Zhanlav, C. Chun, K. Otgondorj, Construction and dynamics of efficient high-order methods for nonlinear systems, Int. J. Comp. Meth., 19 (2022), 2250020. http://dx.doi.org/10.1142/S0219876222500207 doi: 10.1142/S0219876222500207
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1111) PDF downloads(68) Cited by(3)

Figures and Tables

Figures(7)  /  Tables(10)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog