(* ***
example.m

Sample Mathematica file for usage of decsusy3l.m by
Alexander Kurz, Matthias Steinhauser and Nikolai Zerf
*** *)

<<decsusy3l.m;

(num = {MB1->mst2, MB2->mst2, 
	mudec->1000.,
	mt->146.7,
	thet->1.26,Cthetat->Cos[thet],Sthetat->Sin[thet],
	mst1->1022.,mst2->1271.,mgl->1326.,msq->1348.})

(repMassdiff = {Dm2B1gl -> MB1^2 - mgl^2,
		Dm2B1sq -> MB1^2 - msq^2,
		Dm2B1st1 -> MB1^2 - mst1^2,
		Dm2B1st2 -> MB1^2 - mst2^2,
		Dm2B2st1 -> MB2^2 - mst1^2,
		Dm2B2t -> MB2^2 - mt^2
		});

{ZetaAlphas[5,"h1"],ZetaAlphas[5,"h2"],ZetaAlphas[5,"h3"]}/.repMassdiff//.num

ZetaAlphas[5,"exact"]/.repMassdiff//.num
