prec = 160;
eps = 1b-24;
f = log1p(x);
I = [1/sqrt(2)-1; sqrt(2)-1];
P = fpminimax(f, [|3, 4, 5, 6, 7|], [|SG...|], I, floating, absolute, x - 0.5*x^2);
print(P);
print(supnorm(P, f, I, absolute, eps));
print("P = ");
printexpansion(P);
print("E = ");
printsingle(exp(1));
print("LN_2 = ");
printsingle(log(2));
print("SQRT_2 = ");
printsingle(sqrt(2));
print("LOG2_E = ");
printsingle(log2(exp(1)));
print("LOG10_E = ");
printsingle(log10(exp(1)));
quit;