# sgn(x) is sign function;

sgn(x_) := If(x>0,1,If(x<0,-1));