#include "supnormtest.sollya";

mode=relative;
I=[-1/8;1/8];
f=sin(x)/(exp(x)-1);

print("Computing remez polynomial...");
p = remez(f, 15, [-1/2, 1/2], 1, 1b-20);
print("... done.");

quality = 2^(-16);

removableDiscontinuitiesOfF = [|0|];
testSupnorm(p, f, I, quality, mode);