36 def main(n, alpha, beta, meth):
38 roots, weights = gjquadpy.gaussjacobiquadccompat.gauss_jacobi_rec_c(
42 roots, weights = gjquadpy.gaussjacobiquadccompat.gauss_jacobi_gw_c(
45 elif meth ==
"algo665":
46 roots, weights = gjquadpy.gaussjacobiquadccompat.gauss_jacobi_algo665_c(
50 for idx, root
in enumerate(roots):
51 sign =
" " if root >= 0
else ""
52 root_str = f
"{root:23.17E}"
53 weight_str = f
"{weights[idx]:23.17E}"
54 print(f
"Root:{sign} {root_str} Weight: {weight_str}")
int main(int argc, char *argv[])