real function signum(x)
real x
	if      (x>0)  { signum =  1 }
	else if (x<0)  { signum = -1 }
	else           { signum =  0 }
return;end