……这个代码很有可能是机械移植自其他的语言,这个Module用的可谓毫无意义.原本在{}写上变量名的话,是可以让里面的变量局部化的,可是它留空了……对于那段完全可以剃成这样：
Tu[en_, [Phi]n_] := (e := en; [CurlyPhi] := 2*Pi*[Phi]n/4; [Theta] := ArcCos[-(((e - [CurlyEpsilon]0)/v0)/2.0)]; t1 := tu; t1)
当然其实还可以更简单,我懒得想了.
flux := .0嘛,自然就是flux为0咯,这种奇怪的写法显然也是移植的痕迹,在这里改成flux = 0也毫无问题.