Home Exploring homotopy with hyperspherical tracking to find complex roots with application to electrical circuits
Article Open Access

Exploring homotopy with hyperspherical tracking to find complex roots with application to electrical circuits

  • Mario A. Sandoval-Hernandez , Hugo Jimenez-Islas EMAIL logo , Hector Vazquez-Leal , Miriam L. Quemada-Villagómez and María de la Luz Lopez-Gonzalez
Published/Copyright: February 3, 2025

Abstract

In the field of applied sciences, systems are frequently modeled using mathematical frameworks that include systems of nonlinear algebraic equations. Identifying their roots, whether real or complex, is of critical importance. The widespread use of complex numbers in science and engineering highlights the importance of accurately determining the complex roots of equations. This article presents a study in which the complex roots of a system of equations are identified through an approach that utilizes homotopy continuation, with the curve being traced using a hyperspherical path tracking technique. Furthermore, this article details five case studies on electrical networks where this method is applied to solve systems of equations containing imaginary coefficients to find mesh currents. The path tracking shows the behavior of system equation in each case study. Finally, an analysis of the precision of the solutions obtained in these case studies is provided, demonstrating an accuracy of up to 15 SDs in a single iteration during the refinement stage.

MSC 2010: 0101; 65D25; 65H04; 65H10; 68W30

1 Introduction

In the history of complex numbers, a significant contribution was made by the mathematician Niccolò Tartaglia, who shared his secret solution techniques for solving equations with Girolamo Cardano, under the condition that Cardano would not disclose these techniques [1]. However, in 1545, Cardano published “Ars Magna”, a work in which he exposed methods for solving cubic equations [2,3]. This publication incited Tartaglia’s fury, as he accused Cardano of betrayal and dishonesty for failing to fulfill his promise. Nonetheless, Lodovico Ferrari, a young mathematician and Cardano’s assistant at the time, claimed that he was present during the meeting between the two mathematicians and asserted that no oath was taken [1].

Rafael Bombelli was familiar with the works of Girolamo Cardano on cubic equations, having studied his “Ars Magna”. However, he found some concepts therein to be confusing and believed it was possible to elucidate them, making the information more accessible to the general public. Bombelli embarked on the development of a formal algebra to handle expressions of the form a + b 1 [1,4]. In addition, in his 1572 work, Bombelli introduced calculus with negative numbers and established the rules for their addition and multiplication. His most significant contribution to algebra was the unconditional acceptance of the symbol 1 as a legitimate number, marking a pivotal advance in the understanding and acceptance of complex numbers within the field of mathematics [1,4].

After the advancements made by Cardano and Bombelli, other mathematicians independently made significant contributions to the development of complex numbers. René Descartes, for instance, introduced the term “imaginary” in his work “La Géométrie” [5], marking the beginning of specific terminology to describe these numbers. On the other hand, John Wallis provided the first geometric interpretation of complex numbers, using line segments in the plane, which represented an important step toward their visual understanding [1].

Caspar Wessel, in his work “On the Analytical Representation of Direction: An Attempt Applied Chiefly to Solving Plane and Spherical Polygons” [6], demonstrated the geometric representation of an imaginary root of the form a + b i using directed segments (vectors), including the addition and multiplication of these directed segments. This contribution was crucial for the geometric representation of complex numbers and their operations. Furthermore, Jean-Robert Argand proposed interpreting the value i as a rotation of 90 degrees in the coordinate plane, known as the Argand plane. In his essay, he also introduced, for the first time, the concept of modulus to indicate the magnitude of vectors and complex numbers, as well as the conventional notation for vectors with a horizontal arrow above the letters that indicate its endpoints [7].

Leonhard Euler standardized the notation for i to represent 1 and used imaginary numbers to solve equations as well as in trigonometry applications [8,9]. In his work, Euler assigned the term “imaginary” to complex numbers, arguing that they were impossible quantities that only existed in the mind, due to their nontangible nature in the context of real numbers. On the other hand, Abraham de Moivre formulated the theorem bearing his name, which establishes the relationship between complex numbers and trigonometry, thus demonstrating their applicability and relevance in broader mathematical areas [10]. William Rowan Hamilton, in 1833, proposed a conceptualization of complex numbers as ordered pairs ( a , b ) , where a and b are real numbers [11,12]. However, in electrical engineering, the notation for complex numbers underwent a significant change, wherein the imaginary unit is denoted by j . This convention was introduced and standardized by Charles Proteus Steinmetz in his seminal work “Complex Quantities and Their Use in Electrical Engineering” in 1893. Steinmetz’s contribution greatly facilitated the analysis of alternating current (AC) circuits by simplifying complex calculations into straightforward algebraic operations. Consequently, he established the standardized use of phasor representation for complex numbers in educational texts on electrical engineering, employing the lower-case letter j to signify the 90-degree rotation operator in AC system analysis [13].

In this context, considering the historical background, complex numbers are characterized by their arithmetic properties and various forms of representation, among them, the phasor [14,15]. These numbers are situated in the complex plane, where the abscissa axis represents the real part and the ordinate axis signifies the imaginary part [14,15]. Their application spans a wide range, covering several fields from mathematics to various engineering disciplines. A notable example of their use is found in fractals, which are geometric objects characterized by a basic structure that, whether fragmented or irregular, repeats at multiple scales. This concept was introduced by mathematician Benoît Mandelbrot in 1977, using the term “fractal” from the Latin “fractus”, meaning broken or fractured [16]. Among the most prominent fractal sets are the Mandelbrot and Julia sets, as highlighted by Monroy [14,17].

In the field of electricity and its subfields the applications of complex numbers have many applications, including the design of transmission lines in communications, where the Smith Chart is utilized for analyzing load impedances [18]. In signal processing, Fourier analysis employs the Fourier transform, which involves the imaginary expression in the exponential function’s argument, while the Fourier series in their complex form incorporates complex numbers [19]. In the qualitative theory of differential equations, the occurrence of complex roots in the eigenvalues of coefficient matrices plays a critical role in identifying behaviors such as centers and spirals, which can be stable or unstable [20]. In addition, in numerical analysis, complex numbers are utilized as starting points (SPs) in attraction basins to analyze the stability of numerical algorithms. They provide valuable information regarding the numerical instability of optimization methods, such as Newton-Raphson (N-R), the secant method, and Halley’s method, among others [2123]. Various authors have developed numerous analytical techniques using attraction basins. For instance, in [24], algorithms were developed for plotting the basins of attraction of fixed-points for a pair of bivariate rational functions using toroidal attraction basins. Moreover, in optics, complex numbers are employed to apply Fresnel integrals in diffraction phenomena [25]. In electrodynamics, they facilitate the study of the reflection and transmission of electromagnetic waves [26]. In the field of antenna design and behavior, complex numbers are applied in analyzing parameters such as impedance measurements, polarization, near-field and electromagnetic analysis, among others [27]. In electrical circuits, complex numbers are fundamental for calculating capacitive and inductive reactances, enabling the analysis of AC circuits [15]. Similarly, in power electrical systems, including industrial motors, networks, and three-phase circuits, they play a crucial role [28]. Furthermore, complex numbers are essential in the analysis of stability and frequency response [29], as well as in electronics for the analysis and synthesis of active filters [30,31]. In engineering, systems are modeled by equations and systems of differential, and occasionally algebraic, equations. Understanding the solutions to these equations is crucial, as they reveal the behavior of the system under given conditions. Specifically, for algebraic equations, the solutions can sometimes be complex. However, when dealing with equations that have real coefficients, these complex solutions appear in conjugate pairs [32].

There are several methods to determine complex roots in polynomial algebraic equations. For example, [32] details the Graeffe method as a technique for finding all the roots of an equation, whether real or imaginary. The Jenkins-Traub method is a globally convergent numerical algorithm that can find the real and complex roots of polynomials of a single variable [33]. Sandoval-Hernandez et al. [34] proposed a hybrid perturbation-N-R method to solve a nonlinear algebraic equation (NAEs) in circuit analysis, specifically to determine the bias points of a tunnel diode using real solutions. Dubeau and Gnang [35] employed the fixed-point and N-R methods to find solutions to nonlinear equations in the complex plane. In [36], Newton’s method was adapted to polynomials with multiple roots in complex variables, applying it to Julia sets. On the other hand, Sandoval-Hernandez et al. [37] used a complex variable to obtain the general formula for solving quadratic equations. This formula is well known because it facilitates finding the different solutions of a quadratic equation, including solutions with complex numbers. In addition to the known methods for determining roots in equations, such as Regula Falsi or N-R, there are homotopy continuation algorithms that allow for solving systems of NAEs that can include transcendental functions. These are considered global convergence methods because they possess the ability to find solutions to nonlinear problems from almost any SP, theoretically ensuring that a solution will be reached regardless of the initial conditions. However, several authors have investigated the identification of complex roots through the application of homotopy continuation methods (HCMs). For instance, Gdawiec et al. [38] employed homotopy continuation for polynomial equations possessing both real and complex roots as a viable alternative to conventional iterative techniques. This approach was adopted to circumvent the issue of divergence.

In [39], an analysis was conducted on the SP utilizing the differential HCM to determine the conditions under which the convergence of these methods remains independent of the SP. This investigation aimed to facilitate the application of this technique to chemical engineering problems. However, the study also explored the behavior of the homotopic path in a case study where certain roots could not be identified due to the homotopic route encountering a region within the complex domain. Following this, in [40], the implementation of continuation homotopy was conducted to determine the roots of polynomial equations with real coefficients by substituting the complex variable z = x + y i to derive a homotopy for the real part and another for the imaginary part, with applications in chemical engineering. Furthermore, in [3840], the authors have focused solely on finding complex roots for a polynomial equation using continuation homotopy without conducting any error or significant digit (SD) analysis. On the other hand, [40] presented the plotting of homotopic tracking for the real and imaginary parts separately against the homotopic parameter τ . Furthermore, Vazquez-Leal et al. [41] proposed the biparameter homotopy with an automated stop criterion, aiming to achieve the lowest number of iterations necessary to complete the tracing of a symmetrical branch. The multiparametric homotopy utilizes multiple parameters, providing greater flexibility in solving systems of nonlinear equations and reducing computation times. In [42], the use of predictor vectors (without correction vectors) was proposed to find the solution to several systems of NAEs with multiple solutions, employing lower central processing unit (CPU) time compared to standard homotopic tracking algorithms.

Due to its ability to find multiple roots in NAEs, homotopy has broad applications in chemical engineering. For example, chemical process control problems and chemical reactions often involve nonlinear equations with multiple solutions, heat diffusion processes, exothermic reaction, natural convection of air in a 2D square cavity at a steady state, and homotopic methods enable the identification of all possible solutions, allowing the selection of the most appropriate one [4345], solution vectors for an elliptic system of 2D partial differential equations for natural convection [46]. Likewise, homotopy methods have been applied in electronic circuits to find the polarization point in nonlinear circuits [47,48,49], and in electronics to determine the operating point in circuits with transistors [5056]. In addition, they have been used in the design of simulators for electronic circuits [57,58].

Recently, Vazquez-Leal and his research team [59] initiated a new line of research in robotics, emphasizing the importance of understanding the behavior of the homotopic path that models the movement trajectory of a mobile device. In [60], homotopy with hyperspherical tracking was applied for route planning, solving an NAE system that models the obstacles to be avoided. Furthermore, in [61,62], homotopy with hyperspherical tracking was utilized in a collision-free trajectory planning algorithm for robotic arms, addressing an NAE system that models the robot’s workspace.

This paper introduces the implementation of continuation homotopy with hyperspherical tracking [44,63], implemented in Maple using symbolic programming. This approach enables comprehensive tracking, considering both the real and imaginary parts integrally without their separation, aimed at resolving various systems of equations with imaginary coefficients within the context of electrical engineering problems. This method facilitates the exploration and examination of the solutions obtained, as well as the analysis of the behavior of their homotopic paths. Furthermore, it enhances the body of literature [64,65] dedicated to the determination of real roots of systems of equations through the use of continuation homotopy.

This paper is organized as follows: Section 2 introduces a brief historical note about continuation homotopy. Section 3 introduces the continuation homotopy method with hyperspherical tracking. Section 4 presents five case studies; the first involves a system of polynomial equations, while the subsequent case studies focus on solving systems of equations applied to electrical circuits. Section 5 details the numerical simulations and discussions. Finally, Section 6 outlines the conclusions of this work.

2 A brief historical review of HCMs

The formal origins of homotopy trace back to the work of French mathematician Henri Poincaré (1854–1912), widely regarded as the founder of algebraic topology. In his series of publications entitled “Analysis Situs”, Poincaré introduced the foundational concepts that would later develop into the theory of homotopy [66]. To fully contextualize the state of the art during Poincaré’s era, it is essential to understand the problem that sparked his interest in topology.

To understand the foundation of homotopy, it is essential to review the transformations involved. Topology studies the properties of geometric bodies that remain invariant under continuous transformations and their inverses [67]. This concept focuses on specific points that do not change when geometric bodies are contracted, expanded, or deformed [68,69]. In Figure 1, the points on the cube are relocated without generating new points during compression or expansion, illustrating a continuous transformation. This property, where f ( x ) does not introduce changes and maintains continuity, is referred to as a “homeomorphism” [70]. While it is easy to visualize deformation in simple shapes like a cube, the correspondence in more complex bodies is less apparent [71]. This introductory concept provides a basis for understanding homotopy methods as a tool for solving algebraic equations.

Figure 1 
               Schematic diagram of Homotopy between two geometric bodies.
Figure 1

