base_freq = 10;

%declare the number of harmonic frequencies that you want to control, where
%each one is a multiple of the base frequency (e.g., 20Hz).
num_harmonics = 2;

T = 1/base_freq;  %period
%setup figure to plot data so main loop below runs faster
figure(1);
subplot(1,2,1)
hold on
set(gca,'XLim',[0 T])
xlabel('time (s)')
ylabel('acceleration in g or possibly (m/s^2)')
subplot(1,2,2)
hold on
set(gca,'XLim',[0 T])
xlabel('time (s)')
ylabel('voltage (V)')