Rework moments routines for improved numerical stability based on

Rework moments routines for improved numerical stability based on
Michael Kiefte's advice.  Any bugs or remaining numerical problems are
still mine though.

There is now a struct moments1 for use with one-pass moments.  It uses
a provisional means algorithm as an attempt to improve accuracy of
higher moments.  The older struct moments now only handles two-pass
moments.