# sgn(x) is sign function; sgn(x_) := If(x>0,1,If(x<0,-1));