#include <math.h>static double xhypot (const double x, const double y);static double xhypot (const double x, const double y){  double xabs = fabs(x) ;  double yabs = fabs(y) ;  double min, max;  if (xabs < yabs) {    min = xabs ;    max = yabs ;  } else {    min = yabs ;    max = xabs ;  }  if (min == 0)     {      return max ;    }  {    double u = min / max ;    return max * sqrt (1 + u * u) ;  }}