07-Jan-2022 20:33:41

gram_polynomial_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Test gram_polynomial()

gram_abscissas_test:
  Test gram_abscissas()

  x = gram_abscissas(3)

x =

   -0.6667
         0
    0.6667


  x = gram_abscissas(4)

x =

   -0.7500
   -0.2500
    0.2500
    0.7500


  x = gram_abscissas(5)

x =

   -0.8000
   -0.4000
         0
    0.4000
    0.8000


  x = gram_abscissas(10)

x =

   -0.9000
   -0.7000
   -0.5000
   -0.3000
   -0.1000
    0.1000
    0.3000
    0.5000
    0.7000
    0.9000


gram_beta_test:
  Test gram_beta()

       M =          0         1         2         3         4         5
  N =  0:    --------   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000
  N =  1:    --------  --------    0.2500    0.2963    0.3125    0.3200
  N =  2:    --------  --------  --------    0.1481    0.2000    0.2240
  N =  3:    --------  --------  --------  --------    0.1125    0.1646
  N =  4:    --------  --------  --------  --------  --------    0.0914
  N =  5:    --------  --------  --------  --------  --------  --------

gram_polynomial_coefficients_test:
  Test polynomial_coefficients_beta()

Gram polynomial g(0,10) = 
     +1.000000

Gram polynomial g(1,10) = 
     +1.000000 * x

Gram polynomial g(2,10) = 
     +1.000000 * x^2
     -0.330000

Gram polynomial g(3,10) = 
     +1.000000 * x^3
     -0.586000 * x

Gram polynomial g(4,10) = 
     +1.000000 * x^4
     -0.820000 * x^2
     +0.077220

Gram polynomial g(5,10) = 
     +1.000000 * x^5
     -1.033333 * x^3
     +0.202233 * x
  Graphics saved as "gram_0_10.png"
  Graphics saved as "gram_1_10.png"
  Graphics saved as "gram_2_10.png"
  Graphics saved as "gram_3_10.png"
  Graphics saved as "gram_4_10.png"
  Graphics saved as "gram_5_10.png"

A =

    1.0000    0.0000    0.0000    0.0000   -0.0000    0.0000
    0.0000    0.3328    0.0000    0.0000    0.0000    0.0000
    0.0000    0.0000    0.0882    0.0000    0.0000    0.0000
    0.0000    0.0000    0.0000    0.0223    0.0000    0.0000
   -0.0000    0.0000    0.0000    0.0000    0.0055    0.0000
    0.0000    0.0000    0.0000    0.0000    0.0000    0.0013


c =

    1.1751
    1.1035
    0.5367
    0.1761
    0.0436
    0.0087
    0.0014
    0.0002
    0.0000
    0.0000
    0.0000


R8POLY_PRINT_TEST
  R8POLY_PRINT prints an R8POLY.

  The R8POLY: = 
     +9.000000 * x^5
     +0.780000 * x^4
    +56.000000 * x^2
     -3.400000 * x
    +12.000000

  The R8POLY: = 
     +0.780000 * x^4
    +56.000000 * x^2
     -3.400000 * x
    +12.000000

  The R8POLY: = 
     +1.000000 * x^5
     +0.780000 * x^4
    +56.000000 * x^2
    +12.000000

  The R8POLY: = 
  0

gram_polynomial_test():
  Normal end of execution.
07-Jan-2022 20:34:17