The real and imaginary part of the vacuum polarization function
Π(q2) induced by a massive quark is calculated in perturbative
QCD up to order α2s. The method is described and the
results are presented. This extends the calculation by
K\“all\'en and Sabry from two to three loops.