#*****************************************************************************80
#
## cubicspline_test tests cubicspline.
#
#  Licensing:
#
#    Copyright 2016 James P. Howard, II
#
#    The computer code and data files on this web page are distributed under
#    https://opensource.org/licenses/BSD-2-Clause, the BSD-2-Clause license.
#
#  Modified:
#
#    23 February 2020
#
#  Author:
#
#    John Burkardt
#
cat ( date ( ), "\n" )
cat ( "\n" )
cat ( "cubicspline_test\n" )
cat ( " ", version$version.string, "\n" )
cat ( "  cubicspline computes a cubic spline to data.\n" )

source ( "/home/burkardt/public_html/r_src/cubicspline/cubicspline.R" )  

cat ( "\n" )

x <- c ( -2, -1, 0, 1 )
cat ( "\n" )
cat ( "  x coordinates:\n" )
cat ( "\n" )
print ( x )

y <- c ( -1, -2, -1, 2 )
cat ( "\n" )
cat ( "  y coordinates:\n" )
cat ( "\n" )
print ( y )

abcd <- cubicspline ( x, y )

cat ( "\n" )
cat ( "  a vector:\n" )
cat ( "\n" )
print ( abcd$a )

cat ( "\n" )
cat ( "  b vector:\n" )
cat ( "\n" )
print ( abcd$b )

cat ( "\n" )
cat ( "  c vector:\n" )
cat ( "\n" )
print ( abcd$c )

cat ( "\n" )
cat ( "  d vector:\n" )
cat ( "\n" )
print ( abcd$d )
#
#  Terminate.
#
cat ( "\n" )
cat ( "cubicspline_test\n" )
cat ( "  Normal end of execution.\n" )
cat ( date ( ), "\n" )

quit ( )