function y = heaviside_f( x )
  y = 0*x ;
  y( x>=0 ) = 1 ;
end