/* 4 -1 {a,s,l,h}=stats(...), avg, stdev, min, and max of the defined cells
*/

double rf_stats(const Node *n, const Cell *c)
{
  int i;
  double a[4];
  Node *lval = Left(n);

  stats_helper( Right(n), c, a);

  a[1] = sqrt(a[1]);

  for( i = 0; lval && i<4; ++i)
  {
    lval_helper( "stats", lval, c, a[i]);
    lval = lval->next;
  }

  return a[0];
}