#! /usr/bin/gawk -f
BEGIN {
  # Number of samples:
  npt = 1000
  # Mean radius:
  rad = 100
  # Oscillation amplitude
  ampl1 = 0.12
  ampl2 = 0.04
  # Winding number:
  repl = 1
  unit = 0.01
  printf "begin PZLR3Chain.T (format of 97-10-29)\n"
  printf "| output of make-curve\n"
  printf "unit = %f\n", unit
  printf "samples = %d\n", npt
  for (i=0;i< npt;i++)
    { theta = 6.2831853071795864*i/npt;
      phi = 5 * theta;
      tau = 30 * theta;
      radius = rad*(1 + ampl1*cos(phi) + ampl2*sin(tau))
      x = int(radius*cos(repl*theta)/unit+0.5)
      y = int(radius*sin(repl*theta)/unit+0.5)
      printf(" %7d %7d  0\n", x, y)
    }
  printf "end PZLR3Chain.T\n"
  exit 0 
}