Schematic diagram of Homotopy between two geometric bodies.

Poincaré’s contributions to topology were of immense importance, as he introduced algebraic tools that enabled mathematicians to study abstract geometric properties without relying on classical geometric intuition. Homotopy became a cornerstone of this approach, offering an algebraic framework for addressing complex topological problems [70].

One of the first significant practical advances derived from homotopy theory was the fixed-point theorem, formulated by the Dutch mathematician Brouwer [72,73]. This theorem states that any continuous function mapping a closed sphere onto itself must have at least one fixed-point, meaning a point that remains invariant under the function. This theorem forms the philosophical link between Poincaré’s homotopy theory and the development of tools for solving algebraic equations and systems of equations. García and Zangwill further elaborate on the concept introduced by Brouwer [74].

In the context of homotopy theory, Brouwer’s theorem established a direct connection between the abstract concepts of algebraic topology and the solution of concrete problems in NAEs. This theorem marked homotopy’s first major application in developing methods for solving systems of equations, a line of research that would be extensively explored in the following decades [74].

As homotopy theory evolved, new tools were developed to tackle more complex problems. In 1934, French mathematicians Jean Leray and Jules Schauder introduced the concept of “topological degree” in their work “Topologie et équations fonctionnelles” as an extension of homotopy theory for systems of nonlinear equations [75]. The topological degree is an algebraic invariant that measures the “number of solutions” of a nonlinear equation by analyzing the signs of the derivatives or Jacobians of the functions involved.

The topological degree was crucial in formalizing HCMs. These methods involve continuously transforming a system of nonlinear equations into a simpler version with a known solution and then following this transformation until the solution of the original system is obtained. This is accomplished through a homotopic function, which links the solutions of the original system to those of the simplified system via the homotopic variable τ , ranging from τ = 0 to τ = 1 . The topological degree ensures that the number of solutions remains constant throughout the transformation, enabling the tracing of all possible solutions of the system [67,68].

In the 1950s, homotopy methods progressed with the development of corrector-predictor methods, proposed by Ficken in 1951 in his article “The Continuation Method for Functional Equations” [76]. These methods combined solution prediction with iterative correction using the N-R method to ensure convergence. One major improvement was the introduction of tangent predictors to the homotopic curve. In 1961, Haselgrove advanced this approach by applying Newton’s Method as a corrector in systems of nonlinear equations, detailed in his work “The Solution of Non-Linear Equations and of Differential Equations with Two-Point Boundary Conditions”, he details the use of Newton’s n-dimensional algorithm in the solution of systems of n -NAEs with n + 1 variables ( n variables, plus the homotopic parameter τ ) [77]. His contribution provided a general description of a predictor-corrector continuation algorithm, expanding on the earlier work of Davidenko [78], they are as follows:

  1. Horizontal Predictors: Proposed by Lahaye in 1934 [79] and later used by Davidenko in 1953 [78,80], these methods sequentially increase the homotopy parameter. However, they are less efficient when there are returns along the homotopy path, as they require precise control of step size to ensure accuracy.

  2. Differential Predictors: Introduced by Haselgrove in 1961 [77], they use the slope of the homotopic path to predict the next point to correct. Compared to horizontal predictors, they offer predictions closer to the actual path, facilitating convergence on local correction methods.

  3. Polynomial Predictors: Proposed by Rheinboldt in 1975 [81], these methods use polynomial functions, such as Lagrange and Hermite polynomials, to predict the next point on the curve. Although they offer greater accuracy compared to earlier methods, their effectiveness decreases on curves with returns.

  4. Lagrange Polynomials: Implemented by Rheinboldt in 1975 to enhance predictions in continuation processes [81]. Hermite Polynomials: Mentioned by Deufhard in 1979 as a potential option, though no practical applications were provided [82].

  5. Closed Surfaces: The use of closed surfaces as predictors has had limited practical application in homotopy methods and was first introduced by Lyness and McHugh in 1963 [83]. They employed “hypercubes” in numerical integration, comparing their method to traditional techniques like Simpson and Euler-Maclaurin, showing improved precision [84]. After 1963, few publications have reported on closed trajectories as a prediction method for homotopy curves.

  6. Hyperspherical Path Tracking: Jiménez-Islas reviewed the closed surface approach in 1985, conducting various analyses. He was the first to propose hyperspherical tracking ad hyperspherical tracking with variable radius in his 1988 master’s thesis, “Paquete Computacional para la Solución de sistemas de ecuaciones no lineales” [85], at Instituto Tecnológico de Celaya, without knowledge of Ushida’s work [86]. Islas published his results until 1996, applying them to chemical engineering [44].

  7. On the other hand, since 1976, Ushida and Chua worked on solving nonlinear equations in circuit theory using continuation homotopy to address sharp turning points in solution curves [48,49]. In their 1984 article “Tracing Solution Curves of Nonlinear Equations with Sharp Turning Points” [86], they introduced an additional equation for homotopy formulation using arc-length. Previously, this idea had been used in other works [4749,87].

  8. The contribution in this article [86] was to solve the system of n nonlinear equations in ( n + 1 ) variables, directly incorporating the equation for arc-length using the implicit backward-differentiation formula [88,89]. Yamamura later extended the work of Ushida and Chua, unaware of Islas’ contributions [85]. In 1990, Yamamura discussed this work with Inoue from Sanyo Electric Company, and in 1993, he published results using hyperspheres in homotopy with applications in electronic circuits [90].

From 1960s, homotopy methods advanced with the introduction of arc-length parameterization, proposed by Klopfenstein in 1961, which improved solution tracking and mitigated issues related to singularities in the Jacobian [91]. This technique provided greater control over curve tracking. In addition, homotopy methods faced challenges with bifurcations, where solution curves split into multiple branches, signaling changes in solution stability [92]. Alexander and Yorke generalize the HCM in “The Homotopy Continuation Method: Numerically Implementable Topological Procedures” [93], freeing it from its dependence on degree and broadening its application to problems beyond traditional degree arguments, particularly in locating bifurcating sets in nonlinear equations.

Homotopy methods have continued to evolve with the introduction of other techniques to improve their robustness and applicability. For example, Chow et al. [94] presented an approach based on a transversality theorem called “the Parametrized Sard Theorem”. This work gave solutions to problems in which the topological scheme is more complicated, for example, finding the Brouwer fixed-point, finding solutions of two point boundary problems for second-order ordinary differential equations in R n , among others. In [95], the Chow-Yorke algorithm [96] was applied to problems such as the elliptic porous slider, the squeezing of a viscous fluid between parallel plates, the squeezing of a viscous fluid between elliptic plates, and viscous flow between rotating discs with injection on the porous disc, among others. These problems range from fairly simple to extremely difficult, and Newton-type methods either partially or totally failed on all of them. In [97], a method was introduced for solving systems of polynomial equations, specifically when these systems had to be solved repeatedly with different sets of values. This method, referred to as “cheater’s homotopy”, employed a path-tracking technique to find all possible solutions, significantly improving computational efficiency. Unlike traditional methods, which required substantial computational resources based on the overall complexity of the equations, this approach targeted the actual number of solutions, resulting in faster and more practical performance for a variety of applications.

Over time, various homotopy methods have been developed, including Newton homotopy [98], fixed-point homotopy [99,100], affine homotopy [84], Combined fixed-point Newton homotopy [101], probability one homotopy [102], and multiparameter homotopies [103]. In all cases, homotopies find a solution when the homotopic parameter τ = 1 . However, in previous years, obtaining the exact solution at τ = 1 was challenging due to numerical errors in the algorithms, requiring the roots to be calculated experimentally. In his thesis, Islas proposed refining the solution using N-R when the neighborhoods of τ = 1 were reached, before and after this value. His proposal is included in GW-BASIC code in the appendix of his work [85]. However, in 1996, Sosonkina introduced the zero-finding strategy in his work “Note on the End Game in Homotopy Zero Curve Tracking” [104] (without knowledge of Islas’ work). Continuation methods were further improved by the introduction of this strategy, which is crucial for ensuring convergence along the homotopy path, as the process may fail to find solutions without it.

3 Basic concepts of HCMs

Solving systems of NAEs can be challenging using traditional approaches such as the N-R method. The N-R method requires an initial approximation of the solution to begin the process. However, it is well known that for some cases, the N-R method may not converge to a solution from any SP. Often, finding an appropriate initial guess is not straightforward or even feasible. As a result, the reliability of the N-R method for solving all types of problems cannot be assured. An alternative strategy to overcome this limitation involves employing HCM. The HCM methods represent a continuous deformation from a basic problem (easy to solve) to the target problem (difficult to solve), facilitating a smoother pathway to finding solutions.

HCMs are utilized to locate multiple solutions within systems of NAEs ( n equations with n unknowns) of the form

(1) f ( x ) = 0 , where f : R n R n .

The HCMs introduce a perturbation in the system f ( x ) by incorporating a function G ( x ) and a homotopic parameter τ . In this manner, the original system is transformed into a homotopic system expressed as follows:

(2) H i ( x , τ ) = H i ( f ( x ) , τ ) = τ f ( x ) + ( 1 τ ) G ( x ) = 0 ,

where H ( x , τ ) : R n + 1 R n , i = 1 , n , and τ is the homotopy parameter. The function G ( x ) exhibits a known solution or is easily solvable.

Equation (2) represents any homotopy formulation that meets the following conditions:

  1. If τ = 0 and H ( x , 0 ) = G ( x ) = 0 , then a solution for H is known or can be easily found using numerical methods.

  2. If τ = 1 and H ( x , 1 ) = f ( x ) = 0 , then the solution of the original system is obtained.

The deformation process of the system initiates at τ = 0 , and the solutions of interest emerge upon reaching τ = 1 and beyond. Each intersection at τ = 1 constitutes a solution vector for (2). The mapped solutions are presented as a family of curves, γ , as referenced in [39,44,105]. Figure 2 displays four examples of homotopic continuous curves. Here, γ 1 represents a successful case, transitioning from the SP A at τ = 0 to the solution at τ = 1 , whereas the other paths illustrate instances where the HCM method does not succeed. It is noteworthy that the trajectory γ 4 is capable of identifying solutions; however, it remains isolated from the initial point, making it impossible to connect it with γ 1 .

Figure 2 
               Family of homotopy curves.
Figure 2

Family of homotopy curves.

The selection of G ( x ) dictates the formulation of the homotopy. In this study, we employ Newton’s homotopy, defined as G ( x ) = f ( x ) f ( x 0 ) , i.e., the function f ( x ) minus the functions evaluated at x 0 [44]. Thus, by substituting Newton’s homotopy into (2), we obtain

(3) H i ( x , τ ) = f ( x ) ( 1 τ ) G ( x ) , x R n ,

where x 0 is the SP.

The application of HCM requires the implementation of a technique to map the homotopic path [44]. For example, when a hypersphere S 2 with radius r and center O 2 is located on the homotopy path, its contour is required to intersect the homotopy curve at two points, O 1 and O 3 , as shown in Figure 3. The center O 1 is the center of the hypersphere of the previous iteration, O 3 is the center of the hypersphere of the next iteration as depicted in Figure 3.

Figure 3 
               Potential homotopy path.
Figure 3

Potential homotopy path.

3.1 Path tracking of the homotopy trajectory

The formulation of a homotopy system involves n equations and n + 1 variables, that is, n variables of the system plus the homotopy parameter τ [44,60]. Therefore, in the context of the hyperspherical path tracking technique, we must include an additional equation known as the hypersphere, which is given by

(4) S i ( x , τ ) = k = 1 n ( x k c k ) 2 + ( τ c n + 1 ) 2 r 2 = 0 ,

where r is the radius and ( c k , c n + 1 ) is the center of the sphere in the homotopy path. i = 1 , , N are iterations for the hyperspheres S i ( x , τ ) used to trace the homotopic path. By using (3) and (4), we have a system with n + 1 equations given by

(5) H s ( x , τ ) = H 1 ( x , τ ) = 0 , H 2 ( x , τ ) = 0 , H n ( x , τ ) = 0 , S i ( x , τ ) = 0 .

For the case of n variables, the homotopic path γ shown in Figure 3 is formed by the centers O 1 , O 2 , O 3 , , O N and these are calculated with help of equation (4) of the hypersphere S i . To find the coordinates of the centers O i ( c k , c n + 1 ) of the S i must be solve the system (5) (with variables x , τ ) by iterating N times in the path following.

During the tracking algorithm, the center of (4) must be updated at each tracking step and resolved for (5). Subsequently, spherical tracking is conducted, as depicted in Figure 3. Initially, we center the first hypersphere at the SP of the homotopy. Next, a prediction step is performed to approach the next point, and a correction step ensures the path following of the homotopy curve. Finally, the center of the hypersphere is updated at (5), and the process is repeated until the procedure reaches τ = 1 , where interpolation and refinement steps are executed to find a solution for f ( x ) = 0 . Thereafter, the procedure continues until another root is found or the maximum number of iterations is reached. During this process, τ can cross τ = 1 several times from left to right and vice versa.

The predictor-corrector scheme enables the spherical algorithm (SA) to accurately track the homotopy curve γ . Figure 4 illustrates the behavior of the predictor-corrector scheme. The Euler predictor establishes a tangent vector ( x p , τ p ) at the point ( x i , τ i ) along the homotopy path. Subsequently, the predictor point is determined at the intersection of the sphere S i with the tangent vector. The N-R method is then employed as a corrector to ascertain the intersection point ( x i + 1 , τ i + 1 ) between the homotopy path and the sphere of radius r . This yields a point in proximity to the solution of the homotopy system for the subsequent tracking step. The Euler predictor point utilized in SA is expressed as follows:

