# Computing singular points of projective plane algebraic curves by homotopy continuation methods.

1. IntroductionAlgebraic curves are a classic research object in mathematics. The related computation of algebraic curves arises in several applications, including number theoretic problems [1], ancient and modern architectural designs, error-correcting [2], biological shape [3], cryptographic algorithms [4, 5], and computer aided geometric design [6, 7]. Singular points and their multiplicities and characters play an important role in the research of algebraic curves [8]. They help to determine the genus of an algebraic curve. Singular points also show some shape features, such as nodes, self-intersections or cusps of real curves in robot motion planning, computer aided geometry design, and machine vision. And the determination of geometric shape and topology of the real curves depends on the singular points. The computation of singular points is also crucial for tracing curves algorithms [9].

As singular points of an algebraic curve are the solutions of a polynomial system, most of algorithms in [10-15] solve the polynomial system either by Grobner basis method described in [16] or by resultant computation. These methods rely on symbolic algebraic computation that requires exact input of coefficients of algebraic curves. It is easy for them to suffer from the overwhelming coefficient swell. So these methods may be limited to relatively small problem. Furthermore, it is difficult to obtain exact coefficients of algebraic curves due to data error in many fields of science and engineering.

Although we have presented an algorithm to compute singular points of irreducible algebraic curves in [17], the algorithm is almost experimental and related analysis of the algorithm is not provided, for instance, the feasibility and complexity. The aims of this paper are analyzing the feasibility of this algorithm for computing singular points of reducible algebraic curves and proving that the algorithm has polynomial time complexity on the degree of an algebraic curve. We also present the effect of tiny perturbations of coefficients of algebraic curves on singular points and it may be a part of reasons why we compute the singular points numerically.

Our algorithm is totally numerical and can deal with algebraic curves with inexact coefficients. It includes the applications of homotopy continuation methods and the method of calculating the multiple roots of univariate polynomials with inexact coefficients without using multiprecision arithmetic. Homotopy continuation method is a reliably and efficiently numerical method to solve the polynomial systems [18]. Several methods have been presented to compute roots of univariate polynomials, such as Laguerre's method, Jenkins-Traub method, and the QR algorithm with the companion matrix. However, they cannot overcome a barrier of attainable accuracy on an m multiple root [19, 20]. If we have [k.sub.1] digits coefficients accuracy and [k.sub.2] digits machine precision, the attainable accuracy is min[[k.sub.1], [k.sub.2]}/m digits. For example, if the standard double precision of 16 decimal digits is used and the accuracy of coefficients of the polynomial is 15 digits, only 3 correct digits of a root of multiplicity 5 can be obtained by the above methods. A method proposed by Zeng [21] is not subject to the accuracy barrier and a lot of numerical experiments have shown its efficiency and robustness.

The following section presents the basic notions on algebraic curves and singular points. In Section 3, we focus on the algorithm of computing singular points. Section 4 is devoted to the numerical experiments. We conclude this paper in Section 5.

2. Preliminaries

In this section, we introduce some basic notions and results on algebraic curves and singular points. Let C be the complex numbers field and let [P.sup.2](C) be the projective plane over C.

Definition 1 (see [8, 22]). A projective plane algebraic curve over C is defined as the set

C = {(a : b : c) [member of] [P.sup.2] (C) | F (a, b, c) = 0} (1)

for a nonconstant square-free homogeneous polynomial F(x, y, z) [member of] C[x, y, z].

We call F(x,y,z) the defining polynomial of C. The degree of polynomial F(x, y, z) is called the degree of C.

We take the line z = 0 in [P.sup.2](C) as the line at infinity. If the projective plane algebraic curve C is defined by the polynomial F(x, y, z), then the corresponding affine plane algebraic curve C* is defined by the dehomogenization f(x, y) of F(x, y, z),

[C.sup.*] = {(a, b) [member of] [A.sup.2] (C) | f (a, b) = 0}, (2)

where [A.sup.2](C) is the affine plane over C.

Therefore, if

F (x, y, z) = [f.sub.n] (x, y) + [f.sub.n-1] (x, y) z + ... + [f.sub.0] (x, y) [z.sup.n], (3)

then

f(x, y) = [f.sub.n](x, y) + [f.sub.n-1] (x, y) + ... + [f.sub.0](x, y), (4)

where [f.sub.k](x, y) is a homogeneous polynomial of degree k and [f.sub.n](x, y) is nonzero.

