2
Reference: Computer Physics Communications,
54
, 295-305 (1989)