(6) ( x p , τ p ) = ( c x , c τ ) + r v p ,

where ( x p , τ p ) is the predictor point, ( c x , c τ ) is the center of the sphere, r is the radius of the sphere, and v p is the tangent vector.

Figure 4 
                  N-R corrector in spherical tracking.
Figure 4

N-R corrector in spherical tracking.

The N-R method is used to correct the predictor point to a point on the homotopy path γ . The N-R corrector is described by

(7) ( x j + 1 , τ j + 1 ) = ( x j , τ j ) [ J ( x j , τ j ) ] 1 H s ( x j , τ j ) ,

where j = 0 , 1 , 2 , , n iteration, and ( x j , τ j ) is the point updated in each iteration until reaching a final value; this final value can be regarded as the solution. [ J ( x j , τ j ) ] 1 denotes the inverse Jacobian of H s evaluated at the point ( x j , τ j ) for each j th iteration. This iterative method necessitates an appropriate SP to obtain the solution in a minimum number of iterations. Therefore, for the first iteration, j = 0 , ( x j , τ j ) = ( x p , τ p ) . The N-R method is reiterated until it reaches a maximum number of iterations j max or it converges to a solution ( x i + 1 , τ i + 1 ) with a minimal margin of error, as depicted in Figure 4.

3.2 Implementation of homotopy with spherical path tracking

The predictive section of the hyperspherical path tracking methodology presented in [44] forms the core of the strategy. The predictive homotopic path tracking is applied to a system of nonlinear equations derived from a set of n deformed algebraic equations with Newton’s homotopy. By parameterizing equation (3) with parameter p , the homotopic functions H i ( x , τ ) is defined as H i ( x ( p ) , τ ( p ) ) , is defined as follows:

(8) H i ( x ( p ) , τ ( p ) ) = f ( x ( p ) ) ( 1 τ ( p ) ) f ( x 0 ) , x R n .

It is noteworthy that the addition of the parameter p allows visualization regardless of the number of equations that the nonlinear algebraic system, solved by f ( x ) , has. Therefore, the generation of the homotopic path is defined by H 1 ( p , τ ( p ) ) . The homotopic route tracking is done by generating a tangent vector, τ p , which passes through the initial point x 0 , and allows us to predict the point x 1 . The slopes p are defined by calculating the partial derivatives of each homotopic function, H i ( x ( p ) , τ ( p ) ) . We define the matrix H as follows:

(9) H = H 1 x 1 x 1 p H 1 x n x n p H 1 τ τ p H 2 x 1 x 1 p H 2 x n x n p H 2 τ τ p H n x 1 x 1 p H n x n x n p H n τ τ p , H M n , n + 1 R .

By expressing equation (9) as a system of equations in its matrix form, we have

(10) H 1 x 1 H 1 x n H n x 1 H n x n x 1 p x n p = H 1 τ τ p H n τ τ p .

By applying Cramer’s rule in equation (10), we have the solution to the system

x 1 p = H 1 τ τ p H 1 x n H n τ τ p H n x n H ,

(11) x n p = H 1 x 1 H 1 τ τ p H n x 1 H n τ τ p H ,

where H is the determinant of the system. Simplifying using properties of the determinants, each variable is expressed

(12) x 1 p = ( 1 ) ( 1 ) τ p H 1 H , x 2 p = ( 1 ) τ p H 2 H , x k p = ( 1 ) ( 1 ) τ p H k H ,

where H k , k = 1 , 2 , , n is the determinant from which column k has been removed, according to Cramer’s rule. To determine τ p , we must substitute τ p = ( 1 ) H into (12), and by generalizing for any number of variables, x i p , we obtain

x k p = ( 1 ) k H k ,

(13) τ p = ( 1 ) n + 1 H n + 1 .

In this way, by rewriting (13) in the form of a column vector, a general expression for the slopes is obtained

(14) y ˙ k = ( 1 ) k H k x k p τ p , k = 1 , , n + 1 .

Equation (14) expresses the slopes for each of the variables in the homotopic path γ , being essential for calculating each predictor vector. To determine the prediction vectors, we must consider the slopes of the homotopic path calculated with (14) and use Euler’s predictor algorithm. In this case, it is expressed as

(15) x p = x 0 + Δ p d x d p , τ p = τ 0 + Δ p d τ d p ,

where x p and τ p are the prediction values obtained by the Euler predictor, which will be found on the surface of the hypersphere of radius r . To calculate the horizontal advance Δ p , equations (15) are used, adding the equation of the hypersphere. In this way, we have

(16) Δ p = ± r 2 k = 1 n x k p 2 + τ p 2 .

The sign will depend on the direction of tracking to the right or left, that is, in a positive or negative direction. In this way, it is possible to calculate the prediction for x p and τ according to (15). To calculate the corresponding point p to the horizontal advance of the homotopic path, the equation will be used

(17) p = p 0 + Δ p .

Once the prediction vectors have been calculated, they now serve as the new initial values to begin the hyperspherical correction procedure, whose objective is to find the intersection between the hypersphere and the homotopic path, as shown in Figure 4. This procedure is carried out using the N-R method to solve the homotopic equations of the system, described by (3), and with the equation of the hypersphere (4). The N-R method, in two or more dimensions, requires the Jacobian matrix and a vector with the function evaluated at an initial point ( x p , τ p )

(18) ( x , τ ) ( k + 1 ) = ( x , τ ) ( k ) [ J ( H s ( x , τ ) ( k ) ) ] 1 H s ( x , τ ) ( k ) .

However, it is possible to identify in H 1 ( p , τ ) the roots found when the homotopic route path τ = 1 using (17). In fact, the determination of the roots occurs in that place; although, it is challenging for the intersection of the homotopic path to occur precisely at that point. Therefore, in the vicinity of this value, when the step transitions from a lower number ( τ j < 1 ) to a number greater than 1 ( τ j + 1 > 1 ) , the step prior to the find zero strategy (FZS, or refinement) of the roots is linear interpolation using the following formula:

(19) x ^ = x j + ( x j + 1 x j ) ( τ j + 1 τ j ) ( 1 τ j ) .

In order to generate the interpolated start points (ISPs) for FZS, the variables x j + 1 and τ j + 1 correspond to values when ( τ > 1 ) , and x j , τ j for ( τ < 1 ) , x represents each of the variables, and x ^ are the interpolated values for each variable by setting τ = 1 in the interpolation. The values found with (19) will be the ISP values for the N-R method in the FSZ procedure for τ = 1 , evaluated at x ^ .

(20) x ( k + 1 ) = x ( k ) [ J ( f ( x ( k ) ) ) ] 1 f ( x ( k ) ) .

3.3 Flowchart and pseudocode

Figure 5 illustrates the flowchart of the algorithm for homotopy with hyperspherical tracking. Each of the key points of the algorithm is summarized as follows:

  1. The formulation of the system of homotopy equations is initiated by including the hypersphere equation.

  2. Set the number of iterations for hyperspheres and initialize the counter.

  3. Establish the nonsquare matrix H . Calculate the determinants and multiply them by ( 1 ) k . Compute the horizontal advance Δ p .

  4. Proceed to calculate the predictor vector. The symbol refers to the chosen predictor. In this work, the Euler predictor, equation (15).

  5. Apply N-R for homotopic correction. Save in a list H 1 ( x , τ ) and H 1 ( p , τ ( p ) ) . Update the counter.

  6. Calculate τ j and τ j + 1 for the FSZ strategy.

  7. Interpolate and apply N-R. Save the solutions in a list.

  8. If the number of hypersphere iterations has been reached, print the solutions and the homotopic tracking.

Note in the flow diagram that a predictor other than Euler’s can be chosen, such as the Heun, Runge-Kutta, and fourth-order Runge-Kutta predictors, and Ralston [106,107]. In this work, the Euler predictor was implemented.

Figure 5 
                  Flow chart for homotopy with hyperspherical tracking.
Figure 5

Flow chart for homotopy with hyperspherical tracking.

Algorithm 1 shows the pseudocode of the flowchart in Figure 5, where each necessary step for coding in programming languages such as Maple, Matlab, C++, Python, among others, is detailed. The presented pseudocode describes the implementation of the homotopy algorithm for hyperspherical tracking in a simplified manner, highlighting the main steps for its coding.

Algorithm 2 shows the pseudocode for the corrector_N-R subroutine of the flowchart in Figure 5. The subroutine receives as arguments the symbolic J ( H s ( x,τ ) ) , the vector x p , τ . Likewise, Algorithm 2 shows FSZ strategy N-R, FSZ_N-R. This subroutine receives, x ^ , and symbolic J ( f ( x k ) ) .

3.3
3.3

4 Case studies

The mesh analyses presented in the case studies are fundamental in the analysis of circuits within electrical and electronic engineering. Both mesh and node analyses are important in the analysis of electrical circuits [28]. In electrical engineering, the relevance of mesh analysis is vast, being applied, for example, in three-phase circuits and those involving motors [28]. Furthermore, in the analysis of AC circuits, the frequencies of the voltage sources generate reactances in components such as capacitors and inductors. With passive devices such as capacitors and inductors, we obtain impedance, whose unit of measurement is the ohm [15,28,31].

Likewise, these analyses allow the application of various theorems, such as Thevenin’s theorem and the superposition theorem, among others, to simplify circuit analysis. In electronics, mesh analysis is essential for studying circuits with transistors, in both direct current and small signal analyses [108,109]. Transistors are the foundation of electronics and the design of analog integrated circuits, as well as integrated circuits of the transistor-transistor logic (TTL) and complementary metal-oxide-semiconductor (CMOS) families [108110], microcontrollers [111], microprocessors [112], and others.

This section presents five case studies. The first involves a system of two algebraic simultaneous equations with real coefficients. The second and third case studies focus on mesh analysis of AC electrical circuits. In the fourth case study, the voltages of a small power system connected to a motor are determined. The final case study presents a mesh analysis of AC electrical circuits incorporating nonlinear impedances. In all case studies, we use the notation “j” for the imaginary part of the complex numbers, and the radius of the hyperspheres is r = 0.003 . The convention for representing AC is to use the lowercase letter “i” in bold, as it denotes a complex quantity [15,28].

4.1 Case study 1

Solve the system of equations given by

(21) x 3 y + 1 = 0 , 1 2 x y 1 = 0 .

This system of equations can be solved by algebraic or numerical methods. To solve it, in this work, two SPs were proposed to solve the system (21), P 1 ( x , y ) = { 0.0001 + 0.0 j , 0.1 + j } and P 2 = { 0.59 + 0.59 j , 1 + 0.01 j } . Figure 6 displays the graphs of the system (21), where the real root is clearly observed.

Figure 6 
                  Plots for equations for case study 1.
Figure 6

Plots for equations for case study 1.

Table 1 summarizes the roots found with SPs used, the direction of the tracking, right (R) and left (L), the hypersphere iterations used, and CPU time in seconds. For example, the solutions S 1 and S 2 were found using SP P 1 , the first one with tracking R, and the second one with tracking L, using 49 hyperspheres. It is noteworthy that the magnitude of the imaginary part of ( S 1 ) is small, on the order of 1 × 1 0 106 , rendering it effectively real. Solution S 3 was found using SP P 2 with L tracking using 17 hyperspheres, and S 2 was found again during R tracking using 29 hyperspheres. The hypersphere radius used was equal to 0.003. The convention employed in Table 1 of this case study remains consistent across other case studies where the revealed solutions are displayed.

Table 1

Homotopic tracking for case study 1

Figure 7 displays the homotopic tracking to find S 1 and S 2 for SP P 1 , and Figure 8 shows the homotopic tracking for SP P 2 to find S 2 and S 3 . Note that the trace for x presents a greater curlicue due to it being the variable with the greatest exponent. Observe that the axes of the figures are τ , imaginary part, and real part, respectively. At the intersection of the tracking trajectories with the plane τ = 1 are the solutions for x , y , marked with two small spheres.

Figure 7 
                  Hyperspheric homotopic tracking with 
                        
                           
                           
                              
                                 
                                    P
                                 
                                 
                                    1
                                 
                              
                           
                           {{\bf{P}}}_{{\bf{1}}}
                        
                      for case study 1. (a) Right tracking and (b) left tracking.
Figure 7

Hyperspheric homotopic tracking with P 1 for case study 1. (a) Right tracking and (b) left tracking.

Figure 8 
                  Hyperspheric homotopic tracking with 
                        
                           
                           
                              
                                 
                                    P
                                 
                                 
                                    2
                                 
                              
                           
                           {{\bf{P}}}_{{\bf{2}}}
                        
                      for case study 1. (a) Right tracking and (b) left tracking.
Figure 8

Hyperspheric homotopic tracking with P 2 for case study 1. (a) Right tracking and (b) left tracking.

4.2 Case study 2

Find the mesh currents of the polarization point Q for the electric circuit in Figure 9.

Figure 9 
                  Electrical circuit with two meshes for case study 2.
Figure 9

Electrical circuit with two meshes for case study 2.

This circuit was solved in [15] using Cramer’s rule. The loop equations of the circuit in Figure 9 were formulated with Kirchhoff’s law of charge conservation [15]. Therefore, the system of simultaneous equations of the circuit is given by

(22) ( 10 5 j ) i 1 + ( 5 j ) i 2 50 = 0 , ( 5 j ) i 1 + ( 3 j ) i 2 = 0 .

