function T = c4(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/943942704726789760*(6578201978966350292765101973134695-104333096265894924273013815357119360*A+51621866664746315*(4084856293924899386173801214704500736*A^2-348535195340540527230057279096543104*A+8565695844712398123959415174607663)^(1/2))^(1/2)   ;
elseif k(2) == 1
   B =    -1/943942704726789760*(6578201978966350292765101973134695-104333096265894924273013815357119360*A+51621866664746315*(4084856293924899386173801214704500736*A^2-348535195340540527230057279096543104*A+8565695844712398123959415174607663)^(1/2))^(1/2)   ;
elseif k(2) == 2
   B =    1/943942704726789760*(6578201978966350292765101973134695-104333096265894924273013815357119360*A-51621866664746315*(4084856293924899386173801214704500736*A^2-348535195340540527230057279096543104*A+8565695844712398123959415174607663)^(1/2))^(1/2)   ;
elseif k(2) == 3
   B =    -1/943942704726789760*(6578201978966350292765101973134695-104333096265894924273013815357119360*A-51621866664746315*(4084856293924899386173801214704500736*A^2-348535195340540527230057279096543104*A+8565695844712398123959415174607663)^(1/2))^(1/2)   ;
end
if k(3) == 0
   C =    1/89351410301958635671897937392513400*(-719485888242130093256353940677539969891570027678997754363099416126770+2305151627253673582784670512673370366832559533138921360115516547085107200*B^2*A+208156044140633914985290914538025738961264503016805280562472187466444800*B^2+31269097019352092372197477412831687980821178881106509197010982012054720*A)^(1/2)   ;
elseif k(3) == 1
   C =    -1/89351410301958635671897937392513400*(-719485888242130093256353940677539969891570027678997754363099416126770+2305151627253673582784670512673370366832559533138921360115516547085107200*B^2*A+208156044140633914985290914538025738961264503016805280562472187466444800*B^2+31269097019352092372197477412831687980821178881106509197010982012054720*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
-1435980106392218940359240771960832/6310675382586770540051755517*B^3*A+849537003379343267588912480256/901525054655252934293107931*B^3+2372679902038289387483738424832/4507625273276264671465539655*B*A-400730498636670599699941441/643946467610894953066505665*B
28550629156767136137156916936704/128789293522178990613301133*B^3*A+106980218753333908087769055232/901525054655252934293107931*B^3-299034478587756325852282937728/643946467610894953066505665*B*A-11027932291305076949357993182/643946467610894953066505665*B
2115323904252992862744302640431104/6310675382586770540051755517*B^3*A-2966097658595752954212900831232/901525054655252934293107931*B^3-5921574628898871432093780987904/4507625273276264671465539655*B*A+10597628455656338302296447817/643946467610894953066505665*B
-186979829850528138698965777383424/901525054655252934293107931*B^3*A+842123601561936220282819231744/901525054655252934293107931*B^3+425801134328382960092870231424/643946467610894953066505665*B*A+3621661089440393601638626386/643946467610894953066505665*B
-1375934958728041892100054135603200/6310675382586770540051755517*B^3*A+1655009224812962578244002611200/901525054655252934293107931*B^3+746218871718119064629233937920/901525054655252934293107931*B*A-736608162636584156707111975/128789293522178990613301133*B
-201644087266919864654154956800/901525054655252934293107931*B^3*A+198226949880807245016387174400/901525054655252934293107931*B^3+4903919449941244170804888960/128789293522178990613301133*B*A-413175766582616226047882010/128789293522178990613301133*B
902861589305933893468160000/12734509617635853272281*B^3*A-985008835804738216755200/1819215659662264753183*B^3-483267588478451911613440/1819215659662264753183*B*A+234136140684700553290/259887951380323536169*B
180126391418905167020/259887951380323536169*B+46257117644117925167104000/1819215659662264753183*B^3*A-398846784825442095513600/1819215659662264753183*B^3-26284764140739471765760/259887951380323536169*B*A
]';
c = [
C
-202861285139595933482242539520/9723838184289469629031713*C*B^2*A+255031177180610185059202334720/204200601870078862209665973*C*B^2+8181923205157202962226913664/29171514552868408887095139*C*A-51569511277721558405191234/3241279394763156543010571*C
25823854389018593021088726983928971264/505965188304653353600177202909805*C*B^2*A-1479053097086588930614370887924097024/505965188304653353600177202909805*C*B^2-238205820088895895101324814070052224/361403705931895252571555144935575*C*A+47487730935659660656084172986228111/1445614823727581010286220579742300*C
-15166914672728209599154815495180910592/505965188304653353600177202909805*C*B^2*A+15434345759405063028980824647446528/10325820169482721502044432712445*C*B^2+123645145823762581229040630703389472/361403705931895252571555144935575*C*A-1605163254360485542504708353840197/103258201694827215020444327124450*C
-555697116723897274822787714010578944/72280741186379050514311028987115*C*B^2*A+297418189387240418583388922508771328/505965188304653353600177202909805*C*B^2+46345282672389910727022573697113728/361403705931895252571555144935575*C*A-4787464872372220593925592096083921/722807411863790505143110289871150*C
2766278158263488818959139323155316736/505965188304653353600177202909805*C*B^2*A-16872518688256227275793598378221568/72280741186379050514311028987115*C*B^2-19327016237569279174169360082723776/361403705931895252571555144935575*C*A+118725642823765809221962462545488/51629100847413607510222163562225*C
203789540969213177998296017569054720/101193037660930670720035440581961*C*B^2*A-15476028407866203763787504623083520/101193037660930670720035440581961*C*B^2-484987842372609099804394063112704/14456148237275810102862205797423*C*A+96310058290764591816408647011831/57824592949103240411448823189692*C
658775011846140756849608960245760/101193037660930670720035440581961*C*B^2*A-105454161950792133585727499816960/4818716079091936700954068599141*C*B^2-67925162192971635770649801387232/14456148237275810102862205797423*C*A+1091054743708330713322170188557/4130328067793088600817773084978*C
0
0
]';
T = [a;b;c]/sqrt(2);
