R_F(x,y,z):=

R_C(x,y):=

R_J(x,y,z,p):=

R_D(x,y,z):=

R_G(x,y,z):=

e1:=

DF(i,j,a,b):=a SUB i*b SUB j-a SUB j*b SUB i

XF(i,a,b,x):=SQRT(a SUB i+b SUB i*x)

U12A(a,b,x,y):=(SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 2+y*b SUB 2)*SQRT(a SUB 3+~
x*b SUB 3)*SQRT(a SUB 4+x*b SUB 4)+SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 2+x*b S~
UB 2)*SQRT(a SUB 3+y*b SUB 3)*SQRT(a SUB 4+y*b SUB 4))/(x-y)

U13A(a,b,x,y):=(SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 2+x*b SUB 2)*SQRT(a SUB 3+~
y*b SUB 3)*SQRT(a SUB 4+x*b SUB 4)+SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 2+y*b S~
UB 2)*SQRT(a SUB 3+x*b SUB 3)*SQRT(a SUB 4+y*b SUB 4))/(x-y)

U23A(a,b,x,y):=(SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 2+y*b SUB 2)*SQRT(a SUB 3+~
y*b SUB 3)*SQRT(a SUB 4+x*b SUB 4)+SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 2+x*b S~
UB 2)*SQRT(a SUB 3+x*b SUB 3)*SQRT(a SUB 4+y*b SUB 4))/(x-y)

I414_2(v):=[(v SUB 1+v SUB 2)*(v SUB 1+v SUB 3),(v SUB 2+v SUB 1)*(v SUB 2+v ~
SUB 3),(v SUB 3+v SUB 1)*(v SUB 3+v SUB 2)]

I414_1(v):=4*R_F(v SUB 1,v SUB 2,v SUB 3)

I414T(a,b,x,y):=I414_1(I414_2([U12A(a,b,x,y),U13A(a,b,x,y),U23A(a,b,x,y)]))

I414(a,b,x,y):=I414_1(I414_2([u12,u13,u23]))

U15A(a,b,x,y):=SQRT(U12A(a,b,x,y)^2-(a SUB 3*b SUB 1-a SUB 1*b SUB 3)*(a SUB ~
4*b SUB 1-a SUB 1*b SUB 4)*(a SUB 5*b SUB 2-a SUB 2*b SUB 5)/(a SUB 5*b SUB 1~
-a SUB 1*b SUB 5))

U15B(a,b,x,y):=SQRT(u12^2-(a SUB 3*b SUB 1-a SUB 1*b SUB 3)*(a SUB 4*b SUB 1-~
a SUB 1*b SUB 4)*(a SUB 5*b SUB 2-a SUB 2*b SUB 5)/(a SUB 5*b SUB 1-a SUB 1*b~
 SUB 5))

Q15A(a,b,x,y):=SQRT(a SUB 5+x*b SUB 5)*SQRT(a SUB 5+y*b SUB 5)*U15A(a,b,x,y)/~
(SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 1+y*b SUB 1))

Q15B(a,b,x,y):=SQRT(a SUB 5+x*b SUB 5)*SQRT(a SUB 5+y*b SUB 5)*u15/(SQRT(a SU~
B 1+x*b SUB 1)*SQRT(a SUB 1+y*b SUB 1))

S15A(a,b,x,y):=(SQRT(a SUB 2+x*b SUB 2)*SQRT(a SUB 3+x*b SUB 3)*SQRT(a SUB 4+~
x*b SUB 4)*(a SUB 5+y*b SUB 5)/SQRT(a SUB 1+x*b SUB 1)+SQRT(a SUB 2+y*b SUB 2~
)*SQRT(a SUB 3+y*b SUB 3)*SQRT(a SUB 4+y*b SUB 4)*(a SUB 5+x*b SUB 5)/SQRT(a ~
SUB 1+y*b SUB 1))/(x-y)

I417_P(b,c,t):=(b*c+t)/(b+c)

B15A(a,b,u,v,w,z):=(a SUB 5*b SUB 1-a SUB 1*b SUB 5)/((a SUB 2*b SUB 1-a SUB ~
1*b SUB 2)*(a SUB 3*b SUB 1-a SUB 1*b SUB 3)*(a SUB 4*b SUB 1-a SUB 1*b SUB 4~
))*(z+u)*(z+v)*(z+w)