Unlike the system of equations presented in case study 1, the coefficients of (22) in this case study are complex. Applying homotopy with hyperspherical tracking to solve the system of simultaneous equations, a SP P ( i 1 , i 2 ) = { 0.0001 + 0.01 j , 0.1 + j } [A] is employed. Table 2 shows the solutions and homotopic tracking carried out. In both tracking directions, the number of hyperspheres was 87. Similar times were obtained in both path tracking processes.

Table 2

Homotopic tracking for case study 2

Figure 10 displays the homotopic tracking in both directions to determine solutions. Observe the linear behavior of the tracking paths for i 1 and i 2 , which is characteristic of a linear system. In the plane τ = 1 , the solution is indicated with small spheres.

Figure 10 
                  Hyperspheric homotopic tracking with 
                        
                           
                           
                              P
                           
                           {\bf{P}}
                        
                      for case study 2. (a) Right tracking and (b) left tracking.
Figure 10

Hyperspheric homotopic tracking with P for case study 2. (a) Right tracking and (b) left tracking.

4.3 Case study 3

Find the mesh currents of the polarization point Q for the circuit in Figure 11.

Figure 11 
                  Electrical circuit with three meshes for case study 3.
Figure 11

Electrical circuit with three meshes for case study 3.

Likewise, this circuit was solved in [15] using the rule of Cramer. The loop equations for the circuit in Figure 11 are given by

(23) ( 5 + 5 j ) i 1 ( 5 j ) i 2 30 = 0 , 5 j i 1 + ( 8 + 8 j ) i 2 6 i 3 = 0 , 6 i 2 + 10 i 3 + 20 = 0 .

Using the SP P ( i 1 , i 2 , i 3 ) = { j , 1 + 0.02 j , 1 } [A] for the homotopy, the mesh currents obtained are shown in Table 3.

Table 3

Homotopic tracking for case study 3

The behavior of the homotopic path exhibited by i 1 , i 2 , i 3 is linear because the system of equations is linear. Figure 12 shows the homotopic tracking in both directions.

Figure 12 
                  Hyperspheric homotopic tracking with 
                        
                           
                           
                              P
                           
                           {\bf{P}}
                        
                      for case study 3. (a) Right tracking and (b) left tracking.
Figure 12

Hyperspheric homotopic tracking with P for case study 3. (a) Right tracking and (b) left tracking.

4.4 Case study 4

In [28], the problem for the single-line diagram of a small power system shown in Figure 13 (a) was solved. The equivalent reactance diagram, with reactances, is shown in Figure 13 (b). A generator with an electromotive force equal to 125 [V] per unit is connected through a transformer to high-voltage node 3, while a motor with internal voltage equal to 0.481 0.481 j [V] is connected to node 4. Figure 14 illustrates the circuit with the transformation of voltage generators into current sources, presenting the nodes to establish the nodal analysis equations. Determine all node voltages.

Figure 13 
                  The single-line diagram of a small power system. (a) Single-line diagram of the four-bus system and (b) reactance diagram.
Figure 13

The single-line diagram of a small power system. (a) Single-line diagram of the four-bus system and (b) reactance diagram.

Figure 14 
                  Electrical circuit with three meshes for case study 4.
Figure 14

Electrical circuit with three meshes for case study 4.

Applying Kirchhoff’s voltage law, the node voltage equations for the equivalent of Figure 14 is given by (24)

(24) ( 14.5 j ) v 1 + ( 8 j ) v 2 + ( 4 j ) v 3 + ( 2.5 j ) v 4 = 0 , ( 8 j ) v 1 ( 17 j ) v 2 + ( 4 j ) v 3 + ( 5 j ) v 4 = 0 , ( 4 j ) v 1 + ( 4 j ) v 2 ( 8.8 ) v 3 j = 0 , ( 2.5 j ) v 1 + ( 5 j ) v 2 ( 8.3 ) v 3 0.481 ( 0.481 ) j = 0 .

Using the SP P ( v 1 , v 2 , v 3 , v 4 ) = { 0.5 j , 0.1 , 0.21 , 0.4 j } [V] in the homotopy to solve (24), the node voltages are shown in Table 4. The operating point was obtained in both directions of hyperspherical tracking, utilizing 184 hyperspheres. Similar times were obtained in both path tracking processes.

Table 4

Homotopic tracking for case study 4

Figures 15 (a) and (b) display the hyperspherical homotopic tracking with linear behavior for each of the voltages of the system equation (24), showing the unique solution in the plane τ = 1 .

Figure 15 
                  Hyperspheric homotopic tracking with 
                        
                           
                           
                              P
                           
                           {\bf{P}}
                        
                      for case study 4. (a) Right tracking and (b) left tracking.
Figure 15

Hyperspheric homotopic tracking with P for case study 4. (a) Right tracking and (b) left tracking.

4.5 Case study 5

In [113], an electric circuit with a loop containing nonlinear elements was solved. In this case study, we will find the mesh currents for the electric circuit in Figure 16.

Figure 16 
                  Electrical circuit with three meshes for case study 5.
Figure 16

Electrical circuit with three meshes for case study 5.

The formulas for nonlinear impedances [113,114] are given by

(25) X C = 1 ω C 1 + 3 γ i 2 2 ω 2 , X L = ω L 1 3 β i 2 2 ,

where ω is the angular frequency in radians per second [rad/s], C is the capacitance in Farads [F], L is the inductance in Henrys [H], and the reactances X C and X L are given in ( Ω ) (Ohms). The coefficients γ and β are measures of nonlinearity and are positive in the absence of direct current polarization. Details on how these coefficients are determined are presented in [114].

For the mesh analysis of this problem, with nonlinear impedances in this case study, the following values are considered: γ = 1.5 [ rad 2 C 2 ], β = 1.5 [ A 2 ], ω = 62.831 [rad/s], L 1 = 0.03 [H], C 1 = 0.0001 F, L 2 = 0.05 H, and C 2 = 0.0003 [F]. Applying Kirchhoff’s current law for each of the circuit meshes, three equations are obtained. By substituting the values for each one of the elements of the circuit and simplifying algebraically, we have a NAEs composed of three simultaneous equations.

(26) ( 1.88495562 j ( 1 2.250 i 1 2 ) i 1 159.1549407 j ( 1 + 0.0005699316411 i 1 2 ) + 10 ) i 1 10 i 2 5 = 0 , 10 i 1 + ( 24 + 18.9 j ( 1 2.250 i 2 2 ) ) i 2 5 i 3 = 0 , 5 i 2 + ( 11 53.051646913 j ( 1 + 0.00056993164114 i 3 2 ) ) i 3 + 3 = 0 .

To determine the operating points of the nonlinear circuit, two different SPs were proposed: P 1 = { 5 , 0.05 j , 0.05 j } [A] and P 2 = { 10 j , 0.05 j , 0.05 j } [A]. By using these SPs in the routine developed in this work, six operating points were found. Table 5 illustrates the six solutions of (26), along with the hyperspheric tracking. It is noteworthy that the operating point of the electric circuit Q 1 was encountered twice during this tracking, specifically in hyperspheres number 40 and 287, respectively.

Table 5

Homotopic tracking for case study 5

Figure 17 displays the homotopic tracking for the SP P 1 . In Figure 17(a), the homotopic tracking is presented to the right, where four operating points were found. In this manner, operating point Q 1 is identified by the red sphere, Q 2 by the blue sphere, Q 3 by the gold sphere, and Q 4 by the black sphere, all located in the plane τ = 1 . It is important to note that hypersphere iterations 40 and 287 (for Q 1 ) overlap because they represent the same root. Likewise, it should be observed that the homotopic path for i 2 and i 3 tends to form a loop since the trajectory does not change significantly. Figure 17(b) shows the tracking L, where only the operating point Q 5 was found.

Figure 17 
                  Hyperspheric homotopic tracking with 
                        
                           
                           
                              
                                 
                                    P
                                 
                                 
                                    1
                                 
                              
                           
                           {{\bf{P}}}_{1}
                        
                      for case study 5. (a) Right tracking and (b) left tracking.
Figure 17

Hyperspheric homotopic tracking with P 1 for case study 5. (a) Right tracking and (b) left tracking.

Figure 18 presents the homotopic tracking on the left with SP P 2 to find Q 6 . In general, in this case study, the homotopic tracking for all electrical intensities showed a nonlinear behavior.

Figure 18 
                  Hyperspheric homotopic tracking with 
                        
                           
                           
                              
                                 
                                    P
                                 
                                 
                                    2
                                 
                              
                           
                           {{\bf{P}}}_{2}
                        
                      for case study 5.
Figure 18

Hyperspheric homotopic tracking with P 2 for case study 5.

5 Numerical simulation and discussion

The computer utilized for the development of this work was equipped with an Intel Core i7-7700U CPU @ 3.60 GHz with 4 cores and 8 threads, an NVIDIA Corporation GP107 [GeForce GTX 1050 Ti] graphics card, and 64 GB of RAM, operating on Linux Ubuntu version 20.04.6 LTS. For the implementation of the homotopy algorithm with hyperspherical tracking, Maple 2019 software was employed. In addition, in the simulations, the built-in Maple command, fsolve, was used to validate the results.

The analysis in this article focuses on identifying the complex roots of systems of algebraic equations, exploring the SDs found in the real and imaginary parts of each solution, as well as the iterations used. Tables 6, 7, 8, 9, 10 present the analysis of the SDs for the solutions of the case studies discussed in this article, utilizing the SPs provided by the linear interpolation of the homotopy algorithm. The criterion was to achieve more than ten SDs of accuracy in the solutions for the case studies. To validate the SD analysis, the SPs obtained from the interpolation (ISPs) were used in Maple’s fsolve function with the respective system of equations for all case studies to obtain the exact solution with 15 SDs.

Table 6

Significant digits in the solutions for case study 1

Table 7

Significant digits in the solutions for case study 2

Table 8

Significant digits in the solutions for case study 3

Table 9

Significant digits in the solutions for case study 4

The bold values are the interpolated values for each variable, as defined in equation (19).

Table 10

Significant digits in the solutions for case study 5

For case study 1, the real parts of all the solutions had 15 digits of accuracy. For S 1 , the imaginary part of its roots was of the order of 1 × 1 0 388 , so it was considered entirely a real solution. The root S 2 was obtained with the two proposed SPs; it was the solution that required more iterations in its refinement. Similarly, when using P 2 , the imaginary part of y in S 2 had 14 SDs of accuracy. Moreover, the solution S 3 was the root that required the fewest iterations in its refinement.

For case study 2, the analysis of SD was precise. The algebraically obtained solutions for case study 2 are i 1 = 15 4 + 2 5 j [A], i 2 = 2 + 4 j [A], while Table 2 presents their decimal equivalents. Table 7 demonstrates that accuracy was achieved in the first iteration in both directions of hyperspherical tracking. It is noteworthy that the SPs obtained through interpolation showed high accuracy, with some being identical to the refined solution values.

For case study 3, Table 8 shows high accuracy in the SD of all the electric currents obtained in the first refinement iteration for the electric circuit depicted in Figure 11. Only 14 SDs were achieved in the imaginary part of i 3 in the homotopic tracking to the left.

For case study 4, the SDs obtained in the nodal voltages v 1 , , v 4 using Maple’s fsolve function did not reach 15 because the exact numerical solution was achieved. Table 9 illustrates this fact. For the root v 4 , seven SDs were obtained in their real part compared to the eight of the exact numerical solution. This result is highlighted in bold font within the table. It is noteworthy that in the left homotopic trace, only one iteration was required for refinement.

Table 10 likewise demonstrates high accuracy in the SDs of the operating points for the electric circuit in case study 5. There is at least 14 SD in the real parts of the solutions, while the imaginary parts achieve 13 SD, as seen with i 3 for Q 3 . The number of iterations required to refine the solutions Q 1 , , Q 4 ranged between 2 and 3. The ISPs provided by the interpolation also exhibited high accuracy. Nonetheless, the repetition of Q 1 necessitated 11 iterations for its refinement, ultimately achieving 15 SD of accuracy in i 1 , i 2 , i 3 .

It is possible to determine the relative errors of all solutions presented in Tables 610 by employing the formula that associates the relative error with SDs [115,116], given by

(27) SD = log 10 v A v E v E ,

where SD is the SDs, v A is the approximate (measured) value, and v E is the exact value.

Continuation homotopy with hyperspherical tracking can also experience numerical errors, similar to any numerical algorithm, such as round-off errors due to the finite precision of floating-point arithmetic in computers, numerical errors caused by nearly identical equations, convergence failure in N-R, or numerical instabilities [106,107]. These errors can cause the algorithm to lose homotopy path tracking. Furthermore, a study on the sensitivity of the algorithm indicates that it is generally more sensitive to the complexity of f ( x ) than to its dimension, as reported in [117]. Table 11 presents a comparison of the computing times for [ J ( H s ( x , τ ) ) ] 1 and H 1 ( x , τ ) across all the case studies discussed in this work. The table includes the dimension n + 1 because the homotopy parameter τ was added to the system of equations to generate the homotopy. For instance, case studies 3 and 5 involve systems with three equations; however, it can be observed that the evaluation time for [ J ( H s ( x , τ ) ) ] 1 is greater in the latter case. In addition, it is important to note that the CPU times for cases 4 and 5 are nearly identical. The comparison of computing times for all case studies was carried out using 300 hyperspheres with a radius of r = 0.003 .

Table 11

CPU time for case [ J ( H s ( x , τ ) ) ] 1 and H 1 ( x , τ )

