dpythag.c

来自「< C语言数值算法程序大全>>配套程序」· C语言 代码 · 共 13 行

#include <math.h>
#include "nrutil.h"

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