function y=HEAVISIDE(v)
   if (v<0) y=0;
   else y=1;
   endif
endfunction