#======================================================================

#                    P Y T H A G . P L 

#                    doc: Sun Aug  1 10:41:34 1999

#                    dlm: Wed Mar 11 14:53:21 2015

#                    (c) 1999 A.M. Thurnherr

#                    uE-Info: 27 2 NIL 0 0 72 0 2 4 NIL ofnI

#======================================================================

# PYTHAG routine from Numerical Recipes adapted to ANTS

# HISTORY:

#	Aug 01, 1999: - manually converted from c-source

#	Mar 11, 2015: - BUG: syntax errors (code had never been used)

sub pythag($$)

{

	my($a,$b) = @_;							# params

	my($absa,$absb);						# float 

	$absa = abs($a);

	$absb = abs($b);

	return $absa*sqrt(1.0+SQR($absb/$absa))

		if ($absa > $absb);

	return ($absb == 0 ? 0 : $absb*sqrt(1+$absa*$absa/$absb/$absb));

}

1;