31-Dec-2022 21:53:40

gegenbauer_rule_test():
  MATLAB/Octave version 4.2.2
  Test gegenbauer_rule().
31-Dec-2022 21:53:40

gegenbauer_rule():
  MATLAB/Octave version 4.2.2

  Compute a Gauss-Gegenbauer rule for approximating
    Integral ( A <= x <= B ) ((x-A)(B-X))^ALPHA f(x) dx
  of order ORDER.

  The user specifies ORDER, ALPHA, A, B, and FILENAME.

  ORDER is the number of points:
  ALPHA is the exponent:
  A is the left endpoint;
  B is the right endpoint;
  FILENAME is used to generate 3 files:
    filename_w.txt - the weight file
    filename_x.txt - the abscissa file.
    filename_r.txt - the region file.

  ORDER = 4
  ALPHA = 2.000000
  A = -1.000000
  B = 1.000000
  FILENAME = "gegen_o4_a2.0".

  Creating quadrature files.

  "Root" file name is   "gegen_o4_a2.0".

  Weight file will be   "gegen_o4_a2.0_w.txt".
  Abscissa file will be "gegen_o4_a2.0_x.txt".
  Region file will be   "gegen_o4_a2.0_r.txt".

GEGENBAUER_RULE:
  Normal end of execution.

31-Dec-2022 21:53:40

gegenbauer_rule_test():
  Normal end of execution.

31-Dec-2022 21:53:40