function y =l_sat(x)
y=sign(x).*min(1,abs(x));