B15B(a,b):=B15A(a,b,u12,u13,u23,u15)

I417_2(a,b,p):=4/p*R_C(1,1-(a SUB 5*b SUB 2-a SUB 2*b SUB 5)*(a SUB 5*b SUB 3~
-a SUB 3*b SUB 5)*(a SUB 5*b SUB 4-a SUB 4*b SUB 5)/(p^2*(a SUB 5*b SUB 1-a S~
UB 1*b SUB 5)))

I417_1A(a,b,x,y):=I417_2(a,b,I417_P(B15A(a,b,U12A(a,b,x,y),U13A(a,b,x,y),U23A~
(a,b,x,y),U15A(a,b,x,y)),S15A(a,b,x,y)+Q15A(a,b,x,y),(a SUB 5*b SUB 2-a SUB 2~
*b SUB 5)*(a SUB 5*b SUB 3-a SUB 3*b SUB 5)*(a SUB 5*b SUB 4-a SUB 4*b SUB 5)~
/(a SUB 5*b SUB 1-a SUB 1*b SUB 5)))

I417_1(a,b,x,y):=I417_2(a,b,I417_P(b15,s15+q15,(a SUB 5*b SUB 2-a SUB 2*b SUB~
 5)*(a SUB 5*b SUB 3-a SUB 3*b SUB 5)*(a SUB 5*b SUB 4-a SUB 4*b SUB 5)/(a SU~
B 5*b SUB 1-a SUB 1*b SUB 5)))

I417_4(v):=[(v SUB 1+v SUB 2)*(v SUB 1+v SUB 3),(v SUB 2+v SUB 1)*(v SUB 2+v ~
SUB 3),(v SUB 3+v SUB 1)*(v SUB 3+v SUB 2),v SUB 4^2+v SUB 1*v SUB 2+v SUB 1*~
v SUB 3+v SUB 2*v SUB 3]

I417_3(v):=4/3*R_J(v SUB 1,v SUB 2,v SUB 3,v SUB 4)

I417T(a,b,x,y):=(a SUB 2*b SUB 1-a SUB 1*b SUB 2)*(a SUB 3*b SUB 1-a SUB 1*b ~
SUB 3)*(a SUB 4*b SUB 1-a SUB 1*b SUB 4)/(a SUB 5*b SUB 1-a SUB 1*b SUB 5)*I4~
17_3(I417_4([U12A(a,b,x,y),U13A(a,b,x,y),U23A(a,b,x,y),U15A(a,b,x,y)]))+I417_~
1A(a,b,x,y)

I417(a,b,x,y):=(a SUB 2*b SUB 1-a SUB 1*b SUB 2)*(a SUB 3*b SUB 1-a SUB 1*b S~
UB 3)*(a SUB 4*b SUB 1-a SUB 1*b SUB 4)/(a SUB 5*b SUB 1-a SUB 1*b SUB 5)*I41~
7_3(I417_4([u12,u13,u23,u15]))+I417_1(a,b,x,y)

I427T(a,b,x,y):=(b SUB 5*I417T(a,b,x,y)-b SUB 1*I414T(a,b,x,y))/(a SUB 1*b SU~
B 5-a SUB 5*b SUB 1)

I427(a,b,x,y):=(b SUB 5*I417(a,b,x,y)-b SUB 1*I414(a,b,x,y))/(a SUB 1*b SUB 5~
-a SUB 5*b SUB 1)

I428_1(v):=2*R_D(v SUB 2,v SUB 3,v SUB 1)

I428T(a,b,x,y):=(b SUB 1*2*((a SUB 2*b SUB 3-a SUB 3*b SUB 2)*(a SUB 2*b SUB ~
4-a SUB 4*b SUB 2)*I428_1(I414_2([U12A(a,b,x,y),U13A(a,b,x,y),U23A(a,b,x,y)])~
)/3+(x-y)*(SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 2+x*b SUB 2)+SQRT(a SUB 1+x*b S~
UB 1)*SQRT(a SUB 2+y*b SUB 2))/(SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 1+y*b SUB ~
1)*(SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 3+x*b SUB 3)+SQRT(a SUB 1+x*b SUB 1)*S~
QRT(a SUB 3+y*b SUB 3))*(SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 4+x*b SUB 4)+SQRT~
(a SUB 1+x*b SUB 1)*SQRT(a SUB 4+y*b SUB 4))))-b SUB 2*I414T(a,b,x,y))/(a SUB~
 2*b SUB 1-a SUB 1*b SUB 2)

