A Hybrid Symbolic-Numeric Method for Multiple Integration Based on Tensor-Product Series Approximations