function y=sat(x,a)
%Saturation at level pm a y=sat(x,a)
y=min(max(x,-a),a);   