function [h0,h1,h2,h3] = symbal2e
A =    33/80+1/160*151^(1/2)   ;
B =    -1/320*(145520*A-46255)^(1/2)   ;
C =    -1/160*(37120*A-12055)^(1/2)   ;
h0 = [
-A+1/2
-2*A+29/32
2*A-29/32
A
]*sqrt(2);
h0 = [h0; h0(4:-1:1); 0; 0; 0; 0];
h1 = [
9/8*A-141/256
3/8*A-45/256
-13/8*A+203/256
1/8*A-29/256
-5/4*A+85/128
5/4*A-15/128
]*sqrt(2);
h1 = [h1; h1(6:-1:1)];
h2 = [
11520/1073*B*A-5647/1073*B
-160/29*B*A+74/29*B
-1760/111*B*A+854/111*B
B
]*sqrt(2);
h2 = [h2; -h2(4:-1:1); 0; 0; 0; 0];
h3 = [
-188/21*C*A+3679/840*C
-52/21*C*A+1031/840*C
268/21*C*A-5249/840*C
-92/21*C*A+1711/840*C
160/9*C*A-80/9*C
C
]*sqrt(2);
h3 = [h3; -h3(6:-1:1)];
