Base Elemet Box Cat Code
Base Elemet Box Cat Code
In[]:=
ClearAll["Global`*"]
In[]:=
ver3c={X32,X33};B3C=1/16t(m1^4t+m2^4t+t(s+X32)^2-2m2^2(st-(t-2X32)(X32-X33))-2m1^2(m2^2t+st+(t-2X33)(X32-X33))+2stX33-2(2s+t)X32X33+tX33^2);
In[]:=
ureg3c=Product[ver3c[[i]]^rho3[i],{i,1,2}]*B3C^((d-5)/2);
In[]:=
uregnum3c=ureg3c/.{s:>17,t:>23,m1:>37,m2:>43,d:>1/97,rho[1]:>1/113,rho[2]:>1/149,rho[3]:>1/173,rho[4]:>1/197};
In[]:=
subsets3c=Subsets[ver3c]
Out[]=
{{},{X32},{X33},{X32,X33}}
In[]:=
Do[currentvars3c=subsets3c[[k]];remvars3c=Join[Complement[ver3c,currentvars3c],{s,t,d,m1,m2,rho3[1],rho3[2],rho3[3],rho3[4]}];numrules3c=Table[remvars3c[[i]]->1/Prime[23+45*i],{i,Length[remvars3c]}];dlistdo3c=Table[Simplify[D[Log[ureg3c/.numrules3c],currentvars3c[[j]]]],{j,1,Length[currentvars3c]}];SOLdo3c=Solve[dlistdo3c==0,currentvars3c];lenSOLdo3c=Length[SOLdo3c];Print[k,"/",Length[subsets3c]," ",currentvars3c," ",lenSOLdo3c];,{k,2,Length[subsets3c]}];
2/4 {X32} 2
3/4 {X33} 2
4/4 {X32,X33} 4
In[]:=
Do[currentvars3c=Take[ver3c,ii];Print["\nFor set ",currentvars3c," :"];remvars3c=Join[Complement[ver3c,currentvars3c],{s,t,d,m1,m2,rho3[1],rho3[2],rho3[3],rho3[4]}];numrules3c=Table[remvars3c[[i]]->1/Prime[23+45*i],{i,Length[remvars3c]}];Do[currentset3c=subsets3c[[k]];ureg32c=Product[currentset3c[[i]]^rho3[i],{i,1,Length[currentset3c]}]*B3C^((d-5)/2);ureg32c=ureg32c/.numrules3c;dlistdo3c=Table[Simplify[D[Log[ureg32c],currentvars3c[[j]]]],{j,1,Length[currentvars3c]}];SOLdo3c=Solve[dlistdo3c==0,currentvars3c];lenSOLdo3c=Length[SOLdo3c];Print[k,"/",Length[subsets3c]," ",currentset3c," ",lenSOLdo3c];,{k,1,Length[subsets3c]}];,{ii,1,Length[ver3c]}];
For set {X32} :
1/4 {} 1
2/4 {X32} 2
3/4 {X33} 1
4/4 {X32,X33} 2
For set {X32,X33} :
1/4 {} 1
2/4 {X32} 2
3/4 {X33} 2
4/4 {X32,X33} 4