Table 12 presents a comparison of the computation times using differential path tracking and hyperspherical path tracking [34,39] for all the case studies analyzed in this work. The computation times obtained with hyperspherical path tracking are consistently shorter compared to those obtained with differential path tracking across all case studies. In addition, the number of iterations required in this study is lower when using hyperspherical path tracking compared to differential path tracking. It is noteworthy that in all case studies, all solutions were successfully found using hyperspherical path tracking, while differential path tracking only managed to identify some of them.

Table 12

Comparison for solutions using hyperspherical path tracking vs differential path tracking

Furthermore, it is important to note that the inclusion of complex numbers in hyperspherical tracking significantly increases computation times, in contrast to those reported in the literature when real numbers are employed.

Table 13 presents a comparison of solutions obtained using homotopy with hyperspherical path tracking and the Super Ostrowski homotopy continuation method (SOHCM) [118,119]. The table displays the solutions found with Super-Ostrowski HCM in the path tracking, except in case study 4. In general, the CPU times with SOHCM are shorter than this work. However, in some cases, the SP leads to finding a different root. For instance, in this work, solution S 1 was obtained with P 1 using the hyperspherical path tracking, whereas S 3 was found with SOHCM after 280 iterations.

Table 13

Comparison for solutions using hyperspherical path tracking vs. super-Ostrowski HCM

For case study 5, using SOHCM, Q 5 was obtained with P 1 , while employing homotopy with hyperspherical path tracking, the roots were identified in the order Q 1 , Q 2 , Q 3 , Q 4 , and Q 5 . It is important to note that the iterations employed in SOHCM generally exceed 70 to achieve 15 SDs. However, in case study 5, 200 iterations were used to obtain 2 SD.

Likewise, the computational time increases as the complexity of the equation systems grows. Case study 5 requires a computational time of 4.005 [s] due to the system being nonlinear. In contrast, case study 4 does not have a solution. Furthermore, in [118], the structure of the homotopy includes a squared homotopic parameter, and the path tracking does not use any predictor, such as Euler, among others. However, it is possible to obtain τ with higher exponents if G ( x , τ ) , as shown in equation (14) in [118], becomes more complicated. Actually, instead of that, SOHCM uses a simple proportional predictor. The algorithm requires the number of iterations to be set up a priori. The SOHCM algorithm uses the number of iterations as the denominator, with a counter in the numerator for this fraction. The solution is considered found when the numerator equals the denominator, resulting in τ = 1 . For this reason, Table 13 only shows right path tracking in SOHCM. Likewise, SOHCM can only find one solution when τ = 1 because the algorithm terminates at τ = 1 . In this research, equation (14) from [118] was implemented with hyperspherical path tracking, using a hypersphere radius of 0.003; however, divergence was encountered. Nonlinear terms in the structure of any homotopy cause numerical instabilities in path tracking. The results presented in this study were obtained using an Euler predictor with path tracking in both directions.

The path tracking behavior in Figures 10, 12, and 15 for case studies 2, 3 and 5, respectively, exhibits a linear trend because the systems of equations in theses case studies are linear. In contrast, Figures 7 and 8 in case study 1 and Figures 17 and 18 in case study 5 display nonlinear behavior. As the nonlinearity in the system of equations increases, the homotopic tracking demonstrates a more pronounced nonlinear behavior, as seen in Figure 17. However, Figures 19 and 20 show the graph homotopic parameter τ versus advance parameter p used in hyperspherical path tracking. These figures highlight the advantage of utilizing the advance parameter p , as it facilitates the visualization of all solutions via homotopic path tracking. In addition, the graphs depict the system’s behavior. For example, Figures 19 (b)–(d) exhibit linear behavior because the systems of equations are linear for case studies 2, 3, and 4. In contrast, Figures 19 (a) and 20 for case studies 1 and 5, respectively, reveal nonlinear homotopic path tracking, attributed to the nonlinearity of the systems of equations in these cases. It is important to note that Figures 19 and 20 are in good agreement with the figures for hyperspheric homotopic tracking shown in each case study. For instance, Figure 20 (a) shows significant nonlinear behavior in homotopic path tracking for six solutions (five solutions in right tracking, and one solution in left tracking), while Figure 17 (a) depicts nonlinear behavior for each solution.

Figure 19 
               Homotopic advance 
                     
                        
                        
                           p
                        
                        p
                     
                   for case studies 1–4. (a) Homotopic tracking 
                     
                        
                        
                           p
                        
                        p
                     
                   vs 
                     
                        
                        
                           τ
                        
                        \tau 
                     
                   for case study 1, (b) homotopic tracking 
                     
                        
                        
                           p
                        
                        p
                     
                   vs 
                     
                        
                        
                           τ
                        
                        \tau 
                     
                   for case study 2, (c) homotopic tracking 
                     
                        
                        
                           p
                        
                        p
                     
                   vs 
                     
                        
                        
                           τ
                        
                        \tau 
                     
                   for case study 3, and (d) homotopic tracking 
                     
                        
                        
                           p
                        
                        p
                     
                   vs 
                     
                        
                        
                           τ
                        
                        \tau 
                     
                   for case study 4.
Figure 19

Homotopic advance p for case studies 1–4. (a) Homotopic tracking p vs τ for case study 1, (b) homotopic tracking p vs τ for case study 2, (c) homotopic tracking p vs τ for case study 3, and (d) homotopic tracking p vs τ for case study 4.

Figure 20 
               Homotopic advance 
                     
                        
                        
                           p
                        
                        p
                     
                   for case study 5. (a) Homotopic tracking 
                     
                        
                        
                           p
                        
                        p
                     
                   vs 
                     
                        
                        
                           τ
                        
                        \tau 
                     
                   for 
                     
                        
                        
                           
                              
                                 
                                    
                                    P
                                    
                                 
                              
                              
                                 1
                              
                           
                        
                        {\text{P}}_{1}
                     
                   and (b) left tracking 
                     
                        
                        
                           p
                        
                        p
                     
                   vs 
                     
                        
                        
                           τ
                        
                        \tau 
                     
                   for 
                     
                        
                        
                           
                              
                                 
                                    
                                    P
                                    
                                 
                              
                              
                                 2
                              
                           
                        
                        {\text{P}}_{2}
                     
                  .
Figure 20

Homotopic advance p for case study 5. (a) Homotopic tracking p vs τ for P 1 and (b) left tracking p vs τ for P 2 .

For SP P 2 in study case 5, Figure 20 (b) shows only left path tracking with its solution. It is important to note that the right tracking is not shown because a divergence occurred, and consequently, the solution cannot be determined.

Figure 21 shows the attractor basins for case study 1. Both figures were generated using the computer program Maple 21, with ε = 1 × 1 0 8 . The attractor basin for N-R is defined as A = [ 2 , 2 ] × [ 2 , 2 ] , and the attractor basin for Homotopy is A = [ 2.5 , 2.5 ] × [ 2.5 , 2.5 ] . The step size used for all iterations is k = 0.003 . The parameter ε indicates the accuracy of the calculated root, and A represents the area within which we draw the basins of attraction. The iterations and the colors used to draw the attraction basins in Figure 21 are listed in Table 14.

Figure 21 
               Basins of attractions for case study 1. (a) Newton basin and (b) homotopy basin.
Figure 21

Basins of attractions for case study 1. (a) Newton basin and (b) homotopy basin.

Table 14

Iterations in attraction basins for case study 1, Figure 21

In the generated images, each point in the complex region that converges to a specific root is colored to match that root, with darker shades indicating a higher number of iterations and lighter shades fewer iterations. Points that do not converge within the iteration limit are colored gray. Figure 21 (a) shows the basins for a system with roots S 1 , S 2 , and S 3 using the N-R method, taking 50 min to generate, while Figure 21 required 29 h and 14 min. Convergence speed and dynamical behavior are reflected in the color shades and variations, respectively, with uniform regions indicating consistent iteration counts, giving a contour-like appearance.

In the figure, stability occurs when the initial point is located within a basin of attraction that contains a root. Conversely, instability arises when the SP leads to a trajectory that exits the colored zone and diverges to a distant region, potentially converging to a different root. Moreover, the basin boundaries generated by the N-R method intersect, forming the Julia set, a complex boundary where chaotic behavior manifests. It is important to keep in mind that a complex analytical function divides the plane into two distinct regions: the Fatou set, a stable area characterized by relatively simple dynamics, and the Julia set [14,120122], which forms the chaotic boundary of the Fatou set [120,122].

Figure 21 (b) displays the attraction basins for the homotopy method with hyperspherical path tracking for x , exhibiting noticeable symmetry at first glance. Since the basins for y are similar, only those for x are shown. However, significant differences can be observed when compared to the basins in Figure 21 (a). Specifically, the basin boundaries do not extend continuously from a common point. In addition, four regions appear where points fail to converge within the iteration limit; these regions are shown in gray, indicating divergence. Figure 21 (b) also exhibits chaotic behavior similar to Newton’s basins when the initial point lies within the fractal structure. For instance, near the boundaries of the basins for S 1 and S 3 , some SPs converge to S 2 in olive green. However, the basins also suggest a high degree of stability, as the areas in light colors are significantly larger compared to those in dark colors (up to 24 iterations were needed). Consequently, convergence occurs toward the root located within the respective basin.

6 Conclusions

This study introduced a comprehensive exploration of the HCM combined with hyperspherical tracking to identify complex roots in systems of equations with applications to electrical circuits. The research conducted spanned a variety of case studies, ranging from systems of algebraic equations to complex electrical circuits with nonlinear elements. The findings of this study emphasize the efficacy and precision of the HCM with hyperspherical tracking in solving complex systems of equations.

The case studies presented illustrate the method’s capability to accurately determine complex roots, even in scenarios involving nonlinear elements and imaginary coefficients. The precision of the solutions was validated against existing literature and through the use of Maple’s fsolve function, demonstrating the method’s reliability and effectiveness. The analysis of SDs further underscored the method’s high level of accuracy, with most solutions achieving close to or exactly the desired number of SDs.

In addition, the study of basins obtained through homotopy requires an in-depth analysis of their behavior, as well as an examination of the potential existence of Julia sets or other related structures. The basin of attraction obtained through continuation homotopy extends the analysis of these basins while also broadening the study of continuation methods using this technique. On the other hand, it remains to study the causes of divergence of the homotopy in the basin of attraction.

Notably, the study also highlighted the method’s versatility in handling both linear and nonlinear systems, as well as its ability to find multiple operating points in electrical circuits. This versatility makes the HCM with hyperspherical tracking a valuable tool for researchers and practitioners in various fields, including electrical engineering, where complex systems of equations frequently arise.

Furthermore, the research presented contributes to the ongoing development of numerical methods for solving systems of equations. By demonstrating the practical applications of the HCM with hyperspherical tracking, this study adds to the body of knowledge and provides a foundation for future research in this area. The study proposes several avenues for future exploration. First, the homotopy routine with hyperspherical path tracking will be implemented in the C + + language using symbolic computation [123] to reduce computational time for tasks such as Jacobian computation, hyperspherical path tracking. Second, the application of homotopy with hyperspherical tracking is suggested to find complex roots of characteristic polynomials in systems theory. This method could aid in analyzing stability, frequency response, and other crucial aspects [29]. Third, there is a proposal to develop a routine for microcontrollers, aimed at facilitating educational research focused on implementing control systems that require real-time monitoring of stability, frequency response, and other parameters. Moreover, the study suggests investigating the tracking of complex homotopic trajectories applied to robotics, particularly in navigating obstacles in the complex plane. This research could build upon previous works in the field [60,105]. Finally, it is proposed to explore in the future the possibility of using nonlinear interpolation [124126] in FZS to reduce the number of iterations.

In conclusion, the HCM with hyperspherical tracking has proven to be a robust and precise approach for identifying complex roots in systems of equations. Its successful application to a diverse range of case studies underscores its potential as a numerical tool for solving complex problems in science and engineering. Future research could explore the method’s applicability to other complex systems and further optimize its implementation for enhanced efficiency and accuracy.

Acknowledgements

The first author is currently research fellow at Tecnológico Nacional de México en Celaya. He expresses his gratitude for the support received.

  1. Funding information: This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

  2. Author contributions: MASH, HJI, HVL: conceived and designed the experiments; performed the experiments; Analyzed and interpreted the data; contributed reagents, materials, analysis tools or data; wrote the article. MLQV, MLLG: performed the experiments; analyzed and interpreted the data; contributed reagents, materials, analysis tools, or data.

  3. Conflict of interest: Dr. Hector Vazquez-Leal is an Editor of the Open Mathematics journal and was not involved in the review and decision-making process of this article.

  4. Data availability statement: All data analyzed during this study are included in this published article.

References

[1] F. Rivero-Mendoza, Una introducción a los números complejos, Editorial Universidad de los Andes, Mérida, Venezuela, 2001.Search in Google Scholar

[2] H. Cardano, Artis magnae, sive de regulis algebraicis, liber unus, Franco Angeli, Milan, Italy, 2011. Search in Google Scholar

[3] C. R. Maluendas, Una breve historia imaginaria, Univ. Nac. Colomb. 1 (2019), no. 2, 1–9. https://editorial.konradlorenz.edu.co/2019/02/paskin-matematico.html. Search in Google Scholar

[4] R. Bombelli, LaAlgebra, Feltrinelli, Milan, Italy, 1966. Search in Google Scholar

[5] R. Descartes, La Géométrie, A. Herman, Librairie Scientifique, París, 1637. Search in Google Scholar

[6] B. Branner, Caspar Wessel. On the Analytical Representation of Direction: An Attempt Applied Chiefly to Solving Plane and Spherical Polygons, The Royal Danish Academy of Sciences and Letters, Copenhagen, Denmark, 1999. Search in Google Scholar

