Summary: This paper considers the construction of 3-band compactly supported symmetric biorthogonal wavelets. Given an interpolatory refinable function, we present an explicit iterative algorithm for constructing its dual scaling function with desired regularity, and discuss the properties of the constructed scaling function, such as its symmetry and decay. Finally, several design examples are given.