function make_sine_table(filename,N)
% make_sine_table.m 
% Generates N sample points of 1 period of a sine wave
% Creates the header file filename.h
% The argument filename must be a string (e.g. 'myfile')

sine = round(20000*sin(2*pi*[0:N-1]./N)); % 8000/N Hz sine wave

fid = fopen([filename '.h'],'w');         % Open/Create file 'filename'
fprintf(fid,['short sine_table[' num2str(N) ']={']);    

for i=1:N-1
    fprintf(fid,'%d,' ,sine(i));          % Print N-1 points
    if mod(i,8) == 0
        fprintf(fid,'\n');
    end
end
    
fprintf(fid,'%d' ,sine(N));               % Print Nth point
fprintf(fid,'};\n');                      % Print closing bracket
fclose(fid);						      % Close file