New anisotropy function.

Now reproduces CHARMM code exactly; removed unused old version of
the function.

Change-Id: I8d8475d4ae5e991066ad3cb6a2736fe7fc76a090