function bulgaria_population ( )

%*****************************************************************************80
%
%% bulgaria_population() plots the population of Bulgaria over time.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    11 April 2018
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'bulgaria_population():\n' );
  fprintf ( 1, '  MATLAB/Octave version %s\n', version ( ) );
  fprintf ( 1, '  Plot population of Bulgaria over time.\n' );

  data = csvread ( 'bulgaria_population.csv', 1 ) ;

  plot ( data(:,1), data(:,2), 'linewidth', 3 )
  grid on
  xlabel ( '<-- Year -->', 'Fontsize', 16 )
  ylabel ( '<-- Population -->', 'Fontsize', 16 )
  title ( 'Population of Bulgaria Over Time', 'Fontsize', 16 )
%
%  Force the Y-axis to run from 0 to 10,000,000
%
  ylim ( [ 0, 10000000 ] )
%
%  Save graphics in a file.
%
  filename= 'bulgaria_population.png';
  print ( '-dpng', filename );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  Graphics saved as "%s"\n', filename );
%
%  Terminate.
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'bulgaria_population():\n' );
  fprintf ( 1, '  Normal end of execution.\n' );

  return
end

