[SOLVED] Create NxN Haar Matrix

I can't find a definition for generating an un-normalized NxN Haar matrix. So what is the equation?

see: http://en.wikipedia.org/wiki/Haar_wavelet

Thx, Chris

Thx all. Wikipedia gives the 'equation':

$H_{2N} = \begin{bmatrix} H_{N} \otimes [1, 1] \\ I_{N} \otimes [1, -1] \end{bmatrix} $$ :where $$ I_{N} = \begin{bmatrix} 1 & 0 & \dots & 0 \\ 0 & 1 & \dots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \dots & 1 \end{bmatrix}$

I've wrote a recursive solution for generating an un-normalized NxN Haar matrix in octave.

function [h] = haar(n)
h = [1];
if n > 2
    h = haar(n/2);
endif
% calculate upper haar part
h_n = kron(h,[1,1]); 
% calculate lower haar part 
h_i = kron(eye(length(h)),[1,-1]);
% combine parts
h = [h_n; h_i];
endfunction

disp(haar(8));