"Source: Text/Basics/FldRe.text";
"Line: 3756";
"Date: Thu Sep 25 14:47:46 2025";
"Main: Fri Sep 26 14:29:19 2025";
// original file: Text/Basics/FldRe.text, line: 3756
// Example: H26E12 ()
print "Example: H26E12";
ei := GetEchoInput();
SetEchoInput(true);
c := 0.128866348618482136010608730745;
N := 10; 
D := 30;
function f(x) 
   return Sin(x)*Cos(x); 
end function;
A, W := GaussJacobiIntegrationPoints(N,D,-1/3,-1/5);                                          
int_GJ := &+[ W[k] * f(A[k]) : k in [1..N] ];
Abs(int_GJ - c);
h := Real(2^-3); 
N := Ceiling(6/h);
A, W1, W2 := TanhSinhIntegrationPoints(N,h);
int_DE := &+[ W1[k]*W2[k]^(1/3)*(1+A[k])^(1/3-1/5)*f(A[k]) : k in [1..2*N+1] ];
Abs(int_DE - c);
SetEchoInput(ei);
