#*****************************************************************************80
#
## geyser_histogram makes a histogram of geyser data.
#
#  Licensing:
#
#    This code is distributed under the MIT license. 
#
#  Modified:
#
#    06 August 2020
#
#  Author:
#
#    John Burkardt
#
cat ( date ( ), "\n" )
cat ( "\n" )
cat ( "geyser_histogram\n" )
cat ( " ", version$version.string, "\n" )
cat ( "  Create a histogram of the time between eruptions" )
cat ( "  of the Old Faithful geyser." )

filename = "geyser_data.txt"
data <- read.table ( filename, header = FALSE )
print ( data )

filename = "geyser_histogram.png"

png ( filename )
hist ( data[,1], 
  xlab = "Minutes", 
  ylab = "Frequency",
  main = "Time between eruptions of Old Faithful",
  col = "red" )
grid ( )
cat ( '  Graphics saved as "', filename, '"\n' )
#
#  Terminate.
#
cat ( "\n" )
cat ( "geyser_histogram:\n" )
cat ( "  Normal end of execution.\n" )

quit ( )