  function [cVec,sVec] = csVecs(n)
  thetas = (2*pi/n)*(0:n-1)';
  cVec = cos(thetas);
  sVec = sin(thetas);