function s = fsigma(t)

indx1 = find(t<0);
if length(indx1)>0
   s(indx1) = 0;
end

indx2 = find(t>=0);
if length(indx2)>0
   s(indx2) = 1;
end
