amplitudes = 0.25:0.25:1.25;
frequencies = (2:2:10);
t_max = 10;
t_step = 0.01;
    
figure()
hold on
 
for i = 1:length(amplitudes)
   for j = 1:length(frequencies)
      [x_data, y_data] = calculate_sinewave(frequencies(j), ...
                                            amplitudes(i), t_max, t_step);
      plot_sinewave(x_data, y_data, sprintf('freq: %5.2f, ampl: %5.2f',...
                                            frequencies(j), amplitudes(i)))
   end
end
legend('show')