function y=Heaviside(x)
y = (x >= 0);

