The partial derivatives with respect to β j are taken using the computer algebra package Mathematica. For all combinations of angular momenta l i , l j = 0, 1, 2, 3 (up to f-functions) C code for the upper bounds is generated automatically. The integral routines were compared with numerical integration on a grid. The second part of the upper bound is the norm of the polarization operator