I428(a,b,x,y):=(b SUB 1*2*((a SUB 2*b SUB 3-a SUB 3*b SUB 2)*(a SUB 2*b SUB 4~
-a SUB 4*b SUB 2)*I428_1(I414_2([u12,u13,u23]))/3+(x-y)*(SQRT(a SUB 1+y*b SUB~
 1)*SQRT(a SUB 2+x*b SUB 2)+SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 2+y*b SUB 2))/~
(SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 1+y*b SUB 1)*(SQRT(a SUB 1+y*b SUB 1)*SQR~
T(a SUB 3+x*b SUB 3)+SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 3+y*b SUB 3))*(SQRT(a~
 SUB 1+y*b SUB 1)*SQRT(a SUB 4+x*b SUB 4)+SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB ~
4+y*b SUB 4))))-b SUB 2*I414(a,b,x,y))/(a SUB 2*b SUB 1-a SUB 1*b SUB 2)

U10A(a,b,x,y):=SQRT(U12A(a,b,x,y)^2-b SUB 2*(a SUB 3*b SUB 1-a SUB 1*b SUB 3)~
*(a SUB 4*b SUB 1-a SUB 1*b SUB 4)/b SUB 1)

U10T(a,b,x,y):=SQRT(u12^2-b SUB 2*(a SUB 3*b SUB 1-a SUB 1*b SUB 3)*(a SUB 4*~
b SUB 1-a SUB 1*b SUB 4)/b SUB 1)

Q10A(a,b,x,y):=U10A(a,b,x,y)/(SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 1+y*b SUB 1)~
)

Q10B(a,b,x,y):=u10/(SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 1+y*b SUB 1))

S10A(a,b,x,y):=(SQRT(a SUB 2+x*b SUB 2)*SQRT(a SUB 3+x*b SUB 3)*SQRT(a SUB 4+~
x*b SUB 4)/SQRT(a SUB 1+x*b SUB 1)+SQRT(a SUB 2+y*b SUB 2)*SQRT(a SUB 3+y*b S~
UB 3)*SQRT(a SUB 4+y*b SUB 4)/SQRT(a SUB 1+y*b SUB 1))/(x-y)

B10A(a,b,u,v,w,z):=-b SUB 1/((a SUB 1*b SUB 2-a SUB 2*b SUB 1)*(a SUB 3*b SUB~
 1-a SUB 1*b SUB 3)*(a SUB 4*b SUB 1-a SUB 1*b SUB 4))*(z+u)*(z+v)*(z+w)

B10B(a,b):=B10A(a,b,u12,u13,u23,u10)

I429_2(a,b,p):=4/p*R_C(1,1-b SUB 2*b SUB 3*b SUB 4/(p^2*b SUB 1))

I429_1A(a,b,x,y):=I429_2(a,b,I417_P(B10A(a,b,U12A(a,b,x,y),U13A(a,b,x,y),U23A~
(a,b,x,y),U10A(a,b,x,y)),S10A(a,b,x,y)+Q10A(a,b,x,y),b SUB 2*b SUB 3*b SUB 4/~
b SUB 1))

I429_1(a,b,x,y):=I429_2(a,b,I417_P(b10,s10+q10,b SUB 2*b SUB 3*b SUB 4/b SUB ~
1))

I429T(a,b,x,y):=(a SUB 2*b SUB 1-a SUB 1*b SUB 2)*(a SUB 3*b SUB 1-a SUB 1*b ~
SUB 3)*(a SUB 4*b SUB 1-a SUB 1*b SUB 4)/b SUB 1*I417_3(I417_4([U12A(a,b,x,y)~
,U13A(a,b,x,y),U23A(a,b,x,y),U10A(a,b,x,y)]))+I429_1A(a,b,x,y)

