# SU(5) normalized

# To run this file in Maple, first load the file "MapleCode".
# (the file RunAll will load+run all files including thise one).

G := (z*(1+y))/(1-exp(-z*(1+y)))-y*z;

X4:= ( subs(z=E[4],G)*subs(z=U[4]-E[4],G)*subs(z=V[4]-E[4],G)*(3*H+6*L-2*E[1]-2*E[2]-E[3]-E[4])  )/
     ( subs(z=3*H+6*L-2*E[1]-2*E[2]-E[3]-E[4],G)*subs(z=U[4],G)*subs(z=V[4],G)  );
X3:= ( subs(z=E[3],G)*subs(z=U[3]-E[3],G)*subs(z=V[3]-E[3],G)  )/
     ( subs(z=U[3],G)*subs(z=V[3],G)  );
X2 := (  subs(z=E[2],G)*subs(z=U[2]-E[2],G)*subs(z=V[2]-E[2],G)*subs(z=W[2]-E[2],G)  )/
     ( subs(z=U[2],G)*subs(z=V[2],G)*subs(z=W[2],G) );
X1 := (  subs(z=E[1],G)*subs(z=U[1]-E[1],G)*subs(z=V[1]-E[1],G)*subs(z=W[1]-E[1],G)*subs(z=H,G)  )/
     (  subs(z=W[1],G)  );

C := X4*X3*X2*X1;

C4 := BlowUp(E[4], {U[4],V[4]}, C):
C4 := subs(U[4]=H+3*L-E[1]-E[2]-E[3],V[4]=E[2],C4):

C3 := BlowUp(E[3], {U[3],V[3]}, C4):
C3 := subs(U[3]=H+3*L-E[1]-E[2],V[3]=E[1]-E[2],C3):

C2 := BlowUp(E[2], {U[2],V[2],W[2]}, C3):
C2 := subs(U[2]=H+2*L-E[1],V[2]=H+3*L-E[1],W[2]=E[1],C2):

C1 := BlowUp(E[1],{U[1],V[1],W[1]},C2):
C1 := subs(U[1] = H + 2*L, V[1] = H + 3*L, W[1] = S, C1):

Q[SU5] := PiStar(C1, [H, H + 2*L, H + 3*L]):

# Q[SU5] evaluated at y=-1 should yield (12*L+42*L^2*S+12*L^2-35*L*S^2+32*L*S-30*S^2)/(1+L)/(1+S)/(1+6*L-5*S)
# Checks out.