function [h0,h1,h2,h3] = cardbal2

% cardbal2 : cardinal 2-balanced orthogonal multiwavelet filters 
%
% [h0,h1,h2,h3] = cardbal2;
% ShowMwlets(h0,h1,h2,h3);

A = -1/8+1/32*15^(1/2);
% A = -1/8-1/32*15^(1/2);

g = [1/32 A+1/4 15/16 -2*A-1/4 1/32 A];

g0 = g;
g1 = g(6:-1:1).*((-1).^(1:6));

K0 = 3+1;
K1 = 5+1;

h0 = up(g0,2); 
h1 = up(g1,2); 
h2 = -h0; 
h3 = -h1; 
h0(K0) = 1;
h1(K1) = 1;
h2(K0) = 1;
h3(K1) = 1;

r = sqrt(2);
h0 = h0/r;
h1 = h1/r;
h2 = h2/r;
h3 = h3/r;