Totally interpolating biorthogonal multiwavelet systems with approximation order