An algorithm for computing mixed sums of products of Bernoulli polynomials and Euler polynomials. (January 2015)