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

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


A = 3/1280+1/2560*31^(1/2);
% OR
% A = 3/1280-1/2560*31^(1/2);
g = [
2*A+1/512
-A+1/32
-6*A+125/512
3*A+15/16
6*A-125/512
-3*A+1/32
-2*A-1/512
A
]';

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

K0 = 5+1;
K1 = 7+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;


