Approximating the approximant: a numerical code for polynomial compression of discrete integral operators