#*****************************************************************************80
#
## geyser_scatter creates a scatter plot from geyser data.
#
#  Licensing:
#
#    This code is distributed under the MIT license. 
#
#  Modified:
#
#    27 February 2020
#
#  Author:
#
#    John Burkardt
#
cat ( date ( ), "\n" )
cat ( "\n" )
cat ( "geyser_scatter\n" )
cat ( " ", version$version.string, "\n" )
cat ( "  Make a scatter plot of eruption length versus subsequent" )
cat ( "  waiting time at Old Faithful." )

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

filename = "geyser_scatter.png"

png ( filename )
plot ( data[,1], data[,2], 
  xlab = "Eruption Duration (minutes)", 
  ylab = "Pause Duration (minutes)",
  main = "Duration of Eruption / Pause at Old Faithful" )
grid ( )
cat ( '  Graphics saved as "', filename, '"\n' )
#
#  Terminate.
#
cat ( "\n" )
cat ( "geyser_scatter:\n" )
cat ( "  Normal end of execution.\n" )

quit ( )