Thu May 25 21:07:58 2023

jacobi_rule_test():
  Python version: 3.8.10
  Test jacobi_rule().

jacobi_rule():
  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.5
  BETA =  1.5
  A =  -1.0
  B =  1.0
  FILENAME = "jac_o4_a0.5_b1.5".
  abscissas saved as "jac_o4_a0.5_b1.5_x.txt"
  weights saved as "jac_o4_a0.5_b1.5_w.txt"
  region saved as "jac_o4_a0.5_b1.5_r.txt"

jacobi_rule_test():
  Normal end of execution.
Thu May 25 21:07:58 2023