Hadamard product for the Riemann zeta function
What is the best way to compute numerically with Sage the product over the first n nontrivial zeros of the Riemann zeta function?
p(s,n) = product((1-s/rho(k))*exp(s/rho(k)) for k in (1..n))
where zeta(rho(k)) == 0 and Im(rho(k)) != 0.
Wikipedia advises: "To ensure convergence the product should be taken over 'matching pairs' of zeroes, i.e. the factors for a pair of zeroes of the form rho(k) and 1-rho(k) should be combined."