In[]:=
<<X`
In[]:=
eps=0.5mpi=0.139mi=10.86(*+0.005I*)mf=10.3552z=0.4650000000000001sa=2*mpi^2+mf^2+(mi+Ieps)^2-2mz^2Y[s_,mz_]=s+2*mz^2-mf^2-(mi+Ieps)^2-2*mpi^2sigma[s_]=Sqrt[1-(4*mpi^2/s)]κ[s_]=Expand[Sqrt[KallenExpand[Kallenλ[s,(mi+Ieps)^2,mf^2]]]*sigma[s]]log[s_,mz_]=Log[(Y[s,mz]+κ[s])/(Y[s,mz]-κ[s])]+2I*Pi*HeavisideTheta[Re[sa]-s]discC0[s_,mz_]=(-1/Expand[Sqrt[KallenExpand[Kallenλ[s,(mi+Ieps)^2,mf^2]]]])log[s,mz]X[mz_]=Table[(discC0[s,mz]/(s-z)),{s,4*mpi*mpi,0.14,(0.14-(4*mpi*mpi))/20}]/.mz10.47
-2
10