% non-smooth right-hand side

function yDot = fNonSmooth(t,y)

global tau

yDot(1) = rem(t,tau);