home *** CD-ROM | disk | FTP | other *** search
/ Frozen Fish 1: Amiga / FrozenFish-Apr94.iso / bbs / alib / d5xx / d568 / schoonschip.lha / Schoonschip / SDocsAMI.LZH / SDocsAMI / Roblocks.e < prev    next >
Encoding:
Text File  |  1991-06-14  |  19.8 KB  |  493 lines
  1. C RO-parameter. Blocks.
  2.  
  3. BLOCK VERT{2,3,4,5,6}
  4. C 2=THREE-VERTICES WITH DERIVATIVES, 3=WITHOUT DERIVATIVES,
  5. C 4=FOUR-VERTICES,
  6. C 5=GHOST VERTICES WITH DERIVATIVES, 6=GHOST WITHOUT DERIVATIVES.
  7. S EP,TIPI2,N,PM,M,M0,SM,C,S,E=3,LE2
  8. F FK,FP,FQ,FSC,F,VE4,VE3,PROP,LG,FALE
  9. I AL=N,BE=N,GA=N,MU=N,NU=N,MUP=N,NUP=N,
  10.   L1=N,L2=N,L3=N,L4=N,L5=N,L6=N,L7=N,L8=N,L9=N,L0=N
  11. V P,Q,K,P0=z,PL,LE
  12. C A=A, F=FI, G=FP, H=FM, U=WP, U_=WM, W=W0, Z=HIGGS
  13. C X=ACP, Y=CP, O=ACM, P=CM, Q=AY0, R=Y0, S=AYA, T=YA.
  14. C        A B C D E F G H  I J K L M N O P  Q R S T U V W  X Y Z
  15. T QN(K1)=0,0,0,0,0,0,1,-1,0,0,0,0,0,0,1,-1,0,0,0,0,1,0,0,-1,1,0,
  16.  @G_i0,-1,@O_i0,-1,1,0,0,0,0,-1,0,0,1,-1,0
  17. C         O_P_Q_R_S_T_U_ V_W_X_Y_ Z_
  18. C        A B C D E F G H I J K L M N O  P  Q R S  T U V W  X Y Z
  19. T XN(K1)=0,0,0,0,0,0,0,0,0,0,0,0,0,0,-1,1,-1,1,-1,1,0,0,0,-1,1,0,
  20.  @O_i0,1,-1,1,-1,1,-1,0,0,0,1,-1,0
  21. C      O_P_ Q_R_ S_T_ U_V_W_X_Y_ Z_
  22. C QN, QG AND QH ARE THE CHARGES, XN, XG AND XH GHOST NRS.
  23. T CHG(K1)="A,"W,"F,"Z,"U,"U_,"G,"G_,"Y_,"Y,"P_,"P,"R_,"R,"T_,"T
  24. T QG(K1) = 0, 0, 0, 0, 1,-1,  1,-1,-1, 1, 1,-1, 0, 0, 0, 0
  25. T XG(K1) = 0, 0, 0, 0, 0, 0, 0, 0,-1, 1,-1, 1,-1, 1,-1, 1
  26. T CHH(K1)="A,"W,"F,"Z,"U_,"U,"G_,"G,"Y,"Y_,"P,"P_,"R,"R_,"T,"T_
  27. T QH(K1) = 0, 0, 0, 0,-1,  1,-1, 1, 1,-1,-1, 1, 0, 0, 0, 0
  28. T XH(K1) = 0, 0, 0, 0, 0, 0, 0, 0, 1,-1, 1,-1, 1,-1, 1,-1
  29. T TAP:A,F,W,Z
  30. T TFE:P,R,T,Y
  31. T THI: Z=1:1
  32. X DEDE(L1,L2,L3,L4)=2.*D(L1,L2)*D(L3,L4) - D(L1,L3)*D(L2,L4)
  33.                                          - D(L1,L4)*D(L2,L3)
  34. X AA(L1,L2;K): NOM(K,PM)*D(L1,L2)
  35. X FF(L1,L2;K): NOM(K,M0)
  36. X GG_(L1,L2;K): NOM(K,M)
  37. X PP_(L1,L2;K): NOM(K,M)
  38. X RR_(L1,L2;K): NOM(K,M0)
  39. X TT_(L1,L2;K): NOM(K,PM)
  40. X UU_(L1,L2;K): NOM(K,M)*D(L1,L2)
  41. X WW(L1,L2;K): NOM(K,M0)*D(L1,L2)
  42. X YY_(L1,L2;K): NOM(K,M)
  43. X ZZ(L1,L2;K): NOM(K,SM)
  44. G 1
  45. '2'X AGG_(AL,P,BE,Q,GA,K): S*Q(AL)-S*K(AL)
  46. '3'X AGU_(AL,P,BE,Q,GA,K): -I*S*M*D(AL,GA)
  47. '3'X AUG_(AL,P,BE,Q,GA,K): I*S*M*D(AL,BE)
  48. '5'X APP_(AL,P,BE,Q,GA,K): S*K(AL)
  49. '2'X AUU_(AL,P,BE,Q,GA,K): -S*D(AL,BE)*(Q(GA)-P(GA))
  50. '2'       -S*D(AL,GA)*(P(BE)-K(BE))-S*D(BE,GA)*(K(AL)-Q(AL))
  51. '5'X AYY_(AL,P,BE,Q,GA,K): -S*K(AL)
  52. '3'X FFZ(AL,P,BE,Q,GA,K): -0.5*SM**2/M
  53. '2'X FGU_(AL,P,BE,Q,GA,K): -0.5*Q(GA)+0.5*P(GA)
  54. '2'X FUG_(AL,P,BE,Q,GA,K): 0.5*K(BE)-0.5*P(BE)
  55. '6'X FPP_(AL,P,BE,Q,GA,K): -0.5*I*M
  56. '2'X FWZ(AL,P,BE,Q,GA,K): I*(P(BE)-K(BE))/2/C
  57. '6'X FYY_(AL,P,BE,Q,GA,K): 0.5*I*M
  58. '2'X GWG_(AL,P,BE,Q,GA,K): C2MS2*(P(BE)-K(BE))/2/C
  59. '3'X GZG_(AL,P,BE,Q,GA,K): -0.5*SM**2/M
  60. '6'X GPR_(AL,P,BE,Q,GA,K): 0.5*I*M*C**-1
  61. '6'X GRY_(AL,P,BE,Q,GA,K): -0.5*I*C2MS2*M/C
  62. '6'X GTY_(AL,P,BE,Q,GA,K): -I*S*M
  63. '3'X GWU_(AL,P,BE,Q,GA,K): I*S**2*M*D(BE,GA)/C
  64. '2'X GZU_(AL,P,BE,Q,GA,K): 0.5*I*(P(GA)-Q(GA))
  65. '6'X RG_P_(AL,P,BE,Q,GA,K): 0.5*I*C2MS2*M/C
  66. '6'X TG_P_(AL,P,BE,Q,GA,K): I*S*M
  67. '6'X YG_R_(AL,P,BE,Q,GA,K): -0.5*I*M*C**-1
  68. '3'X UWG_(AL,P,BE,Q,GA,K): -I*S**2*M*D(AL,BE)/C
  69. '2'X UZG_(AL,P,BE,Q,GA,K): 0.5*I*(K(AL)-Q(AL))
  70. '5'X PWP_(AL,P,BE,Q,GA,K): C*K(BE)
  71. '6'X PZP_(AL,P,BE,Q,GA,K): -0.5*M
  72. '5'X RP_U_(AL,P,BE,Q,GA,K): -C*Q(GA)
  73. '5'X TP_U_(AL,P,BE,Q,GA,K): -S*Q(GA)
  74. '5'X PUR_(AL,P,BE,Q,GA,K): -C*K(BE)
  75. '5'X PUT_(AL,P,BE,Q,GA,K): -S*K(BE)
  76. '6'X RZR_(AL,P,BE,Q,GA,K): -0.5*M/C**2
  77. '5'X YR_U_(AL,P,BE,Q,GA,K): C*Q(GA)
  78. '5'X RUY_(AL,P,BE,Q,GA,K): C*K(BE)
  79. '5'X YT_U_(AL,P,BE,Q,GA,K): S*Q(GA)
  80. '5'X TUY_(AL,P,BE,Q,GA,K): S*K(BE)
  81. '2'X UWU_(AL,P,GA,K,BE,Q): -C*D(AL,BE)*(Q(GA)-P(GA))
  82. '2'   -C*D(AL,GA)*(P(BE)-K(BE))-C*D(BE,GA)*(K(AL)-Q(AL))
  83. '3'X UZU_(AL,P,BE,Q,GA,K): -M*D(AL,GA)
  84. '3'X WWZ(AL,P,BE,Q,GA,K): -M*D(AL,BE)/C**2
  85. '5'X WYY_(AL,P,BE,Q,GA,K): -C*K(AL)
  86. '6'X YZY_(AL,P,BE,Q,GA,K): -0.5*M
  87. '3'X ZZZ(AL,P,BE,Q,GA,K): -1.5*SM**2/M
  88. G 2
  89. '4'X AAGG_(L1,L2,L3,L4): -2.*S**2*D(L1,L2)
  90. '4'X AAUU_(L1,L2,L3,L4): -S**2*DEDE(L1,L2,L3,L4)
  91. '4'X AFGU_(L1,L2,L3,L4): 0.5*S*D(L1,L4)
  92. '4'X AFUG_(L1,L2,L3,L4): 0.5*S*D(L1,L3)
  93. '4'X AGWG_(L1,L2,L3,L4): -S*C2MS2*D(L1,L3)/C
  94. '4'X AGZU_(L1,L2,L3,L4): -0.5*I*S*D(L1,L4)
  95. '4'X AUZG_(L1,L2,L3,L4): 0.5*I*S*D(L1,L2)
  96. '4'X AUWU_(L1,L2,L3,L4): -S*C*DEDE(L1,L3,L4,L2)
  97. '4'X FFFF(L1,L2,L3,L4): -0.75*SM**2/M**2
  98. '4'X FFGG_(L1,L2,L3,L4): -0.25*SM**2/M**2
  99. '4'X FFUU_(L1,L2,L3,L4): -0.5*D(L3,L4)
  100. '4'X FFWW(L1,L2,L3,L4): -0.5*D(L3,L4)/C**2
  101. '4'X FFZZ(L1,L2,L3,L4): -0.25*SM**2/M**2
  102. '4'X FGWU_(L1,L2,L3,L4): -0.5*S**2*D(L3,L4)/C
  103. '4'X FUWG_(L1,L2,L3,L4): -0.5*S**2*D(L2,L3)/C
  104. '4'X GGG_G_(L1,L2,L3,L4): -0.5*SM**2/M**2
  105. '4'X GUG_U_(L1,L2,L3,L4): -0.5*D(L2,L4)
  106. '4'X GWWG_(L1,L2,L3,L4): -0.5*C2MS2**2*D(L2,L3)/C**2
  107. '4'X GZZG_(L1,L2,L3,L4): -0.25*SM**2/M**2
  108. '4'X GWZU_(L1,L2,L3,L4): 0.5*I*S**2*D(L2,L4)/C
  109. '4'X UWZG_(L1,L2,L3,L4): -0.5*I*S**2*D(L1,L2)/C
  110. '4'X UUU_U_(L1,L2,L3,L4): DEDE(L1,L2,L3,L4)
  111. '4'X UWWU_(L1,L2,L3,L4): -C**2*DEDE(L1,L4,L3,L2)
  112. '4'X UZZU_(L1,L2,L3,L4): -0.5*D(L1,L4)
  113. '4'X WWZZ(L1,L2,L3,L4): -0.5*D(L1,L2)/C**2
  114. '4'X ZZZZ(L1,L2,L3,L4): -0.75*SM**2/M**2
  115. G
  116. ENDBLOCK
  117. BLOCK SVERT{2,3,4}
  118. C 2=TWO-POINT VERTICES, 3=THREE-POINT VERTICES,
  119. C 4=FOUR-POINT VERTICES.
  120. C ************** SPECIAL VERTICES **************************
  121. '3'X AFZ(AL,P,BE,Q,GA,K): I*(K(AL)-Q(AL))*S*TIPI2*ZERO/EP
  122. '3'X AWZ(AL,P,BE,Q,GA,K): D(AL,BE)*2*M*S*TIPI2*ZERO/EP/C
  123. '3'X TZR_(AL,P,BE,Q,GA,K): 0.5*M*S*TIPI2*ZERO/EP/C
  124. '4'X AFFW(L1,L2,L3,L4): D(L1,L4)*S*TIPI2*ZERO/EP/C
  125. '4'X AWZZ(L1,L2,L3,L4): D(L1,L2)*S*TIPI2*ZERO/EP/C
  126. '2'X ANW(AL,P,BE,Q,GA,K): 2*M**2*S*TIPI2/C*D(AL,GA)
  127. '2'      +PDK*(3*C*S+S/6/C)*D(AL,GA)*TIPI2-P(AL)*K(GA)*
  128. '2'       (3*C*S+S/6/C)*TIPI2
  129. '2'X AAN(AL,P,BE,Q,GA,K): 3*S**2*TIPI2*(PDQ*D(AL,BE)-P(AL)*Q(BE))
  130. '2'X AFN(AL,P,BE,Q,GA,K): 2*I*M*S*P(AL)*TIPI2
  131. '2'X FNW(AL,P,BE,Q,GA,K): I*K(GA)*(2*M*C-0.25*M/C**3-49*M/12/C
  132. '2'       -M*E1/C/TIPI2-M*EZF/C/TIPI2)*TIPI2
  133. '2'X FFN(AL,P,BE,Q,GA,K): -0.5*SM**2*ET+PDQ*TIPI2+0.5*PDQ*TIPI2/C**2
  134. '2'   +2*EZF*PDQ
  135. '2'X GNG_(AL,P,BE,Q,GA,K): -0.5*SM**2*ET+PDK*TIPI2+0.5*PDK*TIPI2/C**2
  136. '2'    +2*EZF*PDK
  137. '2'X GNU_(AL,P,BE,Q,GA,K): -I*K(GA)*TIPI2*(25*M/12+0.25*M/C**2+M*E1/TIPI2
  138. '2'   +M*EZF/TIPI2)
  139. '2'X NUG_(AL,P,BE,Q,GA,K): -I*Q(BE)*TIPI2*(25*M/12+0.25*M/C**2+M*E1/TIPI2
  140. '2'   +M*EZF/TIPI2)
  141. '2'X NUU_(AL,P,BE,Q,GA,K): QDK*D(BE,GA)*TIPI2*19/6-Q(BE)*K(GA)*TIPI2*19/6
  142. '2'    -D(BE,GA)*19/6*TIPI2*M**2-2*D(BE,GA)*M**2*E1
  143. '2'X NWW(AL,P,BE,Q,GA,K): (1/3+3*C**2-1/6/C**2)*TIPI2*
  144. '2'   (D(BE,GA)*QDK-Q(BE)*K(GA))
  145. '2'   -D(BE,GA)*TIPI2*M**2*(43/6/C**2-4)-2*D(BE,GA)*M**2*E1/C**2
  146. '2'X NZZ(AL,P,BE,Q,GA,K): QDK*TIPI2+0.5*QDK*TIPI2/C**2
  147. '2'   +(-SM**2-0.5*SM**2/C**2)*TIPI2
  148. '2'   +SM**2*(2*E2-2*E1-0.5*ET)
  149. '2'   -2*EZF*SM**2+2*QDK*EZF
  150. '2'X NPP_(AL,P,BE,Q,GA,K): TIPI2*M**2*(0.25/C**2-25/12)-M**2*E1+QDK*TIPI2
  151. '2'   +M**2*EZF
  152. '2'X NRR_(AL,P,BE,Q,GA,K): TIPI2*M**2*(1+0.25/C**4-37/12/C**2)
  153. '2'    -M**2*E1/C**2+QDK*C**2*TIPI2
  154. '2'   +M**2*EZF/C**2
  155. '2'X NTR_(AL,P,BE,Q,GA,K): TIPI2*S*(M**2/C+QDK*C)
  156. '2'X NRT_(AL,P,BE,Q,GA,K): TIPI2*C*S*QDK
  157. '2'X NTT_(AL,P,BE,Q,GA,K): TIPI2*QDK*S**2
  158. '2'X NYY_(AL,P,BE,Q,GA,K): TIPI2*M**2*(0.25/C**2-25/12)-M**2*E1+QDK*TIPI2
  159. '2'   +M**2*EZF
  160. C **********************************************************
  161. ENDBLOCK
  162. BLOCK CVERT{3,4,5}
  163. C 3=THREE-VERTICES, 4=FOUR-VERTICES, 5=GHOST VERTICES.
  164. '3'Id,AGG_~K: 0.5*TIPI2/C**2+2*EZF
  165. '3'Al,AGU_~K: TIPI2/4/C**2+25*TIPI2/12+E1+EZF
  166. '3'Al,AUG_~K: TIPI2/4/C**2+25*TIPI2/12+E1+EZF
  167. '3'Al,AUU_~K: 7*TIPI2/6
  168. '3'Al,FGU_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  169. '3'Al,FUG_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  170. '3'Al,FWZ~K: TIPI2-2*C**2*TIPI2+0.5*TIPI2/C**2+2*EZF
  171. '3'Al,FFZ~K: -25*TIPI2/12+3*TIPI2/4/C**2+E1-2*E2+3*EZF
  172. '3'Al,GZG_~K: -25*TIPI2/12+3*TIPI2/4/C**2+E1-2*E2+3*EZF
  173. '3'Al,GWG_~K: -0.5*TIPI2/C**2/C2MS2+2*EZF
  174. '3'Id,GWU_~K: E1+7*TIPI2/3/C**2-25*S**2*TIPI2/12/C**2+EZF
  175. '3'Al,GZU_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  176. '3'Al,UWG_~K: E1+7*TIPI2/3/C**2-25*S**2*TIPI2/12/C**2+EZF
  177. '3'Al,UZG_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  178. '3'Al,UWU_~K: 7*TIPI2/6
  179. '3'Al,UZU_~K: TIPI2/12+0.25*TIPI2/C**2+E1+EZF
  180. '3'Al,WWZ~K: -4*C**2*TIPI2+0.25*TIPI2/C**2+49*TIPI2/12+E1+EZF
  181. '3'Al,ZZZ~K: -25*TIPI2/12+0.75*TIPI2/C**2+E1-2*E2+3*EZF
  182. '5'Id,APP_~K: -TIPI2
  183. '5'Al,AYY_~K: -TIPI2
  184. '5'Al,FPP_~K: -TIPI2
  185. '5'Al,FYY_~K: -TIPI2
  186. '5'Al,GPR_~K: -TIPI2
  187. '5'Al,GRY_~K: -C**2*TIPI2/C2MS2
  188. '5'Al,GTY_~K: -0.5*TIPI2
  189. '5'Al,RP_G_~K: -TIPI2*C**2/C2MS2
  190. '5'Al,TG_P_~K: -TIPI2/2
  191. '5'Al,YG_R_~K: -TIPI2
  192. '5'Al,PZP_~K: -TIPI2
  193. '5'Al,PWP_~K: -TIPI2
  194. '5'Al,RP_U_~K: -TIPI2
  195. '5'Id,TP_U_~K: -TIPI2
  196. '5'Al,PUR_~K: -TIPI2
  197. '5'Al,PUT_~K: -TIPI2
  198. '5'Al,RZR_~K: -TIPI2*C**2
  199. '5'Al,YR_U_~K: -TIPI2
  200. '5'Al,RUY_~K: -TIPI2
  201. '5'Al,YT_U_~K: -TIPI2
  202. '5'Al,TUY_~K: -TIPI2
  203. '5'Al,WYY_~K: -TIPI2
  204. '5'Al,YZY_~K: -TIPI2
  205. '4'Id,AAUU_~K: -5*TIPI2/6
  206. '4'Al,AFGU_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  207. '4'Al,AFUG_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  208. '4'Al,AGWG_~K: -TIPI2-0.5*TIPI2/C**2/C2MS2+2*EZF
  209. '4'Al,AGZU_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  210. '4'Al,AUZG_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  211. '4'Al,AUWU_~K: -5*TIPI2/6
  212. '4'Al,AAGG_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  213. '4'Al,FGWU_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  214. '4'Al,FUWG_~K: -TIPI2+0.5*TIPI2/C**2+2*EZF
  215. '4'Id,FFGG_~K: -31*TIPI2/6+TIPI2/C**2-2*E2+4*EZF
  216. '4'Al,FFUU_~K: -3*TIPI2+0.5*TIPI2/C**2+2*EZF
  217. '4'Al,FFWW~K: -4*TIPI2*C**2+0.5*TIPI2/C**2+TIPI2+2*EZF
  218. '4'Al,FFZZ~K: -31*TIPI2/6+TIPI2/C**2-2*E2+4*EZF
  219. '4'Al,FFFF~K: -31*TIPI2/6+TIPI2/C**2-2*E2+4*EZF
  220. '4'Al,GUG_U_~K: -3*TIPI2+0.5*TIPI2/C**2+2*EZF
  221. '4'Al,GWWG_~K: -0.5*(4*C**4+1)*TIPI2/C**2/C2MS2+2*EZF
  222. '4'Id,GZZG_~K: -31*TIPI2/6+TIPI2/C**2-2*E2+4*EZF
  223. '4'Al,GWZU_~K: 0.5*TIPI2/C**2-TIPI2+2*EZF
  224. '4'Al,GGG_G_~K: -31*TIPI2/6+TIPI2/C**2-2*E2+4*EZF
  225. '4'Al,UWZG_~K: 0.5*TIPI2/C**2-TIPI2+2*EZF
  226. '4'Al,UUU_U_~K: -5*TIPI2/6
  227. '4'Al,UZZU_~K: -3*TIPI2+0.5*TIPI2/C**2+2*EZF
  228. '4'Al,UWWU_~K: -5*TIPI2/6
  229. '4'Al,WWZZ~K: TIPI2-4*TIPI2*C**2+0.5*TIPI2/C**2+2*EZF
  230. '4'Al,ZZZZ~K: -31*TIPI2/6+TIPI2/C**2-2*E2+4*EZF
  231. ENDBLOCK
  232. BLOCK SCVERT{3,4,5}
  233. C 3=THREE-VERTICES, 4=FOUR-VERTICES, 5=GHOST VERTICES.
  234. C *************** SPECIAL VERTICES **************************
  235. '3'Id,AFZ~K: EP*ZERO**-1
  236. '3'Al,AWZ~K: EP*ZERO**-1
  237. '5'Id TZR_~K: EP*ZERO**-1
  238. '4'Al,AFFW~K: EP*ZERO**-1
  239. '4'Al,AWZZ~K: EP*ZERO**-1
  240. C ***********************************************************
  241. ENDBLOCK
  242. BLOCK F1t5{}
  243. Id,F(1,M~)=2*I*PI**2*M**2*(1/EP-0.5+0.5*LG(M)
  244.    +EP*(1/4-LG(M)/4+LG(M)*LG(M)/8+PI**2/48))
  245. Al,F(2,M~)=2*I*PI**2*(-1/EP-0.5*LG(M)
  246.     -EP*(PI**2/48+LG(M)*LG(M)/8))
  247. Al,F(3,M~)=I*PI**2/2/M**2+I*EP*PI**2/2/M**2/2*LG(M)
  248. Al,F(4,M~)=I*PI**2/6/M**4
  249.    -1/12*I*EP*PI**2/M**4*(1-LG(M))
  250. Al,F(5,M~)=I*PI**2/M**6*(1/12-EP/16+EP*LG(M)/24)
  251. ENDBLOCK
  252. BLOCK F21F22{}
  253. Id,F221(M~,M~,M~,M~,SM)=PI**4/SM**2*(-PI**2/3
  254.    +2*LG(SM,M)-LG(SM,M)*LG(SM,M))
  255.    +PI**4*M**2/SM**4*(4-4*PI**2/3+12*LG(SM,M)
  256.    -4*LG(SM,M)*LG(SM,M))
  257. Al,F211(SM,SM,SM,SM)=PI**4*(-2/EP**2+(1-2*LG(SM) )/EP
  258.     -0.5-PI**2/12+2*CLPW3+LG(SM)-LG(SM)*LG(SM)  )
  259. Al,F211(SM,SM,M~,SM)=PI**4*(-2/EP**2+(1-2*LG(SM))/EP
  260.    -0.5-PI**2/12+LG(SM)-LG(SM)*LG(SM)-0.5*M**2/SM**2*LG(M,SM)
  261.    +M**2/SM**2 )
  262. Al,F211(M~,M~,SM,SM)=PI**4*(-2/EP**2+(1-2*LG(M))/EP
  263.    +1.5-PI**2/12+LG(M)-LG(M)*LG(M)+LG(SM)+LG(SM,M)*LG(SM,M)/2)
  264. Al,F211(M~,M~,M~,SM)=PI**4*(-2/EP**2+1/EP*(1-2*LG(M)))
  265.    +PI**4*(-0.5+PI**2/12+LG(M)-LG(M)*LG(M)
  266.    +0.5*LG(SM,M)*LG(SM,M)+(M**2*SM**2+3*M**4)/SM**4*LG(SM,M)
  267.    *LG(SM,M) - (2*M**2*SM**2+7*M**4)/SM**4*LG(SM,M)
  268.    +2*(PI**2/6-1)*M**2/SM**2 + (PI**2-11/2)*M**4/SM**4 )
  269. Al,F311(M~,M~,M~,M~,SM)=PI**4/M**2/EP+PI**4/M**2*(-0.5
  270.    +0.5*LG(M)+0.5*LG(SM)+M**2/SM**2*LG(SM,M) )
  271. Al,F211(SM,SM,M~,M0~)=PI**4*(-2/EP**2+(1-2*LG(SM))/EP
  272.    -0.5-PI**2/4+LG(SM)-LG(SM)*LG(SM)
  273.    -(0.5*LG(M,SM)*LG(M0,SM)+PI**2/6)*
  274.    (2*M**2*M0**2/SM**4+(4*M**4*M0**2+4*M**2*M0**4)/SM**6)
  275.    -0.5*LG(M,SM)*(2*M**2/SM**2
  276.      +(M**4+2*M0**2*M**2)/SM**4
  277.      +(2*M**6/3+8*M**4*M0**2+2*M**2*M0**4)/SM**6)
  278.    -0.5*LG(M0,SM)*(2*M0**2/SM**2
  279.      +(M0**4+2*M**2*M0**2)/SM**4
  280.      +(2*M0**6/3+8*M0**4*M**2+2*M0**2*M**4)/SM**6)
  281.    +(M**2+M0**2)/SM**2+(M**4/4+M**2*M0**2+M0**4/4)/SM**4
  282.    +(M**6/9+3*M**4*M0**2+3*M**2*M0**4+M0**6/9)/SM**6 )
  283. Al,F211(M~,M~,M0~,SM)=PI**4*(-2/EP**2+(1-2*LG(M))/EP
  284.   -0.5+PI**2/12+LG(M)-LG(M)*LG(M)+0.5*LG(SM,M)*LG(SM,M)
  285.   +(M0**2*SM**2+2*M**2*M0**2+M0**4)/SM**4*LG(SM,M)*LG(SM,M)
  286.   -(M0**2*SM**2+2*M**2*M0**2+M0**4)/SM**4*LG(SM,M)*LG(M0,M)
  287.   -(M**2*SM**2+M0**2*SM**2+4*M**2*M0**2+5/2*M0**4+0.5*M**4)
  288.   /SM**4*LG(SM,M) + (M0**2*SM**2+M**2*M0**2+5/2*M0**4)/SM**4
  289.   *LG(M0,M) - (M**2+M0**2)/SM**2
  290.   - (1/4*M**4+3*M**2*M0**2+9/4*M0**4)/SM**4
  291.   +PI**2/3*(M0**2*SM**2+2*M**2*M0**2+M0**4)/SM**4  )
  292. ENDBLOCK
  293. BLOCK F21F51{}
  294. Id,F311(M~,M0~,SM)=
  295.   +PI**4*(-1./2.*M**(-2)+M**(-2)*EP**(-1))
  296.   +PI**6*(-1./3.*M0**2*SM**(-4))
  297.   +LG(M)*PI**4*(-1./2.*M**2*SM**(-4)+1./2.*M**(-2)-3./2.*M0**2*SM*
  298.   *(-4)-1./2.*SM**(-2))
  299.   +LG(M)*LG(M0)*PI**4*(-M0**2*SM**(-4))
  300.   +LG(M)*LG(SM)*PI**4*(M0**2*SM**(-4))
  301.   +LG(M0)*PI**4*(-1./2.*M**(-2)*M0**2*SM**(-2)-1./2.*M**(-2)*M0**4
  302.   *SM**(-4)-3./2.*M0**2*SM**(-4))
  303.   +LG(M0)*LG(SM)*PI**4*(M0**2*SM**(-4))
  304.   +LG(SM)*PI**4*(1./2.*M**2*SM**(-4)+1./2.*M**(-2)+1./2.*M**(-2)*M
  305.   0**2*SM**(-2)+1./2.*M**(-2)*M0**4*SM**(-4)+3.*M0**2*SM**(-4)+1./
  306.   2.*SM**(-2))
  307.   +LG(SM)*LG(SM)*PI**4*(-M0**2*SM**(-4))
  308. Al,F221(M~,M0~,SM)=
  309.   +PI**4*(2.*M**2*SM**(-4)+2.*M0**2*SM**(-4))
  310.   +PI**6*(-2./3.*M**2*SM**(-4)-2./3.*M0**2*SM**(-4)-1./3.*SM**(-2)
  311.   )
  312.   +LG(M)*PI**4*(-5.*M**2*SM**(-4)-M0**2*SM**(-4)-SM**(-2))
  313.   +LG(M)*LG(M0)*PI**4*(-2.*M**2*SM**(-4)-2.*M0**2*SM**(-4)-SM**(-2
  314.   ))
  315.   +LG(M)*LG(SM)*PI**4*(2.*M**2*SM**(-4)+2.*M0**2*SM**(-4)+SM**(-2)
  316.   )
  317.   +LG(M0)*PI**4*(-M**2*SM**(-4)-5.*M0**2*SM**(-4)-SM**(-2))
  318.   +LG(M0)*LG(SM)*PI**4*(2.*M**2*SM**(-4)+2.*M0**2*SM**(-4)+SM**(-2
  319.   ))
  320.   +LG(SM)*PI**4*(6.*M**2*SM**(-4)+6.*M0**2*SM**(-4)+2.*SM**(-2))
  321.   +LG(SM)*LG(SM)*PI**4*(-2.*M**2*SM**(-4)-2.*M0**2*SM**(-4)-SM**(-
  322.   2))
  323. Al,F211(M~,M0~,SM)=
  324.   +PI**4*(-1./2.-3.*M**2*M0**2*SM**(-4)-M**2*SM**(-2)-1./4.*M**4*S
  325.   M**(-4)-M0**2*SM**(-2)-9./4.*M0**4*SM**(-4)-2.*EP**(-2)+EP**(-1)
  326.   )
  327.   +PI**6*(1./12.+2./3.*M**2*M0**2*SM**(-4)+1./3.*M0**2*SM**(-2)+1.
  328.   /3.*M0**4*SM**(-4))
  329.   +LG(M)*PI**4*(1.+3.*M**2*M0**2*SM**(-4)+M**2*SM**(-2)+1./2.*M**4
  330.   *SM**(-4)-2.*EP**(-1))
  331.   -1./2.*LG(M)*LG(M)*PI**4
  332.   +LG(M)*LG(M0)*PI**4*(2.*M**2*M0**2*SM**(-4)+M0**2*SM**(-2)+M0**4
  333.   *SM**(-4))
  334.   +LG(M)*LG(SM)*PI**4*(-1.-2.*M**2*M0**2*SM**(-4)-M0**2*SM**(-2)-M
  335.   0**4*SM**(-4))
  336.   +LG(M0)*PI**4*(M**2*M0**2*SM**(-4)+M0**2*SM**(-2)+5./2.*M0**4*SM
  337.   **(-4))
  338.   +LG(M0)*LG(SM)*PI**4*(-2.*M**2*M0**2*SM**(-4)-M0**2*SM**(-2)-M0*
  339.   *4*SM**(-4))
  340.   +LG(SM)*PI**4*(-4.*M**2*M0**2*SM**(-4)-M**2*SM**(-2)-1./2.*M**4*
  341.   SM**(-4)-M0**2*SM**(-2)-5./2.*M0**4*SM**(-4))
  342.   +LG(SM)*LG(SM)*PI**4*(1./2.+2.*M**2*M0**2*SM**(-4)+M0**2*SM**(-2
  343.   )+M0**4*SM**(-4))
  344. Al,F411(M~,M0~,SM)=
  345.   +PI**4*(-1./3.*M**(-4)+1./3.*M**(-4)*EP**(-1)+1./2.*M**(-2)*M0**
  346.   2*SM**(-4)+1./6.*M**(-2)*SM**(-2)+1./6.*SM**(-4))
  347.   +LG(M)*PI**4*(1./6.*M**(-4)+1./6.*SM**(-4))
  348.   +LG(M0)*PI**4*(-1./6.*M**(-4)*M0**2*SM**(-2)-1./6.*M**(-4)*M0**4
  349.   *SM**(-4)+1./3.*M**(-2)*M0**2*SM**(-4))
  350.   +LG(SM)*PI**4*(1./6.*M**(-4)+1./6.*M**(-4)*M0**2*SM**(-2)+1./6.*
  351.   M**(-4)*M0**4*SM**(-4)-1./3.*M**(-2)*M0**2*SM**(-4)-1./6.*SM**(-
  352.   4))
  353. Al,F321(M~,M0~,SM)=
  354.   +PI**4*(1./2.*M**(-2)*M0**2*SM**(-4)+1./2.*M**(-2)*SM**(-2)+3./2
  355.   .*SM**(-4))
  356.   +PI**6*(1./3.*SM**(-4))
  357.   +LG(M)*PI**4*(5./2.*SM**(-4))
  358.   +LG(M)*LG(M0)*PI**4*(SM**(-4))
  359.   +LG(M)*LG(SM)*PI**4*(-SM**(-4))
  360.   +LG(M0)*PI**4*(M**(-2)*M0**2*SM**(-4)+1./2.*M**(-2)*SM**(-2)+3./
  361.   2.*SM**(-4))
  362.   +LG(M0)*LG(SM)*PI**4*(-SM**(-4))
  363.   +LG(SM)*PI**4*(-M**(-2)*M0**2*SM**(-4)-1./2.*M**(-2)*SM**(-2)-4.
  364.   *SM**(-4))
  365.   +LG(SM)*LG(SM)*PI**4*(SM**(-4))
  366. Al,F421(M~,M0~,SM)=
  367.   +PI**4*(1./6.*M**(-4)*M0**2*SM**(-4)+1./6.*M**(-4)*SM**(-2)-5./6
  368.   .*M**(-2)*SM**(-4))
  369.   +LG(M0)*PI**4*(1./3.*M**(-4)*M0**2*SM**(-4)+1./6.*M**(-4)*SM**(-
  370.   2)-1./3.*M**(-2)*SM**(-4))
  371.   +LG(SM)*PI**4*(-1./3.*M**(-4)*M0**2*SM**(-4)-1./6.*M**(-4)*SM**(
  372.   -2)+1./3.*M**(-2)*SM**(-4))
  373. Al,F511(M~,M0~,SM)=
  374.    +PI**4*(-5./24.*M**(-6)+1./6.*M**(-6)*EP**(-1)+1./8.*M**(-4)*M0*
  375.    *2*SM**(-4)+1./24.*M**(-4)*SM**(-2)-1./24.*M**(-2)*SM**(-4))
  376.    +LG(M)*PI**4*(1./12.*M**(-6))
  377.    +LG(M0)*PI**4*(-1./12.*M**(-6)*M0**2*SM**(-2)-1./12.*M**(-6)*M0*
  378.    *4*SM**(-4)+1./12.*M**(-4)*M0**2*SM**(-4))
  379.    +LG(SM)*PI**4*(1./12.*M**(-6)+1./12.*M**(-6)*M0**2*SM**(-2)+1./1
  380.    2.*M**(-6)*M0**4*SM**(-4)-1./12.*M**(-4)*M0**2*SM**(-4))
  381. ENDBLOCK
  382. BLOCK E1EZF{1,2,3,4,5}
  383. C 1=E1, 2=E2, 3=ET, 4=E1F, 5=EZF.
  384. '1'Id,E1=0.5*TIPI2/C**2-31*TIPI2/12 + E1F
  385. '2'Al,E2=-25*TIPI2/12+0.75*TIPI2/C**2-1.5*M**2*TIPI2/SM**2
  386. '2'      -0.75*M**2*TIPI2/C**4/SM**2-0.75*SM**2*TIPI2/M**2
  387. '2'    +TIPI2*EP*SM**2/M**2*(-9/32*PI/W3+25/32-3/8*LG(SM) )
  388. '3'Id,ET=TIPI2*(1.5*M**2/C**4/SM**2+0.25/C**2+0.5
  389. '3'      +3*M**2/SM**2+0.75*SM**2/M**2)
  390. '3'      +EP*TIPI2*(-0.25*M**2/C**4/SM**2-1/8/C**2-0.25
  391. '3'      -0.5*M**2/SM**2-3*SM**2/8/M**2
  392. '3'      +(0.25+1.5*M**2/SM**2)*LG(M)
  393. '3'      +(1/8/C**2+0.75*M**2/C**4/SM**2)*LG(M0)
  394. '3'      +3*SM**2/8/M**2*LG(SM) )
  395. '3'  -EP**2*TIPI2*(PI2*(-M**2/C**4/SM**2/32-1/C**2/192-1/96
  396. '3'    -M**2/SM**2/16-SM**2/M**2/64)
  397. '3'   -M**2/SM**2/8/C**4-1/C**2/16-1/8-M**2/SM**2/4-3*SM**2/M**2/16
  398. '3'   +LG(M)*(1/8+M**2/SM**2/4)
  399. '3'   +LG(M)*LG(M)*(-1/16-3*M**2/SM**2/8)
  400. '3'   +LG(M0)*(M**2/SM**2/C**4/8+1/C**2/16)
  401. '3'   +LG(M0)*LG(M0)*(-3*M**2/C**4/SM**2/16-1/C**2/32)
  402. '3'   +LG(SM)*3*SM**2/M**2/16
  403. '3'   -LG(SM)*LG(SM)*3*SM**2/M**2/32)
  404. '4'Id,E1F=EP*TIPI2*SM**2/M**2*(1/32 + EP*(LG(SM)/64-3/128) )
  405. '5'Al,EZF=EP*TIPI2*SM**2/M**2*(1/32 + EP*(LG(SM)/64-3/128) )
  406. ENDBLOCK
  407. BLOCK PROJ{1,2,3,4}
  408. C 1=PROJ AND RATIONALIZATION P, 2=MAKE C**2=1-S**2,
  409. C 3=MAKE S**2=1-C**2, 4=ALSO RATIONALIZATION K AND Q AND COUNT SM.
  410. '1'Id,LEDLE=D(MU,NU)*D(MUP,NUP)/N*LE2
  411. '1'Id,Dotpr,LE(MU~)=FALE(MU)
  412. '1'Al,LE(MU~)=FALE(MU)
  413. '1'Id,Adiso,FALE(L1~)*FALE(L2~)=D(L1,AL)*D(L2,BE)*FACT
  414. '1'Id,FACT=TRN*(2*D(MU,NU)*D(AL,BE) - N*D(MU,AL)*D(NU,BE)
  415. '1'        -N*D(MU,BE)*D(NU,AL) ) * LE(MUP)*LE(NUP)
  416. '1'       +TRN*(-NP1*D(MU,NU)*D(AL,BE) + D(MU,AL)*D(NU,BE)
  417. '1'             + D(MU,BE)*D(NU,AL) ) * D(MUP,NUP) * LE2
  418. '1'*yep
  419. '4'Id,QDP=0.5*KDK-0.5*PDP-0.5*QDQ
  420. '4'Al,QDK=-0.5*PDP+0.5*KDK+0.5*QDQ
  421. '4'Al,PDK=-0.5*QDQ+0.5*KDK+0.5*PDP
  422. '2'Id,Multi,C**2=1-S**2
  423. '2'Al,C2MS2=1-2*S**2
  424. '3'Id,Multi,S**2=1-C**2
  425. '3'Al,C2MS2=2*C**2-1
  426. '1'*yep
  427. '1'Id,PDP**K1~*FP(L1~,M~)=(-1)**K1*M**K1*M**K1*
  428. '1'  DS(IX1,0,K1+(L1-1-K1)*DT(K1-L1+1),(M**-IX1*M**-IX1*FP(L1-IX1,M)),
  429. '1'          (-(K1-IX1+1)/IX1) )
  430. '4'   +EI*(-1)**K1*(-1)**L1*M**K1*M**K1/M**L1/M**L1*
  431. '4'   DT(K1-L1)*DB(K1-1,K1-L1)* DS(IX2,0,K1-L1,
  432. '4'     (PDP**IX2/M**IX2/M**IX2),(-(K1-L1-IX2+1)/(K1-IX2)) )
  433. '4'Al,KDK**K1~*FK(1,M~)=(-1)**K1*M**K1*M**K1*FK(1,M)
  434. '4'   -EI*(-1)**K1*M**K1*M**K1/M**2*DT(K1-1)*
  435. '4'    DS(IY2,0,K1-1,(KDK**IY2/M**IY2/M**IY2),(-1) )
  436. '4'Al,QDQ**K1~*FQ(L1~,M~)=(-1)**K1*M**K1*M**K1*
  437. '4'  DS(IZ1,0,K1+(L1-1-K1)*DT(K1-L1+1),(M**-IZ1*M**-IZ1*FQ(L1-IZ1,M)),
  438. '4'          (-(K1-IZ1+1)/IZ1) )
  439. '4'   +EI*(-1)**K1*(-1)**L1*M**K1*M**K1/M**L1/M**L1*
  440. '4'   DT(K1-L1)*DB(K1-1,K1-L1)* DS(IZ2,0,K1-L1,
  441. '4'     (QDQ**IZ2/M**IZ2/M**IZ2),(-(K1-L1-IZ2+1)/(K1-IZ2)) )
  442. '4'Id,PDP**2=QDQ**2+KDK**2+2*QDQ*KDK+4/N*QDQ*KDK
  443. '4'Al,PDP=QDQ+KDK
  444. '4'Al,KDK**2=QDQ**2+PDP**2+2*PDP*QDQ+4/N*QDQ*PDP
  445. '4'Al,KDK=QDQ+PDP
  446. '4'Al,QDQ**2=KDK**2+PDP**2+2*PDP*KDK+4/N*PDP*KDK
  447. '4'Al,QDQ=KDK+PDP
  448. '4'C FACTORS WITH PDP, QDQ OR KDK HAVE ONE EI (NOTE EI**2=0).
  449. '4'Id,KDK**K1~*FK(1,M~)=(-1)**K1*M**K1*M**K1*FK(1,M)
  450. '4'Al,PDP**K1~*FP(L1~,M~)=(-1)**K1*M**K1*M**K1*
  451. '4'  DS(IX3,0,K1+(L1-1-K1)*DT(K1-L1+1),(M**-IX3*M**-IX3*FP(L1-IX3,M)),
  452. '4'          (-(K1-IX3+1)/IX3) )
  453. '4'Al,QDQ**K1~*FQ(L1~,M~)=(-1)**K1*M**K1*M**K1*
  454. '4'  DS(IZ3,0,K1+(L1-1-K1)*DT(K1-L1+1),(M**-IZ3*M**-IZ3*FQ(L1-IZ3,M)),
  455. '4'          (-(K1-IZ3+1)/IZ3) )
  456. '4'Id,Adiso,FK(1,M~)*FQ(1,SM)=FK(1,SM)*FQ(1,M)
  457. '4'Al,Adiso,FK(1,M~)*FP(1,SM)=FK(1,SM)*FP(1,M)
  458. '4'Id,Numer,EI,1.
  459. '4'Id,Commu,FK,FP,FQ
  460. '4'Id,FP(3,M~)*FQ(2,S~)=FP(2,S)*FQ(3,M)
  461. '4'Al,FP(3,M~)*FQ(1,S~)=FP(1,S)*FQ(3,M)
  462. '4'Al,FP(2,M~)*FQ(1,S~)=FP(1,S)*FQ(2,M)
  463. '4'*yep
  464. '4'Id,FK(1,C~)*FP(5,M~)*FQ(1,S~)=F511(M,S,C)
  465. '4'Al,FK(1,C~)*FP(4,M~)*FQ(1,S~)=F411(M,S,C)
  466. '4'Id,FK(1,S~)*FP(2,C~)*FQ(4,M~)=F421(M,C,S)
  467. '4'Al,FK(1,S~)*FP(2,C~)*FQ(3,M~)=F321(M,C,S)
  468. '4'Al,FK(1,S~)*FP(2,C~)*FQ(2,M~)=F221(M,C,S)
  469. '4'Al,FK(1,S~)*FP(1,C~)*FQ(4,M~)=F411(M,C,S)
  470. '4'Al,FK(1,S~)*FP(1,C~)*FQ(3,M~)=F311(M,C,S)
  471. '4'Al,FK(1,S~)*FP(1,C~)*FQ(2,M~)=F211(M,C,S)
  472. '4'Al,FK(1,S~)*FP(1,C~)*FQ(1,M~)=F111(C,M,S)
  473. '1'Id,FP(N~,M~)=FSC(N,M)*M**4/M**N/M**N
  474. '4'Al,FQ(N~,M~)=FSC(N,M)*M**4/M**N/M**N
  475. '4'Al,FK(N~,M~)=FSC(N,M)*M**4/M**N/M**N
  476. '4'Id,Count,1,SM,1
  477. '1'*yep
  478. ENDBLOCK
  479. BLOCK NEXP{}
  480. Id,N**2*NP1*TRN =  - 10./9. + 1./18.*EP - 7./324.*EP**2
  481. Al,N*NP1*TRN =  - 5./18. + 1./12.*EP - 17./648.*EP**2
  482. Al,NP1*TRN =  - 5./72. + 11./288.*EP - 167./10368.*EP**2
  483. Al,N**3*TRN =  - 8./9. - 1./162.*EP**2
  484. Al,N**2*TRN =  - 2./9. + 1./18.*EP - 5./324.*EP**2
  485. Al,N*TRN =  - 1./18. + 1./36.*EP - 7./648.*EP**2
  486. Al,TRN =  - 1./72. + 1./96.*EP - 55./10368.*EP**2
  487. Al,N*OMN** - 1 =  - 4/3 + EP/9 - EP**2/27
  488. Al,N = 4 + EP
  489. Al,N**-2 = 1/16 - EP/32 + 3*EP**2/256
  490. Al,N**-1 = 1/4 - EP/16 + EP**2/64
  491. ENDBLOCK
  492. End
  493.