#*****************************************************************************80
#
## tridiagmatrix_test tests tridiagmatrix.
#
#  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:
#
#    09 February 2020
#
#  Author:
#
#    John Burkardt
#
cat ( date ( ), "\n" )
cat ( "\n" )
cat ( "tridiagmatrix_test\n" )
cat ( " ", version$version.string, "\n" )
cat ( "  tridiagmatrix(A) solves a linear system with a tridiagonal matrix.\n" )

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

L <- c ( 4, 2, 3 )
D <- c ( 3, 5, 5, 5 )
U <- c ( 4, 2, 3 )
b <- c ( 20, 28, 18, 18 )
cat ( "\n" )
cat ( "The lower diagonal L:\n" )
cat ( "\n" )
L
cat ( "\n" )
cat ( "The diagonal D:\n" )
cat ( "\n" )
D
cat ( "\n" )
cat ( "The upper diagonal U:\n" )
cat ( "\n" )
U
cat ( "\n" )
cat ( "The right hand side b:\n" )
cat ( "\n" )
b

x <- tridiagmatrix ( L, D, U, b )

cat ( "\n" )
cat ( "The computed solution x:\n" )
cat ( "\n" )
x
#
#  Terminate.
#
cat ( "\n" )
cat ( "tridiagmatrix_test\n" )
cat ( "  Normal end of execution.\n" )
cat ( date ( ), "\n" )

quit ( )