Skip to content
GitLab
Switch branch/tag
Find file Normal viewHistoryPermalink
pythag.c 347 Bytes
Newer Older
Fang Yuedong's avatar
Baseline  
Fang Yuedong committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14 
#include <math.h>
#define NRANSI
#include "nrutil.h"

double pythag(double a, double b)
{
	double absa,absb;
	absa=fabs(a);
	absb=fabs(b);
	if (absa > absb) return absa*sqrt(1.0+SQR(absb/absa));
	else return (absb == 0.0 ? 0.0 : absb*sqrt(1.0+SQR(absa/absb)));
}
#undef NRANSI
/* (C) Copr. 1986-92 Numerical Recipes Software )1!. */