sgn(a)

           function s = sgn(a)
           
           if a > 0
             s = 1;
           elseif a < 0
             s = -1;
           else
             s = 0;
           end;