double ab(double a)
{
  if (a>=0)
    return a;
  else
    return -a;
}