Let f be a squarefree polynomial in x and y. Then f defines an algebraic curve in the plane C^2, and also in the projective plane P^2 by making f homogeneous. This procedure computes the singular points of the curve in the projective plane. The points are given by homogeneous co-ordinates [X,Y,Z].

For each singularity this procedure also computes the multiplicity $m$, the delta invariant delta, and the number of local branches $r$. An ordinary double point is characterized by $m=2,\mathrm{\delta }=1,r=2$. For a cusp one has $m=2,\mathrm{\delta }=1,r=1$. In general $r\le m$ and $\frac{m\left(m-1\right)}{2}\le \mathrm{\delta }$, and both of these are equalities when the singularity is an ordinary $m$-multiple point. The Milnor number equals $2\mathrm{\delta }-r+1$.

The output of this procedure is a set consisting of lists of the following form $\left[\mathrm{point}\,m\,\mathrm{\delta }\,r\right]$.

This procedure computes all singularities up to conjugation. So if a singularity $\left[\mathrm{RootOf}\left({\mathrm{\_Z}}^2-2\right)\,1\,1\right]$ is given in the output, and if $\mathrm{RootOf}\left({\mathrm{\_Z}}^2-2\right)$ does not appear in the input, then $\left[-\mathrm{RootOf}\left({\mathrm{\_Z}}^2-2\right)\,1\,1\right]$ is a singular point as well but will not be given in the output.

The genus of a curve is the number (d-1)*(d-2)/2 - Sum(delta invariants) where $d$ is the degree of the curve. Note that if we apply this formula to compute the genus, then for each singularity we must multiply the delta invariant by the degree of the algebraic extension over which the singularity is defined, because only one singularity of each conjugacy class is given in the output.

$\mathrm{with}\left(\mathrm{algcurves}\right)\:$

$f≔-8y^5-207y{x}^4+180x^5-35y^4-128y^{3}x+621y{x}^3-450x^4+82y^3-521y{x}^2+369x^3-19y^2+135yx-100x^2-28y-7x+8$

$f≔180x^5-207y{x}^4-8y^5-450x^4+621y{x}^3-128y^{3}x-35y^4+369x^3-521y{x}^2+82y^3-100x^2+135yx-19y^2-7x-28y+8$

$\mathrm{degree}\left(f\,\left\{x\,y\right\}\right)$

$5$

$\mathrm{singularities}\left(f\,x\,y\right)$

$\left\{\left[\left[0\,1\,1\right]\,2\,1\,2\right]\,\left[\left[1\,\mathrm{-1}\,1\right]\,2\,1\,2\right]\,\left[\left[1\,0\,1\right]\,2\,1\,2\right]\,\left[\left[\mathrm{RootOf}\left({\mathrm{\_Z}}^2+1\right)\,\mathrm{RootOf}\left({\mathrm{\_Z}}^2+1\right)\,1\right]\,2\,1\,2\right]\right\}$