If I am not mistaken, the connection of multiple zeta values to renormalized Feynman integrals was discovered through computations (by David Broadhurst). Here is a survey of ongoing computer research around this field (from 1996).