u = 2 * pi * sort (rand (1, 2^7));
xy = diag([1; (1 + sqrt(5) / 2)]) * [cos(u); sin(u)];
for p = 1:columns(xy)
  printf("Point(%d) = {%g, %g, 0};\n", p, xy(:, p))
  if p > 1	  
    printf("Line(%d) = {%d, %d};\n", p - 1, p - 1, p)
  end%if	  
end%for
printf("Line(%d) = {%d, 1};\n", p, p)
