function [h0,h1] = coeff1(k)
if k == 0
   A = 1/3+1/4*2^(1/2);
elseif k == 1
   A = 1/3-1/4*2^(1/2);
end
h0 = [1/3 4/3 1/3 0 0 0 0 0]'/sqrt(2);
h1 = [-A+2/3 -1/6 A 1 A -1/6 -A+2/3 0]'/sqrt(2);

