We present a powerful method to compute one-scale master integrals
using differential equations combined with numerical matching.
The latter is achieved with the help of deep expansions around
regular and singular kinematic points.
As an applications we consider massive three loop form factors.
We show that the method can be applied to systems of
differential equations involving several hundred master integrals.