This paper proposes a method for interpolating motion data with multilinear models. A third-order motion tensor is approximated by low-dimensional core tensor using high order tensor approximation. The resulting motion is accurately synthesized by interpolating the principal components using geostatistics. Our method reduces computational and storage cost of motion interpolation.