#! /usr/bin/gawk -f
BEGIN {
  # Number of samples:
  # npt = 4
  npt = 1000
  # Mean radius:
  rad = 100
  # Spike amplitude
  ampl = 100
  # ampl = 10
  unit = 0.01
  printf "begin PZLR3Chain.T (format of 97-10-29)\n"
  printf "| output of make-dirac\n"
  printf "unit = %f\n", unit
  printf "samples = %d\n", npt
  for (i=0;i< npt;i++)
    { theta = 6.2831853071795864*i/npt;
      r = rad;
      if (2*i == npt) { r -= ampl; }
      x = int(r*cos(theta)/unit+0.5)
      y = int(r*sin(theta)/unit+0.5)
      printf(" %7d %7d  0\n", x, y)
    }
  printf "end PZLR3Chain.T\n"
  exit 0 
}