I429(a,b,x,y):=(a SUB 2*b SUB 1-a SUB 1*b SUB 2)*(a SUB 3*b SUB 1-a SUB 1*b S~
UB 3)*(a SUB 4*b SUB 1-a SUB 1*b SUB 4)/b SUB 1*I417_3(I417_4([u12,u13,u23,u1~
0]))+I429_1(a,b,x,y)

I414A(a,b,x,y):=I414([a SUB 1,a SUB 2,a SUB 3,1],[b SUB 1,b SUB 2,b SUB 3,0],~
x,y)

I417A(a,b,x,y):=I417([a SUB 1,a SUB 2,a SUB 3,1,a SUB 4],[b SUB 1,b SUB 2,b S~
UB 3,0,b SUB 4],x,y)

I427A(a,b,x,y):=I427([a SUB 1,a SUB 2,a SUB 3,1,a SUB 4],[b SUB 1,b SUB 2,b S~
UB 3,0,b SUB 4],x,y)

I428A(a,b,x,y):=I428([a SUB 1,a SUB 2,a SUB 3,1],[b SUB 1,b SUB 2,b SUB 3,0],~
x,y)

I429A(a,b,x,y):=I429([a SUB 1,a SUB 2,a SUB 3,1],[b SUB 1,b SUB 2,b SUB 3,0],~
x,y)

IA(x):=

list1998:=[[IA(0),'I414(a,b,x,y)],[IA(e1-e5),'I417(a,b,x,y)],[IA(-e5),'I427(a~
,b,x,y)],[IA(-e1),'I428(a,b,x,y)],[IA(-e2),'I428([a SUB 2,a SUB 1,a SUB 3,a S~
UB 4],[b SUB 2,b SUB 1,b SUB 3,b SUB 4],x,y)],[IA(-e3),'I428([a SUB 3,a SUB 2~
,a SUB 1,a SUB 4],[b SUB 3,b SUB 2,b SUB 1,b SUB 4],x,y)],[IA(-e4),'I428([a S~
UB 4,a SUB 2,a SUB 3,a SUB 1],[b SUB 4,b SUB 2,b SUB 3,b SUB 1],x,y)],[IA(e1)~
,'I429(a,b,x,y)],[IA(e2),'I429([a SUB 2,a SUB 1,a SUB 3,a SUB 4],[b SUB 2,b S~
UB 1,b SUB 3,b SUB 4],x,y)],[IA(e3),'I429([a SUB 3,a SUB 2,a SUB 1,a SUB 4],[~
b SUB 3,b SUB 2,b SUB 1,b SUB 4],x,y)],[IA(e4),'I429([a SUB 4,a SUB 2,a SUB 3~
,a SUB 1],[b SUB 4,b SUB 2,b SUB 3,b SUB 1],x,y)]]

IAS_AUX(v):=SELECT(IF(k SUB 1=v,1,0)=1,k,list1998)

IAS(v):=VECTOR((IAS_AUX(v SUB i_)) SUB 1 SUB 2,i_,1,DIMENSION(v))

list1998a:=[[IA(0),'I414A(a,b,x,y)],[IA(e1-e4),'I417A(a,b,x,y)],[IA(-e4),'I42~
7A(a,b,x,y)],[IA(-e1),'I428A(a,b,x,y)],[IA(-e2),'I428A([a SUB 2,a SUB 1,a SUB~
 3],[b SUB 2,b SUB 1,b SUB 3],x,y)],[IA(-e3),'I428A([a SUB 3,a SUB 2,a SUB 1]~
,[b SUB 3,b SUB 2,b SUB 1],x,y)],[IA(e1),'I429A(a,b,x,y)],[IA(e2),'I429A([a S~
UB 2,a SUB 1,a SUB 3],[b SUB 2,b SUB 1,b SUB 3],x,y)],[IA(e3),'I429A([a SUB 3~
,a SUB 2,a SUB 1],[b SUB 3,b SUB 2,b SUB 1],x,y)]]

IASA_AUX(v):=SELECT(IF(k SUB 1=v,1,0)=1,k,list1998a)

IASA(v):=VECTOR((IASA_AUX(v SUB i_)) SUB 1 SUB 2,i_,1,DIMENSION(v))

