double sq(double a)
{
  return a*a;
}