[7] J. R. Argand, Essai sur une manière de représenter les quantités imaginaires dans, les constructions géométriques, Gauthier-Villars, Imprimeur-Libraire, Paris, France, 1874. Search in Google Scholar

[8] L. Euleri, Introductio in analysin infinitorum, Tomus primus, Lausannae: apud Marcum-Michaelem Bousquet & Socies, Lausanne, Switzerland, 1748, http://hdl.handle.net/10481/31374.Search in Google Scholar

[9] L. Euler, Elements of Algebra, Springer-Verlag, WI, USA, 1972. Search in Google Scholar

[10] M. Hazewinkel, Encyclopaedia of Mathematics, Springer Science & Business Media, New York, USA, 1994. 10.1007/978-94-009-5983-5Search in Google Scholar

[11] S. Gutiérrez, Hamilton: La liberación del álgebra, Suma 49 (2005), no. 49, 95–99. https://revistasuma.fespm.es/sites/revistasuma.fespm.es/IMG/pdf/49/095-099.pdf. Search in Google Scholar

[12] J. M. Sánchez, Historias de Matemáticas: Hamilton y el descubrimiento de los cuaterniones, Pens. Mat. 1 (2011), no. 27. https://dialnet.unirioja.es/servlet/articulo?codigo=3744306. Search in Google Scholar

[13] C. P. Steinmetz, Complex Quatities and their use in electrical engineering, Proc. Int. Electr. Congr. Am. Inst. Electr. Eng. August 21st to 25th American Institute of Electrical Engineers, USA, 1893, pp. 33–75. Search in Google Scholar

[14] A. D. Wunsch, Variable Compleja con Aplicaciones, Pearson Educación, México, 1997. Search in Google Scholar

[15] J. A. Edminister, Circuitos Eléctricos, Serie Schaum, McGraw-Hill, México, 1982. Search in Google Scholar

[16] B. Mandelbrot, The Fractal Geometry of Nature, Free-man and Co., Mexico, 1977. Search in Google Scholar

[17] C. Monroy-Olivares, Curvas Fractales, Alfaomega, México, 2002. Search in Google Scholar

[18] M. J. Neve, Reinterpreting the Smith chart using conformal geometric algebra, IEEE Access 11 (2023), 138827–138838, DOI: https://doi.org/10.1109/ACCESS.2023.3340143. 10.1109/ACCESS.2023.3340143Search in Google Scholar

[19] A. N. Kani, Digital Signal Processing, CBS Publishers and Distributors, New Delhi, India, 2022. Search in Google Scholar

[20] G. A. Muñoz-Fernández and J. B. Seaone-Sepulveda, Fundamentos y Problemas Resueltos de Teoría Cualitativa de Ecuaciones Diferenciales, Ediciones Paraninfo, 2017. Search in Google Scholar

[21] A. L. Granados, Fractal techniques to measure the numerical instability of optimization methods, Comput. Mech. 15 (1995), no. 9, 369–374. Search in Google Scholar

[22] J. M. Gutiérrez, L. J. Hernández-Paricio, M. Marañón-Grandes, and M. T. Rivas-Rodríguez, Influence of the multiplicity of the roots on the basins of attraction of Newton’s method, Numer. Algorithms 66 (2014), 431–455. 10.1007/s11075-013-9742-7Search in Google Scholar

[23] J. M. Gutiérrez, Numerical properties of different root-finding algorithms obtained for approximating continuous Newton’s method, Algorithms 8 (2015), no. 4, 1210–1218, DOI: https://doi.org/10.3390/a8041210. 10.3390/a8041210Search in Google Scholar

[24] L. J. Hernández-Paricio, Bivariate Newton-Raphson method and toroidal attraction basins, Numer. Algorithms 71 (2016), 349–381, DOI: https://doi.org/10.1007/s11075-015-9996-3. 10.1007/s11075-015-9996-3Search in Google Scholar

[25] M. Sandoval-Hernandez, H. Vazquez-Leal, L. Hernandez-Martinez, U. A. Filobello-Nino, V. M. Jimenez-Fernandez, A. L. Herrera-May, et al., Approximation of Fresnel integrals with applications to diffraction problems, Math. Probl. Eng. 2018 (2018), 4031793, DOI: https://doi.org/10.1155/2018/4031793. 10.1155/2018/4031793Search in Google Scholar

[26] M. Maggiore, A Modern Introduction to Classical Electrodynamics, Oxford University Press, Oxford, UK, 2023. 10.1093/oso/9780192867421.001.0001Search in Google Scholar

[27] A. Cardama-Aznar, L. Jofre-Roca, J. M. Rius-Casals, J. Romeu-Robert, S. Blanch-Boris, and M. Ferrando-Bataller, Antenas, Univ. Politècnica de Catalunya, Alfaomega, Espannna, 2004. Search in Google Scholar

[28] J. J. Grainger and W. D. Stevenson Jr, Power System Analysis, McGraw-Hill Series in Electrical and Computer Engineering, Singapore, 1994. Search in Google Scholar

[29] K. Ogata, Ingeniería de Control Moderna, Quinta edición, Pearson Educación, México, 2010. Search in Google Scholar

[30] T. Deliyannis, Y. Sun, and J. K. Fidler, Continuous-time Active Filter Design, CRC Press, Boca Raton, Florida, USA, 2020. 10.1201/9781439821879Search in Google Scholar

[31] L. P. Huelsman, Active and Passive Analog Filter Design: An Introduction, McGraw Hill, New York, 1993. Search in Google Scholar

[32] J. V. Uspensky, Teoría de Ecuaciones, Limusa, México, 1987. Search in Google Scholar

[33] M. A. Jenkins and J. F. Traub, Algorithm 419: zeros of a complex polynomial [C2], Commun. ACM, 15 (1972), no. 2, 97–99, DOI: https://doi.org/10.1145/361254.361262. 10.1145/361254.361262Search in Google Scholar

[34] M. A. Sandoval-Hernandez, O. Alvarez-Gasca, A. D. Contreras-Hernandez, J. E. Pretelin-Canela, B. E. Palma-Grayeb, V. M. Jimenez-Fernandez, et al., Exploring the classic perturbation method for obtaining single and multiple solutions of nonlinear algebraic problems with application to microelectronic circuits, Int. J. Eng. Res. Technol. 8 (2019), no. 9, 636–645. https://www.ijert.org/exploring-the-classic-perturbation-method-for-obtaining-single-and-multiple-solutions-of-nonlinear-algebraic-problems-with-application-to-microelectronic-circuits. Search in Google Scholar

[35] F. Dubeau and C. Gnang, Fixed point and Newton’s methods in the complex plane, J. Complex Anal. 2018 (2018), 1–11, DOI: https://doi.org/10.1155/2018/7289092. 10.1155/2018/7289092Search in Google Scholar

[36] F. Çilingir, Singular perturbations arising in complex Newton’s method, Hacet. J. Math. Stat. 52 (2023), no. 5, 1–7, DOI: https://doi.org/10.15672/hujms.1132257. 10.15672/hujms.1132257Search in Google Scholar

[37] M. Sandoval-Hernandez, H. Vazquez-Leal, U. Filobello-Nino, E. De-Leo-Baquero, A. C. Bielma-Perez, J. C. Vichi-Mendoza, et al., The quadratic equation and its numerical roots, Int. J. Eng. Res. Technol. 10 (2021), no. 6, 301–305. https://www.ijert.org/research/the-quadratic-equation-and-its-numerical-roots-IJERTV10IS060100.pdf. Search in Google Scholar

[38] K. Gdawiec, I. K. Argyros, S. Qureshi, and A. Soomro, An optimal homotopy continuation method: Convergence and visual analysis, J. Comput. Sci. 74 (2023), 1–12, DOI: https://doi.org/10.1016/j.jocs.2023.102166. 10.1016/j.jocs.2023.102166Search in Google Scholar

[39] T. L. Wayburn and J. D. Seader, Homotopy continuation methods for computer-aided process design, Comput. Chem. Eng. 11 (1987), no. 1, 7–25, DOI: https://doi.org/10.1016/0098-1354(87)80002-9. 10.1016/0098-1354(87)80002-9Search in Google Scholar

[40] J. Farhang, J. D. Seader, and S. Khaleghi, Global solution approaches in equilibrium and stability analysis using homotopy continuation in the complex domain, Comput. Chem. Eng. 32 (2008), no. 10, 2333–2345, DOI: https://doi.org/10.1016/j.compchemeng.2007.12.001. 10.1016/j.compchemeng.2007.12.001Search in Google Scholar

[41] H. Vazquez-Leal, et al., Biparameter homotopy-based direct current simulation of multistable circuits, Br. J. Math. Comput. Sci. 2 (2012), no. 3, 137–150, https://journaljamcs.com/index.php/JAMCS/article/view/1240.10.9734/BJMCS/2012/1346Search in Google Scholar

[42] M. D. L. L. López-González, M. L. Quemada-Villagómez, M. G. Martínez-González, J. M. Oliveros-Muñoz, and H. Jiménez-Islas, A novel predictive homotopic path tracking algorithm to solve non-linear algebraic equations, Can. J. Chem. Eng. 101 (2023), no. 6, 3382–3408, DOI: https://doi.org/10.1002/cjce.24694. 10.1002/cjce.24694Search in Google Scholar

[43] J. R. White, Case study: Protein folding using homotopy methods, Scientific Computing with Case Studies, SIAM Press, Philadelphia, USA, 2009, pp. 1–8. https://www.cs.umd.edu/oleary/SCCS/supp/White/Homotopy_paper.pdf. Search in Google Scholar

[44] H. Jimenez-Islas, Sehpe: Programa para la solución de sistemas de ecuaciones no lineales mediante método homotópico con seguimiento hiperesférico, Av. Ing. Quím. 6 (1996), no. 2, 174–179. Search in Google Scholar

[45] J. M. Oliveros-Munoz and H. Jiménez-Islas, Hyperspherical path tracking methodology as correction step in homotopic continuation methods, Chem. Eng. Sci. 97 (2013), no. 28, 413–429, DOI: https://doi.org/10.1016/j.ces.2013.03.053. 10.1016/j.ces.2013.03.053Search in Google Scholar

[46] H. Jiménez-Islas, M. Calderón-Ramírez, G. M. Martínez-González, M. P. Calderón-Álvarado, and J. M. Oliveros-Munoz, Multiple solutions for steady differential equations via hyperspherical path-tracking of homotopy curves, Comput. Math. Appl. 78 (2020), no. 8, 2216–2239, DOI: https://doi.org/10.1016/j.camwa.2019.10.023. 10.1016/j.camwa.2019.10.023Search in Google Scholar

[47] L. O. Chua, Nonlinear Network Analysis-the Parametric Approach, Ph.D. Thesis, University of Illinois, June 1964. https://hdl.handle.net/2142/74281. Search in Google Scholar

[48] L. O., Chua and A. Ushida, A switching-parameter algorithm for finding multiple solutions of nonlinear resistive networks, Int. J. Circuit Theory Appl. 4 (1976), no. 3, 215–237, DOI: https://doi.org/10.1002/cta.4490040302. 10.1002/cta.4490040302Search in Google Scholar

[49] L. O. Chua and N. N. Wang, A new approach to overcome the overflow problem in computer-aided analysis of nonlinear resistive circuits, Int. J. Circuit Theory Appl. 3 (1975), no. 3, 261–284, DOI: https://doi.org/10.1002/cta.4490030305. 10.1002/cta.4490030305Search in Google Scholar

[50] L. Trajkovic, R. C. Melville, and S.-C. Fang, Finding DC operating points of transistor circuits using homotopy methods, 1991 IEEE Int. Symp. Circuits Syst. Singapore 2 (1991), 758–761, DOI: https://doi.org/10.1109/ISCAS.1991.176473. 10.1109/ISCAS.1991.176473Search in Google Scholar

[51] M. Tolikas, L. Trajkovic, and M. D. Ilic, Homotopy methods for solving decoupled power flow equations, 1992 IEEE Int. Symp. Circuits Syst. San Diego, CA, USA 6 (1992), 2833–2839, DOI: https://doi.org/10.1109/ISCAS.1992.230609. 10.1109/ISCAS.1992.230609Search in Google Scholar

[52] R. C. Melville, L. Trajkovic, S. -C. Fang, and L. T. Watson, Artificial parameter homotopy methods for the DC operating point problem, IEEE Trans. Comput.-Aided Des. Integr. Circuits Syst. 12 (1993), no. 6, 861–877, DOI: https://doi.org/10.1109/43.229761. 10.1109/43.229761Search in Google Scholar

[53] A. Dyes, E. Chan, H. Hofmann, W. Horia, and L. Trajkovic, Simple implementations of homotopy algorithms for finding DC solutions of nonlinear circuits, 1999 IEEE Int. Symp. Circuits Syst. Orlando, FL, USA 6 (1999), 290–293, DOI: https://doi.org/10.1109/ISCAS.1999.780152. 10.1109/ISCAS.1999.780152Search in Google Scholar

[54] L. Trajkovic, Homotopy Methods for Computing DC Operating Points, Wiley, 1999. 10.1002/047134608X.W2526Search in Google Scholar

[55] L. Kronenberg, W. Mathis, and L. Trajkovic, Limitations of criteria for testing transistor circuits for multiple DC operating points, Proc. 43rd IEEE Midwest Symp. Circuits Syst. Lansing, MI, USA 1 (2000), 148–151, DOI: https://doi.org/10.1109/MWSCAS.2000.951607. 10.1109/MWSCAS.2000.951607Search in Google Scholar

