prec = 160;
eps = 1b-24;
f = exp(x);
I = [-log(2)/2; log(2)/2];
P = fpminimax(f, [|3, 4, 5, 6, 7|], [|SG...|], I, floating, absolute, 1 + 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("LN_2_INV = ");
printsingle(1/log(2));
quit;