function florida_cvt_geo_test ( )

%*****************************************************************************80
%
%% FLORIDA_CVT_GEO_TEST tests FLORIDA_CVT_GEO.
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    29 June 2016
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'FLORIDA_CVT_GEO_TEST:\n' );
  fprintf ( 1, '  FLORIDA_CVT_GEO can estimate the geometric centroidal\n' );
  fprintf ( 1, '  Voronoi tessellation (CVT) of a set of points\n' );
  fprintf ( 1, '  within the boundaries of Florida,\n' );
  fprintf ( 1, '  expressed as longitudes and latitudes.\n' );
%
%  Set the data.
%
  names = [ ...
  'Pensacola', ...
  'Tallahassee', ...
  'Jacksonville', ...
  'Tampa', ...
  'Gainesville', ...
  'Orlando', ...
  'Miami' ];

  gen_lat = [ ...
   30.43, ...
   30.46, ...
   30.34, ...
   27.97, ...
   29.65, ...
   28.42, ...
   25.78 ];

  gen_lon = [ ...
    -87.20, ...
    -84.25, ...
    -81.66, ...
    -82.48, ...
    -82.32, ...
    -81.30, ...
    -80.21 ];
%
%  Estimate the CVT.
%
  [ gen_lon, gen_lat ] = florida_cvt_geo ( gen_lon, gen_lat );
%
%  Estimate the centroids.
%
  sample_num = 5000;
  seed = 123456789;
  [ cen_lon, cen_lat, seed ] = florida_centroid_geo ( gen_lon, ...
    gen_lat, sample_num, seed );
%
%  Show generators, centroids, and Voronoi region together.
% 
  filename = 'florida_cvt_geo_test.png';
  plot_title = 'Estimated 7-Generator geometric CVT for Florida';
  florida_voronoi_display ( gen_lon, gen_lat, cen_lon, cen_lat, filename, plot_title );

  return
end

