void averages(double u,double v,double *p,double *q)
{
  *p = (u + v) / 2;
  *q = sqrt(u * v);
}