U12Y(a,b,y):=SQRT(b SUB 1)*SQRT(b SUB 2)*SQRT(a SUB 3+y*b SUB 3)*SQRT(a SUB 4~
+y*b SUB 4)+SQRT(b SUB 3)*SQRT(b SUB 4)*SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 2+~
y*b SUB 2)

U13Y(a,b,y):=SQRT(b SUB 1)*SQRT(b SUB 3)*SQRT(a SUB 2+y*b SUB 2)*SQRT(a SUB 4~
+y*b SUB 4)+SQRT(b SUB 2)*SQRT(b SUB 4)*SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 3+~
y*b SUB 3)

U23Y(a,b,y):=SQRT(b SUB 2)*SQRT(b SUB 3)*SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 4~
+y*b SUB 4)+SQRT(b SUB 1)*SQRT(b SUB 4)*SQRT(a SUB 2+y*b SUB 2)*SQRT(a SUB 3+~
y*b SUB 3)

I414Y(a,b,y):=I414_1(I414_2([U12Y(a,b,y),U13Y(a,b,y),U23Y(a,b,y)]))

I48(z):=I414Y(z,[1,1,1,1],0)

I414X(a,b,x):=I414Y(a,-b,-x)

U15Y(a,b,y):=SQRT(U12Y(a,b,y)^2-(a SUB 3*b SUB 1-a SUB 1*b SUB 3)*(a SUB 4*b ~
SUB 1-a SUB 1*b SUB 4)*(a SUB 5*b SUB 2-a SUB 2*b SUB 5)/(a SUB 5*b SUB 1-a S~
UB 1*b SUB 5))

Q15Y(a,b,y):=SQRT(b SUB 5)*SQRT(a SUB 5+y*b SUB 5)/(SQRT(b SUB 1)*SQRT(a SUB ~
1+y*b SUB 1))*U15Y(a,b,y)

S15Y(a,b,y):=(a SUB 5*SQRT(b SUB 2)*SQRT(b SUB 3)*SQRT(b SUB 4)*SQRT(a SUB 1+~
y*b SUB 1)+b SUB 5*(SQRT(b SUB 1)*SQRT(a SUB 2+y*b SUB 2)*SQRT(a SUB 3+y*b SU~
B 3)*SQRT(a SUB 4+y*b SUB 4)+y*SQRT(b SUB 2)*SQRT(b SUB 3)*SQRT(b SUB 4)*SQRT~
(a SUB 1+y*b SUB 1)))/(SQRT(b SUB 1)*SQRT(a SUB 1+y*b SUB 1))

I417Y_1(a,b,y):=I417_2(a,b,I417_P(B15A(a,b,U12Y(a,b,y),U13Y(a,b,y),U23Y(a,b,y~
),U15Y(a,b,y)),S15Y(a,b,y)+Q15Y(a,b,y),(a SUB 5*b SUB 2-a SUB 2*b SUB 5)*(a S~
UB 5*b SUB 3-a SUB 3*b SUB 5)*(a SUB 5*b SUB 4-a SUB 4*b SUB 5)/(a SUB 5*b SU~
B 1-a SUB 1*b SUB 5)))

I417Y(a,b,y):=(a SUB 2*b SUB 1-a SUB 1*b SUB 2)*(a SUB 3*b SUB 1-a SUB 1*b SU~
B 3)*(a SUB 4*b SUB 1-a SUB 1*b SUB 4)/(a SUB 5*b SUB 1-a SUB 1*b SUB 5)*I417~
_3(I417_4([U12Y(a,b,y),U13Y(a,b,y),U23Y(a,b,y),U15Y(a,b,y)]))+I417Y_1(a,b,y)

I410(z):=I417Y(z,[1,1,1,1,1],0)

I417X(a,b,x):=I417Y(a,-b,-x)

I427Y(a,b,y):=(b SUB 5*I417Y(a,b,y)-b SUB 1*I414Y(a,b,y))/(a SUB 1*b SUB 5-a ~
SUB 5*b SUB 1)

I427X(a,b,x):=I427Y(a,-b,-x)

