15-May-2025 06:25:21

gram_polynomial_test():
  MATLAB/Octave version 6.4.0
  Test gram_polynomial()

gram_abscissas_test():
  Test gram_abscissas()

  x = gram_abscissas(3)
x =

  -6.6667e-01
  -5.5511e-17
   6.6667e-01


  x = gram_abscissas(4)
x =

  -0.7500
  -0.2500
   0.2500
   0.7500


  x = gram_abscissas(5)
x =

  -8.0000e-01
  -4.0000e-01
   5.5511e-17
   4.0000e-01
   8.0000e-01


  x = gram_abscissas(10)
x =

  -0.900000
  -0.700000
  -0.500000
  -0.300000
  -0.100000
   0.100000
   0.300000
   0.500000
   0.700000
   0.900000


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   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   0.0000   0.0000   0.0000   0.0055   0.0000
   0.0000   0.0000   0.0000   0.0000   0.0000   0.0013

c =

   1.1751e+00
   1.1035e+00
   5.3666e-01
   1.7611e-01
   4.3589e-02
   8.6573e-03
   1.4355e-03
   2.0427e-04
   2.5454e-05
   2.8211e-06
   2.8152e-07


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.
15-May-2025 06:25:52