16-Aug-2022 10:00:28

r83_np_test():
  MATLAB version
  Test r83_np().

R83_NP_DET_TEST
  R83_NP_DET computes the determinant of a tridiagonal
  matrix factored by R83_NP_FA.

  Matrix order N = 10

  The factored R83 matrix:

  Col:       1             2             3             4             5       
  Row
  ---
    1:             2            -1                                          
    2:          -0.5           1.5            -1                            
    3:                   -0.666667       1.33333            -1              
    4:                                     -0.75          1.25            -1
    5:                                                    -0.8           1.2
    6:                                                             -0.833333

  Col:       6             7             8             9            10       
  Row
  ---
    5:            -1                                                        
    6:       1.16667            -1                                          
    7:     -0.857143       1.14286            -1                            
    8:                      -0.875         1.125            -1              
    9:                                 -0.888889       1.11111            -1
   10:                                                    -0.9           1.1

  R83_NP_DET computes determinant = 11
  Exact determinant = 11

R83_NP_FA_TEST
  R83_NP_FA factors a tridiagonal matrix with no pivoting;

  Matrix order N = 10

  The tridiagonal matrix:

  Col:       1             2             3             4             5       
  Row
  ---
    1:      0.218418      0.829509                                          
    2:      0.956318      0.561695     0.0661187                            
    3:                    0.415307      0.257578      0.043829              
    4:                                  0.109957      0.633966      0.449539
    5:                                               0.0617272      0.401306
    6:                                                              0.754673

  Col:       6             7             8             9            10       
  Row
  ---
    5:      0.797287                                                        
    6:    0.00183837      0.350752                                          
    7:      0.897504     0.0945448      0.859097                            
    8:                   0.0136169      0.840847    0.00751236              
    9:                                  0.123104      0.260303      0.113664
   10:                                                0.912484      0.351629

  Solution:

     1:            1
     2:            2
     3:            3
     4:            4
     5:            5
     6:            6
     7:            7
     8:            8
     9:            9
    10:           10

  Solution to tranposed system:

     1:            1
     2:            2
     3:            3
     4:            4
     5:            5
     6:            6
     7:            7
     8:            8
     9:            9
    10:           10

R83_NP_FS_TEST
  R83_NP_FS factors and solves a tridiagonal
  linear system.

  Matrix order N = 10

  Solution:

     1:            1
     2:            2
     3:            3
     4:            4
     5:            5
     6:            6
     7:            7
     8:            8
     9:            9
    10:           10

R83_NP_ML_TEST
  R83_NP_ML computes A*x or A'*x
  where A has been factored by R83_FA.

  Matrix order N = 10

  A*x and PLU*x:

     1:        1.87744         1.87744
     2:        2.27806         2.27806
     3:        1.77866         1.77866
     4:        5.11343         5.11343
     5:        7.03716         7.03716
     6:        6.23966         6.23966
     7:        12.9196         12.9196
     8:        6.88971         6.88971
     9:         4.4642          4.4642
    10:        11.7286         11.7286

  A'*x and (PLU)'*x

     1:        1.35715         1.35715
     2:        4.39904         4.39904
     3:        6.03519         6.03519
     4:        5.59777         5.59777
     5:        5.81525         5.81525
     6:        6.74571         6.74571
     7:        7.50541         7.50541
     8:        12.9291         12.9291
     9:        14.2607         14.2607
    10:        9.79448         9.79448

R83_NP_SL_TEST
  R83_NP_SL solves a linear system after the tridiagonal
  matrix has been factored by R83_NP_FA.

  Matrix order N = 10

  The tridiagonal matrix:

  Col:       1             2             3             4             5       
  Row
  ---
    1:      0.218418      0.829509                                          
    2:      0.956318      0.561695     0.0661187                            
    3:                    0.415307      0.257578      0.043829              
    4:                                  0.109957      0.633966      0.449539
    5:                                               0.0617272      0.401306
    6:                                                              0.754673

  Col:       6             7             8             9            10       
  Row
  ---
    5:      0.797287                                                        
    6:    0.00183837      0.350752                                          
    7:      0.897504     0.0945448      0.859097                            
    8:                   0.0136169      0.840847    0.00751236              
    9:                                  0.123104      0.260303      0.113664
   10:                                                0.912484      0.351629

  Solution:

     1:            1
     2:            2
     3:            3
     4:            4
     5:            5
     6:            6
     7:            7
     8:            8
     9:            9
    10:           10

  Solution to tranposed system:

     1:            1
     2:            2
     3:            3
     4:            4
     5:            5
     6:            6
     7:            7
     8:            8
     9:            9
    10:           10

r83_np_test():
  Normal end of execution.

16-Aug-2022 10:00:28