function [a,b] = generate_golay(N)
% [a,b] = generate_golay(N)
%
% Generate the Golay codes a and b with length 2^N.
%
% Then write them to disk as golayA.wav and golayB.wav.


% These initial a and b values are Golay
a = [1 1];
b = [1 -1];

% Iterate to create a longer Golay sequence
while (N>1)
    olda = a;
    oldb = b;
    a = [olda oldb];
    b = [olda -oldb];

    N = N - 1;
end

% Guess the sampling rate. It doesn't really matter.
fs = 44100;

% Scaling by 0.9999 suppresses a warning message about clipping.
wavwrite(a*0.9999,fs,'golayA.wav');
wavwrite(b*0.9999,fs,'golayB.wav');