function T = c2(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/421988480*(484926687032655-9535041965805440*A+23077495*(170713507906207744*A^2-11419307696671616*A+169606415611927)^(1/2))^(1/2)   ;
elseif k(2) == 1
   B =    -1/421988480*(484926687032655-9535041965805440*A+23077495*(170713507906207744*A^2-11419307696671616*A+169606415611927)^(1/2))^(1/2)   ;
elseif k(2) == 2
   B =    1/421988480*(484926687032655-9535041965805440*A-23077495*(170713507906207744*A^2-11419307696671616*A+169606415611927)^(1/2))^(1/2)   ;
elseif k(2) == 3
   B =    -1/421988480*(484926687032655-9535041965805440*A-23077495*(170713507906207744*A^2-11419307696671616*A+169606415611927)^(1/2))^(1/2)   ;
end
if k(3) == 0
   C =    1/151764040252182895546080*(-58989258302518894587632769791362130918348910+3848151610865189741151828772397599970503006617600*B^2*A+28744089220769467970154767193597855426991718400*B^2+864670941504079952261124570016181791422293760*A)^(1/2)   ;
elseif k(3) == 1
   C =    -1/151764040252182895546080*(-58989258302518894587632769791362130918348910+3848151610865189741151828772397599970503006617600*B^2*A+28744089220769467970154767193597855426991718400*B^2+864670941504079952261124570016181791422293760*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
-1444927673113000213/2048650651352360901*B-5203811293263046975160320/14340554559466526307*B^3*A+34406404399055785984000/14340554559466526307*B^3+620922719571416355328/2048650651352360901*B*A
-3585147956221854866/682883550450786967*B-77655551518164532917698560/33461293972088561383*B^3*A+90695562480399781478400/4780184853155508769*B^3+1092577070298953913984/4780184853155508769*B*A
-42314393997434519957/3414417752253934835*B+196742688958055917358678016/33461293972088561383*B^3*A-212381366031363128393728/4780184853155508769*B^3-43540914883261331542016/23900924265777543845*B*A
233644751157047517698/10243253256761804505*B-20907444062232700572925952/14340554559466526307*B^3*A+105207945501869765476352/14340554559466526307*B^3+4106270648273229642112/10243253256761804505*B*A
-5765746660505667081/3414417752253934835*B-83306858692484923897413632/33461293972088561383*B^3*A+100175139451639589666816/4780184853155508769*B^3+26408068305383404442112/23900924265777543845*B*A
-21516190290000716998/3414417752253934835*B-711724355641328332177408/4780184853155508769*B^3*A+11594082955585759920128/4780184853155508769*B^3+491634593261260608128/3414417752253934835*B*A
48026588364166922240000/72270613330644841*B^3*A-57319654266299187200/10324373332949263*B^3-3324201852365347840/10324373332949263*B*A-247628478364310/1474910476135609*B
2434971868445016064000/10324373332949263*B^3*A-21776534696776089600/10324373332949263*B^3-190470692623663360/1474910476135609*B*A+1027518993178220/1474910476135609*B
]';
c = [
C
640/7*C*A+2*C
-232783294224520553/11320644997903500*C-219715425923642547503104/27735580244863575*C*B^2*A+270164976282465329152/3962225749266225*C*B^2+65010164133218036864/19811128746331125*C*A
56079110652856139/1886774166317250*C+182187361998735321595904/9245193414954525*C*B^2*A-222448226700805271552/1320741916422075*C*B^2-55993280545558125664/6603709582110375*C*A
-18664617197293/5989759258150*C-7571094443104532430848/616346227663635*C*B^2*A+27075492831456919552/264148383284415*C*B^2+6833777580050637184/1320741916422075*C*A
-5804984165787238/404308749925125*C-8588643845956919361536/3962225749266225*C*B^2*A+1714717150495514624/80861749985025*C*B^2+3298267055802688576/2830161249475875*C*A
11277684428875189/3773548332634500*C+18222516431965875863552/9245193414954525*C*B^2*A-22628539722261733376/1320741916422075*C*B^2-5898334400298783232/6603709582110375*C*A
4388220861995681/1886774166317250*C+919311785378159525888/1320741916422075*C*B^2*A-8362736619799441408/1320741916422075*C*B^2-107002469630003936/314462361052875*C*A
0
0
]';
T = [a;b;c]/sqrt(2);