07-Jan-2022 21:58:05

jacobi_rule_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Test jacobi_rule().
07-Jan-2022 21:58:05

JACOBI_RULE():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2

  Compute a Gauss-Jacobi rule for approximating
    Integral ( A <= x <= B ) (B-x)^alpha (x-A)^beta f(x) dx
  of order ORDER.

  The user specifies ORDER, ALPHA, BETA and FILENAME.

  ORDER is the number of points
  ALPHA is the exponent of (B-x):
  BETA is the exponent of (x-A);
  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 = 0.500000
  BETA = 1.500000
  A = -1.000000
  B = 1.000000
  FILENAME = "jac_o4_a0.5_b1.5".

  Creating quadrature files.

  "Root" file name is   "jac_o4_a0.5_b1.5".

  Weight file will be   "jac_o4_a0.5_b1.5_w.txt".
  Abscissa file will be "jac_o4_a0.5_b1.5_x.txt".
  Region file will be   "jac_o4_a0.5_b1.5_r.txt".

JACOBI_RULE:
  Normal end of execution.

07-Jan-2022 21:58:05

jacobi_rule_test():
  Normal end of execution.

07-Jan-2022 21:58:05