function T = c1(k)
if k(1) == 0
   A =    7/256+1/256*21^(1/2)   ;
elseif k(1) == 1
   A =    7/256-1/256*21^(1/2)   ;
end
if k(2) == 0
   B =    1/188788540945357952*(48466980846931897102088579270655-983245269816062943610766036051840*A+51621866664746315*(362790847485380478974462266064896*A^2-26998021636070769168881055403904*A+480419612082462247915209800863)^(1/2))^(1/2)   ;
elseif k(2) == 1
   B =    -1/188788540945357952*(48466980846931897102088579270655-983245269816062943610766036051840*A+51621866664746315*(362790847485380478974462266064896*A^2-26998021636070769168881055403904*A+480419612082462247915209800863)^(1/2))^(1/2)   ;
elseif k(2) == 2
   B =    1/188788540945357952*(48466980846931897102088579270655-983245269816062943610766036051840*A-51621866664746315*(362790847485380478974462266064896*A^2-26998021636070769168881055403904*A+480419612082462247915209800863)^(1/2))^(1/2)   ;
elseif k(2) == 3
   B =    -1/188788540945357952*(48466980846931897102088579270655-983245269816062943610766036051840*A-51621866664746315*(362790847485380478974462266064896*A^2-26998021636070769168881055403904*A+480419612082462247915209800863)^(1/2))^(1/2)   ;
end
if k(3) == 0
   C =    1/56538280293200*(2168461700908060479940235+175382233365724650317149110272*B^2*A+89936081632916799286902784*B^2-54394341069642946454741120*A)^(1/2)   ;
elseif k(3) == 1
   C =    -1/56538280293200*(2168461700908060479940235+175382233365724650317149110272*B^2*A+89936081632916799286902784*B^2-54394341069642946454741120*A)^(1/2)   ;
end
a = [
A
7*A-5/128
20*A-7/32
28*A-7/16
14*A-7/32
-14*A+35/64
-28*A+35/32
-20*A+7/8
-7*A+11/32
-A+7/128
]';
b = [
B
640/7*B*A+2*B
1218957733/7632057275*B-134444395102470144/267122004625*B^3*A+863072521388032/267122004625*B^3+2168159768064/7632057275*B*A
-21199367274/7632057275*B-6035408946332696576/1869854032375*B^3*A+6977879226466304/267122004625*B^3+9329628464256/53424400925*B*A
-221783756423/22896171825*B-2085565836795314176/801366013875*B^3*A+22620684929302528/801366013875*B^3-12585779605504/22896171825*B*A
-50989306138/7632057275*B+46296366289485037568/1869854032375*B^3*A-55893841716559872/267122004625*B^3-121271887913088/53424400925*B*A
-7522839167816957952/373970806475*B^3*A+8554713706430464/53424400925*B^3+25975933409792/10684880185*B*A+49734719267/1526411455*B
-13267516402/915846873*B-512395464518139904/160273202775*B^3*A+5308410155450368/160273202775*B^3+3105440390528/4579234365*B*A
-10123814/3296785*B+116541632806912/32308493*B^3*A-3477306179584/115387475*B^3-13594541056/23077495*B*A
3165068/3296785*B+29543593148416/23077495*B^3*A-69530370048/6073025*B^3-43492096/173515*B*A
]';
c = [
C
640/7*C*A+2*C
-48314155147039/21153268376100*C-310038519908112269312/185091098290875*C*B^2*A+2509475797778702336/185091098290875*C*B^2+1972605653279872/5288317094025*C*A
-13090099811681/1175181576450*C-1509535884991014633472/143959743115125*C*B^2*A+5709943759822680064/61697032763625*C*B^2+7078984461222496/12339406552725*C*A
50322299443017839280128/1295637688036125*C*B^2*A-62051511994469875712/185091098290875*C*B^2-164664723680272768/37018219658175*C*A+114865485048569/10576634188050*C
32215683414856/5288317094025*C-47749264145970868781056/1295637688036125*C*B^2*A+8325827189830623232/26441585470125*C*B^2+177722210906582336/37018219658175*C*A
-10984834305329/1410217891740*C+648418001289911730176/86375845869075*C*B^2*A-109273797914140672/1762772364675*C*B^2-264788146761984/274209034505*C*A
2665495745237/2115326837610*C+14109335868899065856/5288317094025*C*B^2*A-878967326121291776/37018219658175*C*B^2-449737237435552/1057663418805*C*A
0
0
]';
T = [a;b;c]/sqrt(2);