Numerically evaluating a theta function