Throughout this paper, whenever we speak of a "curve" we mean an "algebraic curve" For brevity, we will usually use the phrase "the curve F(x,y,z) = 0" instead of "the projective plane algebraic curve whose defining polynomial is F(x, y, z)" Of course, a polynomial G = cF, for some nonzero c e C, defines the same curve, so F is unique only up to multiplication by nonzero constants. The above usage is also for the case of affine plane algebraic curve.

Every point (a, b) on [C.sup.*] corresponds to a point (a : b : 1) on C, and every additional point on C is a point at infinity.

In other words, the first two coordinates of additional points are the nontrivial solutions of [f.sub.n](x, y) = 0, with the third coordinates being 0. Thus, the curve C has only finitely many points at infinity.

Definition 2 (see [8, 22]). Let [C.sup.*] be an affine plane algebraic curve over C defined by f(x,y) [member of] C[x, y], and let P = (a,b) [member of] [C.sup.*]. P is of multiplicity r on [C.sup.*] if and only if all derivatives of f(x, y) up to and including the (r - 1)th vanish at P, but at least one rth derivative does not vanish at P.

A point of multiplicity two or more is said to be a singular point, and especially the point of multiplicity two is called a double point. A point of multiplicity one is called a simple point.

It is evident that a necessary and sufficient condition that a point (a, b) of the curve f(x, y) = 0 is singular is that

f (a, b) = 0, [f.sub.x] (a, b) = 0, [f.sub.y] (a, b) = 0. (5)

Definition 3 (see [8]). Let the parametric equations of tangents to the curve f(x, y) = 0 at P = (a, b) be

x = a + [lambda]t, y = b + [mu]t, (6)

and the tangents to the curve f(x, y) = 0 at P = (a, b) of multiplicity r are determined by the ratio A: p and correspond to the roots of

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (7)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], and they are counted with multiplicities equal to the multiplicities of the corresponding roots of this equation.

Definition 4 (see [8,22]). A singular point P of multiplicity r on an affine plane algebraic curve [C.sup.*] is called ordinary if and only if the r tangents to [C.sup.*] at P are distinct and nonordinary otherwise. The property of singular point P of being ordinary or nonordinary is called the character of P.

The criteria for the multiplicity of a point on the projective plane algebraic curve can be characterized as follows.

Proposition 5 (see [8]). Q [member of] [P.sup.2](C) is a point of multiplicity r on the projective plane algebraic curve F(x, y, z) = 0 if and only if all the (r - 1)th derivatives of F(x, y, z), but not all the rth derivatives, vanish at Q.

As a corollary of this proposition, a point Q of the curve F(x, y, z) = 0 being singular can be put in a more convenient form.

Proposition 6 (see [8, 22]). Q [member of] [P.sup.2](C) is a singular point of the projective plane algebraic curve F(x, y, z) = 0 if and only if

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (8)

The following two theorems are of use for analyzing the algorithm of the next section.

Theorem 7 (see [8]). If a curve of degree n with no multiple components has multiplicities [r.sub.i] at points [P.sub.i], then

(n - 1)(n - 2) [greater than or equal to] [summation] [r.sub.i] ([r.sub.i] - l). (9)

