double r_abs(x)
float *x;
{
if(*x >= 0)
	return(*x);
return(- *x);
}