Calculating four-loop tadpoles with one non-zero mass
  An efficient method to calculate tadpole diagrams is proposed. Its capability is demonstrated by analytically evaluating two four-loop tadpole diagrams of current interest in the literature, including their $O(\epsilon)$ terms in $D=4-2\epsilon$ space-time dimensions. 
Bernd A. Kniehl, Anatoly V. Kotikov