// set up default output format and precision of digits after decimal
void setOutputFormat(int digits)
{
  cout.precision(digits);
  cout.setf(ios::fixed | ios::dec | ios::showpoint);
}