v1r23a = {p1r23*u1b1r23, p1r23*u2b1r23, p1r23*u3b1r23}
 
v1r23b = {p1r23*(1 - u2b1r23 - u3b1r23), p1r23*u2b1r23, p1r23*u3b1r23}

v1r23c = {(a23*a3f2*(1 - af23/2)*(2 + n1 + n2 + 2*n1*n2)*
       (1 - u2b1r23 - u3b1r23))/(6*(1 - (1 - a23)*(1 - af23)*(1 - e4))*
       (1 - (1 - e4)*n1*n2*n3)), (a23*a3f2*(1 - af23/2)*
       (2 + n1 + n2 + 2*n1*n2)*u2b1r23)/
      (6*(1 - (1 - a23)*(1 - af23)*(1 - e4))*(1 - (1 - e4)*n1*n2*n3)), 
     (a23*a3f2*(1 - af23/2)*(2 + n1 + n2 + 2*n1*n2)*u3b1r23)/
      (6*(1 - (1 - a23)*(1 - af23)*(1 - e4))*(1 - (1 - e4)*n1*n2*n3))}

v1r23d = {(a23*(4 - a1f2 - a1f3 - a2f1 + 2*(1 - a1f2 - a1f3)*
         (1 - a2f1 - a2f3) - a2f3)*a3f2*(1 - af23/2)*(1 - u2b1r23 - u3b1r23))/
      (6*(1 - (1 - a1f2 - a1f3)*(1 - a2f1 - a2f3)*(1 - a3f1 - a3f2)*(1 - e4))*
       (1 - (1 - a23)*(1 - af23)*(1 - e4))), 
     (a23*(4 - a1f2 - a1f3 - a2f1 + 2*(1 - a1f2 - a1f3)*(1 - a2f1 - a2f3) - 
        a2f3)*a3f2*(1 - af23/2)*u2b1r23)/
      (6*(1 - (1 - a1f2 - a1f3)*(1 - a2f1 - a2f3)*(1 - a3f1 - a3f2)*(1 - e4))*
       (1 - (1 - a23)*(1 - af23)*(1 - e4))), 
     (a23*(4 - a1f2 - a1f3 - a2f1 + 2*(1 - a1f2 - a1f3)*(1 - a2f1 - a2f3) - 
        a2f3)*a3f2*(1 - af23/2)*u3b1r23)/
      (6*(1 - (1 - a1f2 - a1f3)*(1 - a2f1 - a2f3)*(1 - a3f1 - a3f2)*(1 - e4))*
       (1 - (1 - a23)*(1 - af23)*(1 - e4)))}
    
v12r3a = {p12r3*u1b12r3, p12r3*u2b12r3, p12r3*u3b12r3}
 
v12r3b = {p12r3*(1 - u2b12r3 - u3b12r3), p12r3*u2b12r3, p12r3*u3b12r3}
  
v12r3c = {((1 - a12/2)*a2f1*af12*(2 + n1 + n3 + 2*n1*n3)*
       (1 - u2b12r3 - u3b12r3))/(6*(1 - (1 - a21)*(1 - af21)*(1 - e4))*
       (1 - (1 - e4)*n1*n2*n3)), ((1 - a12/2)*a2f1*af12*
       (2 + n1 + n3 + 2*n1*n3)*u2b12r3)/
      (6*(1 - (1 - a21)*(1 - af21)*(1 - e4))*(1 - (1 - e4)*n1*n2*n3)), 
     ((1 - a12/2)*a2f1*af12*(2 + n1 + n3 + 2*n1*n3)*u3b12r3)/
      (6*(1 - (1 - a21)*(1 - af21)*(1 - e4))*(1 - (1 - e4)*n1*n2*n3))}
  
v12r3d = {((1 - a12/2)*a2f1*(4 - a1f2 - a1f3 - a3f1 + 2*(1 - a1f2 - a1f3)*
         (1 - a3f1 - a3f2) - a3f2)*af12*(1 - u2b12r3 - u3b12r3))/
      (6*(1 - (1 - a1f2 - a1f3)*(1 - a2f1 - a2f3)*(1 - a3f1 - a3f2)*(1 - e4))*
       (1 - (1 - a21)*(1 - af21)*(1 - e4))), 
     ((1 - a12/2)*a2f1*(4 - a1f2 - a1f3 - a3f1 + 2*(1 - a1f2 - a1f3)*
         (1 - a3f1 - a3f2) - a3f2)*af12*u2b12r3)/
      (6*(1 - (1 - a1f2 - a1f3)*(1 - a2f1 - a2f3)*(1 - a3f1 - a3f2)*(1 - e4))*
       (1 - (1 - a21)*(1 - af21)*(1 - e4))), 
     ((1 - a12/2)*a2f1*(4 - a1f2 - a1f3 - a3f1 + 2*(1 - a1f2 - a1f3)*
         (1 - a3f1 - a3f2) - a3f2)*af12*u3b12r3)/
      (6*(1 - (1 - a1f2 - a1f3)*(1 - a2f1 - a2f3)*(1 - a3f1 - a3f2)*(1 - e4))*
       (1 - (1 - a21)*(1 - af21)*(1 - e4)))}