#ifndef PYTHAG_H
#define PYTHAG_H
//----------------------------------------------------------------------------
#include <cmath>
#define NRANSI
#include "nrutil.h"
using namespace std;

double pythag(double a, double b)
{
   double absa,absb;
   absa=fabs(a);
   absb=fabs(b);
   if (absa > absb) return absa*sqrt(1.0+pow(absb/absa, 2.0));
   else return (absb == 0.0 ? 0.0 : absb*sqrt(1.0+pow(absa/absb, 2.0)));
}

#undef NRANSI
/* (C) Copr. 1986-92 Numerical Recipes Software ;#. */
//----------------------------------------------------------------------------
#endif
// Fin------------------------------------------------------------------------