[56] L. Trajkovic, Homotopy Algorithms for Finding DC and Steady-State Solutions of Nonlinear Circuits, Electronics and Signal Processing Summer Symposium-LEOS, British Columbia, Canada, 2002. Search in Google Scholar

[57] L. Trajkovic, R. C. Melville, and S.-C. Fang, Improving DC convergence in a circuit simulator using a homotopy method, Proc. IEEE 1991 Custom Integr. Circuits Conf. San Diego, CA, USA, 1991, pp. 8.1/1–8.1/4, DOI: https://doi.org/10.1109/CICC.1991.164032. 10.1109/CICC.1991.164032Search in Google Scholar

[58] L. Trajkovic, E. Fung, and S. Sanders, HomSPICE: simulator with homotopy algorithms for finding DC and steady-state solutions of nonlinear circuits, 1998 IEEE Int. Symp. Circuits Syst. Monterey, CA, USA 6 (1998), 227–231, DOI: https://doi.org/10.1109/ISCAS.1998.705253. 10.1109/ISCAS.1998.705253Search in Google Scholar

[59] H. Vazquez-Leal, A. Marin-Hernandez, Y. Khan, A. Yıldırım, U. Filobello-Nino, R. Castaneda-Sheissa, et al., Exploring collision-free path planning by using homotopy continuation methods, Appl. Math. Comput. 219 (2013), no. 14, 7514–7532, DOI: https://doi.org/10.1016/j.amc.2013.01.038. Search in Google Scholar

[60] G. Diaz-Arango, H. Vazquez-Leal, L. Hernandez-Martinez, M. T. Sanz-Pascual, and M. Sandoval-Hernandez, Homotopy path planning for terrestrial robots using spherical algorithm, IEEE Trans. Autom. Sci. Eng. 15 (2017), no. 2, 567–585, DOI: https://doi.org/10.1109/TASE.2016.2638208. 10.1109/TASE.2016.2638208Search in Google Scholar

[61] G. C. Velez-Lopez, H. Vazquez-Leal, L. Hernandez-Martinez, A. Sarmiento-Reyes, G. Diaz-Arango, and J. Huerta-Chua, et al., A novel collision-free homotopy path planning for planar robotic arms, Sensors 22 (2022), no. 11:4022, 1–27, DOI: https://doi.org/10.3390/s22114022. 10.3390/s22114022Search in Google Scholar PubMed PubMed Central

[62] G. C. Velez-Lopez, L. Hernandez-Martinez, H. Vazquez-Leal, M. A. Sandoval-Hernández, V. M. Jimenez-Fernandez, and M. Gonzalez-Lee, Collision-free path planning applied to multi-degree-of-freedom robotic arms using homotopy methods, IEEE Access 12 (2024), 150702–150718, DOI: https://doi.org/10.1109/ACCESS.2024.3479095. 10.1109/ACCESS.2024.3479095Search in Google Scholar

[63] M. L. Quemada-Villagómez, M. D. L. L. López-González, J. M. Oliveros-Muñoz, and H. Jiménez-Islas, Ensennnanza del método de continuación homotópica con seguimiento hiperesférico para estudiantes de ingeniería, Acta Univ. 32 (2022), 1–43, DOI: https://doi.org/10.15174/au.2022.3358. 10.15174/au.2022.3358Search in Google Scholar

[64] M. A. Sandoval-Hernández, H Vázquez-Leal, J. Huerta-Chua, U. A. Filobello-Nino, and D. Mayorga-Cruz, La didáctica del cálculo integral: el caso de los procedimientos de integración, Rev. Iberoam. Investig. Desarro. Educ. 13 (2022), no. 25, 1–28, DOI: https://doi.org/10.23913/ride.v13i25.1245. 10.23913/ride.v13i25.1245Search in Google Scholar

[65] R. C. Uriza, G. M. Espinosa, and D. R. Gasperini, Análisis del discurso matemático escolar en los libros de texto: una mirada desde la teoría socioepistemológica, Av. Investig. Educ. Mat. 8 (2015), 9–28, DOI: https://doi.org/10.35763/aiem.v1i8.123. 10.35763/aiem.v1i8.123Search in Google Scholar

[66] C. Prieto-De-Castro, Biografía Henri Poincaré, Instituto de Matemáticas, Universidad Nacional Autónoma de México, 1977.Search in Google Scholar

[67] W. S. Massey, Introdución a la Topología Algebraica, Reverté, España, 2006. Search in Google Scholar

[68] M. Macho-Stadler, Homotopía, Publicaciones de divulgación, Universidad del País Vasco-Euskal Herriko Unibertsitatea, 2003. Search in Google Scholar

[69] J. Strom, Modern Classical Homotopy Theory, Graduate Studies in Mathematics, Vol. 127, American Mathematical Society, 2011. 10.1090/gsm/127Search in Google Scholar

[70] S. Lipschutz, Theory and Problems of General Topology, Schaumas Outline Series, Mc Graw-Hill, 1965. Search in Google Scholar

[71] M. Macho-Stadler, Qué es la Topología?, Publicaciones de divulgación, Universidad del País Vasco-Euskal, Herriko Unibertsitatea, Euskadi, 2018, https://divulgacioncientificadecientificos.blogspot.com/2018/08/que-es-la-topologia-marta-macho-stadler.html. Search in Google Scholar

[72] L. E. J. Brouwer, Beweis des Jordanschen Kurvensatzes, Math. Ann. 69 (1910), 169–175, DOI: https://doi.org/10.1007/BF01456867. 10.1007/BF01456867Search in Google Scholar

[73] L. E. J. Brouwer, Über Abbildung von Mannigfaltigkeiten, Math. Ann. Springer 71 (1911), 97–115, DOI: https://doi.org/10.1007/BF01456931. 10.1007/BF01456931Search in Google Scholar

[74] C. B. García and W. I. Zangwill, Pathways to Solution Fixed Point in Equilibria, Prentice Hall, 1981. Search in Google Scholar

[75] J. Leray and J. Schauder, Topologie et équations fonctionnelles, Ann. Sci. Éc. Norm. Supér. 51 (1934), 45–78. http://www.numdam.org/item/10.24033/asens.836.pdf. 10.24033/asens.836Search in Google Scholar

[76] F. A. Ficken, The continuation method for functional equations, Comm. Pure Appl. Math. 4 (1951), no. 4, 435–456, DOI: https://doi.org/10.1002/cpa.3160040405. 10.1002/cpa.3160040405Search in Google Scholar

[77] C. B. Haselgrove, The solution of non-linear equations and of differential equations with two-point boundary conditions, Computer J. Oxford University Press 4 (1961), no. 3, 255–259, DOI: https://doi.org/10.1093/comjnl/4.3.255. 10.1093/comjnl/4.3.255Search in Google Scholar

[78] D. F. Davidenko, On the approximate solution of a system of nonlinear equations, Ukr. Mat. Zh. 5 (1953), no. 2, 196–206. Search in Google Scholar

[79] E. Lahaye, Use Methode de Resolution daune Categorie deeeeeéquations Transcendentes, C.R. Acad. Sci. Paris 198 (1934), 1840–1842. Search in Google Scholar

[80] D. F. Davidenko, On the new method of numerically integrating a system of nonlinear equations, Dokl. Akad. Nauk SSSR 88 (1953), 601–604. Search in Google Scholar

[81] W. C., Rheinboldt, An adaptive continuation process for solving systems of nonlinear equations, Univ. Md. Banach Cent. Publ. 3 (1978), no. 1, 129–142. http://eudml.org/doc/208686. 10.4064/-3-1-129-142Search in Google Scholar

[82] P. Deuflhard, A Step Size Control for Continuation Methods and its Special Application to Multiple Shooting Techniques, Univertitat Hidelberg, Instiut für Angewandte Mathematik, Numerische Mathematik, Springer-Verlag, 1979. 10.1007/BF01399549Search in Google Scholar

[83] J. N. Lyness and B. J. J. McHugh, Integration over multidimensional hypercubes I. A progressive procedure, Comput. J. 6 (1963), no. 3, 264–270, DOI: https://doi.org/10.1093/comjnl/6.3.264. 10.1093/comjnl/6.3.264Search in Google Scholar

[84] E. L. Allgower and K. Georg, Introduction to numerical continuation methods, Classics in Applied Mathematics, Springer-Verlag, Philadelphia, 2003. 10.1137/1.9780898719154Search in Google Scholar

[85] H. Jiménez-Islas, Paquete Computacional Para la Solución de Sistemas de Ecuaciones no Lineales, PhD thesis, Instituto Tecnológico de Celaya, 1988, DOI: https://doi.org/10.13140/RG.2.2.12052.31363, https://acortar.link/cPGvWn. Search in Google Scholar

[86] A. Ushida and L. O. Chua, Tracing solution curves of non-linear equations with sharp turning points, Int. J. Circuit Theory Appl. 12 (1984), no. 1, 1–21, DOI: https://doi.org/10.1002/cta.4490120102. 10.1002/cta.4490120102Search in Google Scholar

[87] Y. Shinohara, A geometric method of numerical solution of nonlinear equation and error estimation by Urabeas proposition, RIMS Kyoto Univ. 5 (1969), no. 1, 1–9, DOI: https://doi.org/10.2977/PRIMS/1195194748. 10.2977/prims/1195194748Search in Google Scholar

[88] R. K. Brayton, F. G. Gustavson, and G. D. Hachtel, A new efficient algorithm for solving differential-algebraic systems using implicit backward differentiation formulas, Proc. IEEE 60 (1972), no. 1, 98–108, DOI: https://doi.org/10.1109/PROC.1972.8562. 10.1109/PROC.1972.8562Search in Google Scholar

[89] C. W. Gear, The automatic integration of ordinary differential equations, Commun. ACM 14 (1971), no. 3, 176–179, DOI: https://doi.org/10.1145/362566.362571. 10.1145/362566.362571Search in Google Scholar

[90] K. Yamamura, Simple algorithms for tracing solution curves, IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 40 (1993), no. 8, 537–541, DOI: https://doi.org/10.1109/81.242328. 10.1109/81.242328Search in Google Scholar

[91] R. W. Klopfenstein, Zeros of nonlinear functions, RCA Lab. N.J. ACM Digit. Libr. 8 (1961), no. 3, 366–373, DOI: https://doi.org/10.1145/321075.321080. 10.1145/321075.321080Search in Google Scholar

[92] E. Allgower and K. Georg, Introduction to Numerical Continuation Methods, Classics in Applied Mathematics, Springer-Verlag, Philadelphia, USA, 1990. 10.1007/978-3-642-61257-2Search in Google Scholar

[93] J. C. Alexander and J. A. Yorke, The homotopy continuation method: numerically implementable topological procedures, Trans. Amer. Math. Soc. 242 (1978), 271–284, DOI: https://doi.org/10.1090/S0002-9947-1978-0478138-5. 10.1090/S0002-9947-1978-0478138-5Search in Google Scholar

[94] S. N. Chow, J. Mallet-Paret, and J. A. Yorke, Finding zeros of maps: homotopy methods that are constructive with probability one, Math. Comput. 32 (1978), 887–899, DOI: https://doi.org/10.1090/S0025-5718-1978-0492046-9. 10.1090/S0025-5718-1978-0492046-9Search in Google Scholar

[95] L. T. Watson, Engineering applications of the Chow-Yorke algorithm, Appl. Math. Comput. 9 (1981), no. 2, 111–133, DOI: https://doi.org/10.1016/0096-3003(81)90010-2. 10.1016/0096-3003(81)90010-2Search in Google Scholar

[96] S. N. Chow, J. Mallet-Paret, and J. A. Yorke, A homotopy method for locating all zeros of a system of polynomials, Lecture Notes Math. 730 (1979). 77–88, DOI: https://doi.org/10.1007/BFb0064312. 10.1007/BFb0064312Search in Google Scholar

[97] T. Y. Li, T. Sauer, and J. A. Yorke, The Cheateras homotopy: an efficient procedure for solving systems of polynomial equations, SIAM J. Numer. Anal. 26 (1989), no. 5, 1241–1251, DOI: https://doi.org/10.1137/0726069. 10.1137/0726069Search in Google Scholar

[98] J. D. Seader, M. Kuno, W. Lin, S. A. Johnson, K. Unsworth, and J. W. Wiskin, Mapped continuation methods for computing all solutions to general systems of nonlinear equations, Comput. Chem. Eng. 14 (1990), no. 1, 71–85, DOI: https://doi.org/10.1016/0098-1354(90)87006-B. 10.1016/0098-1354(90)87006-BSearch in Google Scholar

[99] M. Kuno and J. D. Seader, Computing all real solutions systems of nonlinear equations with a global fixed-point homotopy, Ind. Eng. Chem. Res. 27 (1988), no. 7, 1320–1329, DOI: https://doi.org/10.1021/ie00079a037. 10.1021/ie00079a037Search in Google Scholar

[100] S. H. Choi and N. L. Book, Unreachable roots for global homotopy continuation methods, AIChE J. 37 (1991), no. 7, 1093–1095, DOI: https://doi.org/10.1002/aic.690370713. 10.1002/aic.690370713Search in Google Scholar

[101] S. Khaleghi-Rahimian, F. Jalali, J. D. Seader, and R. E. White, A new homotopy for seeking all real roots of a nonlinear equation, Comput. Chem. Eng. 35 (2011), no. 3 403–411, DOI: https://doi.org/10.1016/j.compchemeng.2010.04.007. 10.1016/j.compchemeng.2010.04.007Search in Google Scholar