I428Y(a,b,y):=(b SUB 1*2*((a SUB 2*b SUB 3-a SUB 3*b SUB 2)*(a SUB 2*b SUB 4-~
a SUB 4*b SUB 2)*I428_1(I414_2([U12Y(a,b,y),U13Y(a,b,y),U23Y(a,b,y)]))/3+(SQR~
T(b SUB 1)*SQRT(a SUB 2+y*b SUB 2)+SQRT(b SUB 2)*SQRT(a SUB 1+y*b SUB 1))*(SQ~
RT(b SUB 1)*SQRT(a SUB 3+y*b SUB 3)-SQRT(b SUB 3)*SQRT(a SUB 1+y*b SUB 1))*(S~
QRT(b SUB 1)*SQRT(a SUB 4+y*b SUB 4)-SQRT(b SUB 4)*SQRT(a SUB 1+y*b SUB 1))/(~
SQRT(b SUB 1)*SQRT(a SUB 1+y*b SUB 1)*(a SUB 3*b SUB 1-a SUB 1*b SUB 3)*(a SU~
B 4*b SUB 1-a SUB 1*b SUB 4)))-b SUB 2*I414Y(a,b,y))/(a SUB 2*b SUB 1-a SUB 1~
*b SUB 2)

I428X(a,b,x):=I428Y(a,-b,-x)

list1998y:=[[IA(0),'I414Y(a,b,y)],[IA(e1-e5),'I417Y(a,b,y)],[IA(-e5),'I427Y(a~
,b,y)],[IA(-e1),'I428Y(a,b,y)],[IA(-e2),'I428Y([a SUB 2,a SUB 1,a SUB 3,a SUB~
 4],[b SUB 2,b SUB 1,b SUB 3,b SUB 4],y)],[IA(-e3),'I428Y([a SUB 3,a SUB 2,a ~
SUB 1,a SUB 4],[b SUB 3,b SUB 2,b SUB 1,b SUB 4],y)],[IA(-e4),'I428Y([a SUB 4~
,a SUB 2,a SUB 3,a SUB 1],[b SUB 4,b SUB 2,b SUB 3,b SUB 1],y)]]

IASY_AUX(v):=SELECT(IF(k SUB 1=v,1,0)=1,k,list1998y)

IASY(v):=VECTOR((IASY_AUX(v SUB i_)) SUB 1 SUB 2,i_,1,DIMENSION(v))

list1998x:=[[IA(0),'I414X(a,b,x)],[IA(e1-e5),'I417X(a,b,x)],[IA(-e5),'I427X(a~
,b,x)],[IA(-e1),'I428X(a,b,x)]]

IASX_AUX(v):=SELECT(IF(k SUB 1=v,1,0)=1,k,list1998x)

IASX(v):=VECTOR((IASX_AUX(v SUB i_)) SUB 1 SUB 2,i_,1,DIMENSION(v))

U12B(a,b,x,y):=(SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 2+x*b SUB 2)+SQRT(a SUB 1+~
y*b SUB 1)*SQRT(a SUB 2+y*b SUB 2))/(x-y)

U10B(a,b,x,y):=(SQRT(a SUB 1+y*b SUB 1)*SQRT(a SUB 2+x*b SUB 2)+SQRT(a SUB 1+~
x*b SUB 1)*SQRT(a SUB 2+y*b SUB 2))/(x-y)

I414B_1(u,v):=4*R_C(1,2*v/(u+v))/(u+v)

I414B(a,b,x,y):=I414B_1(U12B(a,b,x,y),*U10B(a,b,x,y))

U13C(a,b,x,y):=SQRT(a SUB 3+x*b SUB 3)*SQRT(a SUB 3+y*b SUB 3)*(SQRT(a SUB 1+~
y*b SUB 1)*SQRT(a SUB 2+x*b SUB 2)+SQRT(a SUB 1+x*b SUB 1)*SQRT(a SUB 2+y*b S~
UB 2))/(x-y)

I427B_1(x1,x2,x3,y1,y2,y3):=SQRT(x1)*SQRT(x2)*y3+SQRT(y1)*SQRT(y2)*x3

U12C(a,b,x,y):=I427B_1(a SUB 1+x*b SUB 1,a SUB 2+x*b SUB 2,a SUB 3+x*b SUB 3,~
a SUB 1+y*b SUB 1,a SUB 2+y*b SUB 2,a SUB 3+y*b SUB 3)/(x-y)

