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])}];

(*mons0=MonomialList[
(Subscript[p, 3])(Subscript[x, 1])(Subscript[x, 2])/.
{(Subscript[(v:p|x), α_]):>(Subscript[v, α,i])+(Subscript[v, α,j])}
]/.c_Integer*mon_:>mon;*)
mons1=MonomialList[
1+(Subscript[p, 1])(Subscript[x, 1])+(Subscript[p, 2])(Subscript[x, 2])+(Subscript[p, 3])(Subscript[x, 3])+(Subscript[p, 1])(Subscript[p, 1])(Subscript[x, 1])(Subscript[x, 1])+(Subscript[p, 2])(Subscript[p, 2])(Subscript[x, 2])(Subscript[x, 2])+(Subscript[p, 1])(Subscript[p, 2])(Subscript[x, 1])(Subscript[x, 2])+(Subscript[p, 1])(Subscript[p, 3])(Subscript[x, 1])(Subscript[x, 3])+(Subscript[p, 2])(Subscript[p, 3])(Subscript[x, 2])(Subscript[x, 3])+(Subscript[p, 3])(Subscript[p, 3])(Subscript[x, 3])(Subscript[x, 3])+0(Subscript[p, 1])(Subscript[p, 2])(Subscript[x, 3])/.
{(Subscript[(v:p|x), α_]):>(Subscript[v, α,i])+(Subscript[v, α,j])}
]/.c_Integer*mon_:>mon;
(*$A,$B,$C,$D part from DeterminingThePXXandXPPCoefficients.nb file,
sol3b and sol3bα sα3 from UC4A2R42.nb*)
(*Subs={(Subscript[c, _])->1,BB->BB,(Subscript[α, 4])->BB};*)(Subscript[r, 0])[1,i_,j_]:=($A (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 2,i])-((($A+$B (Subscript[T, 2])) (Subscript[p, 3,j]) (Subscript[x, 1,j]) (Subscript[x, 2,i]))/(Subscript[T, 1]))+$B (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 2,j]))
(Subscript[r, 0])[-1,i_,j_]:=((((-$A+$B (Subscript[T, 1])-$B (Subscript[T, 2])) (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 2,i]))/((Subscript[T, 1]^2) (Subscript[T, 2])))+((($A+$B (Subscript[T, 2])) (Subscript[p, 3,j]) (Subscript[x, 1,j]) (Subscript[x, 2,i]))/((Subscript[T, 1]) (Subscript[T, 2])))-(($B (Subscript[p, 3,j]) (Subscript[x, 1,i]) (Subscript[x, 2,j]))/(Subscript[T, 1])))
(Subscript[r, 1])[1,i_,j_]:=(Evaluate[$C (Subscript[p, 1,j]) (Subscript[p, 2,i]) (Subscript[x, 3,i])+$D (Subscript[p, 1,i]) (Subscript[p, 2,j]) (Subscript[x, 3,i])+(-$C-$D) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 3,i])+ Sum[(Subscript[c, z])mons1[[z]],{z,1,Length@mons1}]])
(Subscript[r, 1])[-1,i_,j_]:=
(Evaluate[-(($C (Subscript[p, 1,j]) (Subscript[p, 2,i]) (Subscript[x, 3,i]))/(Subscript[T, 1]))-(($D (Subscript[p, 1,i]) (Subscript[p, 2,j]) (Subscript[x, 3,i]))/(Subscript[T, 2]))+((($D (Subscript[T, 1])+$C (Subscript[T, 2])) (Subscript[p, 1,j]) (Subscript[p, 2,j]) (Subscript[x, 3,i]))/((Subscript[T, 1]) (Subscript[T, 2])))+ Sum[(Subscript[d, z])mons1[[z]],{z,1,Length@mons1}]])
(Subscript[r, 0])[1,4,7]
(Subscript[r, 1])[1,4,7];

(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]]}]

Reidemeister 2b

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

{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