We report on our recent calculation of massive quark form factors using
a semi-numerical approach based on series expansions of the master integrals
around singular and regular kinematic points and numerical matching.
The methods allows to cover the whole kinematic range of negative and
positive values of the virtuality $s$ with at least seven significant digits accuracy.