I427B(a,b,x,y):=I414B_1(U12C(a,b,x,y),U13C(a,b,x,y))

I417B(a,b,x,y):=(I427B(a,b,x,y)*(a SUB 1*b SUB 3-a SUB 3*b SUB 1)+I414B(a,b,x~
,y)*b SUB 1)/b SUB 3

I428B_1(x1,x2,y1,y2):=SQRT(x1)*y1*SQRT(x2)+SQRT(y1)*x1*SQRT(y2)

I428B(a,b,x,y):=2*(x-y)/I428B_1(a SUB 1+x*b SUB 1,a SUB 2+x*b SUB 2,a SUB 1+y~
*b SUB 1,a SUB 2+y*b SUB 2)

list1998b:=[[IA(0),'I414B(a,b,x,y)],[IA(e1-e3),'I417B(a,b,x,y)],[IA(-e3),'I42~
7B(a,b,x,y)],[IA(-e1),'I428B(a,b,x,y)],[IA(-e2),'I428B([a SUB 2,a SUB 1],[b S~
UB 2,b SUB 1],x,y)]]

IASB_AUX(v):=SELECT(IF(k SUB 1=v,1,0)=1,k,list1998b)

IASB(v):=VECTOR((IASB_AUX(v SUB i_)) SUB 1 SUB 2,i_,1,DIMENSION(v))

I414C(a,b,x,y):=2*SQRT(a SUB 1+x*b SUB 1)/b SUB 1-2*SQRT(a SUB 1+y*b SUB 1)/b~
 SUB 1

U12D(a,b,x,y):=((a SUB 2+y*b SUB 2)*SQRT(a SUB 1+x*b SUB 1)+(a SUB 2+x*b SUB ~
2)*SQRT(a SUB 1+y*b SUB 1))/(x-y)

U10D(a,b,x,y):=SQRT(a SUB 2+x*b SUB 2)*SQRT(a SUB 2+y*b SUB 2)*(SQRT(a SUB 1+~
y*b SUB 1)+SQRT(a SUB 1+x*b SUB 1))/(x-y)

I427C(a,b,x,y):=I414B_1(U12D(a,b,x,y),U10D(a,b,x,y))

I417C(a,b,x,y):=(I427C(a,b,x,y)*(a SUB 1*b SUB 2-a SUB 2*b SUB 1)+2*(SQRT(a S~
UB 1+x*b SUB 1)-SQRT(a SUB 1+y*b SUB 1)))/b SUB 2

I428C(a,b,x,y):=2*(1/SQRT(a SUB 1+y*b SUB 1)-1/SQRT(a SUB 1+x*b SUB 1))/b SUB~
 1

list1998c:=[[IA(0),'I414C(a,b,x,y)],[IA(e1-e2),'I417C(a,b,x,y)],[IA(-e2),'I42~
7C(a,b,x,y)],[IA(-e1),'I428C(a,b,x,y)]]

IASC_AUX(v):=SELECT(IF(k SUB 1=v,1,0)=1,k,list1998c)

IASC(v):=VECTOR((IASC_AUX(v SUB i_)) SUB 1 SUB 2,i_,1,DIMENSION(v))

IS(v,h):=IF(h=4,IAS(v),IF(h=3,IASA(v),IF(h=2,IASB(v),IF(h=1,IASC(v),[]))))

I414AY(a,b,y):=I414Y([a SUB 1,a SUB 2,a SUB 3,1],[b SUB 1,b SUB 2,b SUB 3,0],~
y)

I414AX(a,b,x):=I414AY(a,-b,-x)

I427AY(a,b,y):=I427Y([a SUB 1,a SUB 2,a SUB 3,1,a SUB 4],[b SUB 1,b SUB 2,b S~
UB 3,0,b SUB 4],y)

I427AX(a,b,x):=I427AY(a,-b,-x)

I417AY(a,b,y):=I417Y([a SUB 1,a SUB 2,a SUB 3,1,a SUB 4],[b SUB 1,b SUB 2,b S~
UB 3,0,b SUB 4],y)

I417AX(a,b,x):=I417AY(a,-b,-x)

I428AY(a,b,y):=I428Y([a SUB 1,a SUB 2,a SUB 3,1],[b SUB 1,b SUB 2,b SUB 3,0],~
y)

I428AX(a,b,x):=I428AY(a,-b,-x)

list1998ay:=[[IA(0),'I414AY(a,b,y)],[IA(e1-e4),'I417AY(a,b,y)],[IA(-e4),'I427~
AY(a,b,y)],[IA(-e1),'I428AY(a,b,y)],[IA(-e2),'I428AY([a SUB 2,a SUB 1,a SUB 3~
],[b SUB 2,b SUB 1,b SUB 3],y)],[IA(-e3),'I428AY([a SUB 3,a SUB 2,a SUB 1],[b~
 SUB 3,b SUB 2,b SUB 1],y)]]

