procedure EE(m,m1,m2); (m^2+m1^2-m2^2)/(2*m)$ procedure PP(m,m1,m2); sqrt((m+m1+m2)*(m-m1-m2)*(m+m1-m2)*(m-m1+m2))/(2*m)$ procedure JJ(s,m1,m2); sqrt(s)*PP(sqrt(s),m1,m2)$ procedure ph2(m,m1,m2); PP(m,m1,m2)/(4*pi*m)$ procedure Dalitz(m,m1,m2,m3,E1,E2m,E2p); begin scalar a,b,d; a:=2*m*E1-m^2-m1^2; b:=(m-E1)*(a-m2^2+m3^2); d:=sqrt((a+(m2+m3)^2)*(a+(m2-m3)^2))*sqrt(E1^2-m1^2); set(E2m,(b-d)/(2*a)); set(E2p,(b+d)/(2*a)); end$ procedure Mandel(m1,m2,m3,m4,s,tm,tp); begin scalar b,d; b:=s^2-(m1^2+m2^2+m3^2+m4^2)*s+(m1^2-m2^2)*(m3^2-m4^2); d:=sqrt((s-(m1+m2)^2)*(s-(m1-m2)^2) *(s-(m3+m4)^2)*(s-(m3-m4)^2)); set(tm,-(b-d)/(2*s)); set(tp,-(b+d)/(2*s)); end$ end;