This command calculates the bond angle and its standard deviation
for the atoms numbered <a1>, <a2> and <a3>. The angle is measured
at <a2>. Note that no periodic boundaries are applied and it might
be necessary to add +-1 to the coordinates of some of the atoms to
get the desired angle (-> variables).