IASAY_AUX(v):=SELECT(IF(k SUB 1=v,1,0)=1,k,list1998ay)

IASAY(v):=VECTOR((IASAY_AUX(v SUB i_)) SUB 1 SUB 2,i_,1,DIMENSION(v))

ISY(v,h):=IF(h=4,IASY(v),IF(h=3,IASAY(v),[]))

list1998ax:=[[IA(0),'I414AX(a,b,x)],[IA(e1-e4),'I417AX(a,b,x)],[IA(-e4),'I427~
AX(a,b,x)],[IA(-e1),'I428AX(a,b,x)]]

IASAX_AUX(v):=SELECT(IF(k SUB 1=v,1,0)=1,k,list1998ax)

IASAX(v):=VECTOR((IASAX_AUX(v SUB i_)) SUB 1 SUB 2,i_,1,DIMENSION(v))

ISX(v,h):=IF(h=4,IASX(v),IF(h=3,IASAX(v),[]))

IHAT(x):=

ISHAT(v,a,b,x,y):=SUBST(v,[IHAT(0),IHAT(e1-e5),IHAT(-e5),IHAT(-e6),IHAT(-e1),~
IHAT(-e2),IHAT(-e3),IHAT(-e4),IHAT(e1),IHAT(e2),IHAT(e3),IHAT(e4)],['I414(a,b~
,x,y),'I417(a,b,x,y)*b SUB 5/b SUB 1,'I427(a,b,x,y)*b SUB 5,'I427([a SUB 1,a ~
SUB 2,a SUB 3,a SUB 4,a SUB 6],[b SUB 1,b SUB 2,b SUB 3,b SUB 4,b SUB 6],x,y)~
*b SUB 6,'I428(a,b,x,y)*b SUB 1,'I428([a SUB 2,a SUB 1,a SUB 3,a SUB 4],[b SU~
B 2,b SUB 1,b SUB 3,b SUB 4],x,y)*b SUB 2,'I428([a SUB 3,a SUB 2,a SUB 1,a SU~
B 4],[b SUB 3,b SUB 2,b SUB 1,b SUB 4],x,y)*b SUB 3,'I428([a SUB 4,a SUB 2,a ~
SUB 3,a SUB 1],[b SUB 4,b SUB 2,b SUB 3,b SUB 1],x,y)*b SUB 4,'I429(a,b,x,y)/~
b SUB 1,'I429([a SUB 2,a SUB 1,a SUB 3,a SUB 4],[b SUB 2,b SUB 1,b SUB 3,b SU~
B 4],x,y)/b SUB 2,'I429([a SUB 3,a SUB 2,a SUB 1,a SUB 4],[b SUB 3,b SUB 2,b ~
SUB 1,b SUB 4],x,y)/b SUB 3,'I429([a SUB 4,a SUB 2,a SUB 3,a SUB 1],[b SUB 4,~
b SUB 2,b SUB 3,b SUB 1],x,y)/b SUB 4])

ISU(a,b,x,y):=[u12=U12A(a,b,x,y),u13=U13A(a,b,x,y),u23=U23A(a,b,x,y),u15=U15B~
(a,b,x,y),q15=Q15B(a,b,x,y),s15=S15A(a,b,x,y),b15=B15B(a,b),u10=U10T(a,b,x,y)~
,q10=Q10B(a,b,x,y),s10=S10A(a,b,x,y),b10=B10B(a,b)]