Task: Get a parametrization of a plane curve singularity via Puiseux expansion.

Curve C={x3-y2=0} LIB "hnoether.lib"; ring r=0,(x,y),ls; poly f=x^3-y^2; list hn=develop(f,6); param(hn); ==> _[1]=x2 _[2]=x3

t |--> (t2,t3)