function [h0,h1] = symbal2o(k)

% [h0,h1] = symbal2o(k);
%
% Balanced GHM-like multiwavelet scaling filters:
% 2-balanced orthogonal scaling filters based on
% odd length symmetric FIR filters.
%
% the only continuous scaling filters are generated by:
% [h0,h1] = symbal2o([0 0]) 

if k(1) == 0
   A =    8/17+3/34*2^(1/2)   ;
elseif k(1) == 1
   A =    8/17-3/34*2^(1/2)   ;
end
if k(2) == 0
   B =    -1/2*A+1/2+1/204*(10404*A^2-10812*A+2907)^(1/2)   ;
elseif k(2) == 1
   B =    -1/2*A+1/2-1/204*(10404*A^2-10812*A+2907)^(1/2)   ;
end
h0 = [
B+2/3*A-2/3
-2*B-7/6*A+7/6
B
A
]*sqrt(2);
h0 = [h0; h0(3:-1:1); 0; 0; 0; 0; 0];
h1 = [
5/24+17/24*B*A-11/24*B-1/3*A
-35/96-17/12*B*A+11/12*B+7/12*A
17/24*B*A-11/24*B
-1/2*A+1/4
-17/12*B*A-13/12*B-1/3*A+23/24
17/6*B*A+13/6*B+7/6*A-53/48
]*sqrt(2);
h1 = [h1; h1(5:-1:1); 0];