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.