% Step function 
function y = StepFun(n)
y = zeros(size(n));
y (find(n>=eps)) = 1;
end