p(z) = z^3-1
pp(z) = 3*z^2
roots = apply (r -> exp(2*r*I*Pi/3), [0..2])
eps = 1e-7

A(x,y) = {
	if (x || y,
		my (z = x+. + (y+.)*I);
		while (1,
			for (r = 0, #roots-1,
				if (abs(z - roots[r+1])<eps,
					return (r);
				);
			);

			z -= p(z)/pp(z);
		),
		return (-1);
	);
}

{
	m = -1;
	for (d = 0, 140,
		for (k = 0, d,
			print (m++ " " A(d-k, k));
		);
	);
}

quit