prec = 160;
eps = 1b-24;
f = 2^x;
I = [0; 0.5];
P = fpminimax(f, 5, [|1,SG...|], I, floating, absolute);
print(P);
print(supnorm(P, f, I, absolute, eps));
print("P_POW2 = ");
printexpansion(P);
f = 10^x;
I = [0; 0.5];
P = fpminimax(f, 5, [|1,SG...|], I, floating, absolute);
print(P);
print(supnorm(P, f, I, absolute, eps));
print("P_POW10 = ");
printexpansion(P);
quit;