[102] K. Yamamura, A fixed-point homotopy method for solving modified nodal equations, IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 46 (1999), no. 6, 654–665, DOI: https://doi.org/10.1109/81.768822. 10.1109/81.768822Search in Google Scholar

[103] D. M. Wolf and S. R. Sanders, Multiparameter, homotopy methods for finding DC operating points of nonlinear circuits, IEEE Trans. Circuits Syst. I Fundam. Theory Appl. 43 (1996), no. 10, 824–838, DOI: https://doi.org/10.1109/81.538989. 10.1109/81.538989Search in Google Scholar

[104] M. Sosonkina, L. T. Watson, and D. E. Stewart, Note on the end game in homotopy zero curve tracking, ACM Trans. Math. Software 22 (1996), no. 3, 281–287, DOI: https://doi.org/10.1145/232826.232843. 10.1145/232826.232843Search in Google Scholar

[105] H. Vazquez-Leal, A. Marin-Hernandez, Y. Khan, A Yıldırım, U. Filobello-Nino, and R. Castaneda-Sheissa, et al., Exploring collision-free path planning by using homotopy continuation methods, Appl. Math. Comput. Elsevier 219 (2013), no. 14, 7514–7532, DOI: https://doi.org/10.1016/j.amc.2013.01.038. 10.1016/j.amc.2013.01.038Search in Google Scholar

[106] R. Burden and J. Faires, Numerical Analysis, 9th edn., Cengage Learning, Boston, 2011. Search in Google Scholar

[107] S. C. Chapra and R. P. Canale. Métodos Numéricos Para Ingenieros, McGraw-Hill, Mexico, 2011. Search in Google Scholar

[108] A. S. Sedra, K. C. Smith, T .C. Carusone, and V. Gaudet, Microelectronic Circuits, 8th edn., Oxford University Press, Oxford, 2020. Search in Google Scholar

[109] R. Boylestad and L. Nashelsky, Electrónica: teoría de circuitos y dispositivos electrónicos, Pearson Educación, Mexico, 2009. Search in Google Scholar

[110] R. L. Tokheim, Digital Electronics, Glencoe, Glencoe, UK, 1994. Search in Google Scholar

[111] J. M. Angulo-Usategui, B. Garcia-Zapirain, I. Angulo-Martínez, and J. V. Sáez, Microcontroladores avanzados dsPIC: controladores digitales de señales: arquitectura, programación y aplicaciones, Thomson, Mexico, 2006. Search in Google Scholar

[112] B. B. Brey, Los microprocesadores intel: arquitectura, programación e interfaces, Pearson-Prentice-Hall, Mexico, 2007. Search in Google Scholar

[113] J. S. Thomsen and S. Schlesinger, Analysis of nonlinear circuits using impedance concepts, IRE Trans. Circuit Theory 2 (1955), no. 3, 271–278, DOI: https://doi.org/10.1109/TCT.1955.1085250. 10.1109/TCT.1955.1085250Search in Google Scholar

[114] J. S. Thomsen, Graphical analysis of nonlinear circuits using impedance concepts, J. Appl. Phys. AIP Publishing 24 (1953), no. 11, 1379–1382, DOI: https://doi.org/10.1063/1.1721182. 10.1063/1.1721182Search in Google Scholar

[115] M. A. Sandoval-Hernández, H. Vazquez-Leal, U. A. Filobello-Nino, G. J. Morales-AlarcóN, G. C. Velez-LóPez, and R. CastañEda-Sheissa, et al., Using Horner’s algorithm to reduce computing times, Int. J. Eng. Res. Technol. 12 (2023), no. 06, 325–331. https://www.ijert.org/research/using-horners-algorithm-to-reduce-computing-times-IJERTV12IS060143.pdf. Search in Google Scholar

[116] H. Vazquez-Leal, M. A. Sandoval-Hernandez, J. L. Garcia-Gervacio, A. L. Herrera-May, and U. A. Filobello-Nino, PSEM aproximations for both branches of Lambert W with applications, Discrete Dyn. Nat. Soc. 2019 (2019), 267951, DOI: https://doi.org/10.1155/2019/8267951. 10.1155/2019/8267951Search in Google Scholar

[117] L. T. Watson, A globally convergent algorithm for computing fixed-points of C2 maps, Appl. Math. Comput. 5 (1979), no. 4, 297–311, DOI: https://doi.org/10.1016/0096-3003(79)90020-1. 10.1016/0096-3003(79)90020-1Search in Google Scholar

[118] H. Mohamad-Nor, A. Rahman, A. I. Md-Ismail, and A. Abdul-Majid, Super Ostrowski homotopy continuation method for solving polynomial system of equations, Matematika 32 (2016), no. 1, 53–67, DOI: https://doi.org/10.11113/matematika.v32.n1.685. 10.11113/matematika.v32.n1.685Search in Google Scholar

[119] H. Mohamad-Nor, A. I. Md-Ismail, and A. Abdul-Majid, Linear fixed-point function for solving system of polynomial equations, AIP Conf. Proc. 1602 (2014), no. 1, 105–112, DOI: https://doi.org/10.1063/1.4882474. 10.1063/1.4882474Search in Google Scholar

[120] S. Sutherland, An introduction to Julia and Fatou sets, in: Fractals, Wavelets, and their Applications: Contributions from the International Conference and Workshop on Fractals and Wavelets, vol. 92, Springer Proceedings in Mathematics & Statistics, Switzerland, 2014, pp. 37–60. 10.1007/978-3-319-08105-2_3Search in Google Scholar

[121] G. N. Rubiano, Método de Newton, Mathematica y fractales: historia de una página, Bol. Mat. 14 (2007), no. 1, 44–63.Search in Google Scholar

[122] K. Madhu and J. Jayaraman, Higher order methods for nonlinear equations and their basins of attraction, Mathematics 4 (2016), no. 2, 22, DOI: https://doi.org/10.3390/math4020022. 10.3390/math4020022Search in Google Scholar

[123] B. Haible and R. B. Kreckel, CLN, A Class Library for Numbers, GiNaC C++ library, Germany, 1996. https://www.ginac.de/CLN/cln.pdf. Search in Google Scholar

[124] P. R. Bell and J. M. Doughertym, Nonlinear image processing method, IEEE Trans. Nucl. Sci. 25 (1978), no. 2, 928–938, DOI: https://doi.org/10.1109/TNS.1978.4329438. 10.1109/TNS.1978.4329438Search in Google Scholar

[125] H. Daisuke, E. Takaya, M. Kadowaki, Y. Kobayashi, and T. Ueda, Effect of the pixel interpolation method for downsampling medical images on deep learning accuracy, J. Comput. Commun. 9 (2021), no. 11, 150–156, DOI: https://doi.org/10.4236/jcc.2021.911010. 10.4236/jcc.2021.911010Search in Google Scholar

[126] A. Singh and J. Singh, Review and comparative analysis of various image interpolation techniques, 2019 IEEE 2nd Int. Conf. Intell. Comput. Instrument. Control Technol. Kannur, India 1 (2019), 1214–1218, DOI: https://doi.org/10.1109/ICICICT46008.2019.8993258. 10.1109/ICICICT46008.2019.8993258Search in Google Scholar

Received: 2024-04-03
Revised: 2024-11-21
Accepted: 2024-12-03
Published Online: 2025-02-03

© 2025 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Special Issue on Contemporary Developments in Graphs Defined on Algebraic Structures
  2. Forbidden subgraphs of TI-power graphs of finite groups
  3. Finite group with some c#-normal and S-quasinormally embedded subgroups
  4. Classifying cubic symmetric graphs of order 88p and 88p 2
  5. Simplicial complexes defined on groups
  6. Two-sided zero-divisor graphs of orientation-preserving and order-decreasing transformation semigroups
  7. Further results on permanents of Laplacian matrices of trees
  8. Research Articles
  9. Dynamics of particulate emissions in the presence of autonomous vehicles
  10. The regularity of solutions to the Lp Gauss image problem
  11. Exploring homotopy with hyperspherical tracking to find complex roots with application to electrical circuits
  12. The ill-posedness of the (non-)periodic traveling wave solution for the deformed continuous Heisenberg spin equation
  13. Some results on value distribution concerning Hayman's alternative
  14. 𝕮-inverse of graphs and mixed graphs
  15. A note on the global existence and boundedness of an N-dimensional parabolic-elliptic predator-prey system with indirect pursuit-evasion interaction
  16. On a question of permutation groups acting on the power set
  17. Chebyshev polynomials of the first kind and the univariate Lommel function: Integral representations
  18. Blow-up of solutions for Euler-Bernoulli equation with nonlinear time delay
  19. Spectrum boundary domination of semiregularities in Banach algebras
  20. Statistical inference and data analysis of the record-based transmuted Burr X model
  21. A modified predictor–corrector scheme with graded mesh for numerical solutions of nonlinear Ψ-caputo fractional-order systems
  22. Dynamical properties of two-diffusion SIR epidemic model with Markovian switching
  23. Classes of modules closed under projective covers
  24. On the dimension of the algebraic sum of subspaces
  25. Periodic or homoclinic orbit bifurcated from a heteroclinic loop for high-dimensional systems
  26. On tangent bundles of Walker four-manifolds
  27. Regularity of weak solutions to the 3D stationary tropical climate model
  28. A new result for entire functions and their shifts with two shared values
  29. Freely quasiconformal and locally weakly quasisymmetric mappings in metric spaces
  30. On the spectral radius and energy of the degree distance matrix of a connected graph
  31. Solving the quartic by conics
  32. A topology related to implication and upsets on a bounded BCK-algebra
  33. On a subclass of multivalent functions defined by generalized multiplier transformation
  34. Local minimizers for the NLS equation with localized nonlinearity on noncompact metric graphs
  35. Approximate multi-Cauchy mappings on certain groupoids
  36. Multiple solutions for a class of fourth-order elliptic equations with critical growth
  37. A note on weighted measure-theoretic pressure
  38. Majorization-type inequalities for (m, M, ψ)-convex functions with applications
  39. Recurrence for probabilistic extension of Dowling polynomials
  40. A generalized fixed-point theorem for set-valued mappings in b-metric spaces
  41. Unraveling chaos: A topological analysis of simplicial homology groups and their foldings
  42. Global existence and blow-up of solutions to pseudo-parabolic equation for Baouendi-Grushin operator
  43. A characterization of the translational hull of a weakly type B semigroup with E-properties
  44. Some new bounds on resolvent energy of a graph
  45. Carmichael numbers composed of Piatetski-Shapiro primes in Beatty sequences
  46. The number of rational points of some classes of algebraic varieties over finite fields
  47. Singular direction of meromorphic functions with finite logarithmic order
  48. Pullback attractors for a class of second-order delay evolution equations with dispersive and dissipative terms on unbounded domain
  49. Eigenfunctions on an infinite Schrödinger network
  50. Boundedness of fractional sublinear operators on weighted grand Herz-Morrey spaces with variable exponents
  51. On SI2-convergence in T0-spaces
  52. Bubbles clustered inside for almost-critical problems
  53. Classification and irreducibility of a class of integer polynomials
  54. Existence and multiplicity of positive solutions for multiparameter periodic systems
  55. Averaging method in optimal control problems for integro-differential equations
  56. On superstability of derivations in Banach algebras
  57. Investigating the modified UO-iteration process in Banach spaces by a digraph
  58. The evaluation of a definite integral by the method of brackets illustrating its flexibility
  59. Existence of positive periodic solutions for evolution equations with delay in ordered Banach spaces
  60. Tilings, sub-tilings, and spectral sets on p-adic space
  61. The higher mapping cone axiom
  62. Continuity and essential norm of operators defined by infinite tridiagonal matrices in weighted Orlicz and l spaces
  63. A family of commuting contraction semigroups on l 1 ( N ) and l ( N )
  64. Pullback attractor of the 2D non-autonomous magneto-micropolar fluid equations
  65. Maximal function and generalized fractional integral operators on the weighted Orlicz-Lorentz-Morrey spaces
  66. On a nonlinear boundary value problems with impulse action
  67. Normalized ground-states for the Sobolev critical Kirchhoff equation with at least mass critical growth
  68. Decompositions of the extended Selberg class functions
  69. Subharmonic functions and associated measures in ℝn
  70. Some new Fejér type inequalities for (h, g; α - m)-convex functions
  71. The robust isolated calmness of spectral norm regularized convex matrix optimization problems
  72. Multiple positive solutions to a p-Kirchhoff equation with logarithmic terms and concave terms
  73. Joint approximation of analytic functions by the shifts of Hurwitz zeta-functions in short intervals
  74. Green's graphs of a semigroup
  75. Some new Hermite-Hadamard type inequalities for product of strongly h-convex functions on ellipsoids and balls
  76. Infinitely many solutions for a class of Kirchhoff-type equations
  77. On an uncertainty principle for small index subgroups of finite fields
  78. On a generalization of I-regularity
  79. Spin(8, ℂ)-Higgs bundles fixed points through spectral data
  80. Coloring the vertices of a graph with mutual-visibility property
  81. Embedding of lattices and K3-covers of an enriques surface
  82. Algorithm and linear convergence of the H-spectral radius of weakly irreducible quasi-positive tensors
  83. q-Stirling sequence spaces associated with q-Bell numbers
  84. Multiple G-Stratonovich integral in G-expectation space
Downloaded on 5.10.2025 from https://www.degruyterbrill.com/document/doi/10.1515/math-2024-0115/html
Scroll to top button