function hw1q3()
  % Data from running the allTails code in Erlang
  data = [ ...  N   size(allTails(lists:seq(1,N)))
           [  100,   402];
           [  200,   802];
           [  300,  1202];
           [  400,  1602];
           [  500,  2002];
           [  600,  2402];
           [  700,  2802];
           [  800,  3202];
           [  900,  3602];
           [ 1000,  4002] ];
  n = [0,1000];
  lobf = plot(n, 4*n+2, 'k-');
  hold('on'); % plot the data points on top of the line of best fit
  plot(data(:,1), data(:,2), 'b*');
  xlabel('N');
  ylabel('size(allTails(seq(1,N)))');
  legend(lobf, {'line of best fit: $4N+2$'}, 'Interpreter', 'latex', 'FontSize', 12);
  hold('off');
end % hw1q3