Theorem 8 (see [8] Bezout's theorem). Two curves of degrees m and n with no common components have exactly mn intersections.

3. Computation of Singular Points

In this section, a method of solving the overdetermined polynomial systems is presented first. We also outline the algorithm on computing the singular points of projective plane algebraic curves, and afterwards we analyze feasibility and complexity of the algorithm. The last subsection includes the effect of tiny perturbations of coefficients of the curve on singular points.

3.1. Solving Overdetermined Polynomial Systems. The following proposition shows how to reduce an overdetermined polynomial system to a square system that most of the papers and software programs on homotopy continuation method focus on. Let n and N be the number of equations and unknowns, respectively, and n > N.

Proposition 9 (see [23,24]). There are nonempty Zariski open dense sets of parameters [[lambda].sub.ij] [member of] [C.sup.N(n-N)] or [[lambda].sub.ij] [member of] [R.sup.N(n-N)] such

that every isolated solution of

[f.sub.1] ([x.sub.1], ..., [x.sub.N]) = 0, [f.sub.n] ([x.sub.1], ..., [x.sub.N]) = 0 (10)

is an isolated solution of

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (11)

The solutions set of an overdetermined system belongs to that of the corresponding square system, but the converse is not true.

3.2. Algorithm. We summarize the algorithm described in [17] for computing the singular points of irreducible projective plane algebraic curves.

Algorithm 1.

Input. Consider projective plane algebraic curve F(x, y, z) = 0 and corresponding affine plane algebraic curve f(x, y) = 0. Threshold is [member of] > 0.

Output. Consider singular points and their multiplicities and characters of the curve F(x, y, z) = 0.

Step 1 (randomization). By Proposition 9 and choosing the random complex numbers [alpha] and [beta], we add random multiples of the last equation to the first two equations of the polynomial system (5):

f(x, y) + [alpha][f.sub.y] (x, y) = 0, [f.sub.x] (x, y) + [beta][f.sub.y] (x, y) = 0. (12)

Step 2 (computing the singular points (a : b : l)). We solve the polynomial system (12) by the polyhedral homotopy continuation method first and denote the solutions set as {(a : b : 1)}. If [absolute value of f(a, b)] < [epsilon], (a : b : 1) will be a singular point.

Step 3 (computing the singular points (a : b : 0) at infinity). Since [f.sub.n](x, y) is a bivariate homogeneous polynomial of degree n, we obtain the points set {(a : b : 0)} at infinity by solving the corresponding univariate polynomial. By Proposition 6, if Max{([partial derivative]F/[partial derivative]x)(a, b, 0), ([partial derivative]F/[partial derivative]y) (a, b, 0), ([partial derivative]F/[partial derivative]z)(a, b, 0)} < [epsilon], then (a : b : 0) will be a singular point at infinity.

When we determined the multiplicity and character of singular point (a : b : 0) at infinity, (a : b : 0) will be transferred to point (a' : b' : 1) by simple linear transformation.

Step 4 (determining the multiplicity). Evaluating the derivatives of f(x, y) at a singular point from order equal to 2, if the moduli of all derivatives of f(x, y) up to and including the (r - 1)th are less than e, but at least one modulus of rth derivatives is greater than e, then the multiplicity of this singular point will be determined as r.

Step 5 (determining the character). The tangents to the curve f(x, y) = 0 at a singular point of multiplicity r correspond to the roots of (7). Whether the singular point is ordinary or not will be determined if there are multiple roots of (7).

3.3. Remarks on Feasibility of Algorithm. The following theorems and remarks will show the feasibility of Algorithm 1 for computing the singular points of reducible projective plane algebraic curves.

Theorem 10. All the solutions of polynomial system (12) are isolated.

Proof. Since f(x,y) is square-free, f(x,y) and [f.sub.x](x,y) have no nonconstant common divisors. It also holds for f(x, y) and [f.sub.y](x, y). In Step 1 of Algorithm 1, for the random complex numbers [alpha] and [beta], f(x, y) + [alpha][f.sub.y](x, y) and [f.sub.x](x, y) + [beta][f.sub.y](x, y) also have no nonconstant common divisors. It follows that all the solutions of polynomial system (12) are isolated.

Remark 11. The inequality in Theorem 7 shows that an irreducible curve has finitely many singular points. So, the singular point of an irreducible curve is isolated.

Theorem 12. A reducible curve f(x, y) = 0 has only finitely many singular points and they are all isolated points.

Proof. Without loss of generality, we may assume that f(x, y) = [f.sub.1](x, y) [f.sub.2](x, y), where [f.sub.1](x, y) and [f.sub.2](x, y) are irreducible polynomials. Obviously, the singular points of f(x, y) = 0 are the singular points of [f.sub.1](x, y) = 0 (i = 1, 2) and the intersections of [f.sub.1](x, y) = 0 and [f.sub.2](x, y) = 0. By Remark 11 and Theorem 8, the set of singular points of f(x, y) = 0 is finite and the singular points are all isolated.

Remark 13. When solving the polynomial system (12) in Step 2 of Algorithm 1 by homotopy continuation methods, Proposition 9 and Theorem 12 imply that all singular points (a : b : 1) of reducible curve can be determined.

Remark 14. In Step 2 of Algorithm 1, the Jacobian of polynomial system (12) is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (13)

This matrix is singular at a singular point of multiplicity three at least. The numerical techniques of homotopy continuation methods could deal with this well.

Remark 15. As f(x, y) is a nonzero polynomial and has finite degree n, the modulus of some derivative of order less than or equal to n must be greater than the given threshold at a singular point. Hence, Step 4 of Algorithm 1 ends for determining the multiplicities of finitely many singular points.

Remark 16. We solve univariate polynomials in Steps 3 and 5 of Algorithm 1. Numerical singular points make the coefficients of univariate polynomials in Step 5 inexact. We have mentioned the efficiency of computation of multiple roots of inexact univariate polynomials proposed by Zeng in Section 1. Employing Zeng's method, we compute the singular points at infinity accurately in Step 3 and the characters of all singular points can be determined correctly in Step 5.

3.4. Computational Complexity. In this subsection, we analyze the computational complexity of Algorithm 1. The degree of the curve f(x, y) = 0 is denoted by n.

Theorem 17. The complexity of Step 2 in Algorithm 1 is O([n.sup.2]).

Proof. Shub and Smale [25, 26] present that finding an approximate zero of a polynomial system by homotopy continuation methods can be solved in polynomial time on the average and the number of arithmetic operations is bounded by c[N.sup.4], where c is a universal constant and N is the number of variables. By Bezout's Theorem 8, the number of solutions of polynomial system (12) is n(n - 1) at most. We deduce that the complexity of Step 2 is no more than 16cn(n - 1).

Theorem 18. The complexity of Steps 3 and 5 in Algorithm 1 is O([n.sup.3]) and O([n.sup.5]) or less, respectively.

Proof. Pan reports [19] that the complexity of general root-finders of univariate polynomial of degree d is O([d.sup.2]) or less, but Zeng's method [21] reaches high accuracy on multiple roots at the higher computing cost of O([d.sup.3]). It follows that the complexity of Step 3 is O([n.sup.3]). From Steps 2 and 3, we know that the number of singular points of the curve f(x, y) = 0 of degree n is [n.sup.2] at most. Therefore, we conclude that the complexity of Step 5 is O([n.sup.5]) or less.

Theorem 19. The complexity of Step 4 in Algorithm 1 does not exceed O([n.sup.5]).

Proof. The number of all the nonzero derivatives of bivariate polynomial f(x, y) of degree n is not more than n(n + 3)/2. The complexity of evaluation of a polynomial at a point is linear in the degree of the polynomial and the number of variables. Since there are [n.sup.2] points at most in Step 4 of Algorithm 1 and the complexity of linear change of coordinates is polynomial in the number of variables, the complexity of Step 4 does not exceed O([n.sup.5]).

The following theorem holds at once.

Theorem 20. The complexity of Algorithm 1 is polynomial time in the degree n of the projective plane algebraic curve and is O([n.sup.5]).

3.5. Effect of Tiny Perturbations on Singular Points. The authors [27] conclude that random perturbations of coefficients of plane algebraic curves will almost invariably destroy all singular points. The current methods dealing with plane algebraic curves are almost ill-conditioned with respect to tiny perturbations. For example, let us consider the curve [x.sup.3] - [y.sup.2] = 0 with nonordinary double singular point (0,0) and its perturbation [x.sup.3] - [y.sup.2] + [epsilon][x.sup.2] = 0, where [epsilon] > 0. It is natural to compute the approximate singular point out and the original character for relative magnitude [epsilon] > 0. However, the character of singular point (0, 0) is changed by using symbolic algebraic computation, even if [epsilon] > 0 is sufficiently small.

For the curve [x.sup.3] - [y.sup.2] + [epsilon][x.sup.2] = 0, when [epsilon] = 0.0000001, the numerical procedure based on Algorithm 1 outputs a nonordinary double singular point

(-0.00000006666667 + 0.00000000000000i, 0.00000000000000 - 0.000000000000000i). (14)

4. Numerical Experiments

In this section, Algorithm 1 is implemented in Matlab and some examples are presented to show the efficiency of corresponding numerical procedure. There are many freely available homotopy continuation packages. As Lee et al. [28] report that HOM4PS-2.0 is generally the fastest package for solving small to moderately large sparse systems, it is employed to solve polynomial system (12) in Algorithm 1. Many numerical examples are provided in [17] where the curves are irreducible and the coefficients of the curves are inexact. We finish this section with two reducible curves.

Example 1. Consider

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (15)

There are two components [f.sub.1](x, y) = 0 and [f.sub.2](x, y) = 0 for the curve f(x, y) = 0. The curve [f.sub.1](x, y) = 0 has no singular points and its real part is plotted in Figure 1. The curve [f.sub.2](x, y) = 0 has nonordinary double singular point corresponding to (0 : 0: 1) and ordinary double singular point corresponding to (0 : 1 : 0). Figure 2 shows real part of the curve [f.sub.2](x, y) = 0 and the singular point corresponding to (0 : 0 : 1) is marked by a solid square.

The curves [f.sub.1](x, y) = 0 and [f.sub.2](x, y) = 0 have 8 intersections corresponding to [Q.sub.i](i = 1, ..., 8). Their related projective curves intersect [Q.sub.9] = (0,1, 0) at infinity. (0,0) is the simple point of [f.sub.1](x, y) = 0, nonordinary double point of [f.sub.2](x, y) = 0. Therefore, (0, 0) corresponding to [Q.sub.1] = (0:0:1) is nonordinary singular point of multiplicity of 3 of f(x, y) = 0. [Q.sub.9] = (0:1:0) is ordinary double point and simple point of the corresponding projective curves for [f.sub.2](x, y) = 0 and [f.sub.1](x, y) = 0, respectively. The tangent to the latter at [Q.sub.9] coincides with one tangent to the former at [Q.sub.9]. So, [Q.sub.9] = (0 : 1 : 0) is nonordinary singular point of multiplicity of 3 of the projective curve corresponding to f(x, y) = 0.

Table 1 lists the results of our numerical procedure for the projective curve corresponding to f(x, y) = 0 in Example 1. Figure 3 presents the real part of f(x, y) = 0 in Example 1. The solid dot marks the real intersection of [f.sub.1](x, y) = 0 and [f.sub.2](x,y) = 0 corresponding to [Q.sub.5]. The real singular point corresponding to [Q.sub.1] = (0:0:1) is still marked by a solid square.

We use the abbreviations "Mult." and "Ord" for the words "multiplicity" and "ordinary" in our tables, respectively.

Example 2. Consider

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (16)

The curve f(x, y) = 0 has two components [f.sub.1](x, y) = 0 and [f.sub.2](x, y) = 0. There are six singular points corresponding to [P.sub.i](i = 1, ..., 6) for the curve [f.sub.1](x, y) = 0 and they are marked by solid squares in Figure 4 where the real part of the curve [f.sub.1](x, y) = 0 is plotted. Four singular points corresponding to [P.sub.i](i = 7, ..., 10) of the curve [f.sub.2](x, y) = 0 are illustrated by the solid triangles in Figure 5 where the real part of the curve [f.sub.2](x, y) = 0 is presented. Obviously, [P.sub.i](i = 1, ..., 10) are singular points of the projective curve corresponding to the curve f(x, y) = 0.

There are 29 intersections corresponding to [P.sub.i](i = 11, ..., 39) for the curves [f.sub.1](x, y) = 0 and [f.sub.2](x, y) = 0. Their related projective curves intersect [P.sub.40] at infinity. These 30 intersections [P.sub.i] (i = 11, ..., 40) are also singular points of the projective curve corresponding to the curve f(x, y) = 0. So the curve f(x, y) = 0 has 40 singular points.

Figure 6 shows the real part of the curve f(x, y) = 0 whose real singular points are marked by the solid dots, squares, and triangles. The solid dots mark the real intersections of the curves [f.sub.1](x, y) = 0 and [f.sub.2](x, y) = 0. The solid squares and triangles illustrate the singular points of the curves [f.sub.1](x, y) = 0 and [f.sub.2](x, y) = 0, respectively. Table 2 lists the results of our numerical procedure for the projective curve corresponding to f(x, y) = 0 in Example 2.

It takes 0.28628 seconds and 1.26786 seconds for our numerical procedure to obtain the results of Examples 1 and 2, respectively, on the Lenovo PC with Pentium Dual Core, CPU of 2.5 GHZ, and memory of 1.99 GB. We usually choose [10.sup.-6] as the thresholds in Steps 2, 3, and 4 in Algorithm 1.

5. Conclusions

This paper provides an effective algorithm for computing the singular points of projective plane algebraic curves by homotopy continuation methods. The determination of multiplicities relies on the accuracy of singular points. The characters of singular points are determined by Zeng's method for computing multiplicities of roots of inexact univariate polynomials. The precision of coefficients of inexact univariate polynomials in our algorithm depends on the accuracy of singular points. Several numerical examples are presented to illustrate the efficiency of our algorithm. We will discuss the adaptive thresholds of the algorithm in the future.

http://dx.doi.org/10.1155/2014/230847

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This paper is supported by the National Natural Science Foundation of China (no. 11171052).

References

[1] D. Poulakis and E. Voskos, "Solving genus zero Diophantine equations with at most two infinite valuations," Journal of Symbolic Computation, vol. 33, no. 4, pp. 479-491, 2002.

[2] O. Pretzel, Codes and Algebraic Curves, vol. 8, Oxford University Press, 1998.

[3] C. Bajaj, H. Lee, R. Merkert, and V. Pascucci, "NURBS based B-rep models from macromolecules and their properties," in Proceedings of 4rth Symposium on Solid Modeling and Applications, C. Hoffmann and W. Bronsvort, Eds., pp. 217-228, ACM Press, Atlanta, Ga, USA, 1997

[4] J. A. Buchmann, Introduction to Cryptography, Springer, 2001.

[5] N. Koblitz, "Good and bad uses of elliptic curves in cryptography," Moscow Mathematical Journal, vol. 2, no. 4, pp. 693-715, 2002.

[6] T. W. Sederberg, "Applications to computer aided geometric design," in Proceedings of Symposia in Applied Mathematics. Applications of Computational Algebraic Geometry, vol. 53, pp. 67-89, 1998.

[7] G. Farin, J. Hoschek, and M. S. Kim, Handbook of Computer Aided Geometric Design, North Holland, 2002.

[8] R. J. Walker, Algebraic Curves, Springer, New York, NY, USA, 1978.

[9] C. L. Bajaj, C. M. Hoffmann, R. E. Lynch, and J. E. H. Hopcroft, "Tracing surface intersections," Computer Aided Geometric Design, vol. 5, no. 4, pp. 285-307, 1988.

[10] T. W. Sederberg, D. C. Anderson, and R. N. Goldman, "Implicit representation of parametric curves and surfaces," Computer Vision, Graphics, and Image Processing, vol. 28, pp. 72-84, 1984.

[11] S. S. Abhyankar and C. L. Bajaj, "Automatic parameterization of rational curves and surfaces--III. Algebraic plane curves," Computer Aided Geometric Design, vol. 5, no. 4, pp. 309-321, 1988.

[12] J. R. Sendra and F. Winkler, "Symbolic parametrization of curves," Journal of Symbolic Computation, vol. 12, no. 6, pp. 607-631, 1991.

[13] D. Cox, "Curves, surfaces, and syzygies, topics in algebraic geometry and geometric modeling," Contemporary Mathematics, vol. 334, pp. 131-149, 2003.

[14] Y. Sun and J. Yu, "Implicitization of parametric curves via Lagrange interpolation," Computing, vol. 77, no. 4, pp. 379-386, 2006.

[15] T. Sakkalis and R. Farouki, "Singularpoints of algebraic curves," Journal of Symbolic Computation, vol. 9, no. 4, pp. 405-421, 1990.

[16] D. Cox, J. Little, and D. O'Shea, Varieties, and Algorithms: An Introduction to Computational Algebraic Geometry and Commutative Algebra, Springer, New York, NY, USA, 2nd edition, 1997.

[17] Z.-x. Luo, E.-b. Feng, and W.-y. Hu, "Computing numerical singular points of plane algebraic curves," Communications in Mathematical Research, vol. 28, no. 2, pp. 146-158, 2012.

[18] T. Y. Li, "Numerical solution of polynomial systems by homotopy continuation methods," in Handbook of Numerical Analysis, vol. 11, pp. 209-304, North-Holland, Amsterdam, The Netherlands, 2003.

[19] V. Y. Pan, "Solving a polynomial equation: some history and recent progress," SIAM Review, vol. 39, no. 2, pp. 187-220, 1997

[20] M. Igarashi and T. Ypma, "Relationships between order and efficiency of a class of methods for multiple zeros of polynomials," Journal of Computational and Applied Mathematics, vol. 60, no. 1-2, pp. 101-113, 1995.

[21] Z. Zeng, "Computing multiple roots of inexact polynomials," Mathematics of Computation, vol. 74, no. 250, pp. 869-903, 2005.

[22] J. R. Sendra, F. Winkler, and S. Perez-Diazs, Rational Algebraic Cueves, Springer, Berlin, Germany, 2008.

[23] A. J. Sommese and C. W. Wampler, "Numerical algebraic geometry," in The Mathematics of Numerical Analysis (), J. Renegar, M. Shub, and S. Smale, Eds., vol. 32 of Lectures in Applied Mathematics, pp. 749-763, American Mathematical Society, Park City, Utah, USA, 1996.

[24] A. P. Morgan and A. J. Sommese, "Coefficient-parameter polynomial continuation," Applied Mathematics and Computation, vol. 29, no. 2, pp. 123-160, 1989.

[25] M. Shub and S. Smale, "Complexity of Bezout's theorem--V. Polynomial time," Theoretical Computer Science, vol. 133, no. 1, pp. 141-164, 1994.

[26] M. Shub, "Complexity of Bezout's theorem--VI. Geodesics in the condition (number) metric," Foundations of Computational Mathematics, vol. 9, no. 2, pp. 171-178, 2009.

[27] R. T. Farouki and V. T. Rajan, "On the numerical condition of algebraic curves and surfaces--I. Implicit equations," Computer Aided Geometric Design, vol. 5, no. 3, pp. 215-252, 1988.

[28] T. L. Lee, T. Y. Li, and C. H. Tsai, "HOM4PS-2.0: a software package for solving polynomial systems by the polyhedral homotopy continuation method," Computing. Archives for Scientific Computing, vol. 83, no. 2-3, pp. 109-133, 2008.

Zhongxuan Luo, (1,2) Erbao Feng, (1) and Jielin Zhang (1)

(1) School of Mathematical Sciences, Dalian University of Technology, Dalian 116024, China

(2) School of Software, Dalian University of Technology, Dalian 116620, China

Correspondence should be addressed to Zhongxuan Luo; zxluo@dlut.edu.cn

Received 10 April 2014; Accepted 8 May 2014; Published 5 June 2014

Academic Editor: Baodong Zheng

TABLE 1: Results for f(x, y) = 0 in Example 1. x y [Q.sub.1] 0.00000003 + 0.00000005i 0.00000000 + 0.00000000i [Q.sub.2] -0.89868027 + 0.41873939i -0.15046254 + 1.06341745i [Q.sub.3] 0.36804326 - 1.20320368i -0.92239819-0.40178011i [Q.sub.4] 0.36804326 + 1.20320368i -0.92239819 + 0.40178011i [Q.sub.5] -0.49367385 - 0.00000000i 0.28024456 + 0.00000000i [Q.sub.6] -0.89868027 - 0.41873939i -0.15046254 - 1.06341745i [Q.sub.7] 0.77747393 + 0.21532054i -0.56726154 - 0.66500198i [Q.sub.8] 0.77747393 - 0.21532054i -0.56726154 + 0.66500198i [Q.sub.9] 0.00000000 1.00000000 Z Mult. Ord. [Q.sub.1] 1.00000000 3 No [Q.sub.2] 1.00000000 2 Yes [Q.sub.3] 1.00000000 2 Yes [Q.sub.4] 1.00000000 2 Yes [Q.sub.5] 1.00000000 2 Yes [Q.sub.6] 1.00000000 2 Yes [Q.sub.7] 1.00000000 2 Yes [Q.sub.8] 1.00000000 2 Yes [Q.sub.9] 0 3 No TABLE 2: Results for f(x, y) = 0 in Example 2. x y [P.sub.1] 0.00000000 0.00000000 - 0.00000000i [P.sub.2] -0.50000000 - 0.00000000i 0.99999999 + 0.00000000i [P.sub.3] -0.00000000 - 0.00000000i -1.00000000 + 0.00000000i [P.sub.4] 0.99999999 + 0.00000000i -0.50000000 - 0.00000000i [P.sub.5] -1.00000000 + 0.00000000i -0.50000000 + 0.00000000i [P.sub.6] 0.50000000 + 0.00000000i 1.00000000 + 0.00000000i [P.sub.7] 2.41421356 - 0.00000000i 0.41421356 - 0.00000000i [P.sub.8] -2.41421356 + 0.00000000i -0.41421356--0.00000000i [P.sub.9] 0.41421356 - 0.00000000i 2.41421356-0.00000000i [P.sub.10] -0.41421356 - 0.00000000i -2.41421356 + 0.00000000i [P.sub.11] 1.86368244 + 0.00000000i -0.81167781 + 0.00000000i [P.sub.12] -0.56599742 + 0.66867689i 0.95095432--0.12549545i [P.sub.13] -0.58308154 + 0.00000000i -2.34385106 + 0.00000000i [P.sub.14] 1.67570344 + 0.00000000i 1.32326562 + 0.00000000i [P.sub.15] 0.79829233 1.77539598 + 0.00000000i [P.sub.16] 0.09549542 - 0.39009016i -0.97196258 + 0.89462218i [P.sub.17] 0.52445070 + 0.00000000i -2.19040368--0.00000000i [P.sub.18] 0.05975897 - 0.25434442i -0.06366391 - 0.36672723i [P.sub.19] 0.15502368 - 1.39689139i -0.01303583 + 0.19325672i [P.sub.20] 0.51269502 + 0.43085324i -0.51116222 + 0.08768075i [P.sub.21] 0.15502368 + 1.39689139i -0.01303583 - 0.19325672i [P.sub.22] 0.09549542 + 0.39009016i -0.97196258 - 0.89462218i [P.sub.23] 0.54660137 + 0.59723958i -0.36484540 - 0.17346068i [P.sub.24] 1.40455064 + 0.00000000i 1.22984070 + 0.00000000i [P.sub.25] -2.11390876 - 0.00000000i -0.91044961 - 0.00000000i [P.sub.26] -2.39860964 + 0.00000000i -0.45929305 - 0.00000000i [P.sub.27] -0.35998047 + 0.25179313i 0.61198862 - 0.60194230i [P.sub.28] -1.29601698 - 0.00000000i 1.19480441 + 0.00000000i [P.sub.29] 0.05975897 + 0.25434442i -0.06366391 + 0.36672723i [P.sub.30] -0.51418694 - 0.00000000i -2.16372657 - 0.00000000i [P.sub.31] -1.79027097 + 0.00000000i -0.78332475 - 0.00000000i [P.sub.32] -0.82449210 + 0.00000000i 1.84542750 + 0.00000000i [P.sub.33] -0.35998047 - 0.25179313i 0.61198862 + 0.60194230i [P.sub.34] -2.28499755 - 0.00000000i -0.46036727 + 0.00000000i [P.sub.35] 2.24284114 - 0.00000000i -0.46080446 - 0.00000000i [P.sub.36] 0.92084852 + 0.00000000i 2.10459508 + 0.00000000i [P.sub.37] -0.56599742 - 0.66867689i 0.95095432 + 0.12549545i [P.sub.38] 0.54660137 - 0.59723958i -0.36484540 + 0.17346068i [P.sub.39] 0.51269502 - 0.43085324i -0.51116222 - 0.08768075i [P.sub.40] 1.00000000 0.00000000 z Mult. Ord. [P.sub.1] 1.00000000 2 Yes [P.sub.2] 1.00000000 2 Yes [P.sub.3] 1.00000000 2 Yes [P.sub.4] 1.00000000 2 Yes [P.sub.5] 1.00000000 2 Yes [P.sub.6] 1.00000000 2 Yes [P.sub.7] 1.00000000 2 Yes [P.sub.8] 1.00000000 2 Yes [P.sub.9] 1.00000000 2 Yes [P.sub.10] 1.00000000 2 Yes [P.sub.11] 1.00000000 2 Yes [P.sub.12] 1.00000000 2 Yes [P.sub.13] 1.00000000 2 Yes [P.sub.14] 1.00000000 2 Yes [P.sub.15] 1.00000000 2 Yes [P.sub.16] 1.00000000 2 Yes [P.sub.17] 1.00000000 2 Yes [P.sub.18] 1.00000000 2 Yes [P.sub.19] 1.00000000 2 Yes [P.sub.20] 1.00000000 2 Yes [P.sub.21] 1.00000000 2 Yes [P.sub.22] 1.00000000 2 Yes [P.sub.23] 1.00000000 2 Yes [P.sub.24] 1.00000000 2 Yes [P.sub.25] 1.00000000 2 Yes [P.sub.26] 1.00000000 2 Yes [P.sub.27] 1.00000000 2 Yes [P.sub.28] 1.00000000 2 Yes [P.sub.29] 1.00000000 2 Yes [P.sub.30] 1.00000000 2 Yes [P.sub.31] 1.00000000 2 Yes [P.sub.32] 1.00000000 2 Yes [P.sub.33] 1.00000000 2 Yes [P.sub.34] 1.00000000 2 Yes [P.sub.35] 1.00000000 2 Yes [P.sub.36] 1.00000000 2 Yes [P.sub.37] 1.00000000 2 Yes [P.sub.38] 1.00000000 2 Yes [P.sub.39] 1.00000000 2 Yes [P.sub.40] 0.00000000 2 Yes

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Luo, Zhongxuan; Feng, Erbao; Zhang, Jielin |

Publication: | Discrete Dynamics in Nature and Society |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 5773 |

Previous Article: | Pricing scheme of ocean carrier for inbound container storage for assistance of container supply chain finance. |

Next Article: | Incorporating overconfidence into real option decision-making model of metal mineral resources mining project. |

Topics: |