constants <- function ( )

#*****************************************************************************80
#
## constants() stores, and returns constants "g", "c" and "h".
#
#  Discussion:
#
#    Calling [g,c,h]=constants() returns the values of g, c, and h.
#
#    Because the values never change, and don't need to be computed,
#    we use assignment statements here, instead of persistent data.
#
#  Licensing:
#
#    This code is distributed under the MIT license.
#
#  Modified:
#
#    08 May 2021
#
#  Author:
#
#    John Burkardt
#
#  Output:
#
#    real G: gravitational constant m^3/s^2/kg
#
#    real C: light speed, m/s.
#
#    real H: Planck's constant, j s;
#
{
  g <- 6.67384e-11
  c <- 2.99792458e+8
  h <- 6.626070040e-34

  output <- c ( g, c, h ) 

  return ( output )
}