function T = c2(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/22688*(7653655-35892416*A+1418*(640697344*A^2-107915584*A+7525345)^(1/2))^(1/2)   ;
elseif k(2) == 1
   B =    -1/22688*(7653655-35892416*A+1418*(640697344*A^2-107915584*A+7525345)^(1/2))^(1/2)   ;
elseif k(2) == 2
   B =    1/22688*(7653655-35892416*A-1418*(640697344*A^2-107915584*A+7525345)^(1/2))^(1/2)   ;
elseif k(2) == 3
   B =    -1/22688*(7653655-35892416*A-1418*(640697344*A^2-107915584*A+7525345)^(1/2))^(1/2)   ;
end
if k(3) == 0
   C =    1/9667960*(-12992996923725+16062050591440896*B^2*A+2635012142349824*B^2+27224018740800*A)^(1/2)   ;
elseif k(3) == 1
   C =    -1/9667960*(-12992996923725+16062050591440896*B^2*A+2635012142349824*B^2+27224018740800*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 = [
-32246726656/33392625*B^3*A-2199102464/33392625*B^3+14847872/1335705*B*A+1737518/1335705*B
-15008202752/6678525*B^3*A-1356197888/6678525*B^3+1684864/89047*B*A+401472/89047*B
639219335168/634459875*B^3*A+2005256192/634459875*B^3-135543616/8459465*B*A+12755621/8459465*B
-83975200/5075679*B+51901431808/6042475*B^3*A+109927534592/126891975*B^3-40882624/725097*B*A
6075596/1208495*B-192876314624/30212375*B^3*A-18143865856/30212375*B^3+35732544/1208495*B*A
64/5*B*A+16/5*B
B
]';
c = [
-120483807232/43906625*C*B^2*A-23242676224/131719875*C*B^2+152431616/5268795*C*A+14974349/5268795*C
-171827658752/26343975*C*B^2*A-14075998208/26343975*C*B^2+98703168/1756265*C*A+3358028/351253*C
-46798364/1756265*C+2802630852608/131719875*C*B^2*A+210488010752/131719875*C*B^2-349837952/1756265*C*A
3672712/277305*C-11895046144/990375*C*B^2*A-2050269184/2310875*C*B^2+4518592/39615*C*A
C
0
0
]';
T = [a;b;c]/sqrt(2);