function project_02_part01c(fs, rotation_time1, rotation_time2, rotations1, rotations2, number_of_measurements)

%comment out fs, rotation_time, rotations, and number_of_measurements to specify values from command line
fs = 44100;                     %sampling rate
rotation_time1 = 2;             %time it takes in seconds for a full rotation
rotation_time2 = 1;             %time it takes in seconds for a full rotation
rotations1 = 2;                 %number of rotations to perform (use whole numbers only!)
rotations2 = 4;                 %number of rotations to perform (use whole numbers only!)
elev = 0;                       %elevation angle remains 0 for this simulation
azim = 0;                       %azimuth angle of 0 = directly in front of observer; positive degrees rotate clockwise
number_of_measurements = 37;    %37 HRTF measurements are provided, from 0 to 180 degrees, in 5 degree increments

signal1 = project_02_part01a(fs,rotation_time1,rotations1,number_of_measurements);
signal2 = project_02_part01b(fs,rotation_time2,rotations2,number_of_measurements);

%pads shorter signal with zeros to make signals the same length
pad_length = abs(length(signal1)-length(signal2));

if length(signal1) > length(signal2)
    signal3 = signal1 + [signal2;zeros(pad_length,2)];
else
    signal3 = [signal1;zeros(pad_length,2)] + signal2;
end

subplot(211), plot(signal3(:,1),'Color',[1,0.12,0.12]), grid on, axis tight;
title('Left channel'), xlabel('time (samples)'), ylabel('amplitude');
subplot(212), plot(signal3(:,2),'g'), grid on, axis tight;
title('Right channel'), xlabel('time (samples)'), ylabel('amplitude');

%double max(max()) needed, because max() on a matrix returns a row vector
signal3 = signal3/max(max(abs(signal3))); %normalizes input

%print project_02_part01c -dpng -r100;
%wavwrite(signal3,fs,'project_02_part01c');

soundsc(signal3,fs);