function y = hwrect(x);
% y=hwrect(x) - half wave rectify
%

y = (x + abs(x)) ./2;