function T = c3(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/33265657776681412629632*(1143669615207452182802718460155762471113055-23908827267681199855899583689948421795207040*A+9096078298311323765915*(6908888259343036215029286741851562123083776*A^2-631618284815552078709449150283113030705024*A+14456576594627597331222559374660858557503)^(1/2))^(1/2)   ;
elseif k(2) == 1
   B =    -1/33265657776681412629632*(1143669615207452182802718460155762471113055-23908827267681199855899583689948421795207040*A+9096078298311323765915*(6908888259343036215029286741851562123083776*A^2-631618284815552078709449150283113030705024*A+14456576594627597331222559374660858557503)^(1/2))^(1/2)   ;
elseif k(2) == 2
   B =    1/33265657776681412629632*(1143669615207452182802718460155762471113055-23908827267681199855899583689948421795207040*A-9096078298311323765915*(6908888259343036215029286741851562123083776*A^2-631618284815552078709449150283113030705024*A+14456576594627597331222559374660858557503)^(1/2))^(1/2)   ;
elseif k(2) == 3
   B =    -1/33265657776681412629632*(1143669615207452182802718460155762471113055-23908827267681199855899583689948421795207040*A-9096078298311323765915*(6908888259343036215029286741851562123083776*A^2-631618284815552078709449150283113030705024*A+14456576594627597331222559374660858557503)^(1/2))^(1/2)   ;
end
if k(3) == 0
   C =    1/11590518620431467670157652512*(53592079407148875587283225607026587187632119182859470+2023045359359384952536418982901605537197637762732831801344*B^2*A-27501321797362773266246383967373602620152057124061151232*B^2-1098245236631034489366411780709456257747783845534570240*A)^(1/2)   ;
elseif k(3) == 1
   C =    -1/11590518620431467670157652512*(53592079407148875587283225607026587187632119182859470+2023045359359384952536418982901605537197637762732831801344*B^2*A-27501321797362773266246383967373602620152057124061151232*B^2-1098245236631034489366411780709456257747783845534570240*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
-18956205018422329967246311424/25581577378781689622455*B^3*A+17418411432268754298634240/5116315475756337924491*B^3+155596886984220215050752/730902210822333989213*B*A+3051842094295400610291/730902210822333989213*B
-873115129319434154639439691776/179071041651471827357185*B^3*A+182932174570961961466380288/5116315475756337924491*B^3+221878104049971866146688/5116315475756337924491*B*A+5872898477348370983866/730902210822333989213*B
-4685782289473267171492641112064/639539434469542240561375*B^3*A+45405169270061632285321428992/639539434469542240561375*B^3-4011686875652326832050176/18272555270558349730325*B*A-182767201524923517038307/18272555270558349730325*B
76781418653618858854317129990144/4476776041286795683929625*B^3*A-78971258611119614417071128576/639539434469542240561375*B^3-99801438585668073824274304/127907886893908448112275*B*A-798690589469583957366854/18272555270558349730325*B
-315312845161512697218031/18272555270558349730325*B+28915023838932170594967521067008/639539434469542240561375*B^3*A-257311716206668994898343133184/639539434469542240561375*B^3-24498474262857662760488448/18272555270558349730325*B*A
-342817025338670162996315765604352/4476776041286795683929625*B^3*A+402989850255849190560863961088/639539434469542240561375*B^3+411627633195527871875469952/127907886893908448112275*B*A+2080124974149184391695862/18272555270558349730325*B
1444578542889494229549056/72270613330644841*B^3*A-39400353432189528670208/258109333323731575*B^3-34956986594464566272/51621866664746315*B*A-394519325598785518/7374552380678045*B
370056941152943401336832/51621866664746315*B^3*A-15953871393017683820544/258109333323731575*B^3-4042205484101420288/7374552380678045*B*A-33802953410867684/7374552380678045*B
]';
c = [
C
640/7*C*A+2*C
-68085370722427561999148253184/63042288169999389628125*C*B^2*A+475309046045667553242570752/63042288169999389628125*C*B^2+170390902788110731408768/600402744476184663125*C*A-5754673322305010180573/7204832933714215957500*C
-3031985375929723303883019124736/441296017189995727396875*C*B^2*A+515743269419451678153838592/9006041167142769946875*C*B^2+1071249590397984824197216/12608457633999877925625*C*A-16870709428411599024503/3602416466857107978750*C
-709836958538604080727879319552/189126864509998168884375*C*B^2*A+1177049742984526205988241408/27018123501428309840625*C*B^2-4536137420477575937557888/5403624700285661968125*C*A-9581334729016116879421/1543892771510189133750*C
5585349383183699545666729541632/88259203437999145479375*C*B^2*A-2310243731797688630453641216/4202819211333292641875*C*B^2-5549644845092449855223872/2521691526799975585125*C*A+1334372647539495628958/360241646685710797875*C
-36373705643062909006069874819072/441296017189995727396875*C*B^2*A+44250606979065990717171212288/63042288169999389628125*C*B^2+58156727623558316701234432/12608457633999877925625*C*A+107610996602893551059963/7204832933714215957500*C
40835483804163611111935007260672/1323888051569987182190625*C*B^2*A-49286737003139155123975794688/189126864509998168884375*C*B^2-76879243705972924923384032/37825372901999633776875*C*A-15365134451490576681017/1543892771510189133750*C
0
0
]';
T = [a;b;c]/sqrt(2);
