Las funciones de frontera para wavelets en un intervalo finito suelen construirse como combinaciones lineales de funciones de escalado que cruzan la frontera. Un enfoque alternativo se basa en técnicas de álgebra lineal para truncar la matriz infinita de la Transformada Wavelet Discreta a una finita. En este artículo mostramos cómo un algoritmo de Madych para wavelets escalares puede generalizarse a multiwavelets, dada una condición adicional. Luego desarrollamos un nuevo algoritmo que no requiere esta condición adicional. Finalmente, aplicamos resultados de un artículo previo para resolver la no unicidad del algoritmo imponiendo condiciones de regularidad (incluidos órdenes de aproximación) en las funciones de frontera.