Analytic results for planar three-loop integrals for massive form factors

We use the method of differential equations to analytically evaluate all planar three-loop Feynman integrals relevant for form factor calculations involving massive particles. Our results for ninety master integrals at general $q^2$ are expressed in terms of multiple polylogarithms, and results for fiftyone master integrals at the threshold $q^2=4m^2$ are expressed in terms of multiple polylogarithms of argument one, with indices equal to zero or to a sixth root of unity.