#include "supnormtest.sollya";

mode=relative;
I = [-1b-2;1b-2];
f = expm1(x);
p = x * (1 +  x * ( 2097145 * 2^(-22) + x * ( 349527 * 2^(-21) + x * (87609 * 2^(-21) + x * 4369 * 2^(-19))))); 
quality = 2^(-37.6);

testSupnorm(p, f, I, quality, mode);