We consider the solutions of refinement equations written in the form $$ \begin{array}{*{20}c} {{\varphi {\left( x \right)} = {\sum\limits_{\alpha \in \mathbb{Z}^{s} } {a{\left( \alpha \right)}\varphi {\left( {Mx - \alpha } \right)} + g{\left( x \right)}} },}} & {{x \in \mathbb{R}^{s} ,}} \\ \end{array} $$ where the vector of functions ϕ = (ϕ 1, ..., ϕ r ) T is unknown, g is a given vector of compactly supported functions on ℝ s , a is a finitely supported sequence of r × r matrices called the refinement mask, and M is an s × s dilation matrix with m = |detM|. Inhomogeneous refinement equations appear in the construction of multiwavelets and the constructions of wavelets on a finite interval. The cascade algorithm with mask a, g, and dilation M generates a sequence ϕ n , n = 1, 2, ..., by the iterative process $$ \begin{array}{*{20}c} {{\varphi _{n} {\left( x \right)} = {\sum\limits_{\alpha \in \mathbb{Z}^{s} } {a{\left( \alpha \right)}\varphi _{{n - 1}} {\left( {Mx - \alpha } \right)} + g{\left( x \right)}} },}} & {{x \in \mathbb{R}^{s} ,}} \\ \end{array} $$ from a starting vector of function ϕ 0. We characterize the L p -convergence (0 < p < 1) of the cascade algorithm in terms of the p-norm joint spectral radius of a collection of linear operators associated with the refinement mask. We also obtain a smoothness property of the solutions of the refinement equations associated with the homogeneous refinement equation.