"Source: Text/Basics/FldRe.text";
"Line: 3717";
"Date: Thu Sep 25 14:47:46 2025";
"Main: Fri Sep 26 14:29:19 2025";
// original file: Text/Basics/FldRe.text, line: 3717
// Example: H26E11 ()
print "Example: H26E11";
ei := GetEchoInput();
SetEchoInput(true);
N := 20; 
D := 50;
function f(x) 
   return x*Log(1+x); 
end function;
A_GL, W_GL := GaussLegendreIntegrationPoints(N,D);
/* Map abscissas to the interval [0,1] */
A_GL := [ (1/2)*(a+1) : a in A_GL ];
int_GL := (1/2) * &+[ W_GL[k] * f(A_GL[k]) : k in [1..N] ];
Abs(int_GL - 1/4);
A_CC, W_CC := ClenshawCurtisIntegrationPoints(N,D);
A_CC := [ (1/2)*(a+1) : a in A_CC ];
int_CC := (1/2) * &+[ W_CC[k] * f(A_CC[k]) : k in [1..N+1] ];
Abs(int_CC - 1/4);
h := RealField(D)!6/N;
A_DE, W_DE := TanhSinhIntegrationPoints(N,h);
A_DE := [ (1/2)*(a+1) : a in A_DE ];
int_DE := (1/2) * &+[ W_DE[k] * f(A_DE[k]) : k in [1..2*N+1] ];
Abs(int_DE - 1/4);
SetEchoInput(ei);
