with(LinearAlgebra): #spisok[[2],[2],[1,1]]:=Matrix([[4,1,1],[3,2,1]] :R[6,1][[1],[3]]:= DiagonalMatrix(): Racah[6,1][[3],[1],[3]]:=Matrix([[-((q + q^3 + q^5)/(1 + q^2 + q^4 + q^6)), sqrt(q^6 - q^8 - q^20 + q^22)/(q^3*(-1 + q^8))], [sqrt(q^6 - q^8 - q^20 + q^22)/(q^3*(-1 + q^8)), (q*(-1 + q^6))/(-1 + q^8)]]): R[5,2][[1],[3]]:= DiagonalMatrix(): Racah[5,2][[3],[1],[3]]:=Matrix([[q^2/(1 + q^4), -((sqrt(1 + q^4 + q^8)*(-1 + q^4))/(-1 + q^8))], [-((sqrt(1 + q^4 + q^8)*(-1 + q^4))/(-1 + q^8)), -(q^2/(1 + q^4))]]): R[4,3][[1],[3]]:=DiagonalMatrix(): Racah[4,3][[3],[1],[3]]:=Matrix([[-(q^3/(1 + q^2 + q^4 + q^6)), -(sqrt(q^2 - q^8 - q^12 + q^18)/(q - q^9))], [-(sqrt(q^2 - q^8 - q^12 + q^18)/(q - q^9)), q^3/(1 + q^2 + q^4 + q^6)]]): R[6,1][[3],[3]]:=DiagonalMatrix(): Racah[6,1][[3],[3],[1]]:=Matrix([[-(q^3/sqrt(1 + q^2 + q^4 + 2*q^6 + q^8 + q^10 + q^12)), sqrt((q^6 (-1 + q^14))/((-1 + q^8)*(q^6 + q^12)))], [sqrt((q^6 (-1 + q^14))/((-1 + q^8)*(q^6 + q^12))), q^3/sqrt(1 + q^2 + q^4 + 2*q^6 + q^8 + q^10 + q^12)]]): R[5,2][[3],[3]]:=DiagonalMatrix(<-q^3, q^(-1)>): Racah[5,2][[3],[3],[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[4,3][[3],[3]]:= DiagonalMatrix(): Racah[4,3][[3],[3],[1]]:=Matrix([[-((q*sqrt((1 + q^2 + q^4)/(1 + q^4)))/(1 + q^2)), sqrt((-1 + q^10)/((1 + q^2)*(-1 + q^8)))], [sqrt((-1 + q^10)/((1 + q^2)*(-1 + q^8))), sqrt((q^2 (-1 + q^6))/((1 + q^2)*(-1 + q^8)))]]):