SetDirectory["C:\\Users\\T15Roland\\Wiskunde\\Bn\\HigherRank"];
Once[<<KnotTheory`];
 <<Rot.m
<< FormalGaussianIntegration.m

(*The R3 solutions from UC4A242 (written hard coded below the fold):*)

sol3bα={(Subscript[c, 1])->0,(Subscript[c, 2])->0,(Subscript[c, 3])->0,(Subscript[c, 6])->0,(Subscript[c, 7])->0,(Subscript[c, 8])->0,(Subscript[c, 11])->0,(Subscript[c, 12])->0,(Subscript[c, 14])->(-(Subscript[c, 9])+(Subscript[c, 10])-(Subscript[c, 13])+(Subscript[c, 28])+(Subscript[c, 32])+(Subscript[c, 9]) (Subscript[T, 1])-(Subscript[c, 10]) (Subscript[T, 1])+(Subscript[c, 13]) (Subscript[T, 1])-(Subscript[c, 22]) (Subscript[T, 1])+(Subscript[c, 22]) (Subscript[T, 1]^2)-(Subscript[c, 10]) (Subscript[T, 2])+(Subscript[c, 13]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 2])-(Subscript[c, 32]) (Subscript[T, 2])+(Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 13]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 32]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 32]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[α, 4])+(Subscript[T, 1]) (Subscript[α, 4]))/((-1+(Subscript[T, 1])) (-1+(Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))),(Subscript[c, 16])->0,(Subscript[c, 17])->0,(Subscript[c, 18])->0,(Subscript[c, 19])->0,(Subscript[c, 20])->0,(Subscript[c, 21])->0,(Subscript[c, 23])->0,(Subscript[c, 24])->0,(Subscript[c, 25])->(1/2) (-2 (Subscript[c, 4])+2 (Subscript[c, 4]) (Subscript[T, 1])-(Subscript[c, 5]) (Subscript[T, 1])+(Subscript[c, 5]) (Subscript[T, 1]^2)),(Subscript[c, 26])->(1/2) (-2 (Subscript[c, 4])-(Subscript[c, 5])-(Subscript[c, 5]) (Subscript[T, 1])),(Subscript[c, 27])->0,(Subscript[c, 29])->0,(Subscript[c, 30])->(1/(-1+(Subscript[T, 2])))((Subscript[c, 9])+(Subscript[c, 28])-(Subscript[c, 9]) (Subscript[T, 1])-(Subscript[c, 22]) (Subscript[T, 1])-2 (Subscript[c, 9]) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 2])+2 (Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 2]^2)+(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[T, 2]) (Subscript[α, 4])),(Subscript[c, 31])->(((Subscript[c, 9])-(Subscript[c, 9]) (Subscript[T, 1])+(Subscript[α, 4]))/(-1+(Subscript[T, 2]))),(Subscript[c, 33])->0,(Subscript[c, 34])->(1/(-1+(Subscript[T, 2])))(-(Subscript[c, 9])+(Subscript[c, 10])+(Subscript[c, 13])+(Subscript[c, 28])+(Subscript[c, 32])+(Subscript[c, 9]) (Subscript[T, 1])-(Subscript[c, 10]) (Subscript[T, 1])-(Subscript[c, 13]) (Subscript[T, 1])-(Subscript[c, 10]) (Subscript[T, 2])-(Subscript[c, 13]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 2])-(Subscript[c, 32]) (Subscript[T, 2])+(Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 13]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 32]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 32]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[α, 4])+(Subscript[T, 1]) (Subscript[T, 2]) (Subscript[α, 4])),(Subscript[c, 35])->(((Subscript[c, 13])-(Subscript[c, 13]) (Subscript[T, 1])+(Subscript[c, 22]) (Subscript[T, 1])-(Subscript[c, 22]) (Subscript[T, 1]^2)-(Subscript[T, 1]) (Subscript[α, 4]))/(-1+(Subscript[T, 1]) (Subscript[T, 2]))),(Subscript[c, 37])->(((Subscript[c, 9])-(Subscript[c, 10])-(Subscript[c, 28])-2 (Subscript[c, 9]) (Subscript[T, 1])+(Subscript[c, 10]) (Subscript[T, 1])+(Subscript[c, 28]) (Subscript[T, 1])+(Subscript[c, 9]) (Subscript[T, 1]^2)+(Subscript[c, 10]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[α, 3])+(Subscript[T, 2]) (Subscript[α, 3])+(Subscript[α, 4])-(Subscript[T, 1]) (Subscript[α, 4]))/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (-1+(Subscript[T, 2])))),(Subscript[c, 38])->0,(Subscript[c, 39])->(1/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (-1+(Subscript[T, 2]))))(-(Subscript[c, 9])+(Subscript[c, 10])+(Subscript[c, 28])+(Subscript[c, 9]) (Subscript[T, 1])-(Subscript[c, 10]) (Subscript[T, 1])-(Subscript[c, 22]) (Subscript[T, 1])-(Subscript[c, 28]) (Subscript[T, 1])+(Subscript[c, 22]) (Subscript[T, 1]^2)-(Subscript[c, 10]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 2])+(Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)+(Subscript[α, 3])-(Subscript[T, 2]) (Subscript[α, 3])-(Subscript[α, 4])+(Subscript[T, 1]) (Subscript[α, 4])),(Subscript[c, 40])->0,(Subscript[c, 41])->((Subscript[c, 9])-(Subscript[c, 10])-(Subscript[c, 28])-(Subscript[c, 9]) (Subscript[T, 1])+(Subscript[c, 10]) (Subscript[T, 1])+(Subscript[c, 13]) (Subscript[T, 1])-(Subscript[c, 32]) (Subscript[T, 1])-(Subscript[c, 13]) (Subscript[T, 1]^2)+(Subscript[c, 22]) (Subscript[T, 1]^2)-(Subscript[c, 22]) (Subscript[T, 1]^3)+(Subscript[c, 10]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 2])-(Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 13]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 32]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 13]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 32]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]^3) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[c, 32]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[α, 3])+(Subscript[T, 2]) (Subscript[α, 3])+(Subscript[T, 1]) (Subscript[T, 2]) (Subscript[α, 3])-(Subscript[T, 1]) (Subscript[T, 2]^2) (Subscript[α, 3])+(Subscript[α, 4])-(Subscript[T, 1]^2) (Subscript[α, 4])-(Subscript[T, 1]) (Subscript[T, 2]) (Subscript[α, 4])+(Subscript[T, 1]^2) (Subscript[T, 2]) (Subscript[α, 4]))/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (-1+(Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))),(Subscript[c, 42])->0,(Subscript[c, 43])->(1/((-1+(Subscript[T, 1])) (-1+(Subscript[T, 2]))))((Subscript[c, 9])-(Subscript[c, 10])-(Subscript[c, 22])-(Subscript[c, 28])-(Subscript[c, 32])-(Subscript[c, 9]) (Subscript[T, 1])+(Subscript[c, 10]) (Subscript[T, 1])+2 (Subscript[c, 22]) (Subscript[T, 1])-(Subscript[c, 22]) (Subscript[T, 1]^2)+(Subscript[c, 10]) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 2])+(Subscript[c, 32]) (Subscript[T, 2])-(Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])-2 (Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 32]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 32]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[α, 3])+(Subscript[T, 2]) (Subscript[α, 3])+(Subscript[α, 4])-(Subscript[T, 1]) (Subscript[α, 4])),(Subscript[c, 44])->0,(Subscript[c, 45])->((-(Subscript[c, 15])-(Subscript[c, 36]))/(Subscript[T, 1])),(Subscript[c, 46])->0,(Subscript[c, 47])->0,(Subscript[c, 48])->0,(Subscript[c, 51])->0,(Subscript[c, 52])->0,(Subscript[c, 53])->0,(Subscript[c, 55])->(((Subscript[c, 22])+(Subscript[c, 54])-(Subscript[c, 65])-(Subscript[c, 22]) (Subscript[T, 2])-(Subscript[c, 54]) (Subscript[T, 2])+(Subscript[c, 65]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[T, 2]) (Subscript[α, 4]))/((-1+(Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2])))),(Subscript[c, 57])->0,(Subscript[c, 58])->0,(Subscript[c, 59])->-((Subscript[α, 3])/(-1+(Subscript[T, 1]))),(Subscript[c, 60])->0,(Subscript[c, 61])->0,(Subscript[c, 62])->(1/2) (-2 (Subscript[c, 49])+2 (Subscript[c, 49]) (Subscript[T, 2])-(Subscript[c, 50]) (Subscript[T, 2])+(Subscript[c, 50]) (Subscript[T, 2]^2)),(Subscript[c, 63])->(1/2) (-2 (Subscript[c, 49])-(Subscript[c, 50])-(Subscript[c, 50]) (Subscript[T, 2])),(Subscript[c, 64])->0,(Subscript[c, 66])->0,(Subscript[c, 67])->(1/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (-1+(Subscript[T, 2]))))((Subscript[c, 22]) (Subscript[T, 1])-(Subscript[c, 54]) (Subscript[T, 1])-(Subscript[c, 65]) (Subscript[T, 1])-(Subscript[c, 22]) (Subscript[T, 1]^2)+(Subscript[c, 54]) (Subscript[T, 1]^2)+(Subscript[c, 65]) (Subscript[T, 1]^2)+(Subscript[c, 9]) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 2])-2 (Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])+2 (Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])+2 (Subscript[c, 54]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 65]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2])-2 (Subscript[c, 54]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]^3) (Subscript[T, 2])-(Subscript[c, 65]) (Subscript[T, 1]^3) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 2]^2)+(Subscript[c, 28]) (Subscript[T, 2]^2)-(Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 54]) (Subscript[T, 1]) (Subscript[T, 2]^2)+2 (Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)+(Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[c, 28]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)+(Subscript[c, 54]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[c, 65]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[c, 9]) (Subscript[T, 1]^3) (Subscript[T, 2]^2)+(Subscript[c, 10]) (Subscript[T, 1]^3) (Subscript[T, 2]^2)-(Subscript[c, 22]) (Subscript[T, 1]^3) (Subscript[T, 2]^2)+(Subscript[c, 65]) (Subscript[T, 1]^3) (Subscript[T, 2]^2)-(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^3)-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2]^3)+2 (Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^3)+(Subscript[c, 28]) (Subscript[T, 1]^2) (Subscript[T, 2]^3)-(Subscript[c, 10]) (Subscript[T, 1]^3) (Subscript[T, 2]^3)-(Subscript[T, 2]) (Subscript[α, 3])+(Subscript[T, 2]^2) (Subscript[α, 3])+(Subscript[T, 1]) (Subscript[T, 2]^2) (Subscript[α, 3])-(Subscript[T, 1]) (Subscript[T, 2]^3) (Subscript[α, 3])+(Subscript[T, 2]) (Subscript[α, 4])-(Subscript[T, 1]) (Subscript[T, 2]) (Subscript[α, 4])-(Subscript[T, 1]) (Subscript[T, 2]^2) (Subscript[α, 4])+(Subscript[T, 1]^2) (Subscript[T, 2]^2) (Subscript[α, 4])),(Subscript[c, 68])->((Subscript[c, 54]) (Subscript[T, 1])-(Subscript[c, 54]) (Subscript[T, 1]^2)-(Subscript[c, 9]) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 2])+2 (Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])-2 (Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])-2 (Subscript[c, 54]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2])+2 (Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2])+2 (Subscript[c, 54]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 22]) (Subscript[T, 1]^3) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 2]^2)-(Subscript[c, 28]) (Subscript[T, 2]^2)+(Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 54]) (Subscript[T, 1]) (Subscript[T, 2]^2)-2 (Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)+(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-2 (Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)+(Subscript[c, 28]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[c, 54]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)+(Subscript[c, 9]) (Subscript[T, 1]^3) (Subscript[T, 2]^2)-(Subscript[c, 10]) (Subscript[T, 1]^3) (Subscript[T, 2]^2)+(Subscript[c, 22]) (Subscript[T, 1]^3) (Subscript[T, 2]^2)+(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^3)+(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2]^3)-2 (Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^3)-(Subscript[c, 28]) (Subscript[T, 1]^2) (Subscript[T, 2]^3)+(Subscript[c, 10]) (Subscript[T, 1]^3) (Subscript[T, 2]^3)+(Subscript[T, 2]) (Subscript[α, 3])-(Subscript[T, 2]^2) (Subscript[α, 3])-(Subscript[T, 1]) (Subscript[T, 2]^2) (Subscript[α, 3])+(Subscript[T, 1]) (Subscript[T, 2]^3) (Subscript[α, 3])-(Subscript[T, 2]) (Subscript[α, 4])+2 (Subscript[T, 1]) (Subscript[T, 2]) (Subscript[α, 4])-(Subscript[T, 1]^2) (Subscript[T, 2]) (Subscript[α, 4]))/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (-1+(Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))),(Subscript[c, 70])->((Subscript[c, 9])-(Subscript[c, 10])-(Subscript[c, 28])-2 (Subscript[c, 9]) (Subscript[T, 1])+2 (Subscript[c, 10]) (Subscript[T, 1])+(Subscript[c, 22]) (Subscript[T, 1])+(Subscript[c, 28]) (Subscript[T, 1])+(Subscript[c, 54]) (Subscript[T, 1])-(Subscript[c, 65]) (Subscript[T, 1])+(Subscript[c, 9]) (Subscript[T, 1]^2)-(Subscript[c, 10]) (Subscript[T, 1]^2)-(Subscript[c, 22]) (Subscript[T, 1]^2)-(Subscript[c, 54]) (Subscript[T, 1]^2)+(Subscript[c, 65]) (Subscript[T, 1]^2)+(Subscript[c, 10]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 2])-(Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 54]) (Subscript[T, 1]) (Subscript[T, 2])+2 (Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 54]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 65]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 9]) (Subscript[T, 1]^3) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 1]^3) (Subscript[T, 2])-(Subscript[c, 65]) (Subscript[T, 1]^3) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2]^2)+2 (Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)+(Subscript[c, 28]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[c, 10]) (Subscript[T, 1]^3) (Subscript[T, 2]^2)-(Subscript[α, 3])+(Subscript[T, 2]) (Subscript[α, 3])+(Subscript[T, 1]) (Subscript[T, 2]) (Subscript[α, 3])-(Subscript[T, 1]) (Subscript[T, 2]^2) (Subscript[α, 3])+(Subscript[α, 4])-(Subscript[T, 1]) (Subscript[α, 4]))/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (-1+(Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))),(Subscript[c, 71])->0,(Subscript[c, 72])->((-(Subscript[c, 65])+(Subscript[c, 65]) (Subscript[T, 1])+(Subscript[c, 65]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 65]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[α, 3])+(Subscript[T, 2]) (Subscript[α, 3])+(Subscript[T, 2]) (Subscript[α, 4])-(Subscript[T, 1]) (Subscript[T, 2]) (Subscript[α, 4]))/((-1+(Subscript[T, 1])) (-1+(Subscript[T, 2])))),(Subscript[c, 73])->0,(Subscript[c, 74])->((-(Subscript[c, 56])-(Subscript[c, 69]))/(Subscript[T, 2])),(Subscript[c, 75])->0,(Subscript[c, 76])->0,(Subscript[c, 77])->0,(Subscript[c, 80])->0,(Subscript[c, 82])->0,(Subscript[c, 83])->(1/2) (-2 (Subscript[c, 78])+2 (Subscript[c, 78]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 79]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 79]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)),(Subscript[c, 84])->(1/2) (-2 (Subscript[c, 78])-(Subscript[c, 79])-(Subscript[c, 79]) (Subscript[T, 1]) (Subscript[T, 2])),(Subscript[c, 85])->0,(Subscript[c, 87])->(1/((Subscript[T, 1]^2) (-1+(Subscript[T, 2])) (Subscript[T, 2])))(-(Subscript[c, 9])+(Subscript[c, 10])+(Subscript[c, 28])+(Subscript[c, 9]) (Subscript[T, 1])-(Subscript[c, 10]) (Subscript[T, 1])-(Subscript[c, 22]) (Subscript[T, 1])+(Subscript[c, 81]) (Subscript[T, 1])+(Subscript[c, 86]) (Subscript[T, 1])-(Subscript[c, 10]) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 2])+(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 81]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 86]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[α, 3])-(Subscript[T, 2]) (Subscript[α, 3])-(Subscript[α, 4])),(Subscript[α, 1])->(1/(-1+(Subscript[T, 2])))((Subscript[c, 9])-(Subscript[c, 10])-(Subscript[c, 28])-(Subscript[c, 9]) (Subscript[T, 1])+(Subscript[c, 10]) (Subscript[T, 1])+(Subscript[c, 22]) (Subscript[T, 1])-(Subscript[c, 22]) (Subscript[T, 1]^2)+(Subscript[c, 10]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 2])-(Subscript[c, 9]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 9]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 10]) (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 10]) (Subscript[T, 1]^2) (Subscript[T, 2]^2)-(Subscript[α, 3])+2 (Subscript[T, 2]) (Subscript[α, 3])-(Subscript[T, 2]^2) (Subscript[α, 3])+(Subscript[α, 4])-(Subscript[T, 1]) (Subscript[α, 4])),(Subscript[α, 2])->-(Subscript[c, 22]) (Subscript[T, 1])+(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[T, 2]) (Subscript[α, 4])};

sol3b={$A->-$B (Subscript[T, 2])+(($B (Subscript[c, 22]) (Subscript[T, 1]) (-1+(Subscript[T, 2])))/(Subscript[α, 4])),$C->((Subscript[α, 4]) (-(Subscript[c, 28])+(Subscript[c, 22]) (Subscript[T, 1])-(Subscript[c, 22]) (Subscript[T, 1]^2)+(Subscript[c, 28]) (Subscript[T, 2])-(Subscript[c, 22]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[c, 22]) (Subscript[T, 1]^2) (Subscript[T, 2])-(Subscript[c, 28]) (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[c, 9]) (-1+(Subscript[T, 1])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))+(Subscript[c, 10]) (-1+(Subscript[T, 1])) (-1+(Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))+(Subscript[α, 4])-(Subscript[T, 1]) (Subscript[α, 4])))/($B (-1+(Subscript[T, 2])) ((Subscript[c, 22]) (Subscript[T, 1]) (-1+(Subscript[T, 2]))-(Subscript[α, 4]))),$D->((Subscript[α, 4])/$B)};
{sα3}=Solve[{(Subscript[α, 1])(Subscript[α, 4])-(Subscript[α, 2])(Subscript[α, 3])==0}/.sol3bα[[-2;;-1]],{(Subscript[α, 3]),(Subscript[α, 4])}];

r0p/.Sub//CF
r0m/.Sub//CF

r1p/.Sub//CF

r1m/.Sub//CF

(Subscript[r, 0])[1,i_,j_]:=2 (-(Subscript[T, 1])+(Subscript[T, 2])+(Subscript[T, 1]) (Subscript[T, 2])) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 2,i])-2 (-1+(Subscript[T, 2])) (Subscript[p, 3,j]) (Subscript[x, 1,j]) (Subscript[x, 2,i])-2 (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 2,j])(*from r0p*)
(Subscript[r, 0])[-1,i_,j_]:=-((2 (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 2,i]))/(Subscript[T, 1]))+((2 (-1+(Subscript[T, 2])) (Subscript[p, 3,j]) (Subscript[x, 1,j]) (Subscript[x, 2,i]))/(Subscript[T, 2]))+((2 (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 2,j]))/(Subscript[T, 1]))
(Subscript[r, 1])[1,i_,j_]:=-(3/2)-2 (Subscript[p, 1,i]) (Subscript[x, 1,i])-2 (Subscript[p, 1,j]) (Subscript[x, 1,i])-(Subscript[p, 1,i]) (Subscript[p, 1,j]) (Subscript[x, 1,i]^2)-(1/2) (-1+(Subscript[T, 1])) (2+(Subscript[T, 1])) (Subscript[p, 1,j]^2) (Subscript[x, 1,i]^2)+((4 (Subscript[p, 1,j]) (Subscript[x, 1,j]))/(Subscript[T, 1]))-(Subscript[p, 1,i]) (Subscript[p, 1,j]) (Subscript[x, 1,i]) (Subscript[x, 1,j])+(1/2) (3+(Subscript[T, 1])) (Subscript[p, 1,j]^2) (Subscript[x, 1,i]) (Subscript[x, 1,j])-2 (Subscript[p, 2,i]) (Subscript[x, 2,i])+2 (Subscript[p, 2,j]) (Subscript[x, 2,i])-2 (Subscript[p, 1,j]) (Subscript[p, 2,i]) (Subscript[x, 1,i]) (Subscript[x, 2,i])+(((1-3 (Subscript[T, 1])+2 (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2])+9 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 2]^2)+11 (Subscript[T, 1]) (Subscript[T, 2]^2)-18 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-4 (Subscript[T, 1]) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^3)) (Subscript[p, 1,i]) (Subscript[p, 2,j]) (Subscript[x, 1,i]) (Subscript[x, 2,i]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))+(1+4 (Subscript[T, 1])-6 (Subscript[T, 2])+9 (Subscript[T, 1]) (Subscript[T, 2])-9 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 1,i]) (Subscript[x, 2,i])+(((-4+3 (Subscript[T, 1])+16 (Subscript[T, 1]) (Subscript[T, 2])-9 (Subscript[T, 1]^2) (Subscript[T, 2])+8 (Subscript[T, 2]^2)-22 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,i]) (Subscript[x, 1,j]) (Subscript[x, 2,i]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))-(((-1+(Subscript[T, 2])) (2+(Subscript[T, 1])-11 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 1,j]) (Subscript[x, 2,i]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))-2 (Subscript[p, 2,i]) (Subscript[p, 2,j]) (Subscript[x, 2,i]^2)-(1/2) (-1+(Subscript[T, 2])) (4+(Subscript[T, 2])) (Subscript[p, 2,j]^2) (Subscript[x, 2,i]^2)+2 (Subscript[p, 1,i]) (Subscript[p, 2,j]) (Subscript[x, 1,i]) (Subscript[x, 2,j])-(((3 (Subscript[T, 1])-4 (Subscript[T, 2])+7 (Subscript[T, 1]) (Subscript[T, 2])-9 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 1,i]) (Subscript[x, 2,j]))/(-1+2 (Subscript[T, 2])))-(Subscript[p, 2,i]) (Subscript[p, 2,j]) (Subscript[x, 2,i]) (Subscript[x, 2,j])+(1/2) (5+(Subscript[T, 2])) (Subscript[p, 2,j]^2) (Subscript[x, 2,i]) (Subscript[x, 2,j])+((((Subscript[T, 1])+4 (Subscript[T, 2])-11 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,i]) (Subscript[x, 3,i]))/(2 (-1+2 (Subscript[T, 2]))))+(1/2) (Subscript[p, 1,i]) (Subscript[p, 2,j]) (Subscript[x, 3,i])-(((-1+(Subscript[T, 1])+6 (Subscript[T, 2])-11 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 3,i]))/(2 (-1+2 (Subscript[T, 2]))))-2 (Subscript[p, 3,i]) (Subscript[x, 3,i])+(Subscript[p, 3,j]) (Subscript[x, 3,i])+2 (Subscript[p, 1,j]) (Subscript[p, 3,i]) (Subscript[x, 1,i]) (Subscript[x, 3,i])-2 (Subscript[p, 1,i]) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 3,i])-(((-3 (Subscript[T, 1])+4 (Subscript[T, 2])-5 (Subscript[T, 1]) (Subscript[T, 2])+10 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+2 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-4 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 1,j]) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 3,i]))/(-1+2 (Subscript[T, 2])))+(((6-4 (Subscript[T, 1])+(Subscript[T, 1]^2)+4 (Subscript[T, 2])-13 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 3,i]) (Subscript[x, 1,j]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))))+(Subscript[p, 1,i]) (Subscript[p, 3,j]) (Subscript[x, 1,j]) (Subscript[x, 3,i])+(((-3+(Subscript[T, 1])-(Subscript[T, 1]^2)+6 (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2])+12 (Subscript[T, 1]^2) (Subscript[T, 2])-7 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-4 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 1,j]) (Subscript[p, 3,j]) (Subscript[x, 1,j]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))+(Subscript[p, 2,j]) (Subscript[p, 3,i]) (Subscript[x, 2,i]) (Subscript[x, 3,i])-(Subscript[p, 2,i]) (Subscript[p, 3,j]) (Subscript[x, 2,i]) (Subscript[x, 3,i])-(1/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))(-1+(Subscript[T, 1])+5 (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2])+2 (Subscript[T, 2]^2)-20 (Subscript[T, 1]) (Subscript[T, 2]^2)+10 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-8 (Subscript[T, 2]^3)+18 (Subscript[T, 1]) (Subscript[T, 2]^3)+11 (Subscript[T, 1]^2) (Subscript[T, 2]^3)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)+8 (Subscript[T, 1]) (Subscript[T, 2]^4)-22 (Subscript[T, 1]^2) (Subscript[T, 2]^4)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^4)) (Subscript[p, 2,j]) (Subscript[p, 3,j]) (Subscript[x, 2,i]) (Subscript[x, 3,i])-((-3+2 (Subscript[T, 1])-(Subscript[T, 2])+21 (Subscript[T, 1]) (Subscript[T, 2])-11 (Subscript[T, 1]^2) (Subscript[T, 2])+14 (Subscript[T, 2]^2)-44 (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[T, 1]^2) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-8 (Subscript[T, 2]^3)+10 (Subscript[T, 1]) (Subscript[T, 2]^3)+33 (Subscript[T, 1]^2) (Subscript[T, 2]^3)-18 (Subscript[T, 1]^3) (Subscript[T, 2]^3)+8 (Subscript[T, 1]) (Subscript[T, 2]^4)-22 (Subscript[T, 1]^2) (Subscript[T, 2]^4)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^4)) (Subscript[p, 2,j]) (Subscript[p, 3,i]) (Subscript[x, 2,j]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (-1+(Subscript[T, 2])) (-1+2 (Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2])))+((((Subscript[T, 1])+4 (Subscript[T, 2])-11 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 2,i]) (Subscript[p, 3,j]) (Subscript[x, 2,j]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))-(((-1-3 (Subscript[T, 2])+12 (Subscript[T, 1]) (Subscript[T, 2])-(Subscript[T, 1]^2) (Subscript[T, 2])+6 (Subscript[T, 2]^2)-11 (Subscript[T, 1]) (Subscript[T, 2]^2)-7 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-4 (Subscript[T, 1]) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^3)) (Subscript[p, 2,j]) (Subscript[p, 3,j]) (Subscript[x, 2,j]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (-1+(Subscript[T, 2])) (-1+2 (Subscript[T, 2]))))-(Subscript[p, 3,i]) (Subscript[p, 3,j]) (Subscript[x, 3,i]^2)+(1/(2 (-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))(-1+(Subscript[T, 1]) (Subscript[T, 2])) (-2+2 (Subscript[T, 1])+4 (Subscript[T, 2])-7 (Subscript[T, 1]) (Subscript[T, 2])+2 (Subscript[T, 1]^2) (Subscript[T, 2])-2 (Subscript[T, 1]) (Subscript[T, 2]^2)+18 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)+8 (Subscript[T, 1]) (Subscript[T, 2]^3)-22 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 3,j]^2) (Subscript[x, 3,i]^2)+(((-4+9 (Subscript[T, 1])) (Subscript[p, 3,j]) (Subscript[x, 3,j]))/(Subscript[T, 1]))-((-4+3 (Subscript[T, 1])-(Subscript[T, 1]^2)+12 (Subscript[T, 2])-17 (Subscript[T, 1]) (Subscript[T, 2])+12 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+2 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-4 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 1,i]) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 3,j]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2])))-(((2-4 (Subscript[T, 1])+(Subscript[T, 1]^2)) (Subscript[p, 1,j]) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 3,j]))/(-1+(Subscript[T, 1]) (Subscript[T, 2])))+(((-1-(Subscript[T, 2])+(Subscript[T, 1]) (Subscript[T, 2])) (Subscript[p, 2,i]) (Subscript[p, 3,j]) (Subscript[x, 2,i]) (Subscript[x, 3,j]))/((-1+(Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))))+((-1+(Subscript[T, 2])) (1-(Subscript[T, 1])-2 (Subscript[T, 2])+(Subscript[T, 1]^2) (Subscript[T, 2])-8 (Subscript[T, 2]^2)+26 (Subscript[T, 1]) (Subscript[T, 2]^2)-11 (Subscript[T, 1]^2) (Subscript[T, 2]^2)+8 (Subscript[T, 1]) (Subscript[T, 2]^3)-22 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 2,j]) (Subscript[p, 3,j]) (Subscript[x, 2,i]) (Subscript[x, 3,j]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2])))+(((-3+2 (Subscript[T, 1])-2 (Subscript[T, 2])+18 (Subscript[T, 1]) (Subscript[T, 2])-9 (Subscript[T, 1]^2) (Subscript[T, 2])+8 (Subscript[T, 2]^2)-22 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 3,i]) (Subscript[p, 3,j]) (Subscript[x, 3,i]) (Subscript[x, 3,j]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))-(1/(2 (-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))(-5+4 (Subscript[T, 1])+2 (Subscript[T, 2])+11 (Subscript[T, 1]) (Subscript[T, 2])-7 (Subscript[T, 1]^2) (Subscript[T, 2])+8 (Subscript[T, 2]^2)-24 (Subscript[T, 1]) (Subscript[T, 2]^2)+27 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)+8 (Subscript[T, 1]) (Subscript[T, 2]^3)-22 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 3,j]^2) (Subscript[x, 3,i]) (Subscript[x, 3,j])
(Subscript[r, 1])[-1,i_,j_]:=
(3/2)+2 (Subscript[p, 1,i]) (Subscript[x, 1,i])-((2 (-2+(Subscript[T, 1]^2)) (Subscript[p, 1,j]) (Subscript[x, 1,i]))/(Subscript[T, 1]^2))+(Subscript[p, 1,i]) (Subscript[p, 1,j]) (Subscript[x, 1,i]^2)-(((-1+(Subscript[T, 1])) (1+2 (Subscript[T, 1])) (Subscript[p, 1,j]^2) (Subscript[x, 1,i]^2))/(2 (Subscript[T, 1]^2)))-((4 (Subscript[p, 1,j]) (Subscript[x, 1,j]))/(Subscript[T, 1]))+(Subscript[p, 1,i]) (Subscript[p, 1,j]) (Subscript[x, 1,i]) (Subscript[x, 1,j])-(((1+3 (Subscript[T, 1])) (Subscript[p, 1,j]^2) (Subscript[x, 1,i]) (Subscript[x, 1,j]))/(2 (Subscript[T, 1])))+2 (Subscript[p, 2,i]) (Subscript[x, 2,i])-2 (Subscript[p, 2,j]) (Subscript[x, 2,i])-(((-2+3 (Subscript[T, 1])-4 (Subscript[T, 2])+16 (Subscript[T, 1]) (Subscript[T, 2])-9 (Subscript[T, 1]^2) (Subscript[T, 2])+8 (Subscript[T, 2]^2)-22 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,i]) (Subscript[x, 1,i]) (Subscript[x, 2,i]))/((Subscript[T, 1]) (-1+2 (Subscript[T, 2]))))-(((-1-(Subscript[T, 1])+8 (Subscript[T, 2])-10 (Subscript[T, 1]) (Subscript[T, 2])+9 (Subscript[T, 1]^2) (Subscript[T, 2])-8 (Subscript[T, 2]^2)+15 (Subscript[T, 1]) (Subscript[T, 2]^2)-18 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-4 (Subscript[T, 1]) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^3)) (Subscript[p, 1,i]) (Subscript[p, 2,j]) (Subscript[x, 1,i]) (Subscript[x, 2,i]))/((-1+(Subscript[T, 1])) (Subscript[T, 2]) (-1+2 (Subscript[T, 2]))))+(1/((Subscript[T, 1]) (Subscript[T, 2]) (-1+2 (Subscript[T, 2]))))(-4 (Subscript[T, 1])+8 (Subscript[T, 2])-14 (Subscript[T, 1]) (Subscript[T, 2])+18 (Subscript[T, 1]^2) (Subscript[T, 2])-12 (Subscript[T, 2]^2)+35 (Subscript[T, 1]) (Subscript[T, 2]^2)-36 (Subscript[T, 1]^2) (Subscript[T, 2]^2)+8 (Subscript[T, 2]^3)-26 (Subscript[T, 1]) (Subscript[T, 2]^3)+18 (Subscript[T, 1]^2) (Subscript[T, 2]^3)) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 1,i]) (Subscript[x, 2,i])-(((-4+3 (Subscript[T, 1])+16 (Subscript[T, 1]) (Subscript[T, 2])-9 (Subscript[T, 1]^2) (Subscript[T, 2])+8 (Subscript[T, 2]^2)-22 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,i]) (Subscript[x, 1,j]) (Subscript[x, 2,i]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))+(((-1+(Subscript[T, 2])) (2-4 (Subscript[T, 1])+4 (Subscript[T, 2])-13 (Subscript[T, 1]) (Subscript[T, 2])+9 (Subscript[T, 1]^2) (Subscript[T, 2])) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 1,j]) (Subscript[x, 2,i]))/((-1+(Subscript[T, 1])) (Subscript[T, 2])))+(((1+(Subscript[T, 2])) (Subscript[p, 2,i]) (Subscript[p, 2,j]) (Subscript[x, 2,i]^2))/(Subscript[T, 2]))-(((-1+(Subscript[T, 2])) (3+2 (Subscript[T, 2])) (Subscript[p, 2,j]^2) (Subscript[x, 2,i]^2))/(2 (Subscript[T, 2]^2)))-2 (Subscript[p, 1,i]) (Subscript[p, 2,j]) (Subscript[x, 1,i]) (Subscript[x, 2,j])+(((2+(Subscript[T, 1])-8 (Subscript[T, 2])+11 (Subscript[T, 1]) (Subscript[T, 2])-9 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 1,i]) (Subscript[x, 2,j]))/((Subscript[T, 1]) (-1+2 (Subscript[T, 2]))))+(Subscript[p, 2,i]) (Subscript[p, 2,j]) (Subscript[x, 2,i]) (Subscript[x, 2,j])-((3 (1+(Subscript[T, 2])) (Subscript[p, 2,j]^2) (Subscript[x, 2,i]) (Subscript[x, 2,j]))/(2 (Subscript[T, 2])))-((((Subscript[T, 1])+4 (Subscript[T, 2])-11 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 2,i]) (Subscript[x, 3,i]))/(2 (Subscript[T, 1]) (-1+2 (Subscript[T, 2]))))-(((Subscript[p, 1,i]) (Subscript[p, 2,j]) (Subscript[x, 3,i]))/(2 (Subscript[T, 2])))+(((-(Subscript[T, 1])+3 (Subscript[T, 1]) (Subscript[T, 2])+4 (Subscript[T, 2]^2)-11 (Subscript[T, 1]) (Subscript[T, 2]^2)-4 (Subscript[T, 1]) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^3)) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 3,i]))/(2 (Subscript[T, 1]) (Subscript[T, 2]) (-1+2 (Subscript[T, 2]))))+2 (Subscript[p, 3,i]) (Subscript[x, 3,i])-(((-4+9 (Subscript[T, 1])-(Subscript[T, 1]^2)+8 (Subscript[T, 2])-18 (Subscript[T, 1]) (Subscript[T, 2])+(Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+15 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-4 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 3,j]) (Subscript[x, 3,i]))/((Subscript[T, 1]^2) (Subscript[T, 2]) (-1+2 (Subscript[T, 2]))))-(((4-4 (Subscript[T, 1])+(Subscript[T, 1]^2)+4 (Subscript[T, 2])-11 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 3,i]) (Subscript[x, 1,i]) (Subscript[x, 3,i]))/((Subscript[T, 1]) (-1+(Subscript[T, 1]) (Subscript[T, 2]))))+(((-2 (Subscript[T, 1])+4 (Subscript[T, 2])-7 (Subscript[T, 1]) (Subscript[T, 2])+10 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+2 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-4 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 1,i]) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (Subscript[T, 2]) (-1+2 (Subscript[T, 2]))))-(((2-6 (Subscript[T, 1])+(Subscript[T, 1]^2)+(Subscript[T, 1]) (Subscript[T, 2])+8 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+2 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-4 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 1,j]) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 3,i]))/((Subscript[T, 1]^2) (Subscript[T, 2]) (-1+2 (Subscript[T, 2]))))-(((6-4 (Subscript[T, 1])+(Subscript[T, 1]^2)+4 (Subscript[T, 2])-13 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 3,i]) (Subscript[x, 1,j]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))))-(((Subscript[p, 1,i]) (Subscript[p, 3,j]) (Subscript[x, 1,j]) (Subscript[x, 3,i]))/(Subscript[T, 2]))+(((-4+6 (Subscript[T, 1])-2 (Subscript[T, 1]^2)+8 (Subscript[T, 2])-11 (Subscript[T, 1]) (Subscript[T, 2])+4 (Subscript[T, 1]^2) (Subscript[T, 2])+4 (Subscript[T, 2]^2)-11 (Subscript[T, 1]) (Subscript[T, 2]^2)-4 (Subscript[T, 1]) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^3)) (Subscript[p, 1,j]) (Subscript[p, 3,j]) (Subscript[x, 1,j]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (Subscript[T, 2]) (-1+2 (Subscript[T, 2]))))+((-2+(Subscript[T, 1])-3 (Subscript[T, 2])+22 (Subscript[T, 1]) (Subscript[T, 2])-10 (Subscript[T, 1]^2) (Subscript[T, 2])+14 (Subscript[T, 2]^2)-42 (Subscript[T, 1]) (Subscript[T, 2]^2)-(Subscript[T, 1]^2) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-8 (Subscript[T, 2]^3)+10 (Subscript[T, 1]) (Subscript[T, 2]^3)+33 (Subscript[T, 1]^2) (Subscript[T, 2]^3)-18 (Subscript[T, 1]^3) (Subscript[T, 2]^3)+8 (Subscript[T, 1]) (Subscript[T, 2]^4)-22 (Subscript[T, 1]^2) (Subscript[T, 2]^4)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^4)) (Subscript[p, 2,j]) (Subscript[p, 3,i]) (Subscript[x, 2,i]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (Subscript[T, 2]) (-1+2 (Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2])))+((-2+(Subscript[T, 1])) ((Subscript[T, 1])+3 (Subscript[T, 2])-12 (Subscript[T, 1]) (Subscript[T, 2])-6 (Subscript[T, 2]^2)+17 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)+4 (Subscript[T, 2]^3)-3 (Subscript[T, 1]) (Subscript[T, 2]^3)-18 (Subscript[T, 1]^2) (Subscript[T, 2]^3)-4 (Subscript[T, 1]) (Subscript[T, 2]^4)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^4)) (Subscript[p, 2,i]) (Subscript[p, 3,j]) (Subscript[x, 2,i]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (-1+(Subscript[T, 2])) (Subscript[T, 2]) (-1+2 (Subscript[T, 2])))-(1/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (Subscript[T, 2]^2)))(1+(Subscript[T, 1])-(Subscript[T, 1]^2)+5 (Subscript[T, 2])-24 (Subscript[T, 1]) (Subscript[T, 2])+11 (Subscript[T, 1]^2) (Subscript[T, 2])-8 (Subscript[T, 2]^2)+18 (Subscript[T, 1]) (Subscript[T, 2]^2)+11 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)+8 (Subscript[T, 1]) (Subscript[T, 2]^3)-22 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 2,j]) (Subscript[p, 3,j]) (Subscript[x, 2,i]) (Subscript[x, 3,i])+((-3+2 (Subscript[T, 1])-(Subscript[T, 2])+21 (Subscript[T, 1]) (Subscript[T, 2])-11 (Subscript[T, 1]^2) (Subscript[T, 2])+14 (Subscript[T, 2]^2)-44 (Subscript[T, 1]) (Subscript[T, 2]^2)+(Subscript[T, 1]^2) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-8 (Subscript[T, 2]^3)+10 (Subscript[T, 1]) (Subscript[T, 2]^3)+33 (Subscript[T, 1]^2) (Subscript[T, 2]^3)-18 (Subscript[T, 1]^3) (Subscript[T, 2]^3)+8 (Subscript[T, 1]) (Subscript[T, 2]^4)-22 (Subscript[T, 1]^2) (Subscript[T, 2]^4)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^4)) (Subscript[p, 2,j]) (Subscript[p, 3,i]) (Subscript[x, 2,j]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (-1+(Subscript[T, 2])) (-1+2 (Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2])))-((((Subscript[T, 1])+4 (Subscript[T, 2])-11 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 2,i]) (Subscript[p, 3,j]) (Subscript[x, 2,j]) (Subscript[x, 3,i]))/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (-1+2 (Subscript[T, 2]))))-(1/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (Subscript[T, 2]) (-1+2 (Subscript[T, 2]))))(1+6 (Subscript[T, 2])-25 (Subscript[T, 1]) (Subscript[T, 2])+10 (Subscript[T, 1]^2) (Subscript[T, 2])-12 (Subscript[T, 2]^2)+29 (Subscript[T, 1]) (Subscript[T, 2]^2)+11 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)+12 (Subscript[T, 1]) (Subscript[T, 2]^3)-31 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 2,j]) (Subscript[p, 3,j]) (Subscript[x, 2,j]) (Subscript[x, 3,i])-((2-(Subscript[T, 1])+4 (Subscript[T, 2])-23 (Subscript[T, 1]) (Subscript[T, 2])+11 (Subscript[T, 1]^2) (Subscript[T, 2])-8 (Subscript[T, 2]^2)+20 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)+8 (Subscript[T, 1]) (Subscript[T, 2]^3)-22 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 3,i]) (Subscript[p, 3,j]) (Subscript[x, 3,i]^2))/((-1+(Subscript[T, 1])) (Subscript[T, 1]) (Subscript[T, 2]) (-1+2 (Subscript[T, 2])))+((-1+(Subscript[T, 1]) (Subscript[T, 2])) (1+6 (Subscript[T, 2])-28 (Subscript[T, 1]) (Subscript[T, 2])+13 (Subscript[T, 1]^2) (Subscript[T, 2])-8 (Subscript[T, 2]^2)+18 (Subscript[T, 1]) (Subscript[T, 2]^2)+27 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-18 (Subscript[T, 1]^3) (Subscript[T, 2]^2)+16 (Subscript[T, 1]) (Subscript[T, 2]^3)-44 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+18 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 3,j]^2) (Subscript[x, 3,i]^2))/(2 (-1+(Subscript[T, 1])) (Subscript[T, 1]^2) (Subscript[T, 2]^2) (-1+2 (Subscript[T, 2])))-(((-4+9 (Subscript[T, 1])) (Subscript[p, 3,j]) (Subscript[x, 3,j]))/(Subscript[T, 1]))+((-4+3 (Subscript[T, 1])-(Subscript[T, 1]^2)+12 (Subscript[T, 2])-17 (Subscript[T, 1]) (Subscript[T, 2])+12 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+2 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-9 (Subscript[T, 1]^3) (Subscript[T, 2]^2)-4 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+9 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 1,i]) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 3,j]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2])))-(((2+(Subscript[T, 1])-8 (Subscript[T, 2])+11 (Subscript[T, 1]) (Subscript[T, 2])-9 (Subscript[T, 1]^2) (Subscript[T, 2])-4 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 1,j]) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 3,j]))/((Subscript[T, 1]) (-1+2 (Subscript[T, 2]))))-(((-1-(Subscript[T, 2])+(Subscript[T, 1]) (Subscript[T, 2])) (Subscript[p, 2,i]) (Subscript[p, 3,j]) (Subscript[x, 2,i]) (Subscript[x, 3,j]))/((-1+(Subscript[T, 2])) (-1+(Subscript[T, 1]) (Subscript[T, 2]))))-(((-2+(Subscript[T, 1])) (1+2 (Subscript[T, 2])-9 (Subscript[T, 1]) (Subscript[T, 2])-4 (Subscript[T, 2]^2)+9 (Subscript[T, 1]) (Subscript[T, 2]^2)) (Subscript[p, 2,j]) (Subscript[p, 3,j]) (Subscript[x, 2,i]) (Subscript[x, 3,j]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))-(((-3+2 (Subscript[T, 1])-2 (Subscript[T, 2])+18 (Subscript[T, 1]) (Subscript[T, 2])-9 (Subscript[T, 1]^2) (Subscript[T, 2])+8 (Subscript[T, 2]^2)-22 (Subscript[T, 1]) (Subscript[T, 2]^2)+9 (Subscript[T, 1]^2) (Subscript[T, 2]^2)) (Subscript[p, 3,i]) (Subscript[p, 3,j]) (Subscript[x, 3,i]) (Subscript[x, 3,j]))/((-1+(Subscript[T, 1])) (-1+2 (Subscript[T, 2]))))+((1+6 (Subscript[T, 2])-31 (Subscript[T, 1]) (Subscript[T, 2])+15 (Subscript[T, 1]^2) (Subscript[T, 2])-8 (Subscript[T, 2]^2)+16 (Subscript[T, 1]) (Subscript[T, 2]^2)+45 (Subscript[T, 1]^2) (Subscript[T, 2]^2)-27 (Subscript[T, 1]^3) (Subscript[T, 2]^2)+24 (Subscript[T, 1]) (Subscript[T, 2]^3)-66 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+27 (Subscript[T, 1]^3) (Subscript[T, 2]^3)) (Subscript[p, 3,j]^2) (Subscript[x, 3,i]) (Subscript[x, 3,j]))/(2 (-1+(Subscript[T, 1])) (Subscript[T, 1]) (Subscript[T, 2]) (-1+2 (Subscript[T, 2])))

(Subscript[T, 3])=(Subscript[T, 1])(Subscript[T, 2]);
𝒮={(Subscript[x, __]),(Subscript[p, __])};
q[s_,i_,j_]:=Sum[(Subscript[x, α,i]) ((Subscript[p, α,i])-(Subscript[p, α,i+1]))+(Subscript[x, α,j]) ((Subscript[p, α,j])-(Subscript[p, α,j+1]))+(Subscript[x, α,i])((1-(Subscript[T, α]^s))(Subscript[p, α,i+1])+((Subscript[T, α]^s)-1)(Subscript[p, α,j+1])), {α,3}];
(Subscript[γ, 1])[φ_,k_]:= φ(3/2-(Subscript[x, 1,k])(Subscript[p, 1,k])-(Subscript[x, 2,k])(Subscript[p, 2,k])-(Subscript[x, 3,k])(Subscript[p, 3,k]));
ℒ[(Subscript[X, i_,j_])[s_]]:=(Subscript[T, 3]^s)𝔼[-q[s,i,j]+ (Subscript[r, 0])[s,i,j]+ϵ (Subscript[r, 1])[s,i,j]+((O[ϵ])^2)];
ℒ[(Subscript[C, k_])[φ_]]:=(Subscript[T, 3]^φ)𝔼[-(Subscript[x, 1,k])((Subscript[p, 1,k])-(Subscript[p, 1,k+1]))-(Subscript[x, 2,k])((Subscript[p, 2,k])-(Subscript[p, 2,k+1]))-(Subscript[x, 3,k])((Subscript[p, 3,k])-(Subscript[p, 3,k+1]))+ ϵ  (Subscript[γ, 1])[φ,k]+((O[ϵ])^2)];
ℒ[K_]:=(((2π))^(-Features[K][[1]]))CF[ℒ/@Features[K][[2]]];
(Subscript[vs, i_]):=Sequence[(Subscript[p, 1,i]),(Subscript[x, 1,i]),(Subscript[p, 2,i]),(Subscript[x, 2,i]),(Subscript[p, 3,i]),(Subscript[x, 3,i])];
vs[K_] := Union@@Table[{(Subscript[vs, i])},{i,Features[K][[1]]}]

Features[Knot[3,1]]

FKnot31nullm=Features[16,(Subscript[X, 16,14])[1](Subscript[X, 13,11])[1](Subscript[X, 10,8])[1](Subscript[C, 15])[1](Subscript[C, 12])[1](Subscript[C, 9])[1](Subscript[C, 4])[-1] (Subscript[X, 2,6])[-1] (Subscript[X, 5,1])[-1] (Subscript[X, 7,3])[-1]]

Short[ℒ/@FKnot31nullm[[2]],10]

Short[ℒ[Knot[3,1]],10]

vs[Knot[3,1]]

FKnot31nullm[[1]]

∫CF[ℒ/@FKnot31nullm[[2]]]\[DifferentialD] Union@@Table[{(Subscript[vs, i])},{i,FKnot31nullm[[1]]}]

∫ℒ[Knot[3,1]]\[DifferentialD] vs[Knot[3,1]]

-((16384 I (π^14) (Subscript[T, 1]^2) (Subscript[T, 2]^2) 𝔼[ϵSeries[0,(-7+21 (Subscript[T, 1])-29 (Subscript[T, 1]^2)+23 (Subscript[T, 1]^3)-9 (Subscript[T, 1]^4)+35 (Subscript[T, 2])-90 (Subscript[T, 1]) (Subscript[T, 2])+107 (Subscript[T, 1]^2) (Subscript[T, 2])-66 (Subscript[T, 1]^3) (Subscript[T, 2])+9 (Subscript[T, 1]^4) (Subscript[T, 2])+12 (Subscript[T, 1]^5) (Subscript[T, 2])-79 (Subscript[T, 2]^2)+157 (Subscript[T, 1]) (Subscript[T, 2]^2)-139 (Subscript[T, 1]^2) (Subscript[T, 2]^2)+30 (Subscript[T, 1]^3) (Subscript[T, 2]^2)+52 (Subscript[T, 1]^4) (Subscript[T, 2]^2)-38 (Subscript[T, 1]^5) (Subscript[T, 2]^2)-7 (Subscript[T, 1]^6) (Subscript[T, 2]^2)+109 (Subscript[T, 2]^3)-142 (Subscript[T, 1]) (Subscript[T, 2]^3)+56 (Subscript[T, 1]^2) (Subscript[T, 2]^3)+75 (Subscript[T, 1]^3) (Subscript[T, 2]^3)-73 (Subscript[T, 1]^4) (Subscript[T, 2]^3)-28 (Subscript[T, 1]^5) (Subscript[T, 2]^3)+70 (Subscript[T, 1]^6) (Subscript[T, 2]^3)-9 (Subscript[T, 1]^7) (Subscript[T, 2]^3)-91 (Subscript[T, 2]^4)+19 (Subscript[T, 1]) (Subscript[T, 2]^4)+91 (Subscript[T, 1]^2) (Subscript[T, 2]^4)-108 (Subscript[T, 1]^3) (Subscript[T, 2]^4)-83 (Subscript[T, 1]^4) (Subscript[T, 2]^4)+232 (Subscript[T, 1]^5) (Subscript[T, 2]^4)-167 (Subscript[T, 1]^6) (Subscript[T, 2]^4)+12 (Subscript[T, 1]^7) (Subscript[T, 2]^4)+3 (Subscript[T, 1]^8) (Subscript[T, 2]^4)+46 (Subscript[T, 2]^5)+78 (Subscript[T, 1]) (Subscript[T, 2]^5)-97 (Subscript[T, 1]^2) (Subscript[T, 2]^5)-58 (Subscript[T, 1]^3) (Subscript[T, 2]^5)+334 (Subscript[T, 1]^4) (Subscript[T, 2]^5)-376 (Subscript[T, 1]^5) (Subscript[T, 2]^5)+160 (Subscript[T, 1]^6) (Subscript[T, 2]^5)+22 (Subscript[T, 1]^7) (Subscript[T, 2]^5)-9 (Subscript[T, 1]^8) (Subscript[T, 2]^5)-8 (Subscript[T, 2]^6)-84 (Subscript[T, 1]) (Subscript[T, 2]^6)+4 (Subscript[T, 1]^2) (Subscript[T, 2]^6)+196 (Subscript[T, 1]^3) (Subscript[T, 2]^6)-370 (Subscript[T, 1]^4) (Subscript[T, 2]^6)+270 (Subscript[T, 1]^5) (Subscript[T, 2]^6)-36 (Subscript[T, 1]^6) (Subscript[T, 2]^6)-58 (Subscript[T, 1]^7) (Subscript[T, 2]^6)+7 (Subscript[T, 1]^8) (Subscript[T, 2]^6)+24 (Subscript[T, 1]) (Subscript[T, 2]^7)+72 (Subscript[T, 1]^2) (Subscript[T, 2]^7)-170 (Subscript[T, 1]^3) (Subscript[T, 2]^7)+172 (Subscript[T, 1]^4) (Subscript[T, 2]^7)-30 (Subscript[T, 1]^5) (Subscript[T, 2]^7)-85 (Subscript[T, 1]^6) (Subscript[T, 2]^7)+61 (Subscript[T, 1]^7) (Subscript[T, 2]^7)+(Subscript[T, 1]^8) (Subscript[T, 2]^7)-32 (Subscript[T, 1]^2) (Subscript[T, 2]^8)+22 (Subscript[T, 1]^3) (Subscript[T, 2]^8)+30 (Subscript[T, 1]^4) (Subscript[T, 2]^8)-88 (Subscript[T, 1]^5) (Subscript[T, 2]^8)+78 (Subscript[T, 1]^6) (Subscript[T, 2]^8)-21 (Subscript[T, 1]^7) (Subscript[T, 2]^8)-9 (Subscript[T, 1]^8) (Subscript[T, 2]^8)+16 (Subscript[T, 1]^3) (Subscript[T, 2]^9)-30 (Subscript[T, 1]^4) (Subscript[T, 2]^9)+30 (Subscript[T, 1]^5) (Subscript[T, 2]^9)-10 (Subscript[T, 1]^6) (Subscript[T, 2]^9)-6 (Subscript[T, 1]^7) (Subscript[T, 2]^9)+6 (Subscript[T, 1]^8) (Subscript[T, 2]^9))/((((1-(Subscript[T, 1])+(Subscript[T, 1]^2)))^2) (-1+2 (Subscript[T, 2])) (((1-(Subscript[T, 2])+(Subscript[T, 2]^2)))^2) (((1-(Subscript[T, 1]) (Subscript[T, 2])+(Subscript[T, 1]^2) (Subscript[T, 2]^2)))^2))]])/((1-(Subscript[T, 1])+(Subscript[T, 1]^2)) (1-(Subscript[T, 2])+(Subscript[T, 2]^2)) (1-(Subscript[T, 1]) (Subscript[T, 2])+(Subscript[T, 1]^2) (Subscript[T, 2]^2))))

Timing [∫ℒ[Knot[4,1]]\[DifferentialD] vs[Knot[4,1]]]

Features[Knot[3,1]]

And now the mutants!

KTKnot=Knot[11,NonAlternating,42];
ConwayKnot=Knot[11,NonAlternating,34];
Alexander/@{KTKnot,ConwayKnot}

A2Conway=Timing [∫ℒ[ConwayKnot]\[DifferentialD] vs[ConwayKnot]]

A2KT=Timing [∫ℒ[KTKnot]\[DifferentialD] vs[KTKnot]]

Reidemeister 2b

Territory: (Subscript[(x|p), (SuperPlus[(i|j)])]), (Subscript[x, i|j]).S

{lhs2b}=Cases[∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i])+(Subscript[π, α,j])(Subscript[p, α,j]),{α,3}]]ℒ/@((Subscript[X, i,j])[1](Subscript[X, i+1,j+1])[-1])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, j]),(Subscript[vs, i+1]),(Subscript[vs, j+1])}, ϵSeries[ℰ_,ℱ_]:>ℱ,∞]

{rhs2b}=Cases[∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i])+(Subscript[π, α,j])(Subscript[p, α,j]),{α,3}]]ℒ/@((Subscript[C, i])[0](Subscript[C, i+1])[0](Subscript[C, j])[0](Subscript[C, j+1])[0])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, j]),(Subscript[vs, i+1]),(Subscript[vs, j+1])}, ϵSeries[ℰ_,ℱ_]:>ℱ,∞]

eqn2b=CF[CF[CF[lhs2b-rhs2b]/.{$A $C->(Subscript[α, 1]),$A $D->(Subscript[α, 2]),$B $C->(Subscript[α, 3]),$B $D->(Subscript[α, 4])}]/.{$A|$B->0}]

cvs2b=Union@Cases[eqn2b,(Subscript[p, __])|(Subscript[π, __]),∞];

eqns2b=CoefficientRules[eqn2b,cvs2b]/.(_->c_):>(c==0);

vars2b=Union@Cases[eqn2b,(Subscript[d, _]),∞];

eqns2b//Column

{sol2b}=Solve[eqns2b,vars2b]

Cases[sol2b,(Subscript[α, _]),∞]

eqn2b/.sol2b//CF

Verification of Invariance Under Reidemeister 3b

{lhs3b}=Cases[∫𝔼[Sum[(Subscript[π, l,i]) (Subscript[p, l,i])+(Subscript[π, l,j]) (Subscript[p, l,j])+(Subscript[π, l,k]) (Subscript[p, l,k]),{l,3}]]ℒ/@((Subscript[X, i,j])[1](Subscript[X, i+1,k])[1](Subscript[X, j+1,k+1])[1])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, j]),(Subscript[vs, k]),(Subscript[vs, i+1]),(Subscript[vs, j+1]),(Subscript[vs, k+1])}, ϵSeries[ℰ_,ℱ_]:>ℱ,∞]

{rhs3b}=Cases[∫𝔼[Sum[(Subscript[π, l,i]) (Subscript[p, l,i])+(Subscript[π, l,j]) (Subscript[p, l,j])+(Subscript[π, l,k]) (Subscript[p, l,k]),{l,3}]]ℒ/@((Subscript[X, j,k])[1](Subscript[X, i,k+1])[1](Subscript[X, i+1,j+1])[1])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, j]),(Subscript[vs, k]),(Subscript[vs, i+1]),(Subscript[vs, j+1]),(Subscript[vs, k+1])}, ϵSeries[ℰ_,ℱ_]:>ℱ,∞]

eqn3b=CF[CF[CF[lhs3b-rhs3b]/.{$A $C->(Subscript[α, 1]),$A $D->(Subscript[α, 2]),$B $C->(Subscript[α, 3]),$B $D->(Subscript[α, 4])}]/.{$A|$B->0}]

eqn3b/.sol3b/.sol3bα/.sα3//CF

Invariance Under R2c

lhs2c=∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i])+(Subscript[π, α,j])(Subscript[p, α,j]),{α,3}]]ℒ/@((Subscript[X, i+1,j])[1](Subscript[X, i,j+2])[-1](Subscript[C, j+1])[1])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, j]),(Subscript[vs, i+1]),(Subscript[vs, j+1]),(Subscript[vs, j+2])}

rhs2c=∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i])+(Subscript[π, α,j])(Subscript[p, α,j]),{α,3}]]ℒ/@((Subscript[C, i])[0](Subscript[C, i+1])[0](Subscript[C, j])[0](Subscript[C, j+1])[1](Subscript[C, j+2])[0])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, j]),(Subscript[vs, i+1]),(Subscript[vs, j+1]),(Subscript[vs, j+2])}

Cases[Expand[rhs2c],$A $C,∞]

eqn2c=CF[CF[CF[Cases[lhs2c,ϵSeries[ℰ_,ℱ_]:>ℱ,∞]-Cases[rhs2c,ϵSeries[ℰ_,ℱ_]:>ℱ,∞]]/.{$A $C->(Subscript[α, 1]),$A $D->(Subscript[α, 2]),$B $C->(Subscript[α, 3]),$B $D->(Subscript[α, 4])}]/.{$A|$B->0}]

Cases[eqn2c,(Subscript[α, _]),∞]

eqn2c/.sol2b//CF

eqn2cred=eqn2c/.sol2b/.sol3b/.sol3bα/.sα3//CF

cvs2c=Union@Cases[eqn2cred,(Subscript[p, __])|(Subscript[π, __]),∞]

{eqns2c}=CoefficientRules[eqn2cred,cvs2c]/.(_->c_):>(c==0)

vars2b=Union@Cases[eqn2cred,(Subscript[c, _]),∞]

{sol2c}=Solve[eqns2c,{(Subscript[c, 5]),(Subscript[c, 50]),(Subscript[c, 79])}]//Simplify

eqn2c/.sol2b/.sol3b/.sol3bα/.sα3/.sol2c//CF

Invariance Under R1l

{lhs1l}=Cases[∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i]),{α,3}]]ℒ/@((Subscript[X, i+2,i])[1](Subscript[C, i+1])[1])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, i+1]),(Subscript[vs, i+2])},ϵSeries[ℰ_,ℱ_]:>ℱ,∞]

rhs1l=∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i]),{α,3}]]ℒ/@((Subscript[C, i])[0](Subscript[C, i+1])[0](Subscript[C, i+2])[0])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, i+1]),(Subscript[vs, i+2])}

eqn1l=CF[CF[CF[lhs1l]/.{$A $C->(Subscript[α, 1]),$A $D->(Subscript[α, 2]),$B $C->(Subscript[α, 3]),$B $D->(Subscript[α, 4])}]/.{$A|$B->0}];

eqn1lred=eqn1l/.sol2b/.sol3b/.sol3bα/.sα3/.sol2c//CF

{sol1l}=Solve[eqn1lred==0,(Subscript[c, 88])]

Invariance Under R1r

{lhs1r}=Cases[∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i]),{α,3}]]ℒ/@((Subscript[X, i,i+2])[1](Subscript[C, i+1])[-1])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, i+1]),(Subscript[vs, i+2])},ϵSeries[ℰ_,ℱ_]:>ℱ,∞]

rhs1r=∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i]),{α,3}]]ℒ/@((Subscript[C, i])[0](Subscript[C, i+1])[0](Subscript[C, i+2])[0])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, i+1]),(Subscript[vs, i+2])}

eqn1r=CF[CF[CF[lhs1r]/.{$A $C->(Subscript[α, 1]),$A $D->(Subscript[α, 2]),$B $C->(Subscript[α, 3]),$B $D->(Subscript[α, 4])}]/.{$A|$B->0}];

eqn1rred=eqn1r/.sol2b/.sol3b/.sol3bα/.sα3/.sol2c/.sol1l//CF

{sol1r}=Solve[eqn1rred==0,{(Subscript[c, 9])}]//Simplify

eqn1r/.sol2b/.sol3b/.sol3bα/.sα3/.sol2c/.sol1l/.sol1r//CF

Invariance Under Sw

lhssw=∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i])+(Subscript[π, α,j])(Subscript[p, α,j]),{α,3}]]ℒ/@((Subscript[X, i+1,j+1])[1](Subscript[C, i])[-1](Subscript[C, j])[-1](Subscript[C, i+2])[1](Subscript[C, j+2])[1])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, j]),(Subscript[vs, i+1]),(Subscript[vs, j+1]),(Subscript[vs, i+2]),(Subscript[vs, j+2])}

rhssw=∫𝔼[Sum[(Subscript[π, α,i])(Subscript[p, α,i])+(Subscript[π, α,j])(Subscript[p, α,j]),{α,3}]]ℒ/@((Subscript[X, i+1,j+1])[1](Subscript[C, i])[0](Subscript[C, j])[0](Subscript[C, i+2])[0](Subscript[C, j+2])[0])\[DifferentialD]{(Subscript[vs, i]),(Subscript[vs, j]),(Subscript[vs, i+1]),(Subscript[vs, j+1]),(Subscript[vs, i+2]),(Subscript[vs, j+2])}

eqnsw=CF[CF[CF[Cases[lhssw,ϵSeries[ℰ_,ℱ_]:>ℱ,∞]-Cases[rhssw,ϵSeries[ℰ_,ℱ_]:>ℱ,∞]]]]

(*single $A and $B never arise so sw equation holds.*)

eqnsw=CF[CF[CF[Cases[lhssw,ϵSeries[ℰ_,ℱ_]:>ℱ,∞]-Cases[rhssw,ϵSeries[ℰ_,ℱ_]:>ℱ,∞]]/.{$A $C->(Subscript[α, 1]),$A $D->(Subscript[α, 2]),$B $C->(Subscript[α, 3]),$B $D->(Subscript[α, 4])}]/.{$A|$B->0}]

In conclusion: here are the values for the R-matrix we found:

{r0p,r0m}={(Subscript[r, 0])[1,i,j],(Subscript[r, 0])[-1,i,j]}/.sol2b/.sol3b/.sol3bα/.sα3/.sol2c/.sol1l/.sol1r//CF

r1p=CF[CF[(Subscript[r, 1])[1,i,j]/.sol3b/.sol3bα]/.sα3/.sol2c/.sol1l]/.sol1r//CF

r1m=CF[CF[CF[(Subscript[r, 1])[-1,i,j]/.sol2b]/.sol3b/.sol3bα]/.sα3/.sol2c/.sol1l]/.sol1r//CF

Union@Cases[r1m,(Subscript[c, _]),∞]

RandomChoice[{-2,-1,1,2}]

Choices={-2,-1,1,2};
Sub={(Subscript[c, 4])->RandomChoice[Choices],(Subscript[c, 10])->RandomChoice[Choices],(Subscript[c, 13])->RandomChoice[Choices],(Subscript[c, 15])->RandomChoice[Choices],(Subscript[c, 22])->RandomChoice[Choices],(Subscript[c, 28])->RandomChoice[{-2,-1,1,2}],(Subscript[c, 32])->RandomChoice[Choices],(Subscript[c, 36])->RandomChoice[Choices],(Subscript[c, 49])->RandomChoice[Choices],(Subscript[c, 54])->RandomChoice[Choices],(Subscript[c, 56])->RandomChoice[Choices],(Subscript[c, 65])->RandomChoice[Choices],(Subscript[c, 69])->RandomChoice[Choices],(Subscript[c, 78])->RandomChoice[Choices],(Subscript[c, 81])->RandomChoice[Choices],(Subscript[c, 86])->RandomChoice[Choices],(Subscript[α, 4])->RandomChoice[Choices],$B->RandomChoice[Choices]}

r0p/.Sub//CF
r0m/.Sub//CF

r1p/.Sub//CF

r1m/.Sub//CF