wolffd@0: function y = hwr(x)
wolffd@0: % Half-Wave Rectifier
wolffd@0: 
wolffd@0: y = 0.5 * (x + abs(x));