BOKASUN: a fast and precise numerical program to calculate the Master   Integrals of the two-loop sunrise diagrams
  We present the program BOKASUN for fast and precise evaluation of the Master Integrals of the two-loop self-mass sunrise diagram for arbitrary values of the internal masses and the external four-momentum. We use a combination of two methods: a Bernoulli accelerated series expansion and a Runge-Kutta numerical solution of a system of linear differential equations. 
Michele Caffo, Henryk Czyz, Michal Gunia and Ettore Remiddi