Using orthogonal polynomials (i.e., ultraspherical polynomials), the authors construct families of
orthogonal, compactly supported spline multiwavelets of
with various approximation orders. In the case of symmetric or antisymmetric multiscaling functions, this method allows the construction of symmetric or antisymmetric multiwavelets. By restriction on
, these results yield
spline multiwavelet bases of
compactly supported spline multiwavelet basis of
is also sketched, but the formulas become very complicated.