function T = c1(k)
if k(1) == 0
   A =    1/32+1/32*6^(1/2)   ;
elseif k(1) == 1
   A =    1/32-1/32*6^(1/2)   ;
end
if k(2) == 0
   B =    1/8548512*(356033503455-1983237686976*A+534282*(13778706433024*A^2-2844013344064*A+169187913745)^(1/2))^(1/2)   ;
elseif k(2) == 1
   B =    -1/8548512*(356033503455-1983237686976*A+534282*(13778706433024*A^2-2844013344064*A+169187913745)^(1/2))^(1/2)   ;
elseif k(2) == 2
   B =    1/8548512*(356033503455-1983237686976*A-534282*(13778706433024*A^2-2844013344064*A+169187913745)^(1/2))^(1/2)   ;
elseif k(2) == 3
   B =    -1/8548512*(356033503455-1983237686976*A-534282*(13778706433024*A^2-2844013344064*A+169187913745)^(1/2))^(1/2)   ;
end
if k(3) == 0
   C =    1/107768*(20199410+161001308160*B^2*A+24063437312*B^2-352889152*A)^(1/2)   ;
elseif k(3) == 1
   C =    -1/107768*(20199410+161001308160*B^2*A+24063437312*B^2-352889152*A)^(1/2)   ;
end
a = [
-A+1/16
-4*A+5/16
-5*A+5/8
5/8
5*A+5/16
4*A+1/16
A
]';
b = [
-955121664/88625*B^3*A-65135616/88625*B^3+101248/3545*B*A+15262/3545*B
10944/709*B-444530688/17725*B^3*A-40169472/17725*B^3+25728/709*B*A
-2002337/67355*B+128119799808/1683875*B^3*A+10264676352/1683875*B^3-7406528/67355*B*A
75968/67355*B-56959107072/1683875*B^3*A-4130746368/1683875*B^3+647872/67355*B*A
62764/13471*B-2156593152/336775*B^3*A-216050688/336775*B^3+305728/13471*B*A
64/5*B*A+16/5*B
B
]';
c = [
-377147817984/9250625*C*B^2*A-26754717696/9250625*C*B^2+29378944/370025*C*A+5382691/370025*C
754295635968/9250625*C*B^2*A+53509435392/9250625*C*B^2-54021568/370025*C*A-8841252/370025*C
-377147817984/9250625*C*B^2*A-26754717696/9250625*C*B^2+19906304/370025*C*A+1904456/370025*C
64/5*C*A+16/5*C
C
0
0
]';
T = [a;b;c]/sqrt(2);
