Radau IIA transformed coefficients

radau.jl

A = [11//45-7sqrt(6)/360         37//225-169sqrt(6)/1800   -2//225+sqrt(6)/75
     37//225+169sqrt(6)/1800     11//45+7sqrt(6)/360       -2//225-sqrt(6)/75
     4//9-sqrt(6)/36             4//9+sqrt(6)/36               1//9]

T = Float64
T11  = convert(T, 9.1232394870892942792e-02)
T12  = convert(T, -0.14125529502095420843e0)
T13  = convert(T, -3.0029194105147424492e-02)
T21  = convert(T, 0.24171793270710701896e0)
T22  = convert(T, 0.20412935229379993199e0)
T23  = convert(T, 0.38294211275726193779e0)
T31  = convert(T, 0.96604818261509293619e0)
T32 = 1
T33 = 0
TI11 = convert(T, 4.3255798900631553510e0)
TI12 = convert(T, 0.33919925181580986954e0)
TI13 = convert(T, 0.54177053993587487119e0)
TI21 = convert(T, -4.1787185915519047273e0)
TI22 = convert(T, -0.32768282076106238708e0)
TI23 = convert(T, 0.47662355450055045196e0)
TI31 = convert(T, -0.50287263494578687595e0)
TI32 = convert(T, 2.5719269498556054292e0)
TI33 = convert(T, -0.59603920482822492497e0)
T = [T11 T12 T13
T21 T22 T23
T31 T32 T33]
foreach(x->normalize!(x), eachcol(T))
V = eigvecs(inv(A))
(V[:, 1] - conj(V[:, 1])) |> normalize
T

TI = [TI11 TI12 TI13
TI21 TI22 TI23
TI31 TI32 TI33]