with(LinearAlgebra): #spisok[[3],[3],[2,1]]:=[[7,2],[7,1,1],[6,3],[6,2,1],[6,1],[5,4],[5,3,1],[5,2,2],[5,2,1,1],[4,4],[4,3,1],[4,3,1,1]] R[7,2][[3],[3]]:=DiagonalMatrix(): Racah[7,2][[3],[3],[2,1]]:=Matrix([[(q^3*sqrt(1 - q^2 + 2*q^4 - q^6 + q^8))/(1 + q^4 + q^6 + q^10), sqrt(1 + 2*q^4 + q^6 + 2*q^8 + 2*q^10 + 2*q^12 + q^14 + 2*q^16 + q^20)/(1 + q^4 + q^6 + q^10)], [sqrt((1 + q^4 + q^6 + q^8 + q^10 + q^12 + q^16)/(1 + q^4))/(1 + q^6), -(q^3/((1 + q^2)*sqrt((1 + q^4)*(1 - q^2 + q^4))))]]): R[7,1,1][[3],[3]]:=DiagonalMatrix(): Racah[7,1,1][[3],[3],[2,1]]:=Matrix([[q^3/(1 + q^6), sqrt(1 + q^6 + q^12)/(1 + q^6)], [sqrt(1 + q^6 + q^12)/(1 + q^6), -(q^3/(1 + q^6))]]): R[6,3][[3],[3]]:=DiagonalMatrix(<-q^3,q^(-1)>): Racah[6,3][[3],[3],[2,1]]:=Matrix([[q^2/(1 + q^4), sqrt(1 + q^4 + q^8)/(1 + q^4)], [sqrt(1 + q^4 + q^8)/(1 + q^4), -(q^2/(1 + q^4))]]): R[6,2,1][[3],[3]]:=DiagonalMatrix(): Racah[6,2,1][[3],[3],[2,1]]:=Matrix([[q^6/((1 + q^6)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)), (q^3*(1 + q^2 + q^4)^2)/((1 + q^2)*(1 + q^2 + q^4 + q^6 + q^8)*sqrt(1 + q^2 + 2*q^4 + q^6 + q^8)), (q^2*(1 + q^4 + q^8 + q^12))/((1 + q^6)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 3*q^8 + 2*q^10 + 3*q^12 + 2*q^14 + 2*q^16 + q^18 + q^20)), ((1 + q^2)*sqrt(1 + q^4 + q^8 + q^12))/(1 + q^2 + q^4 + q^6 + q^8)], [0, (1 + q^4)*sqrt(((1 + q^2 + q^4)*(1 + q^8))/((1 + q^2 + q^4 + q^6 + q^8)*(1 + q^2 + 2*q^6 + q^10 + q^12))), -((q*(1 - q^2 + q^4))/sqrt(1 + q^2 + 2*q^6 + q^10 + q^12)), -(q^5/sqrt((1 + q^2 + q^4 + q^6 + q^8)*(1 + q^2 + 2*q^6 + q^10 + q^12)))], [(q^2*sqrt((1 + q^2 - q^4 + q^6 + q^8)/(1 + q^4)))/(1 + q^6), (q*(1 - q^4 + q^8)*sqrt(((1 + q^2 + q^4)*(1 + 2*q^2 + q^4 + 2*q^6 + 3*q^8 + 2*q^10 + q^12 + 2*q^14 + q^16))/(1 + q^4)))/((1 + q^2)*(1 + q^2 - q^4 + q^6 + q^8)*(1 + q^2 + q^4 + q^6 + q^8)), ((1 + q^8)*sqrt((1 + q^4 + q^8 + q^12)/(1 + q^2 - q^4 + q^6 + q^8)))/(1 + q^4 + q^6 + q^10), -((q^2*(1 + q^2)*(1 + q^8))/sqrt((1 + q^2 - q^4 + q^6 + q^8)*(1 + q^2 + 2*q^4 + 2*q^6 + 3*q^8 + 2*q^10 + 3*q^12 + 2*q^14 + 2*q^16 + q^18 + q^20)))], [(1 + q^8)*sqrt((1 + q^2 + q^4)/(1 + q^2 + 2*q^4 + 2*q^6 + 3*q^8 + 2*q^10 + 3*q^12 + 2*q^14 + 2*q^16 + q^18 + q^20)), -((q^3*(1 + q^8))/((1 + q^2 + q^4 + q^6 + q^8)*sqrt(1 + q^4 + q^8 + q^12))), -(q^2*sqrt((1 + q^2 + q^4)/(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12))), (q^4*sqrt((1 + q^2 + q^4)/(1 + q^4)))/(1 + q^2 + q^4 + q^6 + q^8)]]): R[5,4][[3],[3]]:=DiagonalMatrix(): Racah[5,4][[3],[3],[2,1]]:=Matrix([[(q + q^3 + q^5)/((1 + q^2)*sqrt(1 + q^2 + 2*q^4 + q^6 + q^8)), sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)/(1 + q^2 + q^4 + q^6)], [(1 + q^2 + q^4 + q^6 + q^8)/((1 + q^2)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)), -((q*sqrt((1 + q^2 + q^4)/(1 + q^4)))/(1 + q^2))]]): R[5,3,1][[3],[3]]:=DiagonalMatrix(<-q^3, q^(-1), q^(-1), -q^(-3)>): Racah[5,3,1][[3],[3],[2,1]]:=Matrix([[-(q^4/((1 + q^2 + q^4)*sqrt(1 - q^2 + 2*q^4 - q^6 + q^8))), -(q^2/sqrt(1 + q^2 + 2*q^4 + q^6 + q^8)), -((q*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/((1 + q^2 + q^4)*sqrt(1 + 2*q^4 + q^6 + 2*q^8 + 2*q^10 + 2*q^12 + q^14 + 2*q^16 + q^20))), (-1 - q^2 - q^4 - q^6 - q^8 - q^10 - q^12)/((1 + q^2 + q^4)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + 2*q^10 + 2*q^12 + q^14 + q^16))], [0, ((1 + q^4)*(1 + q^2 + q^4)*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/sqrt((1 + q^2 + q^4 + q^6 + q^8)*(1 + q^2 + 2*q^4 + q^6 + q^8)*(1 + 3*q^2 + 5*q^4 + 6*q^6 + 8*q^8 + 10*q^10 + 11*q^12 + 10*q^14 + 8*q^16 + 6*q^18 + 5*q^20 + 3*q^22 + q^24)), -((q*(1 - q^2 + q^4)*(1 + q^2 + q^4 + q^6 + q^8))/sqrt(1 + 2*q^2 + 4*q^4 + 6*q^6 + 10*q^8 + 11*q^10 + 14*q^12 + 14*q^14 + 14*q^16 + 11*q^18 + 10*q^20 + 6*q^22 + 4*q^24 + 2*q^26 + q^28)), -(q^6/sqrt(1 + 3*q^2 + 6*q^4 + 9*q^6 + 13*q^8 + 15*q^10 + 16*q^12 + 15*q^14 + 13*q^16 + 9*q^18 + 6*q^20 + 3*q^22 + q^24))], [q/((1 + q^6)*sqrt(((1 + q^4)*(1 + q^2 + q^4))/(1 + q^2 + q^4 + q^6 + 3*q^8 + q^10 + q^12 + q^14 + q^16))), (q - q^5 + q^9)/((1 + q^2)*sqrt((1 + q^4)*(1 + 2*q^2 + 2*q^4 + q^6 + 2*q^8 + 2*q^10 + q^12))), (1 + q^2 + q^4 + q^6 + 2*q^8 + 2*q^10 + 2*q^12 + q^14 + q^16 + q^18 + q^20)/((1 + q^2)*sqrt((1 + q^2 + 2*q^4 + q^6 + q^8)*(1 + 2*q^2 + 2*q^4 + q^6 + 2*q^8 + 2*q^10 + q^12)*(1 + q^4 + q^6 + q^8 + q^10 + q^12 + q^16))), -((q*(1 + q^2)*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/sqrt((1 + q^4)*(1 + q^2 + q^4)*(1 + 3*q^2 + 5*q^4 + 6*q^6 + 8*q^8 + 10*q^10 + 11*q^12 + 10*q^14 + 8*q^16 + 6*q^18 + 5*q^20 + 3*q^22 + q^24)))], [(1 + 2*q^2 + 3*q^4 + 4*q^6 + 5*q^8 + 5*q^10 + 5*q^12 + 4*q^14 + 3*q^16 + 2*q^18 + q^20)/((1 + 2*q^2 + 2*q^4 + q^6)*sqrt(1 + q^2 + 3*q^4 + 4*q^6 + 6*q^8 + 7*q^10 + 9*q^12 + 8*q^14 + 9*q^16 + 7*q^18 + 6*q^20 + 4*q^22 + 3*q^24 + q^26 + q^28)), -((q^2*sqrt((1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/(1 + 2*q^2 + 4*q^4 + 5*q^6 + 6*q^8 + 5*q^10 + 4*q^12 + 2*q^14 + q^16)))/(1 + q^2)), -((q*(1 + q^4)*(1 + q^2 + q^4 + q^6 + q^8))/((1 + q^2 + q^4)*sqrt((1 + q^2)^2*(1 + 2*q^4 + q^6 + 2*q^8 + q^10 + 2*q^12 + q^16)))), (q^2*(1 + q^2)*(1 + q^4))/((1 + q^2 + q^4)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12))]]): R[5,2,2][[3],[3]]:=DiagonalMatrix(<-q^3,q^(-1)>): Racah[5,2,2][[3],[3],[2,1]]:=Matrix([[q^2*sqrt((1 + q^2 + q^4)/(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)), (1 + q^4 + q^6 + q^8 + q^10 + q^12 + q^16)/sqrt((1 + q^4 + q^6 + q^8 + q^12)*(1 + 2*q^4 + q^6 + 2*q^8 + 2*q^10 + 2*q^12 + q^14 + 2*q^16 + q^20))], [(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/sqrt(1 + 2*q^2 + 4*q^4 + 6*q^6 + 8*q^8 + 9*q^10 + 10*q^12 + 9*q^14 + 8*q^16 + 6*q^18 + 4*q^20 + 2*q^22 + q^24), -((q^2*(1 + q^2 + q^4))/sqrt((1 + q^2 + q^4 + q^6 + q^8)*(1 + q^2 + 2*q^4 + q^6 + q^8)))]]): R[5,2,1,1][[3],[3]]:=DiagonalMatrix(<-q^3,q^(-1)>): Racah[5,2,1,1][[3],[3],[2,1]]:=Matrix([[q^2/sqrt(1 + q^2 + 2*q^4 + q^6 + q^8), sqrt((1 + q^2 + q^4 + q^6 + q^8)/(1 + q^2 + 2*q^4 + q^6 + q^8))], [(1 + q^2 + q^4 + q^6 + q^8)/sqrt((1 + q^2 + q^4)*(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)), -(q^2/sqrt((1 + q^4)*(1 + q^2 + q^4)))]]): R[4,4,1][[3],[3]]:=DiagonalMatrix(): Racah[4,4,1][[3],[3],[2,1]]:=Matrix([[sqrt(1 + q^2 + q^4)/(1 + q^2), -(q/(1 + q^2))], [-(q/(1 + q^2)), -(sqrt(1 + q^2 + q^4)/(1 + q^2))]]): R[4,3,2][[3],[3]]:=DiagonalMatrix(): Racah[4,3,2][[3],[3],[2,1]]:=Matrix([[(1 - q^2 + q^4)*sqrt((1 + q^2 + q^4)/(1 + q^4 + q^6 + q^8 + q^12)), -((q*(1 + q^4))/sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12))], [-((q*(1 + q^4))/sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)), -((1 + q^2 + q^4)*sqrt((1 - q^2 + q^4)/(1 + 2*q^2 + 3*q^4 + 3*q^6 + 3*q^8 + 2*q^10 + q^12)))]]): R[4,3,1,1][[3],[3]]:=DiagonalMatrix(): Racah[4,3,1,1][[3],[3],[2,1]]:=Matrix([[q/sqrt(1 + q^2 + q^4), (1 + q^4)/sqrt(1 + q^2 + 2*q^4 + q^6 + q^8)], [(1 + q^4)/sqrt(1 + q^2 + 2*q^4 + q^6 + q^8), -(q/sqrt(1 + q^2 + q^4))]]): R[7,2][[3],[2,1]]:=DiagonalMatrix(): Racah[7,2][[3],[2,1],[3]]:=Matrix([[(q + q^3 + q^5)/(1 + q^2 + q^4 + q^6), sqrt(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/(1 + q^2 + q^4 + q^6)], [sqrt(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/(1 + q^2 + q^4 + q^6), -((q + q^3 + q^5)/(1 + q^2 + q^4 + q^6))]]): R[7,1,1][[3],[2,1]]:=DiagonalMatrix(): Racah[7,1,1][[3],[2,1],[3]]:=Matrix([[-(q^3/(1 + q^6)), sqrt(1 + q^6 + q^12)/(1 + q^6)], [sqrt(1 + q^6 + q^12)/(1 + q^6), q^3/(1 + q^6)]]): R[6,3][[3],[2,1]]:= DiagonalMatrix(): Racah[6,3][[3],[2,1],[3]]:=Matrix([[q^2/(1 + q^4), sqrt(1 + q^4 + q^8)/(1 + q^4)], [sqrt(1 + q^4 + q^8)/(1 + q^4), -(q^2/(1 + q^4))]]): R[6,2,1][[3],[2,1]]:=DiagonalMatrix(): Racah[6,2,1][[3],[2,1],[3]]:=Matrix([[-(q^5/(1 + q^4 + q^6 + q^10)), -((q^4*sqrt((1 + q^2 + q^4)/(1 + q^2 + q^4 + q^6 + q^8)))/(1 + q^2 + q^4 + q^6)), (q*sqrt(1 + q^8))/(1 + q^6), (1 + q^8)/sqrt(1 + q^2 + q^4 + q^6 + 2*q^8 + q^10 + q^12 + q^14 + q^16)], [-((q^4*sqrt((1 + q^2 + q^4)/(1 + q^2 + q^4 + q^6 + q^8)))/(1 + q^2 + q^4 + q^6)), (q + q^3 + 3*q^5 + 3*q^7 + 3*q^9 + q^11 + q^13)/((1 + q^2)*(1 + q^4)*(1 + q^2 + q^4 + q^6 + q^8)), ((1 + q^8)*sqrt((1 + q^2 + q^4)/(1 + q^2 + q^4 + q^6 + 2*q^8 + q^10 + q^12 + q^14 + q^16)))/(1 + q^2), -((q*sqrt(1 + q^2 + q^4 + q^8 + q^10 + q^12))/(1 + q^2 + q^4 + q^6 + q^8))], [(q*sqrt(1 + q^8))/(1 + q^6), sqrt((1 + q^8)*(1 + 2*q^2 + 3*q^4 + 3*q^6 + 3*q^8 + 2*q^10 + q^12))/((1 + q^2)*(1 + q^2 + q^4 + q^6 + q^8)), -((q*(-1 + q^2)^2)/(1 + q^6)), q^2/sqrt(1 + q^2 + q^4 + q^6 + q^8)], [(1 + q^8)/sqrt(1 + q^2 + q^4 + q^6 + 2*q^8 + q^10 + q^12 + q^14 + q^16), -((q*sqrt((1 + q^2 + q^4)*(1 + q^8)))/(1 + q^2 + q^4 + q^6 + q^8)), q^2/sqrt(1 + q^2 + q^4 + q^6 + q^8), -((q^3 + q^5)/(1 + q^2 + q^4 + q^6 + q^8))]]): R[5,4][[3],[2,1]]:=DiagonalMatrix(): Racah[5,4][[3],[2,1],[3]]:=Matrix([[q^3/(1 + q^2 + q^4 + q^6), ((1 + q^2 + q^4)*(1 + q^2 + q^4 + q^6 + q^8))/((1 + q^2)*(1 + q^4)*sqrt(1 + 2*q^2 + 3*q^4 + 3*q^6 + 3*q^8 + 2*q^10 + q^12))], [((1 + q^2 + q^4)*(1 + q^2 + q^4 + q^6 + q^8))/((1 + q^2)*(1 + q^4)*sqrt(1 + 2*q^2 + 3*q^4 + 3*q^6 + 3*q^8 + 2*q^10 + q^12)), -(q^3/(1 + q^2 + q^4 + q^6))]]): R[5,3,1][[3],[2,1]]:= DiagonalMatrix(): Racah[5,3,1][[3],[2,1],[3]]:=Matrix([[-(q^7/(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + 2*q^10 + q^12 + q^14)), -(q^5/((1 + q^2 + q^4 + q^6)*sqrt(1 + q^4 + q^8))), (q^2*sqrt(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/((1 + q^2 + q^4)*(1 + q^6)), (1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/((1 + q^2 + q^4)*sqrt(1 + q^4 + q^6 + q^8 + q^10 + q^12 + q^16))], [-(q^5/((1 + q^2 + q^4 + q^6)*sqrt(1 + q^4 + q^8))), (q + q^3 + 2*q^5 + 3*q^7 + 2*q^9 + q^11 + q^13)/((1 + q^2)*(1 + q^4)*(1 + q^2 + q^4 + q^6 + q^8)), (1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/((1 + q^2)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 3*q^8 + 3*q^10 + 3*q^12 + 2*q^14 + 2*q^16 + q^18 + q^20)), -((q^2*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/((1 + q^2 + q^4 + q^6 + q^8)*sqrt(1 + 2*q^2 + 3*q^4 + 3*q^6 + 3*q^8 + 3*q^10 + 3*q^12 + 2*q^14 + q^16)))], [(q^2*sqrt(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/(1 + q^2 + q^4 + q^6 + q^8 + q^10), (1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/((1 + q^2)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 3*q^8 + 3*q^10 + 3*q^12 + 2*q^14 + 2*q^16 + q^18 + q^20)), -((q*(1 + q^8))/((1 + q^2 + q^4)*(1 + q^6))), q^3/(sqrt(1 - q^2 + q^4)*(1 + q^2 + q^4))], [(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/((1 + q^2 + q^4)*sqrt(1 + q^4 + q^6 + q^8 + q^10 + q^12 + q^16)), -((q^2*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/((1 + q^2 + q^4 + q^6 + q^8)*sqrt(1 + 2*q^2 + 3*q^4 + 3*q^6 + 3*q^8 + 3*q^10 + 3*q^12 + 2*q^14 + q^16))), q^3/(sqrt(1 - q^2 + q^4)*(1 + q^2 + q^4)), -((q^5*(1 + q^2))/((1 + q^2 + q^4)*(1 + q^2 + q^4 + q^6 + q^8)))]]): R[5,2,2][[3],[2,1]]:=DiagonalMatrix(): Racah[5,2,2][[3],[2,1],[3]]:=Matrix([[-((q^3 + q^5)/(1 + q^2 + q^4 + q^6 + q^8)), ((1 + q^2 + q^4)*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/((1 + q^2 + q^4 + q^6 + q^8)*sqrt(1 + 2*q^2 + 3*q^4 + 3*q^6 + 3*q^8 + 3*q^10 + 3*q^12 + 2*q^14 + q^16))], [sqrt((1 + q^2 + q^4)*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/(1 + q^2 + q^4 + q^6 + q^8), (q^3 + q^5)/(1 + q^2 + q^4 + q^6 + q^8)]]): R[5,2,1,1][[3],[2,1]]:=DiagonalMatrix(<1,-q^(-3)>): Racah[5,2,1,1][[3],[2,1],[3]]:=Matrix([[(q + q^3)/(1 + q^2 + q^4), sqrt(1 + q^2 + q^4 + q^6 + q^8)/(1 + q^2 + q^4)], [sqrt(1 + q^2 + q^4 + q^6 + q^8)/(1 + q^2 + q^4), -((q + q^3)/(1 + q^2 + q^4))]]): R[4,4,1][[3],[2,1]]:=DiagonalMatrix(<-q^2,1>): Racah[4,4,1][[3],[2,1],[3]]:=Matrix([[q/(1 + q^2), sqrt(1 + q^2 + q^4)/(1 + q^2)], [sqrt(1 + q^2 + q^4)/(1 + q^2), -(q/(1 + q^2))]]): R[4,3,2][[3],[2,1]]:=DiagonalMatrix(): Racah[4,3,2][[3],[2,1],[3]]:=Matrix([[-(q^4/(1 + q^2 + q^4 + q^6 + q^8)), ((1 + q^2)*sqrt(1 + 2*q^4 + 2*q^8 + q^12))/(1 + q^2 + q^4 + q^6 + q^8)], [((1 + q^2)*sqrt(1 + 2*q^4 + 2*q^8 + q^12))/(1 + q^2 + q^4 + q^6 + q^8), q^4/(1 + q^2 + q^4 + q^6 + q^8)]]): R[4,3,1,1][[3],[2,1]]:=DiagonalMatrix(<1,-q^(-3)>): Racah[4,3,1,1][[3],[2,1],[3]]:=Matrix([[q^2/(1 + q^2 + q^4), ((1 + q^2)*sqrt(1 + q^4))/(1 + q^2 + q^4)], [((1 + q^2)*sqrt(1 + q^4))/(1 + q^2 + q^4), -(q^2/(1 + q^2 + q^4))]]): Racah[7,2][[2,1],[3],[3]]:=Matrix([[q^3/((1 + q^2)*sqrt(1 - q^2 + 2*q^4 - q^6 + q^8)), (1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/((1 + q^2)*sqrt(1 + 2*q^4 + q^6 + 2*q^8 + 2*q^10 + 2*q^12 + q^14 + 2*q^16 + q^20))], [(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/((1 + q^2)*sqrt(1 + 2*q^4 + q^6 + 2*q^8 + 2*q^10 + 2*q^12 + q^14 + 2*q^16 + q^20)), -(q^3/((1 + q^2)*sqrt(1 - q^2 + 2*q^4 - q^6 + q^8)))]]): Racah[7,1,1][[2,1],[3],[3]]:=Matrix([[q^3/(1 + q^6), sqrt(1 + q^6 + q^12)/(1 + q^6)], [sqrt(1 + q^6 + q^12)/(1 + q^6), -(q^3/(1 + q^6))]]): Racah[6,3][[2,1],[3],[3]]:=Matrix([[q^2/(1 + q^4), sqrt(1 + q^4 + q^8)/(1 + q^4)], [sqrt(1 + q^4 + q^8)/(1 + q^4), -(q^2/(1 + q^4))]]): Racah[6,2,1][[2,1],[3],[3]]:=Matrix([[q^6/((1 + q^6)*sqrt((1 + q^4)*(1 + q^2 + q^4 + q^6 + q^8))), (q^2*(1 + q^2 + q^4)*(1 + q^8))/sqrt((1 + q^4)*(1 + 2*q^2 + 2*q^4 + 3*q^6 + 5*q^8 + 5*q^10 + 6*q^12 + 6*q^14 + 6*q^16 + 5*q^18 + 5*q^20 + 3*q^22 + 2*q^24 + 2*q^26 + q^28)), (q^3 + q^9 + q^15)/((1 + q^6)*sqrt((1 + q^4)*(1 + 2*q^2 + 2*q^4 + 3*q^6 + 4*q^8 + 3*q^10 + 4*q^12 + 3*q^14 + 2*q^16 + 2*q^18 + q^20))), (1 + q^8)/sqrt(((1 + q^4)*(1 + q^2 + q^4 + q^6 + 2*q^8 + q^10 + q^12 + q^14 + q^16))/(1 + q^2 + q^4))], [(q^3*(1 + q^2 + q^4)^(3/2))/((1 + q^2)*sqrt(1 + q^4)*(1 + q^2 + q^4 + q^6 + q^8)), -((q^5*(1 + q^8))/sqrt(((1 + q^4)*(1 + 3*q^2 + 5*q^4 + 8*q^6 + 13*q^8 + 17*q^10 + 21*q^12 + 25*q^14 + 28*q^16 + 28*q^18 + 28*q^20 + 25*q^22 + 21*q^24 + 17*q^26 + 13*q^28 + 8*q^30 + 5*q^32 + 3*q^34 + q^36))/(1 + q^2 + q^4))), (1 + 2*q^2 + 2*q^4 + 2*q^6 + 3*q^8 + 2*q^10 + 2*q^12 + 2*q^14 + q^16)/((1 + q^2)*sqrt(((1 + q^4)*(1 + 3*q^2 + 5*q^4 + 8*q^6 + 12*q^8 + 14*q^10 + 16*q^12 + 17*q^14 + 16*q^16 + 14*q^18 + 12*q^20 + 8*q^22 + 5*q^24 + 3*q^26 + q^28))/(1 + q^2 + q^4))), -((q^3*sqrt((1 + q^8)/(1 + q^4)))/(1 + q^2 + q^4 + q^6 + q^8))], [(q^2*(1 + q^4 + q^8 + q^12))/((1 + q^6)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 3*q^8 + 2*q^10 + 3*q^12 + 2*q^14 + 2*q^16 + q^18 + q^20)), (1 + q^2 + q^6 + q^10 + q^12)/sqrt(1 + 2*q^2 + 3*q^4 + 5*q^6 + 6*q^8 + 6*q^10 + 8*q^12 + 6*q^14 + 6*q^16 + 5*q^18 + 3*q^20 + 2*q^22 + q^24), -((q^5*(1 + q^4 + q^8 + q^12))/((1 + q^6)*sqrt(1 + 2*q^2 + 3*q^4 + 5*q^6 + 7*q^8 + 8*q^10 + 11*q^12 + 11*q^14 + 12*q^16 + 11*q^18 + 11*q^20 + 8*q^22 + 7*q^24 + 5*q^26 + 3*q^28 + 2*q^30 + q^32))), -(q^2*sqrt((1 + q^2 + q^4)/(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)))], [((1 + q^2)*sqrt(1 + q^4 + q^8 + q^12))/(1 + q^2 + q^4 + q^6 + q^8), -((q^2*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/sqrt(1 + 3*q^2 + 6*q^4 + 11*q^6 + 17*q^8 + 22*q^10 + 28*q^12 + 31*q^14 + 32*q^16 + 31*q^18 + 28*q^20 + 22*q^22 + 17*q^24 + 11*q^26 + 6*q^28 + 3*q^30 + q^32)), -((q^3*(1 + q^2)*(1 + q^4 + q^8 + q^12))/sqrt(1 + 3*q^2 + 6*q^4 + 11*q^6 + 18*q^8 + 25*q^10 + 34*q^12 + 42*q^14 + 49*q^16 + 53*q^18 + 56*q^20 + 53*q^22 + 49*q^24 + 42*q^26 + 34*q^28 + 25*q^30 + 18*q^32 + 11*q^34 + 6*q^36 + 3*q^38 + q^40)), (q^4*sqrt((1 + q^2 + q^4)/(1 + q^4)))/(1 + q^2 + q^4 + q^6 + q^8)]]): Racah[5,4][[2,1],[3],[3]]:=Matrix([[(q*sqrt((1 + q^2 + q^4)/(1 + q^4)))/(1 + q^2), sqrt((1 + q^2 + q^4 + q^6 + q^8)/(1 + q^4))/(1 + q^2)], [(1 + q^2 + q^4 + q^6 + q^8)/((1 + q^2)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)), -((q*sqrt((1 + q^2 + q^4)/(1 + q^4)))/(1 + q^2))]]): Racah[5,3,1][[2,1],[3],[3]]:=Matrix([[q^4/(sqrt((1 + q^4)*(1 - q^2 + q^4))*(1 + q^2 + q^4)), (q + 2*q^3 + 3*q^5 + 4*q^7 + 5*q^9 + 5*q^11 + 5*q^13 + 4*q^15 + 3*q^17 + 2*q^19 + q^21)/((1 + q^2)*sqrt((1 + q^4)*(1 + 3*q^2 + 6*q^4 + 10*q^6 + 17*q^8 + 24*q^10 + 31*q^12 + 37*q^14 + 42*q^16 + 43*q^18 + 42*q^20 + 37*q^22 + 31*q^24 + 24*q^26 + 17*q^28 + 10*q^30 + 6*q^32 + 3*q^34 + q^36))), (q^2 + q^10)/sqrt((1 + q^4)*(1 + q^2 + q^4)*(1 + q^2 + q^4 + q^6 + 3*q^8 + q^10 + q^12 + q^14 + q^16)), (1 + 2*q^2 + 3*q^4 + 4*q^6 + 5*q^8 + 5*q^10 + 5*q^12 + 4*q^14 + 3*q^16 + 2*q^18 + q^20)/((1 + q^2)*(1 + q^2 + q^4)*sqrt((1 + q^4)*(1 + q^2 + 2*q^4 + 3*q^6 + 4*q^8 + 4*q^10 + 5*q^12 + 4*q^14 + 4*q^16 + 3*q^18 + 2*q^20 + q^22 + q^24)))], [q^2/sqrt((1 + q^4)*(1 + q^2 + q^4)), -((q^3*(1 + q^2 + q^4)*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/((1 + q^2)*sqrt((1 + q^2 + 2*q^4 + q^6 + q^8)*(1 + 4*q^2 + 9*q^4 + 15*q^6 + 23*q^8 + 32*q^10 + 40*q^12 + 45*q^14 + 47*q^16 + 45*q^18 + 40*q^20 + 32*q^22 + 23*q^24 + 15*q^26 + 9*q^28 + 4*q^30 + q^32)))), (1 + q^2 + q^6 + q^8)/sqrt((1 + q^4)*(1 + 2*q^2 + 2*q^4 + q^6 + 2*q^8 + 2*q^10 + q^12)), -((q^2*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/((1 + q^2)*sqrt((1 + q^2 + 2*q^4 + q^6 + q^8)*(1 + 2*q^2 + 3*q^4 + 4*q^6 + 5*q^8 + 5*q^10 + 5*q^12 + 4*q^14 + 3*q^16 + 2*q^18 + q^20))))], [(q + q^3 + q^5 + q^7 + q^9 + q^11 + q^13)/((1 + q^2 + q^4)*sqrt((1 + q^4)*(1 + q^4 + q^6 + q^8 + q^10 + q^12 + q^16))), (1 + q^2 + q^4 + q^6 + 2*q^8 + q^10 + q^12 + q^14 + q^16)/((1 + q^2)*sqrt((1 + q^4)*(1 + 2*q^2 + 3*q^4 + 4*q^6 + 7*q^8 + 7*q^10 + 7*q^12 + 7*q^14 + 7*q^16 + 4*q^18 + 3*q^20 + 2*q^22 + q^24))), -((q^3*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/sqrt((1 + q^4)*(1 + q^2 + q^4)*(1 + 2*q^2 + 3*q^4 + 4*q^6 + 7*q^8 + 8*q^10 + 9*q^12 + 9*q^14 + 9*q^16 + 8*q^18 + 7*q^20 + 4*q^22 + 3*q^24 + 2*q^26 + q^28))), -((q*(1 + q^4)*(1 + q^2 + q^4 + q^6 + q^8))/((1 + q^2)*(1 + q^2 + q^4)*sqrt(1 + 2*q^4 + q^6 + 2*q^8 + q^10 + 2*q^12 + q^16)))], [sqrt((1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/(1 + q^4))/(1 + q^2 + q^4), -((q*(1 - q^2 + q^4)*(1 + 2*q^2 + 2*q^4 + q^6))/sqrt((1 + q^4)*(1 + 3*q^2 + 5*q^4 + 6*q^6 + 8*q^8 + 9*q^10 + 8*q^12 + 6*q^14 + 5*q^16 + 3*q^18 + q^20))), -((q^2*(1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12))/sqrt((1 + q^4)*(1 + q^2 + q^4)*(1 + 3*q^2 + 5*q^4 + 6*q^6 + 8*q^8 + 10*q^10 + 11*q^12 + 10*q^14 + 8*q^16 + 6*q^18 + 5*q^20 + 3*q^22 + q^24))), (q^2*(1 + q^2)*(1 + q^4))/((1 + q^2 + q^4)*sqrt(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12))]]): Racah[5,2,2][[2,1],[3],[3]]:=Matrix([[q^2*sqrt((1 + q^2 + q^4)/(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)), (1 + q^2 + q^4 + q^6 + q^8 + q^10 + q^12)/sqrt(1 + 2*q^2 + 4*q^4 + 6*q^6 + 8*q^8 + 9*q^10 + 10*q^12 + 9*q^14 + 8*q^16 + 6*q^18 + 4*q^20 + 2*q^22 + q^24)], [(1 + q^4 + q^6 + q^8 + q^10 + q^12 + q^16)/sqrt((1 + q^4 + q^6 + q^8 + q^12)*(1 + 2*q^4 + q^6 + 2*q^8 + 2*q^10 + 2*q^12 + q^14 + 2*q^16 + q^20)), -(q^2*sqrt((1 + q^2 + q^4)/(1 + q^2 + 2*q^4 + 2*q^6 + 2*q^8 + q^10 + q^12)))]]): Racah[5,2,1,1][[2,1],[3],[3]]:=Matrix([[q^2/sqrt(1 + q^2 + 2*q^4 + q^6 + q^8), sqrt((1 + q^2 + q^4 + q^6 + q^8)/(1 + q^2 + 2*q^4 + q^6 + q^8))], [sqrt((1 + q^2 + q^4 + q^6 + q^8)/(1 + q^2 + 2*q^4 + q^6 + q^8)), -(q^2/sqrt(1 + q^2 + 2*q^4 + q^6 + q^8))]]): Racah[4,4,1][[2,1],[3],[3]]:=Matrix([[-(q/(1 + q^2)), sqrt(1 + q^2 + q^4)/(1 + q^2)], [sqrt(1 + q^2 + q^4)/(1 + q^2), q/(1 + q^2)]]): Racah[4,3,2][[2,1],[3],[3]]:=Matrix([[q*sqrt((1 + q^4)/(1 + q^2 + q^4 + q^6 + q^8)), sqrt((1 + q^4 + q^8)/(1 + q^2 + q^4 + q^6 + q^8))], [sqrt((1 + q^4 + q^8)/(1 + q^2 + q^4 + q^6 + q^8)), -(q*sqrt((1 + q^4)/(1 + q^2 + q^4 + q^6 + q^8)))]]): Racah[4,3,1,1][[2,1],[3],[3]]:=Matrix([[q/sqrt(1 + q^2 + q^4), sqrt((1 + q^4)/(1 + q^2 + q^4))], [sqrt((1 + q^4)/(1 + q^2 + q^4)), -(q/sqrt(1 + q^2 + q^4))]]):