function p = plaid(n, fM)

t = 2 * pi * fM * linspace(0, 1, n) .^ 2;
s = sin(t);
p = s' * s;