LEO is a FORM code for the calculation of the three-loop single scale bubble type diagrams. The full set of solved recurrence relations is presented in hep-ph/9512442 . The one of application - the calculation of the QCD corrections to the electrowek rho-parameter within SM - is done in hep-ph/9406363 .
The construction of the epsilon-expansion
for diagrams relating to QED/QCD problems
(in particular, BN(0,0,1,1,1,1))
was investigated by David Broadhurst
in Z.Phys.C54 (1992) 599-606
and
hep-th/9604128 .
The finite parts of the three-loop bubble with sixth lines
integrals were presented by David Broadhurst in
hep-th/9803091
First several terms of epsilon expansion for diagrams
D5(1,1,1,1,1,0), D3(0,1,0,1,1,1) and E3(1,1,1,1,1)
and for two-loop bubble with theee equal masses
are evaluated by Fleischer,Kalmykov
in
hep-ph/9910223
The all-order epsilon-expansion of the
two-loop bubble diagram with arbitrary non-zero masses have been
constructed by Davydychev in
hep-ph/9910224 .
Particular case, when one of the mass iz zero can be extracted from
Davydychev,Kalmykov
hep-th/0012189 .
Finite parts
for D3(0,1,0,1,1,1) and E3(1,1,1,1,1)
was evaluated by Chetyrkin,Steinhauser in
hep-ph/9911434
The higher order epsilon-expansion for D4(1,1,1,1,1,1),
D3(0,1,0,1,1,1) and E3(1,1,1,1,1)
have been calculated by Davydychev,Kalmykov in
hep-th/0012189 .
The higher order epsilon-expansion for D5(1,1,1,1,1,0),
have been calculated by Kalmykov in
hep-ph/0503070 .
Independently of analytical calculations,
the high-precision numerical values of higher order coefficients of
epsilon-expansion for the some of three-loop bubble-integrals
have been calculated by S.Laporta in
hep-ph/0102033
and all three-loop master-integrals are evaluated
by Y.Schroder and A. Vuorinen in
hep-ph/0503209 .