half-wave-rectified-sine code snippets in matlab