function r=rmp_fn(t)
r=0.5